USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Interactions”

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.

498 records · Page 28Linked to original sources

Toggling of seismicity by the 1997 Kagoshima earthquake couplet: A demonstration of time-dependent stress transfer

Two M ??? 6 well-recorded strike-slip earthquakes struck just 4 km and 48 days apart in Kagoshima prefecture, Japan, in 1997, providing an opportunity to study earthquake interaction. Aftershocks are abundant where the Coulomb stress is calculated to have been increased by the first event, and they abruptly stop where the stress is dropped by the second event. This ability of the main shocks to toggle seismicity on and off argues that static stress changes play a major role in exciting aftershocks, whereas the dynamic Coulomb stresses, which should only promote seismicity, appear to play a secondary role. If true, the net stress changes from a sequence of earthquakes might be expected to govern the subsequent seismicity distribution. However, adding the stress changes from the two Kagoshima events does not fully capture the ensuing seismicity, such as its rate change, temporal decay, or migration away from the ends of the ruptures. We therefore implement a stress transfer model that incorporates rate/state friction, in which seismicity is treated as a sequence of independent nucleation events that are dependent on the fault slip, slip rate, and elapsed time since the last event. The model reproduces the temporal response of seismicity to successive stress changes, including toggling, decay, and aftershock migration. Nevertheless, the match of observed to predicted seismicity is quite imperfect, due perhaps to inadequate knowledge of several model parameters. However, to demonstrate the potential of this approach, we build a probabilistic forecast of larger earthquakes on the expected rate of small aftershocks, taking advantage of the large statistical sample the small shocks afford. Not surprisingly, such probabilities are highly time- and location-dependent: During the first decade after the main shocks, the seismicity rate and the chance of successive large shocks are about an order of magnitude higher than the background rate and are concentrated exclusively in the stress triggering zones. Copyright 2003 by the American Geophysical Uion.

Journal of Geophysical Research B: Solid Earth↗

Complex magmatic-tectonic interactions during the 2020 Makushin Volcano, Alaska, earthquake swarm

On June 15, 2020, at 21:16 UTC, a locally-felt earthquake of magnitude 4.2 struck Unalaska Island, Alaska, ∼15 km west of the town of Unalaska and the large fishing port of Dutch Harbor. The event was followed by a M4.1 earthquake at 00:34 UTC and several M3+ aftershocks, initiating a prolific sequence with hundreds of earthquakes recorded into late December. The earthquakes all locate about 12 km southeast of the summit of Makushin Volcano at 7 to 10 km depth. To date, no eruptive activity or other surface changes have been observed at the volcano in webcam images, GPS or InSAR. Seismic bursts close to volcanoes are often associated with the onset of unrest that can lead to eruption. However, determining whether seismicity reflects magmatic rather than tectonic stresses is often challenging, although critical for hazard assessments and risk management strategies. To investigate the triggering mechanisms of the recent Makushin seismicity, we integrate information from space-time patterns of the earthquake hypocenters with their fault-plane solutions. We relocate the swarm events using double-difference relocation techniques and a 3D velocity model and find that the earthquakes, although they seem to follow two predominant orientations (NW-SE and SW-NE), do not show clear clustering into preferred alignments. Similarly, we do not observe pronounced migration in time and space. Fault-plane solutions (FPS) for all but one M2.5+ earthquakes have P-axis orientations consistent with subhorizontal NW-SE oriented regional maximum compression, whereas many of the lower-magnitude earthquakes have P-axes perpendicular to regional maximum compression. This provides evidence for the presence of a local stress field likely induced by magma intrusion. Results from Coulomb stress modeling are also consistent with dike inflation modulated by stresses induced by the M4+ earthquakes. The seismic swarm is thus likely linked to a superposition of driving stresses from both magmatic and tectonic processes on pre-existing faults. The case of the 2020 Makushin swarm, with its unusual characteristics, challenges traditional swarm classification schemes and suggests that a reconsideration of the definition of seismic swarms as having the maximum magnitude event in the middle of the swarm is warranted.

Alaska↗

Hysteretic response of suspended-sediment in wildfire affected watersheds of the Pacific Northwest and Southern Rocky Mountains

