Geology Reports⌕ Search

Geology topics

Research about Alaska

Source-linked reports with geographic coverage including Alaska.

At least 217 records · Page 12Linked to original sources

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗

Metal mobilization from thawing permafrost to aquatic ecosystems is driving rusting of Arctic streams

Climate change in the Arctic is altering watershed hydrologic processes and biogeochemistry. Here, we present an emergent threat to Arctic watersheds based on observations from 75 streams in Alaska’s Brooks Range that recently turned orange, reflecting increased loading of iron and toxic metals. Using remote sensing, we constrain the timing of stream discoloration to the last 10 years, a period of rapid warming and snowfall, suggesting impairment is likely due to permafrost thaw. Thawing permafrost can foster chemical weathering of minerals, microbial reduction of soil iron, and groundwater transport of metals to streams. Compared to clear reference streams, orange streams have lower pH, higher turbidity, and higher sulfate, iron, and trace metal concentrations, supporting sulfide mineral weathering as a primary mobilization process. Stream discoloration was associated with dramatic declines in macroinvertebrate diversity and fish abundance. These findings have considerable implications for drinking water supplies and subsistence fisheries in rural Alaska.

Alaska↗

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska↗

Snow avalanches are a primary climate-linked driver of mountain ungulate populations

Snow is a major, climate-sensitive feature of the Earth’s surface and catalyst of fundamentally important ecosystem processes. Understanding how snow influences sentinel species in rapidly changing mountain ecosystems is particularly critical. Whereas effects of snow on food availability, energy expenditure, and predation are well documented, we report how avalanches exert major impacts on an ecologically significant mountain ungulate - the coastal Alaskan mountain goat ( Oreamnos americanus ). Using long-term GPS data and field observations across four populations (421 individuals over 17 years), we show that avalanches caused 23−65% of all mortality, depending on area. Deaths varied seasonally and were directly linked to spatial movement patterns and avalanche terrain use. Population-level avalanche mortality, 61% of which comprised reproductively important prime-aged individuals, averaged 8% annually and exceeded 22% when avalanche conditions were severe. Our findings reveal a widespread but previously undescribed pathway by which snow can elicit major population-level impacts and shape demographic characteristics of slow-growing populations of mountain-adapted animals.

Alaska↗

Automatic identification and quantification of volcanic hotspots in Alaska using HotLINK: The hotspot learning and identification network

An increase in volcanic thermal emissions can indicate subsurface and surface processes that precede, or coincide with, volcanic eruptions. Space-borne infrared sensors can detect hotspots—defined here as localized volcanic thermal emissions—in near-real-time. However, automatic hotspot detection systems are needed to efficiently analyze the large quantities of data produced. While hotspots have been automatically detected for over 20 years with simple thresholding algorithms, new computer vision technologies, such as convolutional neural networks (CNNs), can enable improved detection capabilities. Here we introduce HotLINK: the Hotspot Learning and Identification Network, a CNN trained to detect hotspots with a dataset of −3,800 satellite-based, Visible Infrared Imaging Radiometer Suite (VIIRS) images from Mount Veniaminof and Mount Cleveland volcanoes, Alaska. We find that our model achieves an accuracy of 96% (F1-score 0.92) when evaluated on −1,700 unseen images from the same volcanoes, and 95% (F1-score 0.67) when evaluated on −3,000 images from six additional Alaska volcanoes (Augustine Volcano, Bogoslof Island, Okmok Caldera, Pavlof Volcano, Redoubt Volcano, Shishaldin Volcano). In comparison with an existing threshold-based hotspot detection algorithm, MIROVA (Coppola et al., Geological Society, London, Special Publications, 2016, 426, 181–205), our model detects 22% more hotspots and produces 12% fewer false positives. Additional testing on −700 labeled Moderate Resolution Imaging Spectroradiometer (MODIS) images from Mount Veniaminof demonstrates that our model is applicable to this sensor’s data as well, achieving an accuracy of 98% (F1-score 0.95). We apply HotLINK to 10 years of VIIRS data and 22 years of MODIS data for the eight aforementioned Alaska volcanoes and calculate the radiative power of detected hotspots. From these time series we find that HotLINK accurately characterizes background and eruptive periods, similar to MIROVA, but also detects more subtle warming signals, potentially related to volcanic unrest. We identify three advantages to our model over its predecessors: 1) the ability to detect more subtle volcanic hotspots and produce fewer false positives, especially in daytime images; 2) probabilistic predictions provide a measure of detection confidence; and 3) its transferability, i.e., the successful application to multiple sensors and multiple volcanoes without the need for threshold tuning, suggesting the potential for global application.

