USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Maps”

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

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

At least 901 records · Page 50Linked to original sources

Differential uplift and incision of the Yakima River terraces, central Washington State

The fault-related Yakima folds deform Miocene basalts and younger deposits of the Columbia Plateau in central Washington State. Geodesy implies ~2 mm/yr of NNE directed shortening across the folds, but until now the distribution and rates of Quaternary deformation among individual structures has been unclear. South of Ellensburg, Washington, the Yakima River cuts a ~600 m deep canyon across several Yakima folds, preserving gravel-mantled strath terraces that record progressive bedrock incision and related rock uplift. Here we integrate cosmogenic isochron burial dating of the strath terrace gravels with lidar analysis and field mapping to quantify rates of Quaternary differential incision and rock uplift across two folds transected by the Yakima River: Manastash and Umtanum Ridge. Isochron burial ages from in situ produced 26 Al and 10 Be at seven sites across the folds date episodes of strath terrace formation over the past ~2.9 Ma. Average bedrock incision rates across the Manastash (~88 m/Myr) and Umtanum Ridge (~46 m/Myr) anticlines are roughly 4 to 8 times higher than rates in the intervening syncline (~14 m/Myr) and outside the canyon (~10 m/Myr). These contrasting rates demonstrate differential bedrock incision driven by ongoing Quaternary rock uplift across the folds at rates corresponding to ~0.13 and ~0.06 mm/yr shortening across postulated master faults dipping 30 ± 10°S beneath the Manastash and Umtanum Ridge anticlines, respectively. The reported Quaternary shortening across the anticlines accounts for ~10% of the ~2 mm/yr geodetic budget, suggesting that other Yakima structures actively accommodate the remaining contemporary deformation.

Washington↗

Three-dimensional geophysical mapping of shallow water saturated altered rocks at Mount Baker, Washington: Implications for slope stability

Water-saturated hydrothermal alteration reduces the strength of volcanic edifices, increasing the potential for catastrophic sector collapses that can lead to far traveled and destructive debris flows. Intense hydrothermal alteration significantly lowers the resistivity and magnetization of volcanic rock and therefore hydrothermally altered rocks can be identified with helicopter electromagnetic and magnetic measurements. Geophysical models constrained by rock properties and geologic mapping show that intensely altered rock is restricted to two small (500 m diameter), >150 m thick regions around Sherman Crater and Dorr Fumarole Field at Mount Baker, Washington. This distribution of alteration contrasts with much thicker and widespread alteration encompassing the summits of Mounts Adams and Rainier prior to the 5600 year old Osceola collapse, which is most likely due to extreme erosion and the limited duration of summit magmatism at Mount Baker. In addition, the models suggest that the upper ~300 m of rock contains water which could help to lubricate potential debris flows. Slope stability modeling incorporating the geophysically modeled distribution of alteration and water indicates that the most likely and largest (~0.1 km 3 ) collapses are from the east side of Sherman Crater. Alteration at Dorr Fumarole Field raises the collapse hazard there, but not significantly because of its lower slope angles. Geochemistry and analogs from other volcanoes suggest a model for the edifice hydrothermal system.

Washington↗

Estimation of wave phase speed and nearshore bathymetry from video imagery

A new remote sensing technique based on video image processing has been developed for the estimation of nearshore bathymetry. The shoreward propagation of waves is measured using pixel intensity time series collected at a cross-shore array of locations using remotely operated video cameras. The incident band is identified, and the cross-spectral matrix is calculated for this band. The cross-shore component of wavenumber is found as the gradient in phase of the first complex empirical orthogonal function of this matrix. Water depth is then inferred from linear wave theory's dispersion relationship. Full bathymetry maps may be measured by collecting data in a large array composed of both cross-shore and longshore lines. Data are collected hourly throughout the day, and a stable, daily estimate of bathymetry is calculated from the median of the hourly estimates. The technique was tested using 30 days of hourly data collected at the SandyDuck experiment in Duck, North Carolina, in October 1997. Errors calculated as the difference between estimated depth and ground truth data show a mean bias of −35 cm (rms error = 91 cm). Expressed as a fraction of the true water depth, the mean percent error was 13% (rms error = 34%). Excluding the region of known wave nonlinearities over the bar crest, the accuracy of the technique improved, and the mean (rms) error was −20 cm (75 cm). Additionally, under low-amplitude swells (wave height H ≤ 1 m), the performance of the technique across the entire profile improved to 6% (29%) of the true water depth with a mean (rms) error of −12 cm (71 cm).

