USGS Science⌕ Search

SEARCH · USGS Science

Results for “Arctic”

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

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

At least 1,315 records · Page 73Linked to original sources
USGS2021

Permafrost promotes shallow groundwater flow and warmer headwater streams

The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from eleven headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7 % of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5 % moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

Alaska↗
USGS2022

Watching the Cryosphere thaw: Seismic monitoring of permafrost degradation using distributed acoustic sensing during a controlled heating experiment

Permafrost degradation is rapidly increasing in response to a warming Arctic climate, altering landscapes and damaging critical infrastructure. Solutions for monitoring permafrost thaw dynamics are essential to understand biogeochemical feedbacks as well as to issue warnings for hazardous geotechnical conditions. We investigate the feasibility of permafrost monitoring using permanently installed fiber-optic seismic networks. We conducted a 2-month seismic monitoring campaign during a controlled thaw experiment using a permanent surface orbital vibrator (SOV) and a 2D-array of distributed acoustic sensing (DAS) cables, and observed significant (15%) shear-wave velocity ( V s ) reductions and approximately 2 m depression of the permafrost table beneath the heating zone. These observations were validated by time-lapse horizontal-to-vertical spectral ratio (HVSR) analysis from three co-located broadband seismometers. The combination of SOV and DAS provided unique seismic observations for permafrost monitoring at the field scale, as well as a basis for design and development of early warning systems for permafrost thaw.

Geophysical Research Letters↗
USGS2021

Retreat and regrowth of the Greenland Ice Sheet during the Last Interglacial as simulated by the CESM2-CISM2 coupled climate–ice sheet model

During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model-derived and geological reconstruction estimates of the sea-level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre-industrial conditions shows a much smaller ice-sheet retreat, highlighting the importance of the changes in high-latitude vegetation distribution for amplifying the ice-sheet response.

Greenland Ice Sheet↗
USGS2023

Long-term changes in concentrations and yield of riverine dissolved silicon from the poles to the tropics

Riverine exports of silicon (Si) influence global carbon cycling through the growth of marine diatoms, which account for ∼25% of global primary production. Climate change will likely alter river Si exports in biome-specific ways due to interacting shifts in chemical weathering rates, hydrologic connectivity, and metabolic processes in aquatic and terrestrial systems. Nonetheless, factors driving long-term changes in Si exports remain unexplored at local, regional, and global scales. We evaluated how concentrations and yields of dissolved Si (DSi) changed over the last several decades of rapid climate warming using long-term data sets from 60 rivers and streams spanning the globe (e.g., Antarctic, tropical, temperate, boreal, alpine, Arctic systems). We show that widespread changes in river DSi concentration and yield have occurred, with the most substantial shifts occurring in alpine and polar regions. The magnitude and direction of trends varied within and among biomes, were most strongly associated with differences in land cover, and were often independent of changes in river discharge. These findings indicate that there are likely diverse mechanisms driving change in river Si biogeochemistry that span the land-water interface, which may include glacial melt, changes in terrestrial vegetation, and river productivity. Finally, trends were often stronger in months outside of the growing season, particularly in temperate and boreal systems, demonstrating a potentially important role of shifting seasonality for the flux of Si from rivers. Our results have implications for the timing and magnitude of silica processing in rivers and its delivery to global oceans.

Global Biogeochemical Cycles↗
USGS2023

Characterizing performance of freshwater wetland methane models across time scales at FLUXNET-CH4 sites using wavelet analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

JGR Biogeosciences↗
USGS2023

Long-term, high-resolution permafrost monitoring reveals coupled energy balance and hydrogeologic controls on talik dynamics near Umiujaq (Nunavik, Québec, Canada)

Rising temperatures in the Arctic and subarctic are driving the rapid thaw of permafrost by reducing permafrost cooling, increasing active layer thickness, and promoting talik formation. In this study, the cyrohydrogeology of a permafrost mound located within the discontinuous permafrost zone near Umiujaq (Nunavik, Québec, Canada) is characterized through the analysis of a dataset covering more than two decades of monitoring. This dataset captures a high degree of interannual variability in air temperature and ground thermal conditions, as well as the formation and closure of a supra-permafrost talik. Data indicate that variable saturation and advective heat transport directly contribute to the expansion and contraction of the talik. Data further indicate the presence of two distinct thermo-hydrologic settings resulting from differences in surface conditions, as well as subsurface thermal and flow regimes. The first, found at the top of the mound feature, is characterized by very low moisture contents (< 0.05 m 3 /m 3 ), while the second, found at the side of the mound feature, shows higher annual moisture contents that strongly influence the dynamics of heat and groundwater flow. The data were synthesized into a detailed conceptual model of the cyrohydrogeological dynamics that highlights the important role of hydrogeological characterization and long-term datasets in understanding the effects of groundwater flow on seasonal frost and permafrost dynamics. Specifically, the results presented here show that in the absence of long-term datasets, longer-period transient phenomena such as talik opening and closure may be misrepresented as uni-directional feedback loops, as opposed to highly-dynamic temporary phenomena.

