SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90Linked to original sources

USGS2012

Petrologic evolution of divergent peralkaline magmas from the Silent Canyon caldera complex, southwestern Nevada volcanic field

The Silent Canyon volcanic center consists of a buried Miocene peralkaline caldera complex and outlying peralkaline lava domes. Its location has been corroborated by geophysical data and more than 50 drill holes. Two widespread ash flow sheets, the Tub Spring and overlying Grouse Canyon members of the Miocene Belted Range Tuff, were erupted from the caldera complex and have volumes of 60–100 km 3 and 200 km 3 , respectively. Eruption of the ash flows was preceded by widespread extrusion of precaldera comendite domes and was followed by extrusion of postcollapse peralkaline lavas and tuffs within and outside the caldera complex. Lava flows and tuffs were also deposited between the two major ash flow sheets. Rocks of the Silent Canyon center vary significantly in silica content and peralkalinity. The most mafic rocks are precollapse and postcollapse trachytes (65–69% SiO 2 ). Low-silica comendites (69–73% SiO 2 ) were erupted as the mafic upper part of the chemically zoned Grouse Canyon Member and as postcollapse lavas. The lower part of the Grouse Canyon Member and the underlying rhyolite of Split Ridge are moderately peralkaline comendite (PI is molar ratio Na + K/Al is 1.17–1.26). These comendites have major element characteristics and trace element enrichments approaching those of pantellerites. The Tub Spring Member, by contrast, is a weakly peralkaline chemically unzoned silicic comendite (75–76% SiO 2 ) ash flow tuff. Weakly peralkaline silicic comendites (PI 1.0–1.1) are the most abundant precaldera lavas. Postcollapse lavas range from trachyte to silicic comendite; some have anomalous light rare earth element (LREE) enrichments. Silent Canyon rocks follow a common petrologic evolution from trachyte to low-silica comendite; above 73% SiO 2 , compositions of the moderately peralkaline comendites diverge from those of the weakly peralkaline silicic comendites. These contrasting differentiation paths are shown in the behavior of Fe and other transition metals, Al, Na, K; the trace elements Ba, Zr, Nb; and probably F and Cl. Weakly peralkaline silicic comendites show a LREE/heavy REE crossover in early erupted/late erupted rocks; moderately peralkaline comendites are enriched in all REE. The development of divergent peralkaline magmas, toward both pantelleritic and weakly peralkaline compositions, is unusual in a single volcanic center.

Journal of Geophysical Research Solid Earth↗

USGS2012

State of stress and modern deformation of the northern Basin and Range Province

Constraints on the current stress regime of the actively extending northern Basin and Range province are provided by deformation data (focal mechanisms and fault slip studies), hydraulic fracturing in situ stress measurements, borehole elongation (“breakouts”) analyses, and alignment of young volcanic vents. The integrated data indicate significant variations both in principal stress orientations and magnitudes. An approximately E-W least principal stress direction appears to characterize both the eastern and western margins of the Basin and Range province, whereas in the active interior parts of the province extension occurs in response to a least principal stress oriented NW to N60°W. The contrast in stress orientations between the province boundaries and in the interior suggests that along the margins the least principal stress direction may be locally controlled by the generally northerly trending profound lithospheric discontinuities associated with these margins. Active deformation along the southeastern and western province margins is characterized by a combination of strike-slip and normal faulting. Focal mechanisms along northeastern province margin (Wasatch front) and in central Nevada indicate a combination of normal and oblique-normal faulting. Temporal, regional, and depth-dependent variations in the relative magnitudes of the vertical and maximum horizontal stresses can explain much of the observed variations in deformation styles. However, some depth variation in faulting style inferred from focal mechanisms may be apparent and simply a function of the attitude of fault planes being reactivated. Evidence for significant temporal variation (or multiple cycles of variation) in relative stress magnitude comes from the Sierran front-Basin and Range boundary region where recent earthquakes are predominantly strike slip, whereas the profound relative vertical relief across the Sierra frontal fault zone in the last 9–10 m.y. implies a normal faulting stress regime. Using the best data on stress orientation, relative stress magnitudes are constrained from slip vectors of major earthquakes and young fault displacements. Analysis of well-constrained slip vectors in the Owens Valley, California, area indicate that large temporal variations in the magnitude of the approximately N-S oriented maximum horizontal stress are required to explain dominantly dip-slip and strike-slip offsets on subparallel faults. Similar faulting relations are observed throughout much of the boundary zone between the Basin and Range-Sierra Nevada (including the Walker Lane belt). Along the eastern province margin in the Wasatch front area in Utah, available data suggest that the maximum and minimum horizontal stresses may be approximately equal at depths of <4–5 km. Earthquake focal mechanisms in this area suggest more variability in relative magnitude of the two horizontal stresses with depth. Furthermore, superimposed sets of young fault striae along a segment of the Wasatch fault also indicate temporal variations of relative stress magnitudes. Sources of regional and temporal variations in the stress field may be linked to variable shear tractions applied to the base of the brittle crust related to intrusion, thermally induced flow, and the influence of the San Andreas plate boundary. Although difficult to date accurately, the fault slip data suggest that the temporal variations in relative magnitudes stress may occur on the time scale of both a single major earthquake cycle (1000–5000 years) and multiple earthquake cycles (10,000+ years).

