USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Investigations Map”

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.

862 records · Page 22Linked to original sources
USGS2025

Cenozoic stratigraphy of Colorado

As a successor to previous Colorado stratigraphy charts (MS-53 Colorado Stratigraphic Chart and MS-54 Cretaceous Stratigraphy of Colorado), this Colorado Geological Survey (CGS) publication resulted from a collaboration between the CGS, USGS, and the Denver Museum of Nature and Science (DMNS). The chart was designed to illustrate Cenozoic stratigraphy spanning the state’s many sedimentary basins. It builds upon the work of dozens of colleagues and updates Richard Pearl’s seminal 1974 stratigraphy chart. The chart leverages the community’s stratigraphic work in both the subsurface and outcrop, and depicts new geochronologic constraints for many units. To facilitate comparison of strata to external forcing factors, the chart employs a linear timescale. Each unit’s dominant depositional environment is depicted, as are major mountain building events, erosional events, and regional unconformities.

Colorado↗
USGS2026

Stepovers and beyond: Structural control of The Geysers geothermal system and the broader Clear Lake region

Fault geometry exerts a first-order control on geothermal systems by governing stress localization, fracture development, and permeability, yet in complex fault networks or broader shear zones, the relative influence of individual geometric features is often difficult to resolve. In the northern California Coast Ranges, The Geysers geothermal field is commonly interpreted to occur within a releasing stepover, although no single, clearly defined stepover is identified in published studies. To investigate the structural controls on The Geysers and the broader Clear Lake region, a two-dimensional elastic boundary element model is developed to evaluate spatial patterns of dilational strain associated with progressively more complete fault geometries. Model results show that dilation in the region is not controlled by a single structure but instead reflects the combined effects of multiple interacting fault elements. Three primary controls are identified: (1) opposing bends in the regional strike-slip fault system, including a releasing bend along the Maacama fault; (2) the southern fault tip of the Collayomi fault, which generates a prominent dilational lobe beneath the southern Geysers; and (3) a releasing stepover between the Collayomi fault and the Geyser Peak–Mercuryville–Big Sulphur Creek fault system, inferred to collectively behave as a right-lateral shear zone bounding the western margin of The Geysers. Predicted dilational strain magnitudes are sufficient to localize permeability between faults. These results highlight that incorporating complete fault networks and bedrock geological mapping can enhance geothermal assessments and provide a transferable framework for evaluating structurally controlled permeability in tectonically active regions.

California↗
USGS2024

Joint pilot fish habitat framework

This story map will take you through the process of exploring and testing methods necessary for a higher resolution, seamless fish habitat assessment across both inland and estuarine waters through the lens of our joint pilot assessment Fish habitat assessments attempt to relate past, current, or future landscape conditions to the state of fish species occurrence, distribution, abundance, or community and habitat condition in streams, rivers, or estuaries. Previous fish habitat assessments, such as the National Fish Habitat Assessment, conducted separate and disconnected assessments for inland waters and estuaries. In this project, National Oceanic and Atmospheric Administration ( NOAA ) and U.S. Geological Survey ( USGS ) researchers created a seamless spatial framework to allow assessments that integrate influences on fish habitat from headwaters to the estuary. This effort began when the Chesapeake Bay Program Fish Habitat Action Team expressed interest in a Baywide fish habitat assessment spanning tidal salt, tidal fresh, warm non-tidal, and cold non-tidal waters. However, the complexity of the myriad of implementation details to consider when developing such an assessment necessitated the need for a tributary-specific pilot assessment. To conduct this pilot assessment, a NOAA/USGS joint partnership was formed with cooperation and support from the Chesapeake Bay Agreement and Chesapeake Bay Fish Habitat Action Team (FHAT).

Maryland↗
USGS2026

Software to support remote sensing of river discharge based on critical flow theory