Wildfires can have a profound impact on hydrosedimentary interactions, or the relationship between sediment and runoff, in forested headwater streams. Quantification of sediment-runoff dynamics at the event scale is integral for understanding source areas and transport of suspended-sediment through a watershed following wildfire. Here we used high-frequency turbidity and stream discharge data, coupled with discrete suspended-sediment measurements, in burned and unburned watersheds in the Southern Rocky Mountains and the western Cascades Range to evaluate the response of fine-grained (clay- and silt-sized particles) suspended-sediment. Hysteresis analysis was conducted on estimated suspended-sediment concentrations (using turbidity as a proxy) and streamflow through measurement of the difference in sediment concentration on the rising and falling limbs of the event hydrograph. All burned watersheds exhibited elevated fine suspended-sediment concentrations relative to concentrations found in pre-fire conditions. Changes to hysteretic response vary and may depend on a watershed's sediment connectivity limitations. Results suggest a watershed's inherent hillslope-to-channel (or lateral) connectivity is the primary factor controlling the relative magnitude of event-driven fine sediment fluxes in watersheds affected by wildfire. While wildfire did promote lateral connectivity through activation of hillslope sources, snowmelt, precipitation characteristics and antecedent conditions were more important drivers of hysteretic response than wildfire. For watersheds influenced by annual snowpack, we identified a predominantly clockwise hysteretic response during snowmelt and counterclockwise events during the late spring and summer months. There were also proportionally more counterclockwise events after wildfire in watersheds with high sediment connectivity. Results suggest contrasting wildfire-related sediment risk potential. Rivers in burned watersheds with high sediment connectivity may pose a higher risk to receiving waterbodies, such as larger tributaries or reservoirs, while rivers with low sediment connectivity may experience long-term sediment-related risk within the watershed above the outlet.

Oregon↗

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗

Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults

We argue that key features of thrust earthquake triggering, inhibition, and clustering can be explained by Coulomb stress changes, which we illustrate by a suite of representative models and by detailed examples. Whereas slip on surface-cutting thrust faults drops the stress in most of the adjacent crust, slip on blind thrust faults increases the stress on some nearby zones, particularly above the source fault. Blind thrusts can thus trigger slip on secondary faults at shallow depth and typically produce broadly distributed aftershocks. Short thrust ruptures are particularly efficient at triggering earthquakes of similar size on adjacent thrust faults. We calculate that during a progressive thrust sequence in central California the 1983 M w = 6.7 Coalinga earthquake brought the subsequent 1983 M w = 6.0 Nuñez and 1985 M w = 6.0 Kettleman Hills ruptures 10 bars and 1 bar closer to Coulomb failure. The idealized stress change calculations also reconcile the distribution of seismicity accompanying large subduction events, in agreement with findings of prior investigations. Subduction zone ruptures are calculated to promote normal faulting events in the outer rise and to promote thrust-faulting events on the periphery of the seismic rupture and its downdip extension. These features are evident in aftershocks of the 1957 M w = 9.1 Aleutian and other large subduction earthquakes. We further examine stress changes on the rupture surface imparted by the 1960 M w = 9.5 and 1995 M w = 8.1 Chile earthquakes, for which detailed slip models are available. Calculated Coulomb stress increases of 2–20 bars correspond closely to sites of aftershocks and postseismic slip, whereas aftershocks are absent where the stress drops by more than 10 bars. We also argue that slip on major strike-slip systems modulates the stress acting on nearby thrust and strike-slip faults. We calculate that the 1857 M w = 7.9 Fort Tejon earthquake on the San Andreas fault and subsequent interseismic slip brought the Coalinga fault ∼1 bar closer to failure but inhibited failure elsewhere on the Coast Ranges thrust faults. The 1857 earthquake also promoted failure on the White Wolf reverse fault by 8 bars, which ruptured in the 1952 M w = 7.3 Kern County shock but inhibited slip on the left-lateral Garlock fault, which has not ruptured since 1857. We thus contend that stress transfer exerts a control on the seismicity of thrust faults across a broad spectrum of spatial and temporal scales.

California↗

Assessing the potential hydrological impact of the Gibe III Dam on Lake Turkana water level using multi-source satellite data