North Carolina↗

Fluvial transport and surface enrichment of arsenic in semi-arid mining regions: examples from the Mojave Desert, California

As a result of extensive gold and silver mining in the Mojave Desert, southern California, mine wastes and tailings containing highly elevated arsenic (As) concentrations remain exposed at a number of former mining sites. Decades of weathering and erosion have contributed to the mobilization of As-enriched tailings, which now contaminate surrounding communities. Fluvial transport plays an intermittent yet important and relatively undocumented role in the migration and dispersal of As-contaminated mine wastes in semi-arid climates. Assessing the contribution of fluvial systems to tailings mobilization is critical in order to assess the distribution and long-term exposure potential of tailings in a mining-impacted environment. Extensive sampling, chemical analysis, and geospatial mapping of dry streambed (wash) sediments, tailings piles, alluvial fans, and rainwater runoff at multiple mine sites have aided the development of a conceptual model to explain the fluvial migration of mine wastes in semi-arid climates. Intense and episodic precipitation events mobilize mine wastes downstream and downslope as a series of discrete pulses, causing dispersion both down and lateral to washes with exponential decay behavior as distance from the source increases. Accordingly a quantitative model of arsenic concentrations in wash sediments, represented as a series of overlapping exponential power-law decay curves, results in the acceptable reproducibility of observed arsenic concentration patterns. Such a model can be transferable to other abandoned mine lands as a predictive tool for monitoring the fate and transport of arsenic and related contaminants in similar settings. Effective remediation of contaminated mine wastes in a semi-arid environment requires addressing concurrent changes in the amounts of potential tailings released through fluvial processes and the transport capacity of a wash.

California↗

Partitioning of pyroclasts between ballistic transport and a convective plume: Kīlauea volcano, 19 March 2008

We describe the discrete ballistic and wind-advected products of a small, but exceptionally well-characterized, explosive eruption of wall-rock-derived pyroclasts from Kīlauea volcano on 19 March 2008 and, for the first time, integrate the size distribution of the two subpopulations to reconstruct the true size distribution of a population of pyroclasts as it exited from the vent. Based on thinning and fining relationships, the wind-advected fraction had a mass of 6.1 × 10 5 kg and a thickness half distance of 110 m, placing it at the bottom end of the magnitude and intensity spectra of pyroclastic falls. The ballistic population was mapped, in the field and by using structure-from-motion techniques, to a diameter of > 10–20 cm over an area of ~0.1 km 2 , with an estimated mass of 1 × 10 5 kg. Initial ejection velocities of 50–80 m/s were estimated from inversion of isopleths. The total grain size distribution was estimated by using a mass partitioning of 98% of wind-advected material and 2% of ballistics, resulting in median and sorting values of −1.7 ϕ and 3.1 ϕ . It is markedly broader than those calculated for the products of magmatic explosive eruptions, because the grain size of 19 March 2008 clast population is unrelated to a volcanic fragmentation event and instead was “inherited” from a population of talus clasts that temporary blocked the vent prior to the eruption. Despite a conspicuous near-field presence, the ballistic subpopulation has only a minor influence on the grain size distribution because of its rapid thinning and fining away from source.

Hawaii↗

Interpretation of depositional environment in the Plympton Formation (Permian), Southern Pequop Mountains, Nevada, from physical stratigraphy and a faunule

Field mapping in the southern part of the Pequop Mountains has shown the presence of a major structural high which has profoundly affected the stratigraphy of Pennsylvanian and Lower Permian beds; lesser effects of this high persisted into the later Permian. Within the Park City Group (Permian), the Plymton Formation, overlying the Kaibab Limestone, is divided informally into five parts. The uppermost Plympton has yielded silicified fossils. The faunule studied is characterized by the scaphopod Plagioglypta and by bellerophontacean and neritacean gastropods, suggestive of a middle Permian age. Specimens are worn, and several lines of evidence suggest that this is a high-energy beach deposit associated with a residual high. The interpretation supplements the field evidence of a hiatus between the Plympton and the overlying Gerster Formation. Other occurrences of elements of the faunule in the equivalent Phosphoria rock complex are noted, but only a few were probably deposited in a high-energy beach environment. A new species of Naticopsis that has the color pattern preserved is named and described.