Journal of Geophysical Research Solid Earth↗

USGS2012

Compaction of basin sediments: Modeling based on time-temperature history

Porosity decrease of sandstones, carbonates, and shales during burial is modeled here in terms of time-temperature exposure rather than in terms of depth. Loss of porosity (ϕ) in the subsurface is represented by a power function, ϕ = A ( M ) B , where A and B are constants and M is a measure of integrated time-temperature history. Regression lines of carbonate and sandstone porosity upon Lopatin's time-temperature index of thermal maturity (TTI) generally fall within a rather narrow envelope whose axis is approximated by ϕ = 30(TTI) −0.33 . This equation is useful for regional modeling of carbonate and sandstone compaction. Dependence of porosity upon integrated time-temperature history implies that basin sediments are not in equilibrium but compact through geologic time. Calculations show that subsidence resulting from the loss of porosity with increasing time-temperature exposure (as opposed to deeper burial) can produce a second-stage, passively formed basin containing many hundreds of meters of sediments.

Journal of Geophysical Research Solid Earth↗

USGS2012

Thermal infrared (2.5-13.5 μm) spectroscopic remote sensing of igneous rock types on particulate planetary surfaces

Fundamental molecular vibration bands are significantly diminished by scattering. Thus such bands in spectra of fine particulate regoliths (i.e., dominated by <5-μm particles), or regoliths displaying a similar scale of porosity, are difficult to use for mineralogical or rock type identification. Consequently, other spectral features have been sought that may be more useful in spectroscopic remote sensing of composition. We find that mineralogical information is retained in overtones and combination tones of the fundamental molecular vibrations in the 3.0- to 7.0-μm region, but that relatively few minerals have a sufficiently distinctive band structure to be unambiguously identified with currently available techniques. More significantly, identification of general rock type, as defined by the SCFM chemical index (SCFM = SiO 2 /SiO 2 + CaO + FeO + MgO), is possible using spectral features associated with the principal Christiansen frequency and with a region of relative transparency between the Si-O stretching and bending bands. However, environmental factors may affect the appearance and wavelengths of these features. Finally, prominent absorption bands may result from the presence of relatively small amounts of water, hydroxyl or carbonate, because absorption bands exhibited by these materials in the 2.7- to 4.0-μm region, where silicate spectra are otherwise featureless, increase strongly in spectral contrast with decreasing particle size. Such materials are thus detectable in very small amounts in a particulate regolith composed predominantly of silicate minerals.

Journal of Geophysical Research↗

USGS2012

The competition between thermal contraction and differentiation in the stress history of the Moon

The scarcity of both extension and compression features on the Moon strongly constrains the history of the lunar radius—to variations of less than ±1 km over the past 3.8 Gyr. This limit has traditionally been interpreted as requiring a delicate balance between thermal contraction of the near‐surface and expansion of a substantial cold interior region. Recent theories of lunar origin (e.g., giant impact), in contrast, favor a “hot” initial state. We propose that a reconciliation may be possible by taking account of the volume change Δ V / V | d due to differentiation. We calculate STP densities based on simplified normative mineralogies for a suite of estimates of the bulk lunar composition, of primary lunar basalt, and of the residuum left when the maximum amount of the latter is extracted from the former. Typically Δ V / V | d ≃ 2 to 5%—an expansion equivalent to heating by ∼10 3 K. Provided the timing of differentiation is correct, one might offset the cooling of a magma ocean as much as 630 km deep by differentiation of the remainder of the Moon (which need not start much below the solidus temperature). A large but not impossible amount of gabbroic melt production is implied: ∼100 times the volume of mare basalts known to have been extruded. We do not address the detailed genetic relationship of this melt to the basalts observed on the lunar surface but point out that it need not have reached the surface directly or even have entered the crust in order for the expansion to have occurred. To assess the timing of melt formation, we investigate a simple conductive lunar thermal model which takes account of both Δ V / V | d and thermal contraction. Our initial state is characterized by a central temperature T c and a depth Z 0 above which the material (derived from the magma ocean) is already at the solidus and is not suceptible to volume changes upon further differentiation. We find a range of models satisfying the limits on radius increase and decrease. The hottest has T c = 1210 K, Z 0 = 400 km; without Δ V / V | d , we would need a larger or colder (or both) core, e.g., T c ≲ 700 K for Z 0 = 200–400 km, in agreement with previous investigators. Our modeling thus lends credence to the idea that the Moon could have been initially ≳50% molten (with the remainder relatively close to the solidus) and yet experienced little volume change over the last 3.8 Gyr.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Dynamic rupture modeling with laboratory-derived constitutive relations

