USGS Science⌕ Search

SEARCH · USGS Science

Results for “Energy”

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 1,711 records · Page 95Linked to original sources
USGS2003

Strong ground-motion prediction from Stochastic-dynamic source models

In the absence of sufficient data in the very near source, predictions of the intensity and variability of ground motions from future large earthquakes depend strongly on our ability to develop realistic models of the earthquake source. In this article we simulate near-fault strong ground motion using dynamic source models. We use a boundary integral method to simulate dynamic rupture of earthquakes by specifying dynamic source parameters (fracture energy and stress drop) as spatial random fields. We choose these quantities such that they are consistent with the statistical properties of slip heterogeneity found in finite-source models of past earthquakes. From these rupture models we compute theoretical strong-motion seismograms up to a frequency of 2 Hz for several realizations of a scenario strike-slip Mw 7.0 earthquake and compare empirical response spectra, spectra obtained from our dynamic models, and spectra determined from corresponding kinematic simulations. We find that spatial and temporal variations in slip, slip rise time, and rupture propagation consistent with dynamic rupture models exert a strong influence on near-source ground motion. Our results lead to a feasible approach to specify the variability in the rupture time distribution in kinematic models through a generalization of Andrews' (1976) result relating rupture speed to apparent fracture energy, stress drop, and crack length to 3D dynamic models. This suggests that a simplified representation of dynamic rupture may be obtained to approximate the effects of dynamic rupture without having to do full dynamic simulations.

Bulletin of the Seismological Society of America↗
USGS2004

Observed and simulated ground motions in the San Bernardino basin region for the Hector Mine, California, earthquake

During the MW 7.1 Hector Mine earthquake, peak ground velocities recorded at sites in the central San Bernardino basin region were up to 2 times larger and had significantly longer durations of strong shaking than sites just outside the basin. To better understand the effects of 3D structure on the long-period ground-motion response in this region, we have performed finite-difference simulations for this earthquake. The simulations are numerically accurate for periods of 2 sec and longer and incorporate the detailed spatial and temporal heterogeneity of source rupture, as well as complex 3D basin structure. Here, we analyze three models of the San Bernardino basin: model A (with structural constraints from gravity and seismic reflection data), model F (water well and seismic refraction data), and the Southern California Earthquake Center version 3 model (hydrologic and seismic refraction data). Models A and F are characterized by a gradual increase in sediment thickness toward the south with an abrupt step-up in the basement surface across the San Jacinto fault. The basin structure in the SCEC version 3 model has a nearly uniform sediment thickness of 1 km with little basement topography along the San Jacinto fault. In models A and F, we impose a layered velocity structure within the sediments based on the seismic refraction data and an assumed depth-dependent Vp/Vs ratio. Sediment velocities within the SCEC version 3 model are given by a smoothly varying rule-based function that is calibrated to the seismic refraction measurements. Due to computational limitations, the minimum shear-wave velocity is fixed at 600 m/sec in all of the models. Ground-motion simulations for both models A and F provide a reasonably good match to the amplitude and waveform characteristics of the recorded motions. In these models, surface waves are generated as energy enters the basin through the gradually sloping northern margin. Due to the basement step along the San Jacinto fault, the surface wave energy is confined to the region north of this structure, consistent with the observations. The SCEC version 3 model, lacking the basin geometry complexity present in the other two models, fails to provide a satisfactory match to the characteristics of the observed motions. Our study demonstrates the importance of using detailed and accurate basin geometry for predicting ground motions and also highlights the utility of integrating geological, geophysical, and seismological observations in the development and validation of 3D velocity models.

California↗
USGS2004

Precarious rock and overturned transformer evidence for ground shaking in the Ms 7.7 Kern County earthquake: An analog for disastrous shaking from a major thrust fault in the Los Angeles basin

