USGS Science⌕ Search

SEARCH · USGS Science

Results for “One Earth”

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.

1,430 records · Page 80Linked to original sources

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗

Long-term trends in a Dimictic Lake

The one-dimensional hydrodynamic ice model, DYRESM-WQ-I, was modified to simulate ice cover and thermal structure of dimictic Lake Mendota, Wisconsin, USA, over a continuous 104-year period (1911&ndash;2014). The model results were then used to examine the drivers of changes in ice cover and water temperature, focusing on the responses to shifts in air temperature, wind speed, and water clarity at multiyear timescales. Observations of the drivers include a change in the trend of warming air temperatures from 0.081 &deg;C per decade before 1981 to 0.334 &deg;C per decade thereafter, as well as a shift in mean wind speed from 4.44 m s &minus;1 before 1994 to 3.74 m s &minus;1 thereafter. Observations show that Lake Mendota has experienced significant changes in ice cover: later ice-on date(9.0 days later per century), earlier ice-off date (12.3 days per century), decreasing ice cover duration (21.3 days per century), while model simulations indicate a change in maximum ice thickness (12.7 cm decrease per century). Model simulations also show changes in the lake thermal regime of earlier stratification onset (12.3 days per century), later fall turnover (14.6 days per century), longer stratification duration (26.8 days per century), and decreasing summer hypolimnetic temperatures (&minus;1.4 &deg;C per century). Correlation analysis of lake variables and driving variables revealed ice cover variables, stratification onset, epilimnetic temperature, and hypolimnetic temperature were most closely correlated with air temperature, whereas freeze-over water temperature, hypolimnetic heating, and fall turnover date were more closely correlated with wind speed. Each lake variable (i.e., ice-on and ice-off dates, ice cover duration, maximum ice thickness, freeze-over water temperature, stratification onset, fall turnover date, stratification duration, epilimnion temperature, hypolimnion temperature, and hypolimnetic heating) was averaged for the three periods (1911&ndash;1980, 1981&ndash;1993, and 1994&ndash;2014) delineated by abrupt changes in air temperature and wind speed. Average summer hypolimnetic temperature and fall turnover date exhibit significant differences between the third period and the first two periods. Changes in ice cover (ice-on and ice-off dates, ice cover duration, and maximum ice thickness) exhibit an abrupt change after 1994, which was related in part to the warm El Ni&ntilde;o winter of 1997&ndash;1998. Under-ice water temperature, freeze-over water temperature, hypolimnetic temperature, fall turnover date, and stratification duration demonstrate a significant difference in the third period (1994&ndash;2014), when air temperature was warmest and wind speeds decreased rather abruptly. The trends in ice cover and water temperature demonstrate responses to both long-term and abrupt changes in meteorological conditions that can be complemented with numerical modeling to better understand how these variables will respond in a future climate.

Wisconsin↗

Observations at convergent margins concerning sediment subduction, subduction erosion, and the growth of continental crust