A laboratory-derived state variable friction constitutive relation is used in the numerical simulation of the dynamic growth of an in-plane or mode II shear crack. According to this formulation, originally presented by J. H. Dieterich, frictional resistance varies with the logarithm of the slip rate and with the logarithm of the frictional state variable as identified by A. L. Ruina. Under conditions of steady sliding, the state variable is proportional to (slip rate) −1 . Following suddenly introduced increases in slip rate, the rate and state dependencies combine to produce behavior which resembles slip weakening. When rupture nucleation is artificially forced at fixed rupture velocity, rupture models calculated with the state variable friction in a uniformly distributed initial stress field closely resemble earlier rupture models calculated with a slip weakening fault constitutive relation. Additional rupture models are calculated in which rupture nucleation is achieved naturally, with numerical simulations of the quasi-static response of the fault leading to the onset of unstable, dynamic rupture. When rupture nucleation with the state variable friction law takes place naturally, a large fraction of the fault accelerates before accelerating slip is concentrated in what ultimately becomes the rupture nucleation patch. The state evolution accompanying this accelerating slip leads to higher average rupture speeds or a more rapid rupture acceleration to near P wave rupture speeds. Rupture models are also calculated for the seismological asperity problem, that is, the failure of a highly stressed fault patch surrounded by a region of zero stress drop. Dynamic overshoot of slip into the region of zero stress drop roughly agrees with a simple energy balance analysis; the final size of the rupture is proportional to the square of the size of the high stress patch. Earlier frictional stability analyses have led to the definition of a critical fault patch size for rupture nucleation. This critical patch size is generally different from critical crack lengths determined from crack tip energy balance considerations applied to a simpler slip weakening law. In the model calculations, dynamic rupture does not nucleate if the starting patch size is less than the critical patch size. This is consistent with the frictional stability analyses. Thus these model calculations suggest that dynamic rupture following a state variable friction relation is similar to that following a simpler fault slip weakening law. However, when modeling the full cycle of fault motions, rate-dependent frictional responses included in the state variable formulation are important at low slip rates associated with rupture nucleation. The critical rupture nucleation dimension appropriate for a slip weakening fault does not predict the critical nucleation dimension for a state variable fault.

Journal of Geophysical Research Solid Earth↗

USGS2012

Hydrologic effects on water level changes associated with episodic fault creep near Parkfield, California

As part of the Parkfield, California, earthquake prediction experiment, water level is monitored in a well 460 m from the main trace of the San Andreas fault on Middle Mountain, in the preparation zone of the anticipated Parkfield earthquake. The well configuration allows water level to be monitored in two fluid reservoirs at depths of 85 and 250 m below land surface. During 1987, water level changes were recorded during 12 of the 18 episodes of accelerated fault creep detected by a creep meter spanning the fault trace 750 m northwest of the well. The creep-related water level changes in the shallow reservoir have durations of less than 1 day, whereas in the deeper reservoir the changes persist for as long as 2 months. These data suggest that the transient nature of the water level changes in the shallow interval is due to vertical flow to the water table and is not evidence that creep events propagate past the well. Phase leads of earth tidal constituents in the water level data from the shallow interval relative to the same constituents in the local volume strain tide support the interpretation of significant flow to the water table at periods of 1 day or less. The form of the water level changes in the deep interval is affected by horizontal flow to the well bore. This effect can be removed from the water level records using a theoretical response curve constrained by the phases of earth tidal constituents in the deep interval relative to the local volume strain tide. For the events where the signal in the shallow interval has been large enough to measure, the sizes of the simultaneous water level changes in the two reservoirs are consistent with the same amounts of volume strain occurring at both depths.

Journal of Geophysical Research Solid Earth↗

USGS2012