Water resource management requires accurate observations of streamflow but standard field methods for measuring river discharge ( Q ) are costly and can be hazardous for equipment and personnel. Remote sensing has become a viable alternative, but many image-based techniques require field data for calibration and depth and velocity can seldom be mapped with a single sensor. A new approach based on critical flow theory, in contrast, allows both of these attributes to be inferred from readily available image data. This technique only pertains to sites with standing waves, called undular hydraulic jumps (UHJs), but a recent investigation demonstrated its potential to provide accurate discharge estimates. This paper introduces software designed to facilitate Inferring Q from UHJs Identified in River Images (InQUIRI). The package includes modules for retrieving data from image servers, making the measurements of wavelength and width required to calculate discharge, inferring a representative wavelength from a profile digitized along a wave train, combining multiple estimates to obtain an ensemble median discharge, and assessing accuracy via comparison to gage records from the U.S. Geological Survey. By making these steps easier to implement, InQUIRI enables users to apply the workflow to a variety of UHJ-containing images. Accumulating more case studies, some successful and others less so, would help constrain the range of applicability of the critical flow approach and foster development of refined guidelines for selecting and measuring waves. The software described herein could play an important role in promoting informed use of this new technique for non-contact streamflow measurement.

Arizona, Colorado, New Mexico, Utah↗
USGS2026

Methodology for construction of a three-layer geologic model of the conterminous United States using land surface, top of bedrock, and top of basement

This report describes the methodology used for the construction of a digital three-layer geologic model of the conterminous United States by mapping the altitude of three surfaces: land surface, the top of bedrock, and the top of basement. These surfaces are mapped through the compilation and synthesis of published stratigraphic horizons from numerous topical studies. The mapped surfaces create a three-layer geologic model with three geomaterial-based subdivisions: unconsolidated to weakly consolidated sediment; layered consolidated rock strata that constitute bedrock; and crystalline rocks that are described as “basement,” consisting of either igneous, metamorphic, or highly deformed rocks. The data compilation and synthesis are highly dependent on the definition of the informal terms “bedrock” and “basement,” which may describe different ages or types of rock in different parts of the conterminous United States. This report presents the conceptualization of the three mapped layers, describes the datasets used, and summarizes the decisions made while compiling the three-layer model from the various sources. This digital dataset was created as part of efforts by the U.S. Geological Survey to develop subsurface geologic data in geospatial form as part of a broad directive to develop two-dimensional and three-dimensional geologic information at detailed, national, and continental scales. This digital dataset partly fulfills the goal of the U.S. Geological Survey’s National Cooperative Geologic Mapping Program to construct a national-scale three-dimensional geologic model.

conterminous United States↗
USGS2025

Unlocking ecological insights from sub-seasonal visible-to-shortwave infrared imaging spectroscopy: The SHIFT campaign

We stand at the threshold of a transformative era in Earth observation, marked by space-borne visible-to-shortwave infrared (VSWIR) imaging spectrometers that promise consistent global observations of ecosystem function, phenology, and inter- and intra-annual change. However, the full value of repeat spectroscopy, the information embedded within different temporal scales, and the reliability of existing algorithms across diverse ecosystem types and vegetation phenophases have remained elusive due to the absence of suitable sub-seasonal spectroscopy data. In response, the Surface Biology and Geology (SBG) High-Frequency Time Series (SHIFT) campaign was initiated during late February 2022 in Santa Barbara County, California. SHIFT, designed to support NASA's SBG mission, addressed mission scoping, scientific advancement, applications development, and community building. This ambitious endeavor included weekly Airborne Visible InfraRed Imaging Spectrometer-Next Generation (AVIRIS-NG) imagery acquisitions for 13 weeks (spanning February 24 to May 29, 2022), accompanied by coordinated terrestrial vegetation and coastal aquatic data collection. We describe the rich datasets collected and illustrate how the complex sub-seasonal patterns of change can be linked to biological science and applications, surpassing insights from multispectral observations. Leveraging open-source processing methods and cloud-based analysis tools, the SHIFT campaign showcases the readiness of the scientific community to harness ecological insights from remotely sensed hyperspectral time series. We provide an overview of SHIFT's goals, data collections, preliminary results, and the collaborative efforts of early career scientists committed to unlocking the transformative potential of high-frequency time series data from space-borne VSWIR imaging spectrometers.

Ecosphere↗
USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2025

Stories in stone: The geoheritage significance of U.S. National Park Service stratotypes

