Geology Reports⌕ Search

SEARCH · Geology Reports

Search Geology Reports

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.

47 records · Page 3Linked to original sources
USGS2025

Avak Creek oil occurrence, North Slope, Alaska: Newly discovered oil seep on Native lands, near village of Utqiagvik

An unknown occurrence of oil was detected near Avak Creek on Native lands on the North Slope of Alaska. Determining the source of oil was imperative for allowing stakeholders (Federal, State, and local government agencies and the landowner, an Alaska Native corporation) to make timely and informed decisions and mount a mitigation response, if required. The regional and local geological framework of the Avak Creek site was constructed using seismic surveys, well data, and basin modeling results, to identify local petroleum systems, map structural geometry and faults, define source rock thermal maturity distributions, and infer likely oil-migration pathways. Molecular hydrocarbon fingerprints (biomarkers, diamondoids, compound-specific isotopes) of the oil were compared to those of local and regional oil seeps, exploration well tests, and produced oils. Biomarker acid distributions characterized the history and extent of petroleum biodegradation. Integrating subsurface and geochemical parameters, the oil is interpreted to be a natural seep generated locally, predominantly from the Brookian Lower Cretaceous Hue Shale/gamma-ray zone, rather than an anthropogenic source of pollution. Results highlight sophisticated analytical technologies used to characterize complex, compositionally altered hydrocarbons. Results also advance our understanding of Brookian source rock distribution, subsurface petroleum migration pathways, and Arctic Alaska petroleum systems.

Alaska↗
USGS2025

Deformation mechanisms in quartz veins and shear zones elucidate the origin of gold mineralization at Pogo, Alaska

Pogo is a quartz vein hosted, ca. 8 Moz gold deposit. Although it has similarities to orogenic and magmatic-hydrothermal deposits, its origin remains enigmatic. Observations from surface exposures, underground workings, and drill core provide new constraints on quartz vein origins with implications for mineralization. Abundant, largely barren metamorphic segregation quartz veins are found throughout interior Alaska and share characteristics with mineralized quartz veins at Pogo. Pogo quartz veins show crystal plastic deformation fabrics consistent with greenschist to amphibolite facies metamorphism and a lack of internal shear. In contrast, gold in Pogo quartz veins occurs in brittle microfaults and fractures with sulphides such as arsenopyrite and pyrrhotite. Major and minor cataclastic shear zones also exist at Pogo. Cataclastic shear zones commonly cut the mineralized quartz veins and porphyroclasts are dominated by quartz with the same plastic deformation fabrics found in the major Pogo quartz veins. The porphyroclasts, and the carbonaceous clay-rich matrix they sit in, also contain gold indicating that the shear zones postdate quartz and earliest gold deposition. These observations suggest that competency contrasts between the quartz veins and their phyllosilicate-rich host rocks controlled transient permeability formed by late, preferential brittle deformation events localizing sulphide ± gold deposition.

Alaska↗
USGS2025

Combining scanning electron microscopy, X-ray diffraction, and X-ray fluorescence to characterize shear zones at the Pogo gold deposit, Alaska

This study employs a multi-method analytical approach to characterize the mineralogical, geochemical, and textural properties of fault rocks from the Pogo gold mine in the Yukon-Tanana Upland, central Alaska. Specifically, we examine cataclasites, to document the structural and geochemical evolution of shear zones and their associations with gold mineralization. To investigate the shear zone, we integrate portable X-ray fluorescence (pXRF), scanning electron microscopy-based automated mineralogy (SEM-AM), X-ray diffraction (XRD), and high-resolution micro-X-ray fluorescence (micro-XRF) mapping. These methods collectively provide insights into bulk and trace element chemistry, mineralogical composition, and deformation-related textures across multiple scales. Handheld pXRF enables rapid geochemical screening, guiding SEM-AM and XRD analyses to ensure consistent mineralogical interpretation. X-ray diffraction identifies and quantifies crystalline phases, while SEM-AM produces high-resolution mineral maps, revealing mineral abundances, grain-scale textures, and gold associations. Micro-XRF mapping further refines our understanding by showing visual trace element distributions at sub-millimetre resolution. By integrating these techniques, we improve our understanding of the nature and geochemistry of Pogo shear zones, their role in gold mineralization, and support metallurgical processing strategies. This approach enhances exploration models and resource characterization for structurally complex gold deposits.