Synfolding magnetization in the Jurassic Preuss Sandstone, Wyoming- Idaho-Utah thrust belt

The Jurassic Preuss Sandstone, exposed in five thrust plates of the Wyoming-Idaho-Utah thrust belt, carries directions of remanent magnetization that group most tightly after only partial unfolding. Field, petrographic, and rock magnetic evidence indicates that the carrier of this magnetization is detrital, low-Ti titanomagnetite. The detrital titanomagnetite was remagnetized at low temperatures (75°–150°C) probably completely during folding. Anisotropy of magnetic susceptibility and petrographic observations indicate that the detrital titanomagnetite has been affected by tectonic strain. We suggest that low-temperature remagnetization of the detrital titanomagnetite was either a viscous partial thermoremanent magnetization, the acquisition of which was enhanced by stress, or a piezoremanent magnetization that involved stress-induced movement of domain walls during intracrystalline strain, or was a combination of the two mechanisms. Stress may promote remagnetization at temperatures much lower than predicted by current theoretical models. Other mechanisms, such as acquisition of chemical remanent magnetization during folding, deflection of a prefolding magnetization by internal strain, or combination of components of magnetization with different direction cannot account for the geometry of magnetization in the Preuss. The locus of acquisition of synfolding magnetization in the Preuss migrated in conjunction with deformation in the thrust belt. A model is presented in which synfolding magnetization was acquired during cooling and folding as strata moved up thrust ramps. A lack of reverse-polarity directions remains a puzzling feature of the remanence. The remanent direction is tentatively interpreted to reflect the predominant polarity state during its acquisition over an extended rather than a discrete time period during folding in Late Cretaceous and early Tertiary (?) periods of predominantly normal polarity.

Journal of Geophysical Research Solid Earth↗

USGS2012

Spectral characteristics of chlorites and Mg‐serpentines using high‐resolution reflectance spectroscopy

The present laboratory study using high‐resolution reflectance spectroscopy (0.25–2.7 μm) focuses on two primary phyllosilicate groups, serpentines and chlorites. The results show that it is possible to spectrally distinguish between isochemical end‐members of the Mg‐rich serpentine group (chrysotile, antigorite, and lizardite) and to recognize spectral variations in chlorites as a function of Fe/Mg ratio (∼8–38 wt % Fe). The position and relative strength of the 1.4‐μm absorption feature in the trioctahedral chlorites appear to be correlated to the total iron content and/or the Mg/Si ratio and the loss on ignition values of the sample. Spectral differences in the 2.3‐μm wavelength region can be attributed to differences in lattice environments and are characteristic for specific trioctahedral chlorites. The 1.4‐μm feature in the isochemical Mg‐rich serpentines (total iron content ∼1.5–7.0 wt%) show marked spectral differences, apparently due to structural differences.

Journal of Geophysical Research: Solid Earth↗

USGS2012

Teleseismically recorded seismicity before and after the May 7, 1986, Andreanof Islands, Alaska, earthquake

The May 7, 1986, Andreanof Islands earthquake ( M w 8.0) is the largest event to have occurred in that section of the Aleutian arc since the March 9, 1957, Aleutian Islands earthquake ( M w 8.6). Teleseismically well-recorded earthquakes in the region of the 1986 earthquake are relocated with a plate model and with careful attention to the focal depths. The data set is nearly complete for m b ≥4.7 between longitudes 172°W and 179°W for the period 1964 through April 1987 and provides a detailed description of the space-time history of moderate-size earthquakes in the region for that period. Additional insight is provided by source parameters which have been systematically determined for M w ≥5 earthquakes that occurred in the region since 1977 and by a modeling study of the spatial distribution of moment release on the mainshock fault plane. A technically significant component of oblique convergence in the central Aleutian arc results in its breakup into clockwiserotating and westward translating blocks. The western part of the Andreanof block is distinct from and stronger than flanking regions. The greater strength of this block segment and strong coupling along the main thrust zone result in the accumulation of high levels of shear stress, which give rise to great earthquakes near its eastern boundary. The occurrence of the 1986 rupture only 29 years after the 1957 earthquake may indicate that in the central Aleutians M 7+ earthquakes ordinarily do represent the predominant mode of strain release. Segmentation of the main thrust zone into upper and lower planes is supported by spatial and temporal patterns of seismicity and by focal mechanism data. This and other lines of evidence indicate a downdip increase in fault strength and possibly in heterogeneity within the main thrust zone in subduction zone environments. Aftershocks of the 1986 earthquake were bounded to the west by the Andreanof block boundary (Adak Canyon) and to the east by an aseismic segment of the main thrust zone near the subducted extension of the Amlia fracture zone. The aftershock distribution was bounded to the south by a forearc shear zone and to the north by the base of the main thrust zone in the Hawley Ridge segment and by the downdip edge of the upper plane of the main thrust zone in the eastern segment. Aftershocks which occurred near the volcanic line at shallow crustal depths in the upper plate were triggered by the mainshock and manifest a partial decoupling of oblique slip in this region along a west-striking right-lateral fault with low shear strength. Aftershock clustering along the main thrust zone was very similar to the distribution of prior seismicity, suggesting a continuation of long-term processes and the existence of areas with distinct mechanical properties. Interconnecting regions of low seismicity during both the premainshock and aftershock periods coincided with areas of major moment release during the mainshock. Seismicity data prior to the 1986 mainshock gave few clues about the location of the mainshock nucleation point, the mainshock size, and its time of occurrence. However, a large part of the mainshock moment release did coincide with a zone of seismic quiescence monitored by the Adak local seismograph network.