Nunavik, Quebec↗
USGS2024

Controls on stable methane isotope signatures in northern peatlands and potential shifts in signatures under permafrost thaw scenarios

Northern peatlands are a globally significant source of methane (CH 4 ), and emissions are projected to increase due to warming and permafrost loss. Understanding the microbial mechanisms behind patterns in CH 4 production in these systems will be key to predicting annual emissions changes, with stable carbon isotopes (δ 13 C-CH 4 ) being a powerful tool for characterizing these drivers. Given that δ 13 C signatures of CH 4 are used in top-down atmospheric inversion models to partition sources, our ability to model CH 4 production pathways and associated δ 13 C-CH 4 signatures in peatland types impacted by a changing climate is critical. We sought to characterize the role of environmental conditions, including both hydrologic and vegetation patterns associated with permafrost thaw, on δ 13 C-CH 4 signatures from a diverse set of high-latitude peatlands. We measured porewater and emitted CH 4 stable isotopes, pH, and vegetation composition from five boreal-Arctic peatlands. Porewater δ 13 C-CH 4 was strongly associated with peatland type, with δ 13 C enriched values obtained from more minerotrophic fens (-61.2 ± 9.1‰) compared to permafrost-free bogs (-74.1 ± 9.4‰) and raised permafrost bogs (-81.6 ± 11.5‰). Variation in porewater δ 13 C-CH 4 was best explained by sedge cover, CH 4 concentration, and the interactive effect of peatland type and pH ( r 2 = 0.50, p < 0.001). Emitted δ 13 C-CH 4 varied greatly but was positively correlated with porewater δ 13 C-CH 4 , suggesting that porewater data can be used to predict changing emissions signatures from these systems. We calculated a weighted mean mixed atmospheric CH 4 signature for northern peatlands of -65.3 ± 7‰ and show that this signature is more sensitive to landscape drying (4 to 10 % depletion in δ 13 C) than wetting (1.5 to 5% enrichment in δ 13 C) under permafrost thaw scenarios. Our results suggest northern peatland δ 13 C-CH 4 signatures are likely to shift in the future which has important implications for source partitioning in atmospheric inversion models.

Journal of Geophysical Research: Biogeosciences↗
USGS2025

Organic matter composition versus microbial source: Controls on carbon loss from fen wetland and permafrost soils

Wetland and permafrost soils contain some of Earth's largest reservoirs of organic carbon, and these stores are threatened by rapid warming across the Arctic. Nearly half of northern wetlands are affected by permafrost. As these ecosystems warm, the cycling of dissolved organic matter (DOM) and the opportunities for microbial degradation are changing. This is particularly evident as the relationship between wetland and permafrost DOM dynamics evolves, especially with the introduction of permafrost-derived DOM into wetland environments. Thus, understanding the interplay of DOM composition and microbial communities from wetlands and permafrost is critical to predicting the impact of released carbon on global carbon cycling. As little is understood about the interactions between wetland active layer and permafrost-derived sources as they intermingle, we conducted experimental bioincubations of mixtures of DOM and microbial communities from two fen wetland depths (shallow: 0–15 cm, and deep: 15–30 cm) and two ages of permafrost soil (Holocene and Pleistocene). We found that the source of microbial inoculum was not a significant driver of dissolved organic carbon (DOC) degradation across treatments; rather, DOM source and specifically, DOM molecular composition, controlled the rate of DOC loss over 100 days of bioincubations. DOC loss across all treatments was negatively correlated with modified aromaticity index, O/C, and the relative abundance of condensed aromatic and polyphenolic formula, and positively correlated with H/C and the relative abundance of aliphatic and peptide-like formula. Pleistocene permafrost-derived DOC exhibited ∼70% loss during the bioincubation driven by its initial molecular-level composition, highlighting its high bioavailability irrespective of microbial source.