Alaska↗

2020 Volcanic activity in Alaska—Summary of events and response of the Alaska Volcano Observatory

The Alaska Volcano Observatory responded to eruptions, volcanic unrest or suspected unrest, increased seismicity, and other significant activity at nine volcanic centers in Alaska in 2020. The most notable volcanic activity in 2020 was an eruption of Shishaldin Volcano, which produced lava flows, lahars, and ash. Mount Cleveland had one small ash-producing eruption in June but was quiet thereafter. Other activity documented in 2020 consisted of elevated seismicity at the volcanoes Mount Veniaminof, Pavlof Volcano, Makushin Volcano, Atka volcanic complex (Korovin Volcano), Great Sitkin Volcano, and Semisopochnoi Island. Finally, the resuspension of ash deposited during the 1912 Novarupta-Katmai eruption was documented on three occasions.

Alaska↗

2021 Volcanic activity in Alaska and the Commonwealth of the Northern Mariana Islands—Summary of events and response of the Alaska Volcano Observatory

In 2021, the Alaska Volcano Observatory responded to eruptions, volcanic unrest or suspected unrest, increased seismicity, and other significant activity at 15 volcanic centers in Alaska and the Commonwealth of the Northern Mariana Islands. Eruptive activity in Alaska consisted of repeated small, ash-producing, phreatomagmatic explosions from Mount Young on Semisopochnoi Island; an explosion at Great Sitkin Volcano followed by the eruption of a thick lava flow that filled and overflowed the summit crater; weak explosive activity and the eruption of small, channelized flows at Pavlof Volcano; and a short-lived eruption at Mount Veniaminof that produced ash emissions from an intracaldera cone, as well as lava flows confined to a melt pit in the ice mantling the cone’s flank. Mount Cleveland had a period of unrest, but no eruptive activity took place there. Anomalous seismicity was also detected at Atka volcanic complex, Mount Gareloi, and Davidof volcano. New warm springs opened and deposited mud at the summit and north base of Shrub mud volcano. Other activity of note in Alaska consisted of large ice and rock avalanches at Iliamna Volcano and Mount Spurr, ash resuspension events at Mount Katmai and Aniakchak Crater, and anomalous deformation at Mount Okmok that was consistent with a shallow intrusion of magma. In the Commonwealth of the Northern Marianas Islands, a brief, ash-producing eruption occurred at Mount Pagan.

Alaska↗

A high-resolution, daily hindcast (1990-2021) of Alaskan river discharge and temperature from coupled and optimized physical models

Water quality and freshwater ecosystems are affected by river discharge and temperature. Models are frequently used to estimate river temperature on large spatial and temporal scales due to limited observations of discharge and temperature. In this study, we use physically based river routing and temperature models to simulate daily discharge and river temperature for rivers in 138 basins in Alaska, including the entire Yukon River basin, from 1990–2021. The river temperature model was optimized for ice free months using a surrogate-based model optimization method, improving model performance at uncalibrated river gages. A common statistical model relating local air and water temperature was used as a benchmark. The physically based river temperature model exhibited superior performance compared to the benchmark statistical model after optimization, suggesting river temperature model optimization could become more routine. The river temperature model demonstrated high sensitivity to air temperature and model parameterization, and lower sensitivity to discharge. Validation of the models showed a Kling-Gupta Efficiency of 0.46 for daily river discharge and a root mean square error of 2.04°C for daily river temperature, improving on the non-optimized physical model and the benchmark statistical model, which had root mean square errors of 3.24 and 2.97°C, respectively. The simulation shows that rivers in northern Alaska have higher maximum summer temperatures and more variability than rivers in the Central and Southern regions. Furthermore, this framework can be readily adapted for use across models and regions.