Journal of Geophysical Research Solid Earth↗

USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗

USGS2012

Active faulting and deformation of the Coalinga anticline as interpreted from three-dimensional velocity structure and seismicity

This work gives a clear picture of the geometry of aftershock seismicity in a large thrust earthquake. Interpretation of hypocenters and fault plane solutions, from the 1983 Coalinga, Coast Range California, earthquake sequence, in combination with the three-dimensional velocity structure shows that the active faulting beneath the fold primarily consists of a set of southwest dipping thrusts uplifting blocks of higher-velocity material. Above the main listric blind thrust there is a conjugate fault, steeply northeast dipping, that provides the western limit of the aftershocks within the Coalinga Anticline and that corresponds in location and spatial extent with the adjacent Pleasant Valley syncline. The character of the seismicity varies with the degree of previous deformation on each section of the anticline. Where the previous uplift was largest, the shallow seismicity shows secondary faulting on either side of the fold with orientations that correspond to the preexisting geologic structure. Diffuse seismicity characterizes the area with the least previous deformation. The mainshock rupture terminated where the fold trend was no longer uniform but had competing north and west trending features. The upward extent of the mainshock rupture ended at the approximate boundary between Franciscan and Great Valley Sequence rocks. Above that depth the main thrust appears to splay into a steeper segment and a near-horizontal segment. Thus the extent of rupture area is limited by the area of uniform structural orientation and by the variation in the type of material. With the three-dimensional velocity model each individual hypocenter moved slightly (0–2 km) in accord with the details of the surrounding velocity structure, so that secondary features in the seismicity pattern are more detailed than with a local one-dimensional model and station corrections. The overall character of the fault plane solutions was not altered by the three-dimensional model, but the more accurate ray paths did result in distinct changes. In particular, the mainshock has a fault plane dipping 30° southwest instead of the 23° obtained with the one-dimensional model.

Journal of Geophysical Research Solid Earth↗

USGS2012

Tectonic setting of the Yukon-Koyukuk basin and its borderlands, western Alaska

The Yukon-Koyukuk basin of western Alaska is composed of an arcuate belt of Jurassic and Early Cretaceous subduction-related volcanic and plutonic rocks (Koyukuk terrane) flanked by deep subbasins filled with mid-Cretaceous terrigenous sedimentary rocks. The basin is bordered on three sides by metamorphosed Proterozoic and Paleozoic continental rocks (Seward, Arctic Alaska, and Ruby terranes) and is separated from the metamorphic borderlands by a narrow, highly tectonized belt of oceanic crust and mantle rocks (composite Angayucham-Tozitna terrane). The oceanic and mantle rocks, which dip inward beneath the basin and are thrust outward onto the borderlands, are divided into three separate thrust panels: (1) a structurally lowest panel (Slate Creek) composed of phyllite and metagraywacke of probable Devonian age, (2) a middle panel (Narvak) composed of imbricated basalt, chert, and gabbro of Devonian to Early Jurassic age, and (3) a structurally highest panel (Kanuti) composed of gabbro and peridotite of probable Middle and Late Jurassic age. The three thrust panels appear to represent a reversely stacked sequence that progresses from continental slope deposits in the lower panel to cumulus and mantle peridotites in the upper. Metamorphic mineral K-Ar ages from garnet amphibolite on the sole of the upper panel suggest that the upper panel was emplaced on the middle panel in the Middle to Late Jurassic. Subsequent accretion of all three panels to the continental rocks of the borderlands occurred in the latest Jurassic and Early Cretaceous, synchronous with arc volcanism within the basin. Arc volcanism waned and the accretionary phase ended in the middle of Early Cretaceous time. Uplift and erosion of the metamorphic borderlands and the obducted oceanic rocks began in late Early Cretaceous and was accompanied by the rapid filling of two flanking subbasins with turbiditic sediments. In the latest Early Cretaceous and early Late Cretaceous, shallow marine and nonmarine conglomerates were deposited around the margins of the basin, and a prograding delta was built out from the southeast margin of the basin across the turbiditic subbasins and the remnant volcanic arc. In the Late Cretaceous, western Alaska was subjected to strong east-west compression which severely deformed both the Yukon-Koyukuk basin and the borderlands. Several widespread magmatic episodes in the mid- and Late Cretaceous and in early Tertiary transgress the basin boundaries and stitch together the accreted arc and oceanic terranes and the continental borderlands.