The U.S. National Park Service stratotype inventory has systematically documented hundreds of published stratotypes across the country that represent a quintessential component of America’s geoheritage and possess significant scientific, educational, cultural, historic, and aesthetic values. As valuable geologic reference standards, stratotypes record intervals of Earth history relating to our nation’s unique geologic evolution, biological progression, cultural traditions, historical narratives, and stunning landforms. However, the general lack of awareness and knowledge regarding the scope, significance, distribution, and scale of stratotypes can result in these resources being overlooked. Consequently, stratotypes may inadvertently experience accidental alteration, damage, or destruction through both natural and anthropogenic processes. As with many American geoheritage sites, stratotypes lack specific recognition through legislation, policies, or other management strategies to ensure their preservation and protection. Through education and outreach, we hope to inform, increase awareness, and engage the public about these internationally and nationally significant geologic reference sites, unite various geoheritage communities, and promote their preservation for future generations.

Geoheritage↗
USGS2026

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration↗
USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗
USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗
USGS2026

Widespread landslide activity in an extreme wet season and implications for regional sediment management, eastern San Francisco Bay area, California

Watershed sediment production is expected to increase in a warmer future with more extreme rain, with cascading effects throughout drainage and sediment-transport networks. This study investigated landscape-scale sediment movement in the eastern San Francisco Bay area, California, USA, during the extreme 2016–2017 wet season that brought major rainfall, landslides, and flooding. Mapping 8,928 landslides across a 1,050-km 2 study area revealed new sediment yield of 510–956 t/km 2 , equivalent to denudation of 193–361 mm/ky. These results correspond closely to long-term denudation rates in the northern and central California Coast Ranges, indicating that mass wasting in very wet years dominates long-term sediment mobilization. However, due to long residence times in drainage networks, the 2017 landslides contributed at most ∼1%–2% of the estimated locally derived fluvial sediment transport to San Francisco Bay. Although the amount of sediment mobilized did not threaten municipal water supplies, small rangeland impoundments in this mixed-use landscape lost storage capacity to new sedimentation. Considering regional sediment supply and demand, even the exceptionally large sediment delivery in an extreme wet year cannot meet the need for sediment to accrete tidal wetlands in the bay. To keep pace with rising sea levels, this abnormally high terrestrial sediment input would need to occur in 50 of the next 75 years, an unlikely occurrence due to the prevalence of recent drought years. Shoreline protection and restoration in the bay would need additional sources of sediment, such as through management of dredged sediment through beneficial-reuse programs.

California↗
USGS2025

Ordovician stratigraphy, structure, and karst of the Falling Spring Valley, Alleghany County, Virginia, USA

This one-day trip highlights new findings on a preliminary bedrock geologic map that shows results from ongoing geologic mapping in the Falling Spring Valley of Alleghany County, Virginia, USA, which is the southern end of the larger Warm Springs Valley, an elongated anticlinal valley rimmed by Ordovician and Silurian siliciclastic rocks, and which is famous for its thermal springs. This mapping includes stratigraphic, structural, and karst field and lab research focused on the Ordovician strata exposed in the area, the oldest of which is the dolomitic upper part of the Beekmantown Formation (Lower Ordovician, Darriwilian), and the youngest of which is the Juniata Formation (Upper Ordovician, Katian), a sequence of siliciclastic redbeds. Warm Springs Valley is the location of the only known caves in the eastern United States—three at present—with thermal waters flowing in some of their passages. Stops on the trip will highlight key details from mapping efforts, primarily within the structurally deformed Ordovician carbonate sequence that is exposed in the core and limbs of the anticline, as well as the associated karst features that are developed in those carbonate rocks, including results of recent dye traces and water temperature monitoring that have improved our understanding of the karst hydrogeologic systems developed in these strata.

Virginia↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2023

Overview of the Cenozoic geology of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia

The Harrat Rahat volcanic field, located in the west-central part of the Kingdom of Saudi Arabia, is one of the larger Cenozoic harrats among the more than 17 harrats situated upon the Arabia Plate. The map plate contained herein shows, at a scale of 1:100,000, the mapped volcanic geology of northern Harrat Rahat, which consists of the northernmost one-fifth of Harrat Rahat. Northern Harrat Rahat has an area of about 3,340 square kilometers (km 2 ), of which 2,567 km 2 is covered by Harrat Rahat volcanic rocks and deposits, and it encompasses more than 900 exposed vents (that is, craters, cryptodomes, maars, and scoria cones), 289 of which are isolated by younger volcanic deposits and have not been correlated with the 234 volcanic rock units distinguished by geologic mapping. Northern Harrat Rahat is of special interest owing to its proximity to the city of Al Madinah al Munawwarah, which sits within, and is continuing to expand southward over, the north end of the volcanic field. Al Madinah is home to an expanding population, currently at more than 2 million residents, together with the intermittent addition of approximately 3 million pilgrims during Hajj and Umrah (religious visitations). The center of Al Madinah is less than 8 km from lava flows of the only confirmed historically documented eruption, which occurred in 1256 C.E. (654 A.H.). Earlier prehistoric lava flows also encroached into the area of the present-day city limits, as demonstrated by volcanic rocks exposed widely throughout the city in roadcuts, parks, and excavations for new buildings, although no evidence has been found of any other than the 1256 C.E. lava having reached that area during times of human habitation. Geologic mapping was undertaken by the U.S. Geological Survey in collaboration with the Saudi Geological Survey. The features of primary interest within the map area are scoria cones, lava flows, lava domes, craters, and pyroclastic deposits from the Quaternary, which have compositions of basalt, hawaiite, mugearite, benmoreite, and trachyte. The geologic mapping was published by Downs and others (2019) on a single sheet at 1:75,000 scale and two detailed sheets at 1:25,000 scale, accompanied by detailed explanations of the geology. This report presents the geology on a single sheet at 1:100,000 scale and provides condensed geologic explanations for the convenience of readers of this volume. Some minor errors of nomenclature and lava source regions that do not change fundamental interpretations have been corrected herein.

Professional Paper↗
USGS2024

Complex landslide patterns explained by local intra-unit variability of stratigraphy and structure: Case study in the Tyee Formation, Oregon, USA

Lithology and geologic structure are important controls on landslide susceptibility and are incorporated into many regional landslide hazard models. Typically, metrics for mapped geologic units are used as model input variables and a single set of values for material strength are assumed, regardless of spatial heterogeneities that may exist within a map unit. Here we describe how differences in bedding thickness, grain size, inferred uniaxial compressive strength, and bedding dip control the inherent susceptibility of slopes to deep-seated failure within a single mapped geologic unit - the Tyee Formation of Oregon, USA. The Tyee, which covers over 15,000 km2 and underlies much of the Oregon Coast Range, comprises gently folded alternating beds of sandstone and siltstone deposited as turbidites, forming a 2-km thick Eocene submarine fan which has been uplifted and exhumed through the Cenozoic. Deep-seated landslides are widespread in the Tyee, but form a complex spatial pattern such that landslide density ranges from 0 to 24% of the total landscape area. These slides are often extensive and sufficiently deep to reduce local hillslope gradients, resulting in a strong negative correlation between landslide density and mean local slope. Mean annual precipitation and predicted strong ground motions from Cascadia earthquake scenarios also fail to explain the spatial distribution of deep-seated landslides. Consequently, landslide stability models, which are strongly influenced by landscape slope, pore-water pressure, and seismic acceleration, yield landslide susceptibility maps which are broadly anti-correlated with mapped deep-seated landslide density. Through a multivariable linear regression model, we show that much of the variance in deep-seated landslide density can be explained by variability of intra-unit stratigraphic and structural characteristics, which we measure at 128 sites across two study areas totaling ∼3000 km2. Our results suggest bedding dip is only weakly correlated to landslide density, but strongly influences landslide failure style. Subtle increases in bedding dip, even in the gently folded Tyee Formation, result in a substantially higher likelihood of a landslide being cataclinal, or parallel to bedding. Overall, we find a slight majority of landslides fail within these cataclinal slopes, and that these landslides tend to be larger than non-cataclinal landslides. We also show that the lithological and structural properties that influence landslide susceptibility are distinct for these two populations of landslides. Our results demonstrate how localized, intra-unit, geologic variability can exert strong control on landslide susceptibility and failure style. This suggests that in some locations, landslide hazard models could be significantly improved by incorporating detailed, spatially variable, geologic properties rather than relying solely on generalized geologic map units.

Oregon↗