Alaska↗

Versatile modeling of deformation (VMOD) inversion framework: Application to 20 years of observations at Westdahl Volcano and Fisher Caldera, Alaska, US

We developed an open source, extensible Python-based framework, that we call the Versatile Modeling of Deformation (VMOD), for forward and inverse modeling of crustal deformation sources. VMOD abstracts from specific source model implementations, data types and inversion methods. We implement the most common geodetic source models which can be combined to model and analyze multi-source deformation. VMOD supports Global Navigation Satellite System (GNSS), InSAR, electronic distance measurement, Leveling and tilt data. To infer source characteristics from observations, VMOD implements non-linear least squares and Markov Chain Monte-Carlo Bayesian inversions, including joint inversions using different sources of data. VMOD's structure allows for easy integration of new geodetic models, data types, and inversion strategies. We benchmark the forward models against other published results and the inversion approaches against other implementations. We apply VMOD to analyze deformation at Unimak Island, Alaska, observed with continuous and campaign GNSS, and ascending and descending InSAR time series generated from Sentinel-1 satellite radar acquisitions. These data show an inflation pattern at Westdahl volcano and subsidence at Fisher Caldera. We use VMOD to test a range of source models by jointly inverting the GNSS and InSAR data sets. Our final model simultaneously constrains the parameters of two sources. Our results reveal a depressurizing spheroid under Fisher Caldera ∼4–6 km deep, contracting at a rate of ∼2–3 Mm 3 /yr, and a pressurizing spherical source underneath Westdahl volcano ∼6–8 km deep, inflating at ∼5 Mm 3 /yr. This and past applications of VMOD to volcanic unrest benefit from an extensible framework which supports jointly inversions of data sets for parameters of easily composable multi-source models.

Alaska↗

Current and projected flood exposure for Alaska coastal communities

Globally, coastal communities experience flood hazards that are projected to worsen from climate change and sea level rise. The 100-year floodplain or record flood are commonly used to identify risk areas for planning purposes. Remote communities often lack measured flood elevations and require innovative approaches to estimate flood elevations. This study employs observation-based methods to estimate the record flood elevation in Alaska communities and compares results to elevation models, infrastructure locations, and sea level rise projections. In 46 analyzed communities, 22% of structures are located within the record floodplain. With sea level rise projections, this estimate increases to 30–37% of structures by 2100 if structures remain in the same location. Flood exposure is highest in western Alaska. Sea level rise projections suggest northern Alaska will see similar flood exposure levels by 2100 as currently experienced in western Alaska. This evaluation of record flood height, category, and history can be incorporated into hazard planning documents, providing more context for coastal flood exposure than previously existed for Alaska. This basic flood exposure method is transferable to other areas with similar mapping challenges. Identifying current and projected hazardous zones is essential to avoid unintentional development in floodplains and improve long-term safety.

Alaska↗

Canopy cover and microtopography control precipitation-enhanced thaw of ecosystem-protected permafrost