Journal of Geophysical Research Solid Earth↗

USGS2012

Geologic implications of topographic, gravity, and aeromagnetic data in the northern Yukon-Koyukuk province and its borderlands, Alaska

The northern Yukon-Koyukuk province is characterized by low elevation and high Bouguer gravity and aeromagnetic anomalies in contrast to the adjacent Brooks Range and Ruby geanticline. Using newly compiled digital topographic, gravity, and aeromagnetic maps, I have divided the province into three geophysical domains. The Koyukuk domain, which is nearly equivalent to the Koyukuk lithotectonic terrane, is a horseshoe-shaped area, open to the south, of low topography, high gravity, and high-amplitude magnetic anomalies caused by an intraoceanic magmatic arc. The Angayucham and Kanuti domains are geophysical subdivisions of the Angayucham lithotectonic terrane that occur along the northern and southeastern margins of the Yukon-Koyukuk province, where oceanic rocks have been thrust over continental rocks of the Brooks Range and Ruby geanticline. Basalt of the Angayucham domain causes strong gravity highs and weak magnetic highs. The Kanuti domain is distinguished from the Angayucham domain by intense magnetic highs caused by cumulus mafic and ultramafic plutonic rocks, abundant ultramafic mantle tectonites, and magnetic syenite and monzonite. Long-wavelength, low-intensity magnetic highs and undulating gravity anomalies indicate an undulating basement surface of varied lithology beneath the Kobuk-Koyukuk and Lower Yukon basins. Modeling of gravity and magnetic anomalies shows that oceanic rocks of the Angayucham and Kanuti domains dip inward beneath the Kobuk-Koyukuk basin. The modeling supports, but does not prove, the hypothesis that the crust of the Kobuk-Koyukuk basin is 32–35 km thick, consisting of a tectonically thickened section of Cretaceous volcanic and sedimentary rocks and older oceanic crust. Plutons of the Brooks Range and the southern Ruby geanticline are nonmagnetic, ilmenite series, S-type granites that cause magnetic lows. Plutons of the northern Ruby geanticline are variable in their magnetic properties and cause both highs and lows. Plutons of both the eastern and western Yukon-Koyukuk province are variable in their magnetic expression but commonly cause magnetic lows in contrast to andesite.

Journal of Geophysical Research Solid Earth↗

USGS2012

Igneous history of the Koyukuk terrane, western Alaska: Constraints on the origin, evolution, and ultimate collision of an accreted island arc terrane

The Koyukuk terrane of western Alaska consists of volcanic, volcaniclastic, and plutonic rocks which range from Late Paleozoic to Early Cretaceous in age. The terrane crops out in a U-shaped belt which is roughly paralleled by outer belts of ultramafic rocks, oceanic plate basalts and cherts, and retrograded blueschist facies rocks of continental protolith. These rocks have been interpreted as components of a volcanic arc terrane that collided with the North American continental margin in Early Cretaceous time. The Koyukuk terrane consists of four time-stratigraphic units: (1) pre-Middle Jurassic basalts, (2) Middle and Late Jurassic granitic rocks, (3) lower Lower Cretaceous volcanic rocks, and (4) upper Lower Cretaceous volcanic rocks. Limited chemical data from the basalts of unit 1 indicate that they were erupted in a nonarc tectonic environment, possibly in an oceanic island or back arc setting. Units 2, 3, and 4 have the characteristics of subduction-related volcanic rocks (i.e., depleted Nb and Ta and enriched alkaline elements, relative to the light rare earth elements). Unit 3 contains tholeiitic, calc-alkaline, and alkaline rocks with chondrite-normalized rare earth element patterns that range from flat (La N /Yb N = 1) to highly light rare earth element enriched (La N /Yb N > 15). The highly alkaline or shoshonitic lavas were erupted toward the end of unit 3 time (Valanginian) during the final stages of arc-continent collision. These alkaline lavas could have been derived by very small degrees of partial melting of a similar source to that of the earlier arc lavas. Unit 4 lavas are also alkaline or shoshonitic, but their incompatible element composition indicates that they were derived from a different source than that of the earlier arc lavas. These late alkaline lavas are chemically similar to crosscutting mid-Cretaceous plutons whose isotopic compositions (Arth et al., this issue ( a )) suggest derivation by partial melting of distinctly older subcontinental lithosphere. We speculate that the parental magmas of unit 4 lavas may also have been derived by partial melting of this subcontinental mantle which was underthrust beneath the Koyukuk arc terrane during the final stage of arc-continent collision.