Nevada↗

Constraining the relative importance of raindrop- and flow-driven sediment transport mechanisms in postwildfire environments and implications for recovery time scales

Mountain watersheds recently burned by wildfire often experience greater amounts of runoff and increased rates of sediment transport relative to similar unburned areas. Given the sedimentation and debris flow threats caused by increases in erosion, more work is needed to better understand the physical mechanisms responsible for the observed increase in sediment transport in burned environments and the time scale over which a heightened geomorphic response can be expected. In this study, we quantified the relative importance of different hillslope erosion mechanisms during two postwildfire rainstorms at a drainage basin in Southern California by combining terrestrial laser scanner-derived maps of topographic change, field measurements, and numerical modeling of overland flow and sediment transport. Numerous debris flows were initiated by runoff at our study area during a long-duration storm of relatively modest intensity. Despite the presence of a well-developed rill network, numerical model results suggest that the majority of eroded hillslope sediment during this long-duration rainstorm was transported by raindrop-induced sediment transport processes, highlighting the importance of raindrop-driven processes in supplying channels with potential debris flow material. We also used the numerical model to explore relationships between postwildfire storm characteristics, vegetation cover, soil infiltration capacity, and the total volume of eroded sediment from a synthetic hillslope for different end-member erosion regimes. This study adds to our understanding of sediment transport in steep, postwildfire landscapes and shows how data from field monitoring can be combined with numerical modeling of sediment transport to isolate the processes leading to increased erosion in burned areas.

California↗

Geothermal gradients in the conterminous United States

Geothermal gradients from published temperature/depth measurements in drill holes generally deeper than 600 m are used to construct a temperature gradient map of the conterminous United States. The broadly contoured map displays 284 temperature gradients that are applicable to a depth of 2 km. In terms of the number of contoured areas and the fraction of data points having a value not within a contour interval (outliers), the temperature gradient data and associated deep heat flow data have similar measures of contourability. Areally, most of the United States is contoured from 15° to 35°C/km. The eastern United States is generally cooler (average 25°C/km) than the western United States (average 34°C/km), in accordance with broad heat flow trends. Differences between the temperature gradient and heat flow maps are caused by areal differences in rock thermal conductivities. The effect of conductivity on gradients is particularly apparent in the eastern United States where heat flow is relatively constant over large areas. Gradients are elevated where thick, low-conductivity, sedimentary deposits occur such as in the Atlantic Coastal Plain province and in basins in the Allegheny Plateau and the Great Plains provinces. No clear gradient pattern emerges where both heat flow and conductivity vary widely, such as in the northern Basin and Range and Rocky Mountain provinces. Using the temperature gradients determined in this study and associated heat flow values, derived thermal conductivities are calculated for the depth range of a few hundred meters to 2 km; the average conductivity is 2.5±0.8 W/m °C. Some areas show little variation in derived thermal conductivity, while others show a wide range. The Atlantic Coastal Plain, Appalachian Plateaus, Superior Upland, Gulf Coastal Plain, Columbia Plateaus and Cascade Mountains, and Sierra Nevada provinces have restricted ranges of conductivities, which supports the concept of regional conductivities.

Journal of Geophysical Research Solid Earth↗

Possible refugia in the Alexander Archipelago of southeastern Alaska during the late Wisconsin glaciation

The interpretation of the extent of late Wisconsin glaciation in southeastern Alaska has varied between geologists and biologists. Maps and reports of the region prepared by geologists commonly indicated that late Wisconsin ice extended as a large uniform front west to the edge of the continental shelf. However, the distribution of plants and animals in the region has led many biologists to suggest that there may have been ice-free areas that served as refugia during the late Wisconsin. Based on analyses of aerial photographs, topographic maps, and bathymetric charts, in conjunction with a review of previous literature and reconnaissance fieldwork throughout the region, this study presents data supporting a limited ice extent in the Alexander Archipelago during the late Wisconsin and identifies possible ice-free areas that may have served as refugia. These areas include (1) the Fairweather Ground, (2) the Herbert Graves Island area, (3) the western coast of southern Baranof Island and adjacent continental shelf, (4) Coronation Island and the adjacent continental shelf, (5) the Warren Island area, (6) the continental shelf from west of Heceta Island to Forrester Island in the south, (7) parts of the west coast of southern Dall Island, and (8) lowland areas in southern Prince of Wales Island. The identification of these possible refugia has bearing on the recolonization of the Alexander Archipelago, as they could have served as centers of biotic dispersal upon regional deglaciation and as stepping stones for early humans with a maritime tradition entering the western hemisphere from Asia. ?? 2007 NRC Canada.