Northern high-latitudes are projected to get warmer and wetter, which will affect rates of permafrost thaw and mechanisms by which thaw occurs. To better understand the impact of rain, as well as other factors such as snow depth, canopy cover, and microtopography, we instrumented a degrading permafrost plateau in south-central Alaska with high-resolution soil temperature sensors. The site contains ecosystem-protected permafrost, which persists in unfavorable climates due to favorable ecologic conditions. Our study (2020–2022) captured three of the snowiest years and three of the four wettest years since the site was first studied in 2015. Average thaw rates along an across-site transect increased nine-fold from 6 ± 5 cm yr −1 (2015–2020) to 56 ± 12 cm yr −1 (2020–2022). This thaw was not uniform. Hummock locations, residing on topographic high points with relatively dense canopy, experienced only 8 ± 9 cm yr −1 of thaw, on average. Hollows, topographic low points with low canopy cover, and transition locations, which had canopy cover and elevation between hummocks and hollows, thawed 44 ± 6 cm yr −1 and 39 ± 13 cm yr −1 , respectively. Mechanisms of thaw differed between these locations. Hollows had high warm-season soil moisture, which increased thermal conductivity, and deep cold-season snow coverage, which insulated soil. Transition locations thawed primarily due to thermal energy transported through subsurface taliks during individual rain events. Most increases in depth to permafrost occurred below the ∼45 cm thickness seasonally frozen layer, and therefore, expanded existing site taliks. Results highlight the importance of canopy cover and microtopography in controlling soil thermal inputs, the ability of subsurface runoff from individual rain events to trigger warming and thaw, and the acceleration of thaw caused by consecutive wet and snowy years. As northern high-latitudes become warmer and wetter, and weather events become more extreme, the importance of these controls on soil warming and thaw is likely to increase.

Alaska↗

A comparison of contemporary and historical hydrology and water quality in the foothills and coastal plain of the Arctic National Wildlife Refuge, Arctic Slope, northern Alaska

The Arctic National Wildlife Refuge is a unique landscape in northern Alaska with limited water resources, substantial biodiversity of rare and threatened species, as well as oil and gas resources. The region has unique hydrology related to perennial springs, and the formation of large aufeis fields—sheets of ice that grow in the river channels where water reaches the surface in the winter and freezes. This work aims to update our understanding of water resources and water quality in the springs, streams, rivers, and lakes of this region, returning to sites sampled by the U.S. Geological Survey in the 1970s. We resampled eight streams, four springs, and six lakes for hydrological metrics, water quality, and macroinvertebrates, and recalculated flood-frequency metrics for rivers using updated data and modern techniques. Aufeis field melt rates were also assessed for the past several decades. Although the available data preclude trend determinations in most cases, our analysis and comparison to the historical sampling indicates an increase in dissolved ions for streams and springs, faster and earlier aufeis melt, and similar macroinvertebrate populations.

Alaska↗

Insights into glendonite formation from the upper Oligocene Sagavanirktok Formation, North Slope, Alaska

The type locality for the upper Oligocene Nuwok Member of the Sagavanirktok Formation (Carter Creek, North Slope, Alaska, USA) contains abundant occurrence of glendonite, a pseudomorph after the calcium carbonate mineral ikaite, which typically forms in the shallow subsurface of cold marine sediments. The region during the time of Nuwok Member deposition was located at a high latitude, similar to today, and the study site is characterized by sands and silty muds interpreted here to have been deposited in coastal and shelfal marine environments. Isotopic (Sr) and biostratigraphic (foraminifera) evidence presented here refine the depositional age of the outcrop to approximately 24 Ma. Glendonites occur in two basic forms: radial clusters, commonly centered around a single larger primary crystal ( approx. 10 cm; Type A) and larger single blades generally without accessory crystals (approx. 15–25 cm; Type B). Microscopic examination revealed a sequence of multiple types of replacive calcite that formed as a direct result of ikaite transformation: Type 1 rhombohedral crystals characterized by microporous and inclusion-rich cores and concentric zones, Type 2A, composed of clear calcite that overgrew and augmented Type 1 crystals, and inclusion-rich, microcrystalline Type 2B, which formed a matrix surrounding the rhombs and commonly dominates the outer rims of glendonite specimens. Type 3 calcite precipitated as fibrous, botryoidal epitaxial cement atop previous phases and is not ikaite-derived. These phases are distributed in similar ways in all examined specimens and are consistent with several previously described glendonite occurrences around the world, despite differing diagenetic and geologic histories. Stable isotope evidence (δ13C and δ18O) suggests sourcing of glendonite carbon from both organic and methanogenic sources. Glendonites of the Nuwok Member can therefore assist in the determination of a more comprehensive ikaite transformation model, improving our understanding of glendonite formation and the sedimentological and environmental context of their occurrence. Oligocene glendonites are uncommon globally; the well-preserved occurrence described here can allow future studies to better reconstruct Arctic environmental conditions and paleoclimates during this time.