At ocean margins where two plates converge, the oceanic plate sinks or is subducted beneath an upper one topped by a layer of terrestrial crust. This crust is constructed of continental or island arc material. The subduction process either builds juvenile masses of terrestrial crust through arc volcanism or new areas of crust through the piling up of accretionary masses (prisms) of sedimentary deposits and fragments of thicker crustal bodies scraped off the subducting lower plate. At convergent margins, terrestrial material can also bypass the accretionary prism as a result of sediment subduction, and terrestrial matter can be removed from the upper plate by processes of subduction erosion. Sediment subduction occurs where sediment remains attached to the subducting oceanic plate and underthrusts the seaward position of the upper plate's resistive buttress (backstop) of consolidated sediment and rock. Sediment subduction occurs at two types of convergent margins: type 1 margins where accretionary prisms form and type 2 margins where little net accretion takes place. At type 2 margins (∼19,000 km in global length), effectively all incoming sediment is subducted beneath the massif of basement or framework rocks forming the landward trench slope. At accreting or type 1 margins, sediment subduction begins at the seaward position of an active buttress of consolidated accretionary material that accumulated in front of a starting or core buttress of framework rocks. Where small-to-medium-sized prisms have formed (∼16,300 km), approximately 20% of the incoming sediment is skimmed off a detachment surface or decollement and frontally accreted to the active buttress. The remaining 80% subducts beneath the buttress and may either underplate older parts of the frontal body or bypass the prism entirely and underthrust the leading edge of the margin's rock framework. At margins bordered by large prisms (∼8,200 km), roughly 70% of the incoming trench floor section is subducted beneath the frontal accretionary body and its active buttress. In rounded figures the contemporary rate of solid-volume sediment subduction at convergent ocean margins (∼43,500 km) is calculated to be 1.5 km³/yr. Correcting type 1 margins for high rates of terrigenous seafloor sedimentation during the past 30 m.y. or so sets the long-term rate of sediment subduction at 1.0 km³/yr. The bulk of the subducted material is derived directly or indirectly from continental denudation. Interstitial water currently expulsed from accreted and deeply subducted sediment and recycled to the ocean basins is estimated at 0.9 km³/yr. The thinning and truncation caused by subduction erosion of the margin's framework rock and overlying sedimentary deposits have been demonstrated at many convergent margins but only off northern Japan, central Peru, and northern Chile has sufficient information been collected to determine average or long-term rates, which range from 25 to 50 km³/m.y. per kilometer of margin. A conservative long-term rate applicable to many sectors of convergent margins is 30 km³/km/m.y. If applied to the length of type 2 margins, subduction erosion removes and transports approximately 0.6 km³/yr of upper plate material to greater depths. At various places, subduction erosion also affects sectors of type 1 margins bordered by small- to medium-sized accretionary prisms (for example, Japan and Peru), thus increasing the global rate by possibly 0.5 km³/yr to a total of 1.1 km³/yr. Little information is available to assess subduction erosion at margins bordered by large accretionary prisms. Mass balance calculations allow assessments to be made of the amount of subducted sediment that bypasses the prism and underthrusts the margin's rock framework. This subcrustally subducted sediment is estimated at 0.7 km³/yr. Combined with the range of terrestrial matter removed from the margin's rock framework by subduction erosion, the global volume of subcrustally subducted material is estimated to range from 1.3 to 1.8 km³/yr. Subcrustally subducted material is either returned to the terrestrial crust by arc-related igneous processes or crustal underplating or is lost from the crust by mantle absorption. Geochemical and isotopic data support the notion that upper mantle melting returns only a small percent of the subducted material to the terrestrial crust as arc igneous rocks. Limited areal exposures of terrestrial rocks metamorphosed at deep (>20–30 km) subcrustal pressures and temperatures imply that only a small fraction of subducted material is reattached via deep crustal underplating. Possibly, therefore much of the subducted terrestrial material is recycled to the mantle at a rate near 1.6 km³/yr, which is effectively equivalent to the commonly estimated rate at which the mantle adds juvenile igneous material to the Earth's layer of terrestrial rock.

Alaska↗

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↗

Wildland fire ash: Production, composition and eco-hydro-geomorphic effects