Alaska↗
USGS2025

Toward co-designed Earth System Models: Reflecting end-user priorities in local applications from a modeler's perspective

Earth System Models (ESM) are crucial for quantifying climate impacts across Earth's interconnected systems and supporting science-based adaptation and mitigation. However, not including end-users, especially decision-makers representing communities vulnerable to climate change, can limit model utility, increase epistemic risks, and lead to information misuse in decision-making. While the ESM community increasingly values broad community engagement, end-users may not initially perceive models as useful for local planning. Co-designing models with end-users fosters two-way learning: users better understand models and their outputs, while modelers gain insights into fine-scale local processes like monitoring practices and management priorities. Higher-level co-design can lead to more customized, priority-driven, and useful modeling products. Despite these benefits, modelers often struggle to initiate meaningful partnerships with local communities. Therefore, this paper explores model co-design from the perspective of modelers. This study presents two case studies where modelers and social scientists collaborated with Indigenous communities' decision-makers to reflect their priorities in model design and application. In the Arctic Rivers Project, high-resolution climate and hydrology data sets for Alaska were developed with guidance from an Indigenous Advisory Council, using optimized, coupled land-atmosphere models. In the Mid-Klamath Project, we partnered with the Karuk Tribe's Department of Natural Resources to assess climate change and prescribed burning impacts on terrestrial hydrology in the Klamath River Basin. Drawing from these studies, we introduce a four-level framework: (a) Co-design Configuration; (b) Model Tuning; (c) Incorporate Contextual Knowledge; (d) Co-develop New Model Functions. We aim to help researchers consider and compare co-design across diverse modeling projects systematically and coherently.

Alaska, California, Oregon↗
USGS2025

Interacting sea-level rise, sea-ice loss, storm flooding, erosion, and permafrost thaw threaten ecosystems, wildlife, and communities on the Yukon-Kuskokwim Delta

The Yukon-Kuskokwim Delta has the largest intertidal wetland in North America, is a globally critical breeding area for waterbirds, and is home to the largest regional indigenous population in the Arctic. Here, coastal tundra ecosystems, wildlife, and indigenous communities are highly vulnerable to sea-ice loss in the Bering Sea, sea-level rise, storm flooding, erosion, and collapsing ground from permafrost thaw caused by climate warming. These drivers interact in non-linear ways to increase flooding, salinization, and sedimentation, and thus, alter ecosystem trajectories and broader landscape evolution. Rapid changes in these factors over decadal time scales are highly likely to cause transformative shifts in coastal ecosystems across roughly 70% of the outer delta this century. We project saline and brackish ecotypes on the active delta floodplain with frequent sedimentation will maintain dynamic equilibrium with sea-level rise and flooding, slightly brackish ecotypes on the inactive floodplain with infrequent flooding and low sedimentation rates will be vulnerable to increased flooding and likely transition to more saline and brackish ecotypes, and fresh lacustrine and lowland ecotypes on the abandoned floodplain with permafrost plateaus will be vulnerable to thermokarst, salinization and flooding that will shift them toward brackish ecosystems. This will greatly affect bird nesting and foraging habitats, with both winners and losers. Already, some Yup'ik communities are facing relocation of their low-lying villages. The societal challenges and consequences of adapting to these changing landscapes are enormous and will require a huge societal effort.

Alaska↗
USGS2026

A sea ice entrapment event in the southern Chukchi Sea: Analysis and prediction

Amplified Arctic warming is reducing sea ice cover, which is driving an increase in geopolitical interest in the region as it offers the possibility of reducing shipping times between Asia, Europe, and eastern North America, at the risk of increased ice hazards. Here, we examine the case of the Norseman II research ship that was trapped by sea ice in the southern Chukchi Sea for 14 days in June 2024. This is the first study of its sort in this region. We show that anomalously thick and extensive sea ice was present north of the region prior to the event, and that strong northerly winds in early June advected this ice southward, trapping the ship. Later in June, southerly winds advected ice northwards away from the ship, helping to free it. We further show that the event was forecastable.

Alaska↗
USGS2012

Is the extent of glaciation limited by marine gas-hydrates?