Precariously balanced rocks and overturned transformers in the vicinity of the White Wolf fault provide constraints on ground motion during the 1952 Ms 7.7 Kern County earthquake, a possible analog for an anticipated large earthquake in the Los Angeles basin (Shaw et al., 2002; Dolan et al., 2003). On the northeast part of the fault preliminary estimates of ground motion on the footwall give peak accelerations considerably lower than predicted by standard regression curves. On the other hand, on the hanging-wall, there is evidence of intense ground shattering and lack of precarious rocks, consistent with the intense hanging-wall accelerations suggested by foam-rubber modeling, numerical modeling, and observations from previous thrust fault earthquakes. There is clear evidence of the effects of rupture directivity in ground motions on the hanging-wall side of the fault (from both precarious rocks and numerical simulations). On the southwest part of the fault, which is covered by sediments, the thrust fault did not reach the surface ("blind" thrust). Overturned and damaged transformers indicate significant transfer of energy from the hanging wall to the footwall, an effect that may not be as effective when the rupture reaches the surface (is not "blind"). Transformers near the up-dip projection of the fault tip have been damaged or overturned on both the hanging-wall and footwall sides of the fault. The transfer of energy is confirmed in a numerical lattice model and could play an important role in a similar situation in Los Angeles. We suggest that the results of this study can provide important information for estimating the effects of a large thrust fault rupture in the Los Angeles basin, specially given the fact that there is so little instrumental data from large thrust fault earthquakes.

Bulletin of the Seismological Society of America↗
USGS2008

Simulated ground motion in Santa Clara Valley, California, and vicinity from M≥6.7 scenario earthquakes

Models of the Santa Clara Valley (SCV) 3D velocity structure and 3D finite-difference software are used to predict ground motions from scenario earthquakes on the San Andreas (SAF), Monte Vista/Shannon, South Hayward, and Calaveras faults. Twenty different scenario ruptures are considered that explore different source models with alternative hypocenters, fault dimensions, and rupture velocities and three different velocity models. Ground motion from the full wave field up to 1 Hz is exhibited as maps of peak horizontal velocity and pseudospectral acceleration at periods of 1, 3, and 5 sec. Basin edge effects and amplification in sedimentary basins of the SCV are observed that exhibit effects from shallow sediments with relatively low shear-wave velocity (330 m/sec). Scenario earthquakes have been simulated for events with the following magnitudes: (1) M 6.8–7.4 Calaveras sources, (2) M 6.7–6.9 South Hayward sources, (3) M 6.7 Monte Vista/Shannon sources, and (4) M 7.1–7.2 Peninsula segment of the SAF sources. Ground motions are strongly influenced by source parameters such as rupture velocity, rise time, maximum depth of rupture, hypocenter, and source directivity. Cenozoic basins also exert a strong influence on ground motion. For example, the Evergreen Basin on the northeastern side of the SCV is especially responsive to 3–5-sec energy from most scenario earthquakes. The Cupertino Basin on the southwestern edge of the SCV tends to be highly excited by many Peninsula and Monte Vista fault scenarios. Sites over the interior of the Evergreen Basin can have long-duration coda that reflect the trapping of seismic energy within this basin. Plausible scenarios produce predominantly 5-sec wave trains with greater than 30 cm/sec sustained ground-motion amplitude with greater than 30 sec duration within the Evergreen Basin.

California↗
USGS2016

Characterizing potentially induced earthquake rate changes in the Brawley Seismic Zone, southern California

The Brawley seismic zone (BSZ), in the Salton trough of southern California, has a history of earthquake swarms and geothermal energy exploitation. Some earthquake rate changes may have been induced by fluid extraction and injection activity at local geothermal fields, particularly at the North Brawley Geothermal Field (NBGF) and at the Salton Sea Geothermal Field (SSGF). We explore this issue by examining earthquake rate changes and interevent distance distributions in these fields. In Oklahoma and Arkansas, where considerable wastewater injection occurs, increases in background seismicity rate and aftershock productivity and decreases in interevent distance were indicative of fluid‐injection‐induced seismicity. Here, we test if similar changes occur that may be associated with fluid injection and extraction in geothermal areas. We use stochastic epidemic‐type aftershock sequence models to detect changes in the underlying seismogenic processes, shown by statistically significant changes in the model parameters. The most robust model changes in the SSGF roughly occur when large changes in net fluid production occur, but a similar correlation is not seen in the NBGF. Also, although both background seismicity rate and aftershock productivity increased for fluid‐injection‐induced earthquake rate changes in Oklahoma and Arkansas, the background rate increases significantly in the BSZ only, roughly corresponding with net fluid production rate increases. Moreover, in both fields the interevent spacing does not change significantly during active energy projects. This suggests that, although geothermal field activities in a tectonically active region may not significantly change the physics of earthquake interactions, earthquake rates may still be driven by fluid injection or extraction rates, particularly in the SSGF.

California↗
USGS2016

Seismic envelope-based detection and location of ground-coupled airwaves from volcanoes in Alaska