Canadian Journal of Earth Sciences↗

Why the 2002 Denali fault rupture propagated onto the Totschunda fault: implications for fault branching and seismic hazards

The propagation of the rupture of the M w 7.9 Denali fault earthquake from the central Denali fault onto the Totschunda fault has provided a basis for dynamic models of fault branching in which the angle of the regional or local prestress relative to the orientation of the main fault and branch plays a principal role in determining which fault branch is taken. GeoEarthScope LiDAR and paleoseismic data allow us to map the structure of the Denali-Totschunda fault intersection and evaluate controls of fault branching from a geological perspective. LiDAR data reveal the Denali-Totschunda fault intersection is structurally simple with the two faults directly connected. At the branch point, 227.2 km east of the 2002 epicenter, the 2002 rupture diverges southeast to become the Totschunda fault. We use paleoseismic data to propose that differences in the accumulated strain on each fault segment, which express differences in the elapsed time since the most recent event, was one important control of the branching direction. We suggest that data on event history, slip rate, paleo offsets, fault geometry and structure, and connectivity, especially on high slip rate-short recurrence interval faults, can be used to assess the likelihood of branching and its direction. Analysis of the Denali-Totschunda fault intersection has implications for evaluating the potential for a rupture to propagate across other types of fault intersections and for characterizing sources of future large earthquakes.

Alaska↗

Detrital zircons and the magmatic history of Viti Levu, Fiji

We integrate the existing detrital zircon data from multiple modern river sediment samples on Viti Levu, Fiji, with the most current available geological and topographic mapping of the respective river drainage basins to compare detrital populations with potential bedrock sources. The temporal and spatial variations in zircon geochemistry supplement what is known from igneous rocks and confirm the petrological differences between plutonic and volcanic rocks from the Eocene to early Oligocene (Yavuna age, >30 Ma), middle Oligocene to middle Miocene (Wainimala age, 30–12.5 Ma), late Miocene (Colo age, 12.5–6.5 Ma) and latest Miocene (Namosi age, 6.5–5 Ma). The >30 Ma Yavuna-age zircons are restricted to areas that drain the previously mapped Yavuna Group. The 30–12.5 Ma zircons are found across central Viti Levu from west to east, and the 30–15 Ma zircons have distinctively low U/Yb and high Dy/Yb ratios. They are the best radiometric evidence of widespread early to middle Miocene arc magmatism in Fiji that was relatively U-poor. Peak deconvolution of the Colo age zircons from individual basins suggests the following ages for undated or poorly dated plutons from central Viti Levu. The large Mavuvu pluton is probably composed of multiple intrusions in the 12–10 Ma range, the Waiqa pluton is probably ca 10 Ma, and the Noikoro pluton is probably ca 9 Ma. There are zircons from unknown plutonic or volcanic sources between 8 and 7 Ma in western Viti Levu that have distinct Eu/Eu* ratios. We attribute the highest U/Yb ratios in some Colo age zircons to crustal anatexis. Namosi-age zircons are abundant in the Medrausucu Group and can be found scattered across Viti Levu.

Viti Levu↗

Dynamic ice-wallow relief of northern Alaska's nearshore