Lake Turkana, the largest desert lake in the world, is fed by ungauged or poorly gauged river systems. To meet the demand of electricity in the East African region, Ethiopia is currently building the Gibe III hydroelectric dam on the Omo River, which supplies more than 80% of the inflows to Lake Turkana. On completion, the Gibe III dam will be the tallest dam in Africa with a height of 241 m. However, the nature of interactions and potential impacts of regulated inflows to Lake Turkana are not well understood due to its remote location and unavailability of reliable in-situ datasets. In this study, we used 12 years (1998–2009) of existing multi-source satellite and model-assimilated global weather data. We use calibrated multi-source satellite data-driven water balance model for Lake Turkana that takes into account model routed runoff, lake/reservoir evapotranspiration, direct rain on lakes/reservoirs and releases from the dam to compute lake water levels. The model evaluates the impact of Gibe III dam using three different approaches such as (a historical approach, a knowledge-based approach, and a nonparametric bootstrap resampling approach) to generate rainfall-runoff scenarios. All the approaches provided comparable and consistent results. Model results indicated that the hydrological impact of the dam on Lake Turkana would vary with the magnitude and distribution of rainfall post-dam commencement. On average, the reservoir would take up to 8–10 months, after commencement, to reach a minimum operation level of 201 m depth of water. During the dam filling period, the lake level would drop up to 2 m (95% confidence) compared to the lake level modelled without the dam. The lake level variability caused by regulated inflows after the dam commissioning were found to be within the natural variability of the lake of 4.8 m. Moreover, modelling results indicated that the hydrological impact of the Gibe III dam would depend on the initial lake level at the time of dam commencement. Areas along the Lake Turkana shoreline that are vulnerable to fluctuations in lake levels were also identified. This study demonstrates the effectiveness of using existing multi-source satellite data in a basic modeling framework to assess the potential hydrological impact of an upstream dam on a terminal downstream lake. The results obtained from this study could also be used to evaluate alternate dam-filling scenarios and assess the potential impact of the dam on Lake Turkana under different operational strategies.

Hydrology and Earth System Sciences Discussions↗

Stress interaction between subduction earthquakes and forearc strike-slip faults: Modeling and application to the northern Caribbean plate boundary

Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from "unclamping" of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high. Copyright 2004 by the American Geophysical Union.

Puerto Rico↗

Geologic and hydrologic characterization and evaluation of the Basin and Range Province relative to the disposal of high-level radioactive waste: Part I, Introduction and guidelines

The U.S. Geological Survey's program for geologic and hydrologic evaluation of physiographic provinces to identify areas potentially suitable for locating repository sites for disposal of high-level nuclear wastes was announced to the Governors of the eight States in the Basin and Range Province on May 5, 1981. Representatives of Arizona, California, Idaho, New Mexico, Nevada, Oregon, Texas, and Utah, were invited to cooperate with the Federal Government in the evaluation process. Each Governor was requested to nominate an Earth scientist to represent the State in a province working group composed of State and U.S. Geological Survey representatives. This report, Part I of a three-part report, provides the background, introduction and scope of the study. This part also includes a discussion of geologic and hydrologic guidelines that will be used in the evaluation process and illustrates geohydrologic environments and the effect of individual factors in providing multiple natural barriers to radionuclide migration. Part II is a reconnaissance characterization of the geologic and hydrologic factors to be used in the initial screening of the Basin and Range Province. Part III will be the initial evaluation of the Province and will identify regions that appear suitable for further study. The plan for study of the Province includes a stepwise screening process by which successively smaller land units are considered in increasing detail. Each step involves characterization of the geology and hydrology and selection of subunits for more intensive characterization. Selection of subunits for further study is by evaluation of geologic and hydrologic conditions following a set of guidelines. By representation on the Province Working Group, the States participate in a consultation and review role in: (1) Establishing geologic and hydrologic guidelines, and (2) characterizing and evaluating the Province. The States also participate in compilation of geologic and hydrologic data used in characterizing the Province. The current (1983) needs for a high-level radioactive waste repository include: (1) Disposal in a mined repository; (2) retrievability of the waste for as much as 50 years; and (3) confidence of isolation of the waste from the accessible environment. Isolation of the waste needs to be assured using geologic and hydrologic conditions that: (1) Minimize risk of inadvertent future intrusions by man; (2) minimize the possibility of disturbance by processes that would expose the waste or increase its mobility; and (3) provide a system of natural barriers to the migration of waste by ground water. The guidelines adopted by the Province Working Group are designed to provide a standard with which these conditions can be compared. The guidelines can be grouped into four principal categories: (1) Potential host media, (2) ground-water conditions, (3) tectonic conditions, and. (4) occurrence of natural resources. Ideally the host medium constitutes the first natural barrier to migration of radionculides. The host medium ideally should be a rock type that prevents or retards dissolution and transport of radionuclides. Rocks in both the saturated and unsaturated zones may have desirable characteristics for host media. Rocks-other than the host-in the ground-water flow path from the repository ideally should be major barriers to radionuclide migration. Confining beds of low permeability might be present to retard the rate of flow between more permeable beds. Additionally, sorption of radionuclides by materials such as clays and zeolites in the flow path can further retard the flow of radionuclides by several orders of magnitude. Tectonic conditions in an area should not present a probable cause for exhumation or increased mobility of radioactive waste. Natural resources are a factor for consideration because of the problem of future human intrusion and exposure to radioactivity in the quest for minerals, oil, gas, water, and geothermal resources. The ultimate evaluation of the suitability of a geohydrologic environment for developing a mined repository needs to assess all geologic and hydrologic characteristics and their interaction in providing confidence that a geohydrologic environment will effectively isolate radionuclides from human access. Several hypothetical settings with typical geohydrologic conditions in the Basin and Range Province are used to illustrate the effect of multiple barriers in the isolation of radionuclides.