Fire transforms fuels (i.e. biomass, necromass, soil organic matter) into materials with different chemical and physical properties. One of these materials is ash, which is the particulate residue remaining or deposited on the ground that consists of mineral materials and charred organic components. The quantity and characteristics of ash produced during a wildland fire depend mainly on (1) the total burned fuel (i.e. fuel load), (2) fuel type and (3) its combustion completeness. For a given fuel load and type, a higher combustion completeness will reduce the ash organic carbon content, increasing the relative mineral content, and hence reducing total mass of ash produced. The homogeneity and thickness of the ash layer can vary substantially in space and time and reported average thicknesses range from close to 0 to 50 mm. Ash is a highly mobile material that, after its deposition, may be incorporated into the soil profile, redistributed or removed from a burned site within days or weeks by wind and water erosion to surface depressions, footslopes, streams, lakes, reservoirs and, potentially, into marine deposits. Research on the composition, properties and effects of ash on the burned ecosystem has been conducted on material collected in the field after wildland and prescribed fires as well as on material produced in the laboratory. At low combustion completeness (typically T < 450 °C), ash is organic-rich, with organic carbon as the main component. At high combustion completeness (T > 450 °C), most organic carbon is volatized and the remaining mineral ash has elevated pH when in solution. It is composed mainly of calcium, magnesium, sodium, potassium, silicon and phosphorous in the form of inorganic carbonates, whereas at T > 580 °C the most common forms are oxides. Ash produced under lower combustion completeness is usually darker, coarser, and less dense and has a higher saturated hydraulic conductivity than ash with higher combustion completeness, although physical reactions with CO 2 and when moistened produce further changes in ash characteristics. As a new material present after a wildland fire, ash can have profound effects on ecosystems. It affects biogeochemical cycles, including the C cycle, not only within the burned area, but also globally. Ash incorporated into the soil increases temporarily soil pH and nutrient pools and changes physical properties such as albedo, soil texture and hydraulic properties including water repellency. Ash modifies soil hydrologic behavior by creating a two-layer system: the soil and the ash layer, which can function in different ways depending on (1) ash depth and type, (2) soil type and (3) rainfall characteristics. Key parameters are the ash's water holding capacity, hydraulic conductivity and its potential to clog soil pores. Runoff from burned areas carries soluble nutrients contained in ash, which can lead to problems for potable water supplies. Ash deposition also stimulates soil microbial activity and vegetation growth. Further work is needed to (1) standardize methods for investigating ash and its effects on the ecosystem, (2) characterize ash properties for specific ecosystems and wildland fire types, (3) determine the effects of ash on human and ecosystem health, especially when transported by wind or water, (4) investigate ash's controls on water and soil losses at slope and catchment scales, (5) examine its role in the C cycle, and (6) study its redistribution and fate in the environment.

Earth-Science Reviews↗

Tectonic history of the Syria Planum province of Mars

We attribute most of the development of extensive fractures in the Tharsis region to discrete tectonic provinces within the region, rather than to Tharsis as a single entity. One of these provinces is in Syria Planum. Faults and collapse structures in the Syria Planum tectonic province on Mars are grouped into 13 sets based on relative age, areal distribution, and morphology. According to superposition and fault crosscutting relations and crater counts we designate six distinct episodes of tectonic activity in the following chronologic order: stage I is an early to late Noachian deformation forming mostly east-west structures (fault set IA); some large volcanoes also formed. Faults were produced possibly by flexural uplift. Also, arcuate, north trending grabens (set IB) indicate that faulting is transitional to the next stage. Stage II is late Noachian to early Hesperian radial faulting centered in Syria Planum, possibly due to isostatic uplift in late Noachian (set IIA) to early Hesperian (set IIB) time. Stage III is early to late Hesperian faulting tangential to Syria Planum that was related to local centers of uplift (sets III1-III3) on the periphery of Syria Planum. Stage IV is a late Hesperian graben formation that was circumferential to Syria Planum (set IVA), caused either by collapse associated with eruption of magma or by flexure of the lithosphere due to volcanic loading. In association with volcanism, minor faulting occurred, producing an oval pattern of faults in southwest Syria Planum (set IVBl) and a radial pattern south of the planum (set IVB2) that apparently rejuvenated buried stage II faults. Stage V is a late Hesperian to early Amazonian development of grabens and troughs of Noctis Labyrinthus and western Valles Marineris (set V) that was probably instigated by local uplift; exposure of groundwater or ground ice zones may have produced further collapse and trough enlargement. Stage VI is early Amazonian northwest trending faulting in Noctis Fossae (set VIl), perhaps due to Tharsis Monies-centered tectonism, and north northwest normal faulting along the eastern side of the Claritas rise (set VI2) that was due to tectonic subsidence. The duration of tectonic activity in the Syria Planum province was perhaps 2–3 b.y. Photoclinometric topographic profiles across 132 grabens and fault scarps show that Syria Planum grabens have widths (average of 2.5 km, and most range from 1 to 6 km) similar to lunar grabens, but the Martian grabens have slightly higher side walls (average about 132 m) and gentler wall slopes (average of 9° and range of 2°–25°) than lunar grabens (93 m high and 18° slopes). Scarp degradation on Mars has progressed through quakes, impact shaking, and dry slope processes; the lower slopes may be due to Mars’ higher gravity. Estimates of the amount of extension for individual grabens range from 20 to 350 m; most estimates of the thickness of the faulted layer range from 0.5 to 4.5 km (average is 1.5 km). This thickness range corresponds closely to the 0.8-to 3.6-km range in depth for pits, troughs, and canyons in Noctis Labyrinthus and along the walls of Valles Marineris. We propose that the predominant 1-to 1.5-km values obtained for both the thickness of the faulted layer and the depths of the pits, troughs, and theater heads of the canyons reflect the initial depth to the water table in this region, as governed by the depth to the base of ground ice. Maximum depths for these features may indicate lowered groundwater table depths and the base of ejecta material.