Contour maps with 0.5-m depth interval were prepared for a small area seaward of Reindeer Island, a barrier island in the Beaufort Sea, Alaska, by repeated surveys with very accurate navigation and very close trackline spacing. The maps reveal numerous closed depressions and mounds about 50 to 100 m in diameter and 2 to 3 m in relief, presumably related to grounded ice floes common in the area year round. Some of the features were obliterated over the course of three seasons while new ones formed. Although the depressions resemble kettles, they are formed by very different mechanisms. We believe that these bedforms represent erosion and deposition caused by: a) intensified flow around stationary ice floes serving as obstacles and b) pulsating currents generated by vertical oscillations or rocking motions of grounded floes in a seaway. Because sediment transport occurs around the ice, not where it directly touches the sea floor, the depressions are much larger than the base of the acting floes. Ice-wallow bedforms, although not found everywhere, are characteristic of arctic nearshore regions with non-cohesive sediments, and most likely occur in other ice-stressed coastal environments in differing degrees. The bedforms studied here are highly active and must be considered in planning nearshore construction activities.

Journal of Sedimentary Petrology↗

Stratiform chromitite at Campo Formoso, Bahia, Brazil

Serpentinized ultramafic rocks of pre-Minas age (Precambrian) crop out in an elongate belt 1 km wide and 15 km long on the west side of the Serra de Jacobina near Campo Formoso, Brazil. The original stratiform ultramafic body has been so metamorphosed, faulted, and eroded that only narrow dismembered segments occur along the thrust-faulted west side of the Serra de Jacobina. The most economically important segment is that near Campo Formoso where chromitite layers can be followed for about 8 km near the base of the intrusion. Three to five chromitite layers of minable thicknesses have been mapped, the thickest being 1.4 to 1.8 m. The aggregate thickness of chromitite is about 5 m and the average grade is about 43.5 percent Cr 2 0 3 . The Cr:Fe ratio ranges from 1.6:1 in weathered friable chromite to 3:1 in massive unweathered chromitite; in cleaned chromite concentrates the ratio is as much as 3.3:1. Estimates of indicated chromite ore reserves range from 16 to 20 million metric tons, of which only an estimated 6 million tons can profitably be mined. Factors preventing full exploitation of the chromite include (1) the local presence of as much as 7 m of overlying colluvium, (2) the extremely faulted character of the chromitite layers, (3) the presence of an overlying quartzite thrust plate very near the line of chromitite outcrops, and (4) the steep dips of the chromitite layers, 15° to 40° along the Pedrinhas-Limoeiro trends and 50° to 65° along the Cascabulhos trend.

Bahia↗

Growth trishear model and its application to the Gilbertown graben system, southwest Alabama

Fault-propagation folding associated with an upward propagating fault in the Gilbertown graben system is revealed by well-based 3-D subsurface mapping and dipmeter analysis. The fold is developed in the Selma chalk, which is an oil reservoir along the southern margin of the graben. Area-depth-strain analysis suggests that the Cretaceous strata were growth units, the Jurassic strata were pregrowth units, and the graben system is detached in the Louann Salt. The growth trishear model has been applied in this paper to study the evolution and kinematics of extensional fault-propagation folding. Models indicate that the propagation to slip (p/s) ratio of the underlying fault plays an important role in governing the geometry of the resulting extensional fault-propagation fold. With a greater p/s ratio, the fold is more localized in the vicinity of the propagating fault. The extensional fault-propagation fold in the Gilbertown graben is modeled by both a compactional and a non-compactional growth trishear model. Both models predict a similar geometry of the extensional fault-propagation fold. The trishear model with compaction best predicts the fold geometry. ?? 2008 Elsevier Ltd. All rights reserved.

Journal of Structural Geology↗

Fault-related clay authigenesis along the Moab Fault: Implications for calculations of fault rock composition and mechanical and hydrologic fault zone properties

The presence of clays in fault rocks influences both the mechanical and hydrologic properties of clay-bearing faults, and therefore it is critical to understand the origin of clays in fault rocks and their distributions is of great importance for defining fundamental properties of faults in the shallow crust. Field mapping shows that layers of clay gouge and shale smear are common along the Moab Fault, from exposures with throws ranging from 10 to ???1000 m. Elemental analyses of four locations along the Moab Fault show that fault rocks are enriched in clays at R191 and Bartlett Wash, but that this clay enrichment occurred at different times and was associated with different fluids. Fault rocks at Corral and Courthouse Canyons show little difference in elemental composition from adjacent protolith, suggesting that formation of fault rocks at those locations is governed by mechanical processes. Friction tests show that these authigenic clays result in fault zone weakening, and potentially influence the style of failure along the fault (seismogenic vs. aseismic) and potentially influence the amount of fluid loss associated with coseismic dilation. Scanning electron microscopy shows that authigenesis promotes that continuity of slip surfaces, thereby enhancing seal capacity. The occurrence of the authigenesis, and its influence on the sealing properties of faults, highlights the importance of determining the processes that control this phenomenon. ?? 2010 Elsevier Ltd.