Arizona, California, Idaho, Nevada, New Mexico, Or↗

Igneous rocks in the Fish Creek Mountains and environs, Battle Mountain area, north-central Nevada: A microcosm of Cenozoic igneous activity in the northern Great Basin, Basin and Range Province, USA

The Great Basin of the western United States, the northern component of the Basin and Range Province, is a region of Cenozoic lithospheric extension with multiple periods and types of igneous activity. The composition and volume of Cenozoic magmas reflect a complex interaction between mantle-derived magmas and highly diverse crust, where both mantle sources and magmatic processes were modulated by tectonic environment. The Fish Creek Mountains in north-central Nevada underwent multiple igneous events ranging from ca. 40 Ma to 1 Ma that span all of the complex Cenozoic tectono-magmatic episodes of the Great Basin. The Fish Creek Mountains, therefore, is an ideal location to evaluate the different sources and processes involved in magma generation. Many plutons were emplaced in the region between about 40 and 38 Ma, several of which host base and precious metal deposits. Between 36 and 33 Ma, lava fields and calderas of the 37–19 Ma Ignimbrite Flare-up were emplaced. Both these and the preceding plutons resulted from southwestward rollback of the Farallon plate beneath North America during by far the most voluminous phase of Cenozoic magmatism. The lavas range from rare basalt and basaltic andesite to andesite, dacite, and rhyolite, have continental arc-like incompatible element patterns, and high initial 87 Sr/ 86 Sr and low εNd that require a metasomatized lithospheric mantle source combined with minor crustal component. Ignimbrites of the 34.4 Ma Cove Mine (trachydacite to rhyolite) and 34.0 Ma Caetano calderas (rhyolite to high-silica rhyolite) are abundantly porphyritic, include hydrous phases, were largely derived from partial melts of crustal rocks, but likely include 20–30% of a mantle-derived component. Igneous activity ceased in the region as the rollback-arc migrated to the southwest, but at 24.9 Ma a new caldera formed in the southern Fish Creek Mountains that was filled by ignimbrites of the Fish Creek Mountains Tuff. Intracaldera rhyolite ignimbrites range from aphyric, pumice-rich deposits at the base to progressively more quartz-feldspar phyric ignimbrites at higher levels; all flow units lack hydrous phases. No contemporaneous mafic or intermediate igneous activity accompanied caldera formation, but initial 87 Sr/ 86 Sr values in the Fish Creek Mountains tuffs are lower than in the Caetano Tuff, suggesting a greater mantle contribution to the 24.9 Ma ignimbrites. After another hiatus in igneous activity, the region was intruded and overlain by basalt to rhyolite dykes and lavas of the northern Nevada rift between 16.8 and 15.1 Ma. The primarily tholeiitic igneous suite is of the same age, chemistry, and isotopic composition as the Grande Ronde Formation of the Columbia River flood basalts, and evolved members (trachydacite and rhyolite) are crustally contaminated. The youngest northern Nevada rift lava is an alkali olivine basalt with isotopic affinity to basalts of the eastern Snake River Plain. After 10 Ma of quiescence, the region was locally covered by mafic lava flows with high-alumina olivine tholeiite compositions, represented by the 5.4 Ma Pumpernickel Valley flows. Their mid-ocean ridge-like incompatible element compositions indicate a depleted mantle source for the lavas, but radiogenic isotopic compositions indicate that the lavas of this region include a significant contribution from a mafic to ultramafic, high- 87 Sr/ 86 Sr source. The final igneous event in the Fish Creek Mountains region, the 4.0 to 1.0 Ma Buffalo Valley volcanic field, includes flows and spatter cones of transitional to alkalic basalt that are divided into two geochemical groups with identical isotopic compositions. They represent variable, low percent partial melts of the asthenosphere at different depths, yielding different rare earth element characteristics. Similar to the Lunar Crater volcanic field, the Buffalo Valley rocks represent a rare case where the lithosphere in the central Great Basin is now thin enough to allow melting of the underlying asthenosphere. Cenozoic magmatism in the northern Great Basin exhibits several transitions in magma sources and tectonic setting with time. Magmatism began as pre-extension, subduction-related, primarily lithospherically-derived magmas emplaced on/in tectonically-thickened crust. The onset of extension was partially driven by impingement of the Yellowstone plume that resulted in emplacement of rift-related volcanic and intrusive rocks in the northern Nevada rift, followed by the eruption of extension-related HAOT lavas along the northwest margin of the Great Basin. Finally, lithospheric thinning allowed for partial melting of the asthenosphere and eruption of alkaline basaltic lavas.