Journal of Geophysical Research Solid Earth↗

USGS2012

The Kanuti ophiolite, Alaska

The Kanuti ophiolite is a mafic-ultramafic thrust sheet of probable Jurassic age, formerly considered to be the upper part of the Yukon-Koyukuk ophiolite belt (Angayucham terrane). It is here called the Kanuti ophiolite after the Kanuti River region on the southeastern flank of the Yukon-Koyukuk Basin. The thrust sheet crops out discontinuously for a distance of more than 900 km along the northern and southeastern margins of the basin. It is probably correlative with similar ophiolite thrust sheets to the north in the western Brooks Range and to the south in the Ruby geanticline. Technically, the ophiolite is considered to be the Kanuti thrust panel of the Angayucham-Totzitna terrane. The Kanuti consistently overlies another extensive thrust sheet, consisting mostly of pillow basalt and radiolarian chert of Devonian to Jurassic age (Narvak thrust panel). This sheet is thrust over a third sheet consisting of probable Devonian phyllite and metagraywacke, which is in turn thrust over older metamorphic rocks (Slate Creek thrust panel). The Kanuti ophiolite is a partial ophiolite that consists of a lower residual mantle suite and an upper magmatic suite, but dikes, extrusives, and sediments are absent. The residual mantle suite is composed of harzburgite and dunite with refractory mineral compositions. The harzburgite is attributed to partial melting and extraction of basaltic magma; residual dunite is attributed to partial melting or to reaction of orthopyroxene out of harzburgite in contact with ascending melt diapirs. The magmatic suite consists of layered ultramafic and gabbroic rocks, containing minerals having limited iron enrichment. The absence of large volumes of magmatic rocks intermediate in composition between cumulus ultramafics and evolved gabbros favors periodic introduction of magma, rather than closed system fractional crystallization. The ultramafic rocks of both the residual mantle and magmatic suites are tectonites, which have undergone high-temperature deformation involving isoclinal folding on all scales and related syntectonic recrystallization. The olivine fabric is consistent with the glide system {0kl} [100], which has been produced experimentally at 800°–1190°C at 20 kbar. Olivine Z axes and subparallel isoclinal fold axes have consistent, northeast trends throughout the Kanuti region (>100 km NE-SW) and may be close to the original upper mantle flow direction, despite later low-angle thrust faulting. The order of crystallization in the cumulus ultramafic rocks of the magmatic suite is olivine, clinopyroxene, plagioclase, and orthopyroxene. The high Mg numbers of clinopyroxene (0.85–0.93) coexisting with olivine suggest that the cumulus ultramafic rocks crystallized at relatively high pressures (>10 kbar). The effects of parental magma composition cannot be evaluated, but the small difference in Mg numbers of coexisting olivine and clinopyroxene in the cumulus ultramafic rocks and in residual harzburgite suggests that regardless of absolute pressures, the pressure difference between the melting that produced the basalt magma and the initial fractional crystallization of the magma is small. Because of the limited range in rock types in the ophiolite, the tectonic environment cannot be interpreted unambiguously. However, the structural and petrological data are best reconciled with an origin in a volcanic arc tectonic setting.

Journal of Geophysical Research Solid Earth↗

USGS2012

Pillow basalts of the Angayucham terrane: Oceanic plateau and island crust accreted to the Brooks Range