Volcanic explosions and other infrasonic sources frequently produce acoustic waves that are recorded by seismometers. Here we explore multiple techniques to detect, locate, and characterize ground‐coupled airwaves (GCA) on volcano seismic networks in Alaska. GCA waveforms are typically incoherent between stations, thus we use envelope‐based techniques in our analyses. For distant sources and planar waves, we use f ‐ k beamforming to estimate back azimuth and trace velocity parameters. For spherical waves originating within the network, we use two related time difference of arrival (TDOA) methods to detect and localize the source. We investigate a modified envelope function to enhance the signal‐to‐noise ratio and emphasize both high energies and energy contrasts within a spectrogram. We apply these methods to recent eruptions from Cleveland, Veniaminof, and Pavlof Volcanoes, Alaska. Array processing of GCA from Cleveland Volcano on 4 May 2013 produces robust detection and wave characterization. Our modified envelopes substantially improve the short‐term average/long‐term average ratios, enhancing explosion detection. We detect GCA within both the Veniaminof and Pavlof networks from the 2007 and 2013–2014 activity, indicating repeated volcanic explosions. Event clustering and forward modeling suggests that high‐resolution localization is possible for GCA on typical volcano seismic networks. These results indicate that GCA can be used to help detect, locate, characterize, and monitor volcanic eruptions, particularly in difficult‐to‐monitor regions. We have implemented these GCA detection algorithms into our operational volcano‐monitoring algorithms at the Alaska Volcano Observatory.

Alaska↗
USGS2022

High‐frequency rupture processes of the 2014 Mw 8.2 Iquique and 2015 Mw 8.3 Illapel, Chile, earthquakes determined from strong‐motion recordings

Strong‐motion recordings of the 2014 M w "> M w Mw 8.2 Iquique and 2015 M w "> M w Mw 8.3 Illapel, Chile, earthquakes were analyzed to determine rupture propagation and the location, timing, and strength of subevents that produce most of the high‐frequency (≥1 Hz) ground motions. A moving window,cross‐correlation analysis of recordings from a local dense array, band‐pass filtered at 1 Hz, directly shows that the Iquique earthquake ruptured to the southeast over a distance of about 60 km. Array analysis of lower frequency energy (0.03–0.1 Hz) indicates that it occurred updip of the high‐frequency rupture. A methodology was developed for inverting the envelopes of acceleration records (1–5 Hz) to map high‐frequency source factors on the rupture zone and was applied to the two earthquakes. Waveforms of M w "> M w Mw 6 earthquakes were used as empirical Green’s functions in the inversions. High‐frequency subevents within the two M w "> M w Mw 8 earthquakes were located at depths ranging from 25 to 55 km and mostly occurred downdip of the peak slip of these earthquakes. Fourier spectral ratios of the Iquique mainshock with respect to M w "> M w Mw 5–6 aftershocks were fit to determine their stress drops. The stress drops were roughly constant from M w "> M w Mw 5 to 8 at 10–20 MPa. A compound rupture model is described in which subevents occur in areas of spatially heterogeneous strength and stress on the rupture, and produce the high‐frequency radiated energy of the overall earthquake, but are not located in the areas of peak slip. The stress drop of the overall earthquake is shown to equal the root mean square stress drop of subevents averaged over the rupture area.

Iquique, Illapel↗
USGS2024

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California↗
USGS1977

Reflection and refraction of type-II S waves in elastic and anelastic media

The general theory of viscoelasticity, which accounts for elastic as well as anelastic linear behavior of materials, predicts that two types of S waves propagate in anelastic earth materials. The particle motion for an inhomogeneous plane S wave of type I is elliptical in the plane defined by the directions of propagation and attenuation, while the particle motion for an inhomogeneous plane S wave of type II is linear perpendicular to this plane. The general theory predicts that an S -wave incident upon a plane boundary perpendicular to the plane defined by the directions of propagation and attenuation generates S waves only of the same type. General characteristics of the type-II S waves reflected and refracted at plane anelastic boundaries are: (a) velocities and maximum attenuations which depend on the angle of incidence and frequency, (b) maximum energy flow at a different velocity and in a different direction than phase propagation, (c) energy flow across the boundary due to interaction of the incident and reflected waves. The general theory predicts these characteristics for the waves whenever a plane type-II S wave interacts with a plane anelastic boundary such as a soil-bedrock, crust-mantle, or core-mantle interface. None of these characteristics are predicted for the plane SH waves described by elasticity theory.

Bulletin of Seismological Society of America↗
USGS1990

The 1989 earthquake swarm beneath Mammoth Mountain, California: An initial look at the 4 May through 30 September activity