Alaska↗
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2024

Simulating present and future groundwater/surface-water interactions and stream temperatures in Beaver Creek, Kenai Peninsula, Alaska

In many places, coldwater ecosystems are facing increasing pressure from anthropogenic warming. This study examined stream temperatures and the water balance in the Beaver Creek watershed on the Kenai Peninsula in south-central Alaska—an area that is experiencing rapid warming. Low-gradient streams near the Kenai coast provide important spawning and rearing habitat for salmon but may be especially vulnerable to rising temperatures, because of long residence times, inflows from abundant riparian wetlands, and reliance on groundwater discharge that may also warm, or decrease in volume with rising evapotranspiration. In recent decades, observed maximum 7-day temperatures have consistently exceeded statistical (regression-based) projections. Here we simulate total streamflows and temperatures with a physics-based model that links the Soil Water Balance, MODFLOW 6 and SNTEMP simulation codes on a 7-day timestep. The model is based on existing data and groundwater levels, instream flows, and stream temperatures collected during 2019–23. Future climate scenarios were developed for 2023–50 from downscaled climate projections. Results indicate that groundwater discharge is about 64 percent of the total streamflow during the months of May through September. Total streamflow and groundwater discharge are expected to remain similar to current conditions through 2050. Stream temperatures are expected to rise; by midcentury, near the Beaver Creek mouth the model predicts 34 to 63 additional days per year with average weekly temperatures above 13 degrees Celsius, 14 to 81 additional days with average weekly temperatures above 15 degrees Celsius, and routine exceedances of 20 degrees Celsius during the warmest periods. Projected stream temperatures vary spatially. Areas of high groundwater inflows in the lower main stem and some tributaries may be most resilient to warming air temperatures during dry conditions. During storm events, groundwater-dominated tributaries may have the coolest stream temperatures.

Alaska↗
USGS2024

Constraining the earthquake recording threshold of intraslab earthquakes with turbidites in southcentral Alaska’s lakes and fjords

Strong ground motion from intraslab earthquakes, which do not produce primary paleoseismic evidence, may initiate gravity-driven turbidity flows in subaqueous basins. The resulting deposits (turbidites) can provide a paleoseismic proxy if the conditions that initiate these flows are known. To better constrain the initiating conditions, we use two recent intraslab earthquakes in southcentral Alaska, the M w 7.1 30 November 2018 Anchorage earthquake and the M w 7.1 24 January 2016 Iniskin earthquake, as calibration events. Through a multilake investigation, we document the occurrence, or the absence, of earthquake-generated turbidity flows from these two events. Both earthquakes are recorded by centimeter-scale turbidites that can be differentiated from climatically generated deposits, as well as other seismic sources based on deposit thickness, sedimentological properties, and deposit age. We show that a Modified Mercalli Intensity (MMI) of ∼V–V1/2 is the minimum shaking intensity required to generate localized sediment remobilization from deltaic slopes, and an MMI of ∼V1/2 is required to produce a deposit of sufficient thickness that a seismic origin can be confidently assigned. The documentation of seismically generated deposits in quick succession (∼2 years) with diagnostic features highlights the utility of using recent earthquakes as calibration events to investigate the subaqueous response to strong ground motion.

Alaska↗
USGS2024

Oblique contraction along the fastest ocean-continent transform plate boundary focuses rock uplift west of the Fairweather fault, southeast Alaska

Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (>25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes >20° to the west of plate-boundary motion.

Alaska↗
USGS2024

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska↗
USGS2024

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2023

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska↗