Journal of Geophysical Research Solid Earth↗

Modification of wave-cut and faulting-controlled landforms

From a casual observation that the form of degraded fault scarps resembles the error function, this investigation proceeds through an elementary diffusion equation representation of landform evolution to the application of the resulting equations to the modern topography of scarplike landforms. The morphologic observations can be analyzed either in the form of one or more cross-strike elevation profiles or in the form of the slope-offset plot, a point plot of maximum scarp slope versus scarp offset. Working with either or both of these data representations for nine geologic structures, which range in age from 3 to 400 ka B.P. and in offset from 1 to 50 m, we apply analytical solutions for the vertical initial value scarp, the vertical continuous offset scarp, and the finite slope, initial value scarp. The model calculations are intrinsically ambiguous, yielding as the final answer only the product κ t (in the case of the initial value problem) or the product κ A −1 (in the case of the repeated faulting problem); here t is the age of a single scarp-forming event, 2 A is the vertical slip rate, and κ is the “mass diffusivity.” A single profile across three sea cliffs along the Santa Cruz, California, coast is analyzed as three separate initial value problems. A reasonably constrained age for the sea cliff standing above the Highway 1 platform returns κ = 11 GKG (1 GKG = 1 m 2 /ka). With this κ, we can date the two older sea cliffs. In fact, we do the converse: age estimates for these two older sea cliffs based on a uniform rate of uplift both yield the same κ as for the lower sea cliff. We treat a single profile of the Raymond fault in Pasadena/San Marino in terms of the repeated faulting problem; for it the uplift rate of R. Crook and others yields κ = 16 GKG. The very substantial preexisting offset across the Raymond fault must have been buried/leveled some 230 ka B.P., when the modern topography began to form. Our analysis of the Lake Bonneville shoreline scarps reveals a dependence of κ t on 2a, suggestive of nonlinear modification processes. This appearance is treated with the finite slope initial value scarp model to determine κ=1.1 GKG for the Lake Bonneville shoreline scarps. The suggestion of M. N. Machette that approximately 100,000-year-old, meter-high scarps are “unobservable” in weakly consolidated alluvial terranes of the Basin and Range and Rio Grande Rift Valley provinces can be formulated as κ ≳ 1 GKG. The coincidence between this inequality and the Lake Bonneville shoreline κ is striking, and it suggests that the value of κ = 1 GKG may be generally applicable, as a good first approximation, to the modification of alluvial terranes within the semiarid regions of the western United States. The Lake Bonneville shoreline κ is the basis for dating four sets of fault scarps in west-central Utah. The Drum Mountains fault scarps can be modeled in several different circumstances, but the most likely interpretation is that these fault scarps formed as the result of a single episode of normal faulting 3.6 to 5.7 ka B.P. The younger age is associated with quite low initial slope angles (25°). The other three sets of fault scarps show no evidence for finite initial value slopes. Fault scarps along the eastern base of the Fish Springs Range are very young, 3 ka B.P. We estimate the age of fault scarps along the western flank of the Oquirrh Mountains to be 32 ka B.P., which meets the weak geologic constraint that they be older than the Lake Bonneville shoreline. Fault scarps along the northeastern margin of the Sheeprock Mountains are even older, 53 ka B.P. An intriguing consequence of our single-event analysis of these scarps is that an 11.5-m offset occurred in a single earthquake.

Journal of Geophysical Research Solid Earth↗

The potential influence of changing climate on the persistence of salmonids of the inland west