Mammoth Mountain is a 50,000- to 200,000-year-old cumulovolcano standing on the southwestern rim of Long Valley in eastern California. On 4 May 1989, two M = 1 earthquakes beneath the south flank of the mountain marked the onset of a swarm that has continued for more than 6 months. In addition to its longevity, noteworthy aspects of this persistent swarm include (1) an exponential-like increase in the rate of activity through the first month; (2) a vertically oriented, planar distribution of hypocenters at depths between 6 and 9 km with a north-northeast strike (roughly perpendicular to the average T-axis orientation for the swarm earthquakes); (3) recurring spasmodic bursts (rapid-fire sequences of similar-sized earthquakes with overlapping coda) and occasional earthquakes with enhanced low-frequency energy; (4) a uniform temporal distribution of the four largest (M ≈ 3) events over the first 4 months of the swarm with a cumulative seismic moment for the entire sequence through 30 September corresponding to a single M ≈ 4 earthquake; (5) a b-value of 1.2; and (6) submicrostrain perturbations on the nearby borehole dilatometer, the first of which led the onset of swarm activity by more than 2 weeks. These aspects of the swarm, together with its location along the southern extension of the youthful Mono-Inyo volcanic chain, which last erupted 500 to 600 years ago, point to a magmatic source for the modest but persistent influx of strain energy into the crust beneath Mammoth Mountain.

Bulletin of the Seismological Society of America↗
USGS1992

On the characteristics of local geology and their influence on ground motions generated by the Loma Prieta earthquake in the San Francisco Bay region, California

Strong ground motions recorded at 34 sites in the San Francisco Bay region from the Loma Prieta earthquake show marked variations in characteristics dependent on crustal structure and local geological conditions. Peak horizontal acceleration and velocity inferred for sites underlain by “rock” generally occur on the transverse component of motion. They are consistently greater with lower attenuation rates than the corresponding mean value predicted by empirical curves based on previous strong-motion data. Theoretical amplitude distributions and synthetic seismograms calculated for 10-layer models suggest that “bedrock” motions were elevated due in part to the wide-angle reflection of S energy from the base of a relatively thin (25 km) continental crust in the region. Characteristics of geologic and geotechnical units as currently mapped for the San Francisco Bay region show that average ratios of peak horizontal acceleration, velocity and displacement increase with decreasing mean shear-wave velocity. Ratios of peak acceleration for sites on “soil” (alluvium, fill/Bay mud) are statistically larger than those for sites on “hard rock” (sandstone, shale, Franciscan Complex). Spectral ratios establish the existence of predominant site periods with peak amplifications near 15 for potentially damaging levels of ground motion at some sites underlain by alluvium and fill/bay mud. Average spectral amplifications inferred for vertical and the mean horizontal motion are, respectively, (1,1) for sites on the Franciscan Complex (KJf), (1.4, 1.5) for sites on Mesozoic and Tertiary rocks (TMzs), (2.1, 2.0) for sites on the Santa Clara Formation (QTs), (2.3, 2.9) for sites on alluvium (Qal), and (2.1, 4.0) for sites on fill/Bay mud (Qaf/Qhbm). These mean values are not statistically different at the 5% significance level from those inferred from previous low-strain data. Analyses suggest that soil amplification and reflected crustal shear energy were major contributors to levels of ground motion sufficient to cause damage to vulnerable structures at distances near 100 km in the cities of San Francisco and Oakland.

California↗
USGS2008

Multidecadal climate-induced variability in microseisms

Microseisms are the most ubiquitous continuous seismic signals on Earth at periods between approximately 5 and 25 s (Peterson 1993; Kedar and Webb 2005). They arise from atmospheric energy converted to (primarily) Rayleigh waves via the intermediary of wind-driven oceanic swell and occupy a period band that is uninfluenced by common anthropogenic and wind-coupled noise processes on land (Wilson et al. 2002; de la Torre et al. 2005). "Primary" microseisms (near 8-s period) are generated in shallow water by breaking waves near the shore and/or the nonlinear interaction of the ocean wave pressure signal with the sloping sea floor (Hasselmann 1963). Secondary microseisms occur at half of the primary period and are especially strongly radiated in source regions where opposing wave components interfere (Longuett-Higgins 1950; Tanimoto 2007), which principally occurs due to the interaction of incident swell and reflected/scattered wave energy from coasts (Bromirski and Duennebier 2002; Bromirski, Duennebier, and Stephen 2005). Coastal regions having a narrow shelf with irregular and rocky coastlines are known to be especially efficient at radiating secondary microseisms (Bromirski, Duennebier, and Stephen 2005; Shulte-Pelkum et al. 2004). The secondary microseism is globally dominant, and its amplitudes proportional to the square of the standing wave height (Longuett-Higgins 1950), which amplifies its sensitivity to large swell events (Astiz and Creager 1994; Webb 2006).