Alaska↗

Quantifying spatiotemporal variation of nearshore forage fish schools with aerial surveys in Prince William Sound, Alaska

Objective Changes in abundance and distribution of schooling forage fish, such as the Pacific Sand Lance Ammodytes hexapterus and Pacific Herring Clupea pallasii , can be difficult to document using traditional boat-based methods, especially in the shallow, nearshore habitats frequented by these species. In contrast, nearshore fish schools are easily observed and quantified from aircraft when light and sea conditions are favorable. We used aerial shoreline surveys to assess interannual variability in the distribution and abundance of schooling forage fish in Prince William Sound, Alaska, during the summers of 2010 and 2012–2022. Methods During the surveys, aerial observers classified fish schools by their size, species, and (in some cases) age-class. All observations were georeferenced along the flight path, converted to estimated surface area (m 2 ) based on school diameter, and standardized by effort (shoreline kilometers surveyed). Result Pacific Herring were widely distributed, and school densities varied annually; there were several spikes in school density of up to 54.38 m 2 /km interspersed among years of lower average densities (7.73–25.57 m 2 /km). In contrast, Pacific Sand Lance were usually limited in their distribution to a few predictable locations. School density in these consistent areas varied across years, from a high of 50.98 m 2 /km in 2010 to a low of 0.15 m 2 /km in 2017. We validated 88 schools during aerial surveys conducted in 2014–2016 and 2019–2022, of which 76 (86%) were correctly identified to species. Conclusion Here, we provide indices of Pacific Herring and Pacific Sand Lance school density over time in shallow nearshore coastal areas of Prince William Sound, Alaska. These indices were generated from aerial surveys, which offer an effective alternative to boat-based surveys for tracking forage fish schools when they occur in shallow and nearshore coastal habitats.

Alaska↗

Database and time series of nearshore waves along the Alaskan coast from the United States-Canada border to the Bering Sea

Alaska’s Arctic coast has some of the highest coastal erosion rates in the world, primarily driven by permafrost thaw and increasing wave energy. In the Arctic, a warming climate is driving sea ice cover to decrease in space and time. A lack of long-term observational wave data along Alaska’s coast challenges the ability of engineers, scientists, and planners to study and address threats and effects from wave-driven erosion and flooding. To overcome the lack of available observational wave data in the nearshore in this study by the U.S. Geological Survey, waves were downscaled with the Simulating WAves Nearshore numerical wave model (SWAN) for the hindcast period of 1979 to 2019 from the United States-Canada border to the Bering Sea utilizing nine model domains. For each domain, the model was forced at the open boundary with 2,500 representative “sea states,” which are likely combinations of significant wave heights, mean wave periods, mean wave directions, and wind speeds and directions. The sea states were obtained from the European Centre for Medium-Range Weather Forecasts “ERA5” dataset for reanalysis of winds and waves using a multivariant maximum-dissimilarity algorithm. The SWAN runs created a downscaled wave database at each grid point, which was used to reconstruct the 40-year time series in the nearshore along the 5- and 10-meter isobaths at locations approximately 400 m apart and corresponding to transects spaced approximately 50 m alongshore, as developed for USGS shoreline-change assessments. Reconstructed time series were compared to observations to validate the numerical model and the downscaled wave database method and showed overall good agreements.

Alaska↗