Methane may have been released to the atmosphere during the Quaternary from Arctic shelf gas-hydrates as a result of thermal decomposition caused by climatic warming and rising sea-level; this release of methane (a greenhouse gas) may represent a positive feedback on global warming [Revelle, 1983; Kvenvolden, 1988a; Nisbet, 1990]. We consider the response to sea-level changes by the immense amount of gas-hydrate that exists in continental rise sediments, and suggest that the reverse situation may apply&mdash;that release of methane trapped in the deep-sea sediments as gas-hydrates may provide a negative feedback to advancing glaciation. Methane is likely to be released from deep-sea gas-hydrates as sea-level falls because methane gas-hydrates decompose with pressure decrease. Methane would be released to sediment pore space at shallow sub-bottom depths (100's of meters beneath the seafloor, commonly at water depths of 500 to 4,000 m) producing zones of markedly decreased sediment strength, leading to slumping [Carpenter, 1981; Kayen, 1988] and abrupt release of the gas. Methane is likely to be released to the atmosphere in spikes that become larger and more frequent as glaciation progresses. Because addition of methane to the atmosphere warms the planet, this process provides a negative feedback to glaciation, and could trigger deglaciation.

Geophysical Research Letters↗
USGS2012

The Beaufort Sea continental shelf as a seasonal source of atmospheric methane

Methane concentrations in the Beaufort Sea under the winter ice canopy offshore from northern Alaska are 3 to 28 times greater than they are in late summer when the ice is absent in a similar region offshore from northern Canada where methane is in approximate equilibrium with the atmosphere. These observations suggest that methane concentrates in the water under the sea‐ice cover during winter and ventilates rapidly in late summer as the ice melts and retreats. Conditions similar to those on the Beaufort Sea shelf likely exist on the much larger Siberian shelf, making the Arctic Ocean margin a possible seasonal, high‐latitude, marine source of about 0.1 Tg yr −1 atmospheric methane. The small addition of methane likely contributes to the late‐summer increase in atmospheric methane that is observed each year particularly in the northern hemisphere.

Alaska↗
USGS1998

Greenland Sea Odden sea ice feature: Intra-annual and interannual variability

The “Odden” is a large sea ice feature that forms in the east Greenland Sea that may protrude eastward to 5°E from the main sea ice pack (at about 8°W) between 73° and 77°N. It generally forms at the beginning of the winter season and can cover 300,000 km 2 . Throughout the winter the outer edge of the Odden may advance and retreat by several hundred kilometers on timescales of a few days to weeks. Satellite passive microwave observations from 1978 through 1995 provide a continuous record of the spatial and temporal variations of this extremely dynamic phenomenon. Aircraft synthetic aperture radar, satellite passive microwave, and ship observations in the Odden show that the Odden consists of new ice types, rather than older ice types advected eastward from the main pack. The 17-year record shows both strong interannual and intra-annual variations in Odden extent and temporal behavior. For example, in 1983 the Odden was weak, in 1984 the Odden did not occur, and in 1985 the Odden returned late in the season. An analysis of the ice area and extent time series derived from the satellite passive microwave observations along with meteorological data from the International Arctic Buoy Program (IABP) determined the meteorological forcing associated with Odden growth, maintenance, and decay. The key meteorological parameters that are related to the rapid ice formation and decay associated with the Odden are, in order of importance, air temperature, wind speed, and wind direction. Oceanographic parameters must play an important role in controlling Odden formation, but it is not yet possible to quantify this role because of a lack of long-term oceanographic observations.

Journal of Geophysical Research C: Oceans↗
USGS2012

Early Cretaceous paleolatitude of the Yukon-Koyukuk province, Alaska

We report results from a paleomagnetic study of the Yukon-Koyukuk province, a key region for reconstructing the Cretaceous paleogeography of northern Alaska and northeastern Asia. The province lies between the displaced continental fragment of Arctic Alaska and the accreted terranes of southern Alaska. Although Lower Cretaceous volcanogenic rocks of the Yukon-Koyukuk province preserve primary remanent magnetizations, the primary component is masked by a pervasive overprint. The over-printing magnetization, which has dual polarity, was acquired after a Late Cretaceous orogeny left the region highly deformed. By fitting planes to tilt-corrected thermal demagnetization data and solving for the best intersection of the corresponding remagnetization circles, we separated the primary component, presumably a vestige of the original thermoremanent magnetization, from the overprinting component of the natural remanence. The primary magnetization corresponds to a mean paleopole of Early Cretaceous age at 60.2°N, 159°E. The mean paleolatitude is 69°N with the 95% confidence ellipse ranging from 60° to 75°N. This analysis was successful because the tilt corrections of the sampled strata had sufficient variation to ensure a strong intersection of the remagnetization circles. We infer that Lower Cretaceous rocks of the Brooks Range and the Yukon-Koyukuk province were once separated by a wide ocean basin that was closed by subduction, probably during the Cretaceous. The northern limb of the Angayucham terrane, which lies between the Brooks Range and the Yukon-Koyukuk province, may contain remnants of the now-closed ocean basin.

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

The nature of the crust in the Yukon-Koyukuk province as inferred from the chemical and isotopic composition of five Late Cretaceous to Early Tertiary volcanic fields in western Alaska

Late Cretaceous and early Tertiary volcanic and plutonic rocks in western Alaska comprise a vast magmatic province extending from the Alaska Range north to the Arctic Circle, south to Bristol Bay, and west to the Bering Sea Shelf. The chemical and isotopic composition of five of these Late Cretaceous to early Tertiary volcanic fields in the north central part of this province were studied to determine if Paleozoic or older continental crust underlies the Yukon-Koyukuk province. Three of the fields, the Blackburn Hills, Yukon River, and Kanuti, occur within the Yukon-Koyukuk province and two, the Sischu and Nowitna, overlie bordering Precambrian and Paleozoic metamorphic terranes to the southeast. High initial 87 Sr/ 86 Sr of 0.7075–0.7079 and moderate initial 143 Nd/ 144 Nd of 0.51244–0.51247 of rhyolite, dacite, and high-silica andesite of the Sischu volcanic field indicate that the magmas have interacted with the underlying Paleozoic or older continental crust. The relatively limited variation of isotopic (initial 87 Sr/ 86 Sr = 0.7044–0.7051; initial 143 Nd/ 144 Nd = 0.51256–0.51257) and elemental compositions of andesites from the Nowitna field can be accounted for by assimilation of small amounts of Paleozoic or older continental crust during crystal fractionation of andesite parent magmas at crustal levels. The Blackburn Hills field, which consists of medium-K basalt, andesite, and rhyolite intruded by a small granitic pluton, has a large range in initial 87 Sr/ 86 Sr and initial 143 Nd/ 144 Nd that plot in the field for 60 Ma mantle, from near mid-ocean ridge basalts to near “bulk-earth” compositions (initial 87 Sr/ 86 Sr = 0.7033–0.7052; initial 143 Nd/ 144 Nd = 0.51253–0.51290). Andesites and basalts from the Blackburn Hills are divided into two group on the basis of rare earth element (REE) and isotopic composition. Isotopic variation in the more primitive group 1 is best explained by assimilation of the lower crust of the Jurassic to Early Cretaceous Koyukuk terrane by mantle-derived basalts during crystal fractionation, though part of the isotopic variation may be due to metasomatism of an oceanic island basalt type mantle source by fluids derived from subducted sediments. Group 2 andesites from the Blackburn Hills have lower heavy REE abundances and more enriched isotopic compositions. These group 2 andesites and dacites from the Kanuti field, which have ( 87 Sr/ 86 Sr) i = 0.7043–0.7048 and ( 143 Nd/ 144 Nd) i = 0.51248–0.51267, appear to have formed by partial melting of the lower crust of the Koyukuk terrane. The Yukon River field consists of basalt, andesite, dacite, and rhyolite having ( 87 Sr/ 86 Sr) i = 0.7037–0.7051 and ( 143 Nd/ 144 Nd) i = 0.51266–0.51280; its isotopic composition does not require the presence of Paleozoic or older continental crust under the volcanic field and may have formed by interaction between mantle-derived melts and the oceanic Angayucham/Tozitna or island arc Koyukuk terrane. Most of the intrusive rocks and rhyolite domes from the Blackburn Hills volcanic field have ( 87 Sr/ 86 Sr) i = 0.7038–0.7041 and dacites from the Kanuti volcanic field have ( 87 Sr/ 86 Sr) i = 0.7043–0.7048. Thus little or no old continental crust was involved in the genesis of the Late Cretaceous and early Tertiary rocks and therefore probably does not extend beneath this part of the Yukon-Koyukuk province. However, the ultimate source of the small volumes of enriched shoshonitic andesite ( 87 Sr/ 86 Sr = 0.7075, 143 Nd/ 144 Nd = 0.5125) erupted at 118 Ma in the Yukon-Koyukuk province may be continental lithosphere, which may have been thrust under this part of the Yukon-Koyukuk province during arc-continent collision in the Early Cretaceous.

Journal of Geophysical Research Solid Earth↗