Seismological Research Letters↗
USGS2002

How often do fishes "run on empty"?

We used a large data set of African, Neotropical, and North American fishes to examine the frequency with which fishes have empty stomachs ( n species = 254; n individuals = 36 875). Mean percentage of empty stomachs was low across all fishes (16.2 ± 1.2%) but varied from 0% to 79.4% among individual species. Nocturnal fishes had empty stomachs more frequently than diurnal fishes. Trophic classification was strongly associated with the percentage of empty stomachs, a pattern also revealed from an intraspecific analysis. Fishes appear to adjust their feeding intervals relative to the energy density, conversion efficiency, and particle size of their food. Piscivorous fishes seem to be the only trophic group that regularly experience long periods of empty stomachs, with species that consume prey whole and those that provide extended parental care having the highest proportions of empty stomachs. Activity patterns and life histories of some piscivorous species probably have evolved in partial response to energetic benefits of large, energy-rich food resources.

Ecology↗
USGS2008

Thermal and energetic constraints on ectotherm abundance: A global test using lizards

Population densities of birds and mammals have been shown to decrease with body mass at approximately the same rate as metabolic rates increase, indicating that energetic needs constrain endotherm population densities. In ectotherms, the exponential increase of metabolic rate with body temperature suggests that environmental temperature may additionally constrain population densities. Here we test simple bioenergetic models for an ecologically important group of ectothermic vertebrates by examining 483 lizard populations. We find that lizard population densities decrease as a power law of body mass with a slope approximately inverse to the slope of the relationship between metabolic rates and body mass. Energy availability should limit population densities. As predicted, environmental productivity has a positive effect on lizard density, strengthening the relationship between lizard density and body mass. In contrast, the effect of environmental temperature is at most weak due to behavioral thermoregulation, thermal evolution, or the temperature dependence of ectotherm performance. Our results provide initial insights into how energy needs and availability differentially constrain ectotherm and endotherm density across broad spatial scales. ?? 2008 by the Ecological Society of America.

Ecology↗
USGS2009

Changes in vegetation in northern Alaska under scenarios of climate change, 2003-2100: Implications for climate feedbacks

Assessing potential future changes in arctic and boreal plant species productivity, ecosystem composition, and canopy complexity is essential for understanding environmental responses under expected altered climate forcing. We examined potential changes in the dominant plant functional types (PFTs) of the sedge tundra, shrub tundra, and boreal forest ecosystems in ecotonal northern Alaska, USA, for the years 2003–2100. We compared energy feedbacks associated with increases in biomass to energy feedbacks associated with changes in the duration of the snow-free season. We based our simulations on nine input climate scenarios from the Intergovernmental Panel on Climate Change (IPCC) and a new version of the Terrestrial Ecosystem Model (TEM) that incorporates biogeochemistry, vegetation dynamics for multiple PFTs (e.g., trees, shrubs, grasses, sedges, mosses), multiple vegetation pools, and soil thermal regimes. We found mean increases in net primary productivity (NPP) in all PFTs. Most notably, birch ( Betula spp.) in the shrub tundra showed increases that were at least three times larger than any other PFT. Increases in NPP were positively related to increases in growing-season length in the sedge tundra, but PFTs in boreal forest and shrub tundra showed a significant response to changes in light availability as well as growing-season length. Significant NPP responses to changes in vegetation uptake of nitrogen by PFT indicated that some PFTs were better competitors for nitrogen than other PFTs. While NPP increased, heterotrophic respiration ( R H ) also increased, resulting in decreases or no change in net ecosystem carbon uptake. Greater aboveground biomass from increased NPP produced a decrease in summer albedo, greater regional heat absorption (0.34 ± 0.23 W·m −2 ·10 yr −1 [mean ± SD]), and a positive feedback to climate warming. However, the decrease in albedo due to a shorter snow season (−5.1 ± 1.6 d/10 yr) resulted in much greater regional heat absorption (3.3 ± 1.24 W·m −2 ·10 yr −1 ) than that associated with increases in vegetation. Through quantifying feedbacks associated with changes in vegetation and those associated with changes in the snow season length, we can reach a more integrated understanding of the manner in which climate change may impact interactions between high-latitude ecosystems and the climate system.