Characterization of the structural–stratigraphic and reservoir controls on the occurrence of gas hydrates in the Eileen Gas Hydrate Trend, Alaska North Slope

One of the most studied permafrost-associated gas hydrate accumulations in Arctic Alaska is the Eileen Gas Hydrate Trend. This study provides a detailed re-examination of the Eileen Gas Hydrate Trend with a focus on the gas hydrate accumulation in the western part of the Prudhoe Bay Unit. This integrated analysis of downhole well log data and published geophysical data has provided new insight on structural, stratigraphic, and reservoir controls on the occurrence of gas hydrates in the Eileen Gas Hydrate Trend. This study revealed the relatively complex nature of the gas hydrate occurrences in the Eileen Gas Hydrate Trend, with gas hydrates present in a series of coarsening upward, laterally pervasive, mostly fine-grained sand beds exhibiting high gas hydrate saturations. Most of the gas hydrate-bearing reservoirs in the Eileen Gas Hydrate Trend are laterally segmented into distinct northwest- to southeast-trending fault blocks, occur in a combination of structural–stratigraphic traps, and are only partially hydrate filled with distinct down-dip water contacts. These findings suggest that the traditional parts of a petroleum system (i.e., reservoir, gas source, gas migration, and geologic timing of the system formation) also control the occurrence of gas hydrates in the Eileen Gas Hydrate Trend.

Alaska↗

Satellite interferometry landslide detection and preliminary tsunamigenic plausibility assessment in Prince William Sound, southcentral Alaska

Regional mapping of actively deforming landslides, including measurements of landslide velocity, is integral for hazard assessments in paraglacial environments. These inventories are also critical for describing the potential impacts that the warming effects of climate change have on slope instability in mountainous and cryospheric terrain. The objective of this study is to identify slow-moving landslides in the Prince William Sound region, southcentral Alaska, United States, which has had rapid deglaciation since the mid-1800s, and assess their tsunamigenic plausibility. We use an automated time series persistent scatterer interferometric synthetic aperture radar processing method with 7 years of Sentinel-1 data (2016–22) to identify 43 slow-moving slopes with average velocities ranging from approximately 0.2 to 21 millimeters per year. Landslide presence is confirmed using aerial imagery and previous landslide inventory records. We assess the tsunamigenic plausibility of the landslides using empirically derived estimates of landslide mobility based on modeled landslide volumes. Of the identified landslides, our preliminary analysis suggests that 11 have tsunamigenic potential if they were to fail rapidly and catastrophically. Although our estimate of tsunamigenic plausibility is preliminary and can be refined with additional observations and analyses, it can be used to prioritize ongoing and future hazard assessment, surveillance, and research efforts.

Alaska↗

Potential impacts of an autumn oil spill on polar bears summering on land in northern Alaska

Demand for oil and natural gas continues to increase, leading to the development of remote regions where it is riskier to operate. Many of these regions have had limited development, so understanding potential impacts to wildlife could inform management decisions. In 2017, the United States passed legislation allowing oil and gas development in the coastal plain of the Arctic National Wildlife Refuge in northeastern Alaska. This area has received limited industrial development and is an important region for polar bears that use the coastline as a travel corridor in autumn. We sought to understand how an autumn near-shore oil spill in the Refuge could affect polar bears. We simulated oil spills from shallow sub-sea pipelines at 3 locations along the coastline of the Refuge and allowed spills to discharge 4800 barrels of oil per day for 6 days, and tracked oil for 50 days. We interacted the trajectories with simulated polar bear movements to estimate how many bears might be exposed. Oil spread quickly along the coastline and during some weeks exposed an average of 10 bears (95 % CI: 0–37) to lethal levels of oil, and 60 (95 % CI: 9–120) to sub-lethal levels. Our results suggest a significant number of polar bears could become oiled from a spill in the region and require decontamination. Therefore, future mitigation strategies could include careful siting of future development and additional capacity to capture and decontaminate bears in the event of an oil spill.

Alaska↗