The Earth's climate warmed steadily during the 20th century, and mean annual air temperatures are estimated to have increased by 0.6°C (Intergovernmental Panel on Climate Change, 2007). Although many cycles of warming and cooling have occurred in the past, the most recent warming period is unique in its rate and magnitude of change (Siegenthaler and others, 2005) and in its association with anthropogenic emissions of greenhouse gases (Intergovernmental Panel on Climate Change , 2007). The climate in the western United States warmed in concert with the global trend but at an accelerated rate (+0.8°C during the 20th century; Saunders and others, 2008). The region could also prove especially sensitive to future changes because the relatively small human population is growing rapidly, as are demands on limited water supplies. Regional hydrological patterns are dominated by seasonal snow accumulation at upper elevations. Most of the region is relatively dry, and both terrestrial and aquatic ecosystems are strongly constrained b y water availability (Barnett and others, 2008; Brown and others, 2008). Stream environments are dynamic and climatically extreme, and salmonid fishes are the dominant elements of the native biodiversity (McPhail and Lindsey, 1986; Waples and others, 2008). Salmonids have broad economic and ecologic importance, but a century of intensive water resource development, nonnative fish stocking, and land use has significantly reduced many populations and several taxa are now protected under the Endangered Species Act (Thurow and others, 1997; Trotter, 2008). Because salmonids require relatively pristine, cold water environments and are often isolated in headwater habitats, members of this group may be especially vulnerable to the effects of a warming climate (Keleher and Rahel, 1996; Rieman and others, 2007; Williams and others, 2009). Warming during the 20th century drove a series of environmental trends that have profound implications for many aspects of salmonid habitat, including disturbance regimes such as wildfire, and unfavorable changes to thermal and hydrologic properties of aquatic systems. Warmer air temperatures have been associated with decreased winter snow accumulations, have accelerated snowmelt, and have advanced the timing of peak runoff by several days to weeks across most of western North America (Stewart and others, 2005; Barnett and others, 2008). Less snow and earlier runoff decrease aquifer recharge, make less water available for groundwater inputs to streams, and are contributing to widespread decreases in summer low flows (Stewart and others, 2005; Rood and others, 2008; Luce and Holden 2009). Interannual variability in stream flow is increasing, as is the persistence of multi-year extreme conditions (McCabe and others, 2004; Pagano and Garen 2005). In many areas of western North America, flood risks have increased in association with warmer temperatures during the 20th century (Hamlet and Lettenmaier, 2005). Streams where midwinter temperatures are near freezing have proven especially sensitive to increased flooding because of associated transitional hydrological patterns (mixtures of rainfall and snowmelt) and propensity for occasional rain-on-snow events to rapidly melt winter snowpack and generate large floods (Hamlet and Lettenmaier, 2005). Stream temperatures in many areas are increasing (Peterson and Kitchell, 2001; Morrison and others, 2002; Bartholow, 2005; Kaushal and others, 2010), due to both air temperature increases and reduced summer flows that make streams more sensitive to warmer air temperatures (Isaak and others, 2010). In recent decades, wildfires have become more common across much of the western United States during periods of more frequent droughts (Westerling and others, 2006; Hoerling and Eischeid, 2007), and local stream temperature can increase in postfire environments (Gresswell, 1999; Dunham and others, 2007). Fire-related temperature increase within streams is commonly a transient phenomenon, lasting only until riparian vegetation has recovered (Gresswell, 1999); however, ongoing climate change could preclude recovery to higher stature, prefire vegetation types in some areas (McKenzie and others, 2004; van Mantgem and Stephenson, 2007), resulting in a loss of critical riparian shading. Additionally, when wildfires occur in steep mountain topographies, the vegetation that stabilize s soils on hillslopes is often killed and landslides become more prevalent (Gresswell, 1999). Landslides int o stream channels form debris flows composed of sediment slurries and dead trees that can scour channels to bedrock and further exacerbate stream heating, delay recovery of riparian areas, or extirpate fish populations (Gresswell, 1999; May and Gresswell, 2003; Dunham and others, 2007). Changes in stream environments will shift habitat distributions, sometimes unpredictably, in both time and space for many salmonid fishes. Water temperature fundamentally influences aquatic ecosystem health because distribution, reproduction, fitness, and survival of ectothermic organisms are inextricably linked to the thermal regime of the environment. Historically, research has focused on defining lethal thermal limits of salmonids (Eaton and others, 1995; Selong and others, 2001; Todd and others, 2008); however, water temperature is known to be important in biological processes at a variety of spatial scales and levels of biological organization (Rahel and Olden, 2008; McCullough and others, 2009). For instance, trout are affected directly by water temperature through feeding, metabolism, and growth rates, and indirectly by factors such as prey availability and species interactions (Wehrly and others, 2007; Rahel and Olden, 2008). Where cold water temperatures currently limit habitat suitability and distributions of some species (for example, at the highest and most northerly distributional extents; Nakano and others, 1996; Coleman and Fausch, 2007), a warming climate may gradually increase the quality and extent of suitable habitat. Over time, previously constrained populations are expected to expand into these new habitats and increase in number. Some evidence suggests this may already be happening in Alaska, where streams in recently deglaciated areas are being colonized by emigrants from nearby salmon and char populations (Milner and others, 2000). Unfortunately, many of the sensitive salmonid species that are often the focus of western managers are unlikely to benefit from future water temperature increases. Warmer stream temperatures will facilitate invasion by nonnative species that are broadly established in downstream areas into upstream areas where they will compete with native species (Rieman and others, 2006; Rahel and Olden, 2008; Fausch and others, 2009). In other cases, warmer stream temperatures will render thermally suitable habitats unsuitable in downstream areas and effect net losses of habitat because upstream distributions are often constrained by streams that are too small or steep (Hari and others, 2006; Isaak and others, 2010). Both scenarios are realistic for fish species like bull trout ( Salvelinus confluentus ) (Rieman and others, 2006; Rieman and others, 2007), the various subspecies of cutthroat trout ( Oncorhynchus clarkii ) (Williams and others, 2009), Gila trout ( Oncorhynchus gilae gilae ) (Kennedy and others, 2008), and Apache trout ( Oncorhynchus gilae apache ) (Rinne and Minckley, 1985; Carmichael and others, 1993). As native species are increasingly confined to smaller and more isolated habitats by a gradually warming climate, the effects of wildfires (whether related to lethal changes in water quality during a fire, channel debris flows, or chronic postfire warming ) could have greater proportional effects on remaining habitats (for example, Brown and others, 2001; Rieman and others, 2007). If these changes were accompanied by additional hydrologic alterations associated with changes to the magnitude, frequency, duration, timing, and rate of change of discharge patterns (Jager and others, 1999; Henderson and others, 2000), populations may begin to lose some of their historic resilience and become ever more susceptible to local extirpations. As dramatic and extensive as climatic and environmental trends are for salmonid habitats, global climate models (GCMs) project that many of these trends will continue and even accelerate until at least the middle of the 21st century (Intergovernmental Panel on Climate Change, 2007). Current projections suggest mean annual air temperatures will increase by an additional 1–3°C, and early indications are that climate trajectory is at the higher end of this range (Pittock, 2006; Raupach and others, 2007). Although predicted changes vary considerably, even the most conservative estimates suggest a warming rate that will be twice that observed during the 20th century. Projections for the midcentury are most certainly due to the effects of greenhouse gases already emitted or predicted in the short term, uncertainties of the effects of longer-term greenhouse gas emissions, short-term climate cycles, and process errors associated with climate models (Cox and Stephenson, 2007). Projections of changes in total precipitation are less certain than those for air temperatures, but most GCMs project relatively small changes in the Northwest, with the exception of slightly drier summer periods (Mote and others, 2008; Karl and others, 2009). In the Southwest, however, significant decreases (such as 15–30 percent ) are projected during most periods of the year, and this area is one of the few for which Intergovernmental Panel on Climate Change (2007) precipitation projections have a high level of certainty (Hoerling and Eischeid, 2007; Karl and others, 2009). Clearly, managers of native salmonids in the wester n United States should consider adjusting management strategies to accommodate a warmer and possibly drier future (Williams and others, 2009). Tools are needed to forecast where important changes may occur and how conservation efforts should be prioritized. In this Open-File Report, we document our initial efforts in this regard for 10 species and subspecies of inland trout and Montana Arctic grayling ( Thymallus arcticus ) across the western United States.

Arizona, California, Colorado, Idaho, Montana, Nev↗