Journal of Structural Geology↗

The wind-forced response on a buoyant coastal current: Observations of the western Gulf of Maine plume

The Freshwater plume in the western Gulf of Maine is being studied as part of an interdisciplinary investigation of the physical transport of a toxic alga. A field program was conducted in the springs of 1993 and 1994 to map the spatial and temporal patterns of salinity, currents and algal toxicity. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman drift dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to as little as 10 km. Using a simple slab model, we find qualitative agreement between the observed variations of plume width and those predicted by Ekman theory for short time scales of integration. Near surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate a dominant Ekman balance within the plume. A significant correlation between alongshore currents and winds suggests that interfacial drag may be important, although inclusion of a Raleigh drag term does not significantly improve the alongshore momentum balance.

Gulf of Maine↗

Rapid, remote assessment of Hurricane Matthew impacts using four-dimensional structure-from-motion photogrammetry

Timely assessment of coastal landforms and structures after storms is important for evaluating storm impacts, aiding emergency response and restoration, and initializing and assessing morphological models. Four-dimensional multiview photogrammetry, also known as structure from motion (4D SfM), provides a method for generating three-dimensional reconstructions of landscapes at two times (before and after events) using only photos and existing information for ground control points. Here, these techniques were applied using National Oceanic and Atmospheric Administration (NOAA)-obtained oblique aerial photos taken before (2015) and immediately after Hurricane Matthew (2016) to assess coastal changes near Matanzas, Florida. This work demonstrated that 3D digital elevation models can be constructed within 48 hours of postevent photo collection without on-site ground control measurements. One advantage of timely SfM elevation-change assessments is that they avoid confusion of storm impacts with changes that occur after the event but before LIDAR surveys can be performed. The accuracy and precision of the 4D SfM maps were assessed a posteriori using the first-available LIDAR data, which were collected more than a month after the hurricane, and 11 independent ground-truth survey points measured a week after the hurricane. Horizontal coordinates of the 4D SfM reconstruction were biased by an average of 0.79 m (0.83 m root-mean-square difference; RMSD) compared with the ground-truth points, but vertical elevations were more accurate. They were biased from the LIDAR by −0.09 to −0.25 m, with ∼0.20 m RMSD from both the LIDAR data and five ground-truth points with good vertical positioning and 0.25 m RMSD from LIDAR data along a 60-m stretch of pavement. This level of precision was sufficient to quantify geomorphological change that was often in excess of 1 m. The methodology is conducive for rapid assessment of changes along short stretches (tens of kilometers) of coast with modest resources and could be scaled up for larger regions.

Florida↗

Mapping habitat suitability rangewide for headwater stream-associated torrent salamanders

In the moist, coniferous forests of Oregon and Washington, a dense network of headwater streams supports high levels of amphibian species endemism and diversity, including the torrent salamanders ( Rhyacotriton spp.). The Columbia and Cascade torrent salamanders ( R. kezeri and R. cascadae , respectively) are species of greatest conservation need, but current information is insufficient for an adequate status assessment. We aimed to determine the current distributional extent, occurrence, and abundance of both species across their respective ranges, and to characterize watershed-scale habitat suitability. We conducted rangewide surveys for both species and used a novel, robust model development and refinement process to characterize habitat suitability based on climatic, topographic, forest-structure, hydrological, and water-balance drivers of their distributions. Results supported interspecific differences in rangewide occurrence and abundance patterns, as well as ecological associations that provide insights into potentially divergent patterns of resilience. Topographic, water-balance, and streamflow metrics were top predictors of Columbia torrent occurrences, whereas metrics of seasonal temperature and atmospheric moisture were top predictors of Cascade torrent occurrences. Our findings can inform regional conservation planning efforts to identify likely climate refugia and habitat-connectivity pathways for gene flow and watershed-scale ecological resilience.

Oregon, Washington↗