Alaska↗
USGS2012

Optimum swimming pathways of fish spawning migrations in rivers

Fishes that swim upstream in rivers to spawn must navigate complex fluvial velocity fields to arrive at their ultimate locations. One hypothesis with substantial implications is that fish traverse pathways that minimize their energy expenditure during migration. Here we present the methodological and theoretical developments necessary to test this and similar hypotheses. First, a cost function is derived for upstream migration that relates work done by a fish to swimming drag. The energetic cost scales with the cube of a fish's relative velocity integrated along its path. By normalizing to the energy requirements of holding a position in the slowest waters at the path's origin, a cost function is derived that depends only on the physical environment and not on specifics of individual fish. Then, as an example, we demonstrate the analysis of a migration pathway of a telemetrically tracked pallid sturgeon ( Scaphirhynchus albus ) in the Missouri River (USA). The actual pathway cost is lower than 10 5 random paths through the surveyed reach and is consistent with the optimization hypothesis. The implication—subject to more extensive validation—is that reproductive success in managed rivers could be increased through manipulation of reservoir releases or channel morphology to increase abundance of lower-cost migration pathways.

Missouri River↗
USGS2016

Multi-scale connectivity and graph theory highlight critical areas for conservation under climate change

Conservation planning and biodiversity management require information on landscape connectivity across a range of spatial scales from individual home ranges to large regions. Reduction in landscape connectivity due changes in land-use or development is expected to act synergistically with alterations to habitat mosaic configuration arising from climate change. We illustrate a multi-scale connectivity framework to aid habitat conservation prioritization in the context of changing land use and climate. Our approach, which builds upon the strengths of multiple landscape connectivity methods including graph theory, circuit theory and least-cost path analysis, is here applied to the conservation planning requirements of the Mohave ground squirrel. The distribution of this California threatened species, as for numerous other desert species, overlaps with the proposed placement of several utility-scale renewable energy developments in the American Southwest. Our approach uses information derived at three spatial scales to forecast potential changes in habitat connectivity under various scenarios of energy development and climate change. By disentangling the potential effects of habitat loss and fragmentation across multiple scales, we identify priority conservation areas for both core habitat and critical corridor or stepping stone habitats. This approach is a first step toward applying graph theory to analyze habitat connectivity for species with continuously-distributed habitat, and should be applicable across a broad range of taxa.

Ecological Applications↗
USGS2013

Detecting channel riparian vegetation response to best-management-practices implementation in ephemeral streams with the use of spot high-resolution visible imagery

Heavily grazed riparian areas are commonly subject to channel incision, a lower water table, and reduced vegetation, resulting in sediment delivery above normal regimes. Riparian and in-channel vegetation functions as a roughness element and dissipates flow energy, maintaining stable channel geometry. Ash Creek, a tributary of the Bad River in western South Dakota contains a high proportion of incised channels, remnants of historically high grazing pressure. Best management practices (BMP), including off-stream watering sources and cross fencing, were implemented throughout the Bad River watershed during an Environmental Protection Agency (EPA) 319 effort to address high sediment loads. We monitored prairie cordgrass (Spartina pectinata Link) establishment within stream channels for 16 yr following BMP implementation. Photos were used to group stream reaches (n = 103) subjectively into three classes; absent (estimated < 5% cover; n = 64), present (estimated 5–40% cover; n = 23), and dense (estimated > 40% cover; n = 16) based on the relative amount of prairie cordgrass during 2010 assessments of ephemeral channels. Reaches containing drainage areas of 0.54 to 692 ha were delineated with the use of 2010 National Agriculture Imagery Program (NAIP) imagery. Normalized difference vegetation index (NDVI) values were extracted from 5 to 39 sample points proportional to reach length using a series of Satellite Pour l'Observation de la Terre (SPOT) satellite imagery. Normalized NDVI (nNDVI) of 2 152 sample points were determined from pre- and post-BMP images. Mean nNDVI values for each reach ranged from 0.33 to 1.77. ANOVA revealed significant increase in nNDVI in locations classified as present prairie cordgrass cover following BMP implementation. Establishment of prairie cordgrass following BMP implementation was successfully detected remotely. Riparian vegetation such as prairie cordgrass adds channel roughness that reduces the flow energy responsible for channel degradation.

South Dakota↗