Nevada↗

Coherence of Mach fronts during heterogeneous supershear earthquake rupture propagation: Simulations and comparison with observations

We study how heterogeneous rupture propagation affects the coherence of shear and Rayleigh Mach wavefronts radiated by supershear earthquakes. We address this question using numerical simulations of ruptures on a planar, vertical strike-slip fault embedded in a three-dimensional, homogeneous, linear elastic half-space. Ruptures propagate spontaneously in accordance with a linear slip-weakening friction law through both homogeneous and heterogeneous initial shear stress fields. In the 3-D homogeneous case, rupture fronts are curved owing to interactions with the free surface and the finite fault width; however, this curvature does not greatly diminish the coherence of Mach fronts relative to cases in which the rupture front is constrained to be straight, as studied by Dunham and Bhat (2008a). Introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear wave speed. Calculations of the Fourier amplitude spectra (FAS) of ground velocity time histories corroborate the kinematic results of Bizzarri and Spudich (2008a): (1) The ground motion of a supershear rupture is richer in high frequency with respect to a subshear one. (2) When a Mach pulse is present, its high frequency content overwhelms that arising from stress heterogeneity. Present numerical experiments indicate that a Mach pulse causes approximately an ω −1.7 high frequency falloff in the FAS of ground displacement. Moreover, within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources, our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the Mach fronts, ground motion at stations experiencing Mach pulses should be richer in high frequencies compared to stations without Mach pulses. In contrast to the foregoing theoretical results, we find no average elevation of 5%-damped absolute response spectral accelerations (SA) in the period band 0.05–0.4 s observed at stations that presumably experienced Mach pulses during the 1979 Imperial Valley, 1999 Kocaeli, and 2002 Denali Fault earthquakes compared to SA observed at non-Mach pulse stations in the same earthquakes. A 20% amplification of short period SA is seen only at a few of the Imperial Valley stations closest to the fault. This lack of elevated SA suggests that either Mach pulses in real earthquakes are even more incoherent that in our simulations or that Mach pulses are vulnerable to attenuation through nonlinear soil response. In any case, this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20% close to the fault to account for Mach pulses, provided that the existing data are adequately representative of ground motions from supershear earthquakes.

Journal of Geophysical Research B: Solid Earth↗

Paleozoic sedimentary rocks in the Red Dog Zn-Pb-Ag district and vicinity, western Brooks Range, Alaska: provenance, deposition, and metallogenic significance