The Angayucham Mountains (north margin of the Yukon-Koyukuk province) are made up of an imbricate stack of four to eight east-west trending, steeply dipping, fault slabs composed of Paleozoic (Devonian to Mississippean), Middle to Late Triassic, and Early Jurassic oceanic upper crustal rocks (pillow basalt, subordinate diabase, basaltic tuff, and radiolarian chert). Field relations and geochemical characteristics of the basaltic rocks suggest that the fault slabs were derived from an oceanic plateau or island setting and were emplaced onto the Brooks Range continental margin. The basalts are variably metamorphosed to prehnite-pumpellyite and low-greenschist facies. Major element analyses suggest that many are hypersthene-normative olivine tholeiites. Classification based on immobile trace elements confirms the tholeiitic character of most of the basalts but suggests that some had primary compositions transitional to alkali basalt. Although field and petrographic features of the basalts are similar, trace element characteristics allow definition of geographically distinct suites. A central outcrop belt along the crest of the mountains is made up of basalt with relatively flat rare earth element (REE) patterns. This belt is flanked to the north and south by LREE (light rare earth element)-enriched basalts. Radiolarian and conodont ages from interpillow and interlayered chert and limestone indicate that the central belt of basalts is Triassic in age, the southern belt is Jurassic in age, and the northern belt contains a mixture of Paleozoic and Mesozoic ages. Data for most of the basalts cluster in the “within-plate basalt” fields of trace element discriminant diagrams; none have trace-element characteristics of island arc basalt. The Triassic and Jurassic basalts are geochemically most akin to modern oceanic plateau and island basalts. Field evidence also favors an oceanic plateau or island setting. The great composite thickness of pillow basalt probably resulted from obduction faulting, but the lack of fault slabs of gabbro or peridotite suggests that obduction faults did not penetrate below oceanic layer 2, a likely occurrence if layer 2 were anomalously thick, as in the vicinity of an oceanic island. The presence of basaltic tuff interbeds indicates proximity to an explosive basaltic eruptive center. The juxtaposition of submarine basalts of differing chemical affinity and age, adjacent to higher-grade Paleozoic metamorphic rocks of the Brooks Range to the north, may be explained by obduction of internally complex (thickened) oceanic crust formed in an ocean plateau setting. Emplacement and rotation of thrust plates to steep attitudes occurred during accretion of the Brooks Range passive margin, probably beginning in the Late to Middle Jurassic.

Journal of Geophysical Research Solid Earth↗

USGS2012

Isotopic and trace element variations in the Ruby Batholith, Alaska, and the nature of the deep crust beneath the Ruby and Angayucham Terranes

Thirty-six samples from plutons of the Ruby batholith of central Alaska were collected and analyzed for 22 trace elements, and many were analyzed for the isotopic compositions of Sr, Nd, O, and Pb in order to delimit the processes that produced the diversity of granodioritic to granitic compositions, to deduce the nature of the source of magmas at about 110 Ma, and to characterize the deep crust beneath the Ruby and Angayucham terranes. Plutons of the batholith show a substantial range in initial 87 Sr/ 86 Sr (SIR) of 0.7055–0.7235 and a general decrease from southwest to northeast. Initial 143 Nd/ 144 Nd (NIR) have a range of 0.51150–0.51232 and generally increase from southwest to northeast. The δ 18 O values for most whole rocks have a range of +8.4 to +11.8 and an average of +10.3&permil;. Rb, Cs, U, and Th show large ranges of concentration, generally increase as SiO 2 increases, and are higher in southwest than in northeast plutons. Sr, Ba, Zr, Hf, Ta, Sc, Cr, Co, and Zr show large ranges of concentration and generally decrease as SiO 2 increases. Rare earth elements (REE) show fractionated patterns and negative Eu anomalies. REE concentrations and anomalies are larger in the southwest than in the northeast plutons. Uniformity of SIR and NIR in Sithylemenkat and Jim River plutons suggests a strong role for fractional crystallization or melting of uniform magma sources at depth. Isotopic variability in Melozitna, Ray Mountains, Hot Springs, and Kanuti plutons suggests complex magmatic processes such as magma mixing and assimilation, probably combined with fractional crystallization, or melting of a complex source at depth. The large variations in SIR and NIR in the batholith require a variation in source materials at depth. The southwestern plutons probably had dominantly siliceous sources composed of metamorphosed Proterozoic and Paleozoic upper crustal rocks. The northeastern plutons probably had Paleozoic sources that were mixtures of siliceous and intermediate to mafic crustal rocks. The inferred sources could well have been the higher-metamorphic-grade lithologic equivalents of the exposed Proterozoic(?) to Paleozoic schists, orthogneisses, and metavolcanic rocks of Ruby terrane, the silicic portions of which are quite radiogenic. The deeper crustal sources that gave rise to most of the batholithic magmas are inferred to be similar under both the Ruby metamorphic terrane and the Angayucham ophiolitic terrane.

Alaska↗