Geochemical analyses of Paleozoic sedimentary rocks in the western Brooks Range reveal a complex evolutionary history for strata surrounding the large Zn-Pb-Ag deposits of the Red Dog district. Data for major elements, trace elements, and rare earth elements (REE) were obtained on 220 samples of unaltered and unmineralized siliciclastic rocks from the Upper Devonian-Lower Mississippian Endicott Group (Hunt Fork Shale, Noatak Sandstone, Kanayut Conglomerate, Kayak Shale), the overlying Carboniferous Lisburne Group (Kuna Formation, unnamed drowned shelf facies), and the Pennsylvanian-Permian Siksikpuk Formation. Major base metal sulfide deposits of the region are present only in the Kuna Formation, which in the Red Dog district comprises siliceous black shale and black chert, minor limestone (calcareous radiolarite), and sparse lithic turbidite and bedded siliceous rock. Gray and rare black shales of the Kayak Shale and common black shales of the Kuna Formation are anomalously low in iron (avg Fe/Ti = 6.25 and 6.34, respectively) relative to other Paleozoic shales in the region (9.58-10.6) and to average shales worldwide (10.1-10.5). In contrast, the bedded siliceous rocks contain appreciable hematite (avg Fe/Ti = 35.0) and high U/Ti and REE/Ti ratios that are interpreted to reflect low amounts of detrital material and a major Fe-rich eolian component. Geochemical data (e.g., MnO <0.01 wt %; avg Cr = 317 ppm), sizes of framboidal pyrite grains, and limited bioturbation suggest anoxic and denitrifying depositional conditions for most black shales of the Kuna Formation; low Mo/Ti ratios argue against euxinic (sulfate-reducing) conditions. Organic-rich black shales of the Kuna Formation with up to 8.4 wt percent C organic and gray to black shales of the Kayak Shale with up to 4.1 wt percent C organic typically have only sparse pyrite (<1 wt % S) and very low iron-limited S/C ratios (mostly <0.2). Immobile element plots (e.g., Th-Zr/10-Sc) suggest that source terranes for all of the formations were dominated by one or more felsic-rich continental arcs; a small proportion of recycled sediments is present locally. A minor mafic igneous component also occurs in several shales of the Kuna and Siksikpuk Formations. High average values for the chemical index of alteration [Al 2 O 3 /(Al 2 O 3 + CaO + Na 2 O + K 2 O)] &infin; 100 for shales of the Endicott Group (76.4-81.5) imply moderate to intense chemical weathering in source areas of these sediments. A lower average for black shales of the Kuna Formation (73.7) does not require such weathering. Textural and geochemical data record the effects of diagenetic and/or hydrothermal fluid flow in some of the Paleozoic rocks. Mobility of P, F, U, and light REE is documented in black shales of the Kuna Formation by phosphate replacements of carbonate clasts and of matrix material surrounding the clasts. A relatively low average Ce/Ce* value of 0.73 for P-poor black shales of the Kuna Formation ( < 0.05 wt % P 2 O 5 ) and a similar Ce/Ce* value of 0.78 for a siderite concretion in Kayak Shale suggest that these diagenetic fluids were oxidizing. Many shales of the Kuna Formation have high (K 2 O &infin; 100)/(K 2 O + Al 2 O 3 ) ratios of 21.0 to 23.0, which contrast with low ratios of generally <18.0 for shales of the underlying Endicott Group. The high ratios in shales of the Kuna Formation reflect preferential reaction of smectite to illite during the Jurassic-Cretaceous Brookian orogeny, owing to high silica activities in pore fluids that were generated by the dissolution of abundant biogenic silica. The distribution and composition of Paleozoic strata in the western Brooks Range may have played a fundamental role in Zn-Pb mineralization of the Red Dog district. In our model, deposition and early lithification of biogenic chert and bedded siliceous rocks in the upper part of the Kuna Formation served as a regional hydrologic seal, acting as a cap rock to heat and hydrothermal fluids during Late Mississippian base-metal mineralization. Equally important was the iron-poor composition of black shales of the Kuna Formation (i.e., low Fe/Ti ratios), which limited synsedimentary pyrite formation in precursor sediments, resulting in significant H 2 S production in pore waters through the interaction of aqueous sulfate with abundant organic matter. This H 2 S may have been critical to the subsurface deposition of the huge quantities of Zn and Pb in the district. On the basis of this model, we propose that low Fe/Ti and S/C ratios in black shale sequences are potential basin-scale exploration guides for giant sediment-hosted, stratiform Zn-Pb-Ag deposits.

Alaska↗

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↗