USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interfaces”

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 433 records · Page 24Linked to original sources
USGS2017

Development of a coastal drought index using salinity data

A critical aspect of the uniqueness of coastal drought is the effects on the salinity dynamics of creeks, rivers, and estuaries. The location of the freshwater–saltwater interface along the coast is an important factor in the ecological and socioeconomic dynamics of coastal communities. Salinity is a critical response variable that integrates hydrologic and coastal dynamics including sea level, tides, winds, precipitation, streamflow, and tropical storms. The position of the interface determines the composition of freshwater and saltwater aquatic communities as well as the freshwater availability for water intakes. Many definitions of drought have been proposed, with most describing a decline in precipitation having negative impacts on the water supply. Indices have been developed incorporating data such as rainfall, streamflow, soil moisture, and groundwater levels. These water-availability drought indices were developed for upland areas and may not be ideal for characterizing coastal drought. The availability of real-time and historical salinity datasets provides an opportunity for the development of a salinity-based coastal drought index. An approach similar to the standardized precipitation index (SPI) was modified and applied to salinity data obtained from sites in South Carolina and Georgia. Using the SPI approach, the index becomes a coastal salinity index (CSI) that characterizes coastal salinity conditions with respect to drought periods of higher-saline conditions and wet periods of higher-freshwater conditions. Evaluation of the CSI indicates that it provides additional coastal response information as compared to the SPI and the Palmer hydrologic drought index, and the CSI can be used for different estuary types and for comparison of conditions along coastlines.

Bulletin of the American Meteorological Society↗
USGS2021

Spectral damping scaling factors for horizontal components of ground motions from subduction earthquakes using NGA-Subduction data

This article develops global models of damping scaling factors (DSFs) for subduction zone earthquakes that are functions of the damping ratio, spectral period, earthquake magnitude, and distance. The Next Generation Attenuation for subduction earthquakes (NGA-Sub) project has developed the largest uniformly processed database of recorded ground motions to date from seven subduction regions: Alaska, Cascadia, Central America and Mexico, South America, Japan, Taiwan, and New Zealand. NGA-Sub used this database to develop new ground motion models (GMMs) at a reference 5% damping ratio. We worked with the NGA-Sub project team to develop an extended database that includes pseudo-spectral accelerations (PSA) for 11 damping ratios between 0.5% and 30%. We use this database to develop parametric models of DSF for both interface and intraslab subduction earthquakes that can be used to adjust any subduction GMM from a reference 5% damping ratio to other damping ratios. The DSF is strongly influenced by the response spectral shape and the duration of motion; therefore, in addition to the damping ratio, the median DSF model uses spectral period, magnitude, and distance as surrogate predictor variables to capture the effects of the spectral shape and the duration of motion. We also develop parametric models for the standard deviation of DSF. The models presented in this article are for the RotD50 horizontal component of PSA and are compared with the models for shallow crustal earthquakes in active tectonic regions. Some noticeable differences arise from the considerably longer duration of interface records for very large magnitude events and the enriched high-frequency content of intraslab records, compared with shallow crustal earthquakes. Regional differences are discussed by comparing the proposed global models with the data from each subduction region along with recommendations on the applicability of the models.

Earthquake Spectra↗
USGS2021

NGA-subduction global ground motion models with regional adjustment factors

We develop semi-empirical ground motion models (GMMs) for peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral accelerations for periods from 0.01 to 10 s, for the median orientation-independent horizontal component of subduction earthquake ground motion. The GMMs are applicable to interface and intraslab subduction earthquakes in Japan, Taiwan, Mexico, Central America, South America, Alaska, the Aleutian Islands, and Cascadia. The GMMs are developed using a combination of data inspection, data regression with respect to physics-informed functions, ground-motion simulations, and geometrical constraints for certain model components. The GMMs capture observed differences in source and path effects for interface and intraslab events, conditioned on moment magnitude, rupture distance, and hypocentral depth. Site effect and aleatory variability models are shared between event types. Regionalized GMM components include the model constant (that controls ground motion amplitude), anelastic attenuation, magnitude-scaling break point, linear site response, and sediment depth terms. We develop models for the aleatory between-event variability ( τ ) "> ( τ ) , within-event variability ( ϕ ) "> ( ϕ ) , single-station within-event variability ( ϕ SS ) "> ( ϕ SS ) , and site-to-site variability ( ϕ S 2 S ) "> ( ϕ S 2 S ) . Ergodic analyses should use the median GMM and aleatory variability computed using the between-event and within-event variability models. An analysis incorporating non-ergodic site response should use the median GMM at the reference shear-wave velocity condition, a site-specific site response model, and aleatory variability computed using the between-event and single-station within-event variability models. Epistemic uncertainty in the median model is represented by standard deviations on the regional model constants, which facilitates scaled-backbone representations of model uncertainty in hazard analyses.

Earthquake Spectra↗
USGS2016

A decision framework for identifying models to estimate forest ecosystem services gains from restoration

Restoring degraded forests and agricultural lands has become a global conservation priority. A growing number of tools can quantify ecosystem service tradeoffs associated with forest restoration. This evolving “tools landscape” presents a dilemma: more tools are available, but selecting appropriate tools has become more challenging. We present a Restoration Ecosystem Service Tool Selector (RESTS) framework that describes key characteristics of 13 ecosystem service assessment tools. Analysts enter information about their decision context, services to be analyzed, and desired outputs. Tools are filtered and presented based on five evaluative criteria: scalability, cost, time requirements, handling of uncertainty, and applicability to benefit-cost analysis. RESTS uses a spreadsheet interface but a web-based interface is planned. Given the rapid evolution of ecosystem services science, RESTS provides an adaptable framework to guide forest restoration decision makers toward tools that can help quantify ecosystem services in support of restoration.

Forest Ecosystems↗
USGS1994

Large-explosive source, wide-recording aperture, seismic profiling on the Columbia Plateau, Washington

Clear subsurface seismic images have been obtained at low cost on the Columbia Plateau, Washington. The Columbia Plateau is perhaps the most notorious of all 'bad-data' areas because large impedance contrasts in surface flood basalts severely degrade the seismic wavefield. This degradation was mitigated in this study via a large-explosive source, wide-recording aperture shooting method.The shooting method emphasizes the wide-angle portion of the wavefield, where Fermat's principle guarantees reverberation will not interfere with the seismic manifestations of crucial geologic interfaces. The basalt diving wave, normally discarded in standard common midpoint (CMP) seismic profiling, can be used to image basalt velocity structure via traveltime inversion. Maximum depth-penetration of the diving wave tightly constrains basalt-sediment interface depth. An arrival observed only at shot-receiver offsets greater than 15 km can be used to determine the velocity and geometry of basement via simultaneous inversion.The results from this study suggest that previous geologic hypotheses and hydrocarbon play concepts for the Columbia Plateau may have been in error.

Geophysics↗
USGS1994

An image of the Columbia Plateau from inversion of high‐resolution seismic data

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.

Geophysics↗
USGS2007

Seismoelectric imaging of shallow targets

We have undertaken a series of controlled field experiments to develop seismoelectric experimental methods for near-surface applications and to improve our understanding of seismoelectric phenomena. In a set of off-line geometry surveys (source separated from the receiver line), we place seismic sources and electrode array receivers on opposite sides of a man-made target (two sand-filled trenches) to record separately two previously documented seismoelectric modes: (1) the electromagnetic interface response signal created at the target and (2) the coseismic electric fields located within a compressional seismic wave. With the seismic source point in the center of a linear electrode array, we identify the previously undocumented seismoelectric direct field, and the Lorentz field of the metal hammer plate moving in the earth's magnetic field. We place the seismic source in the center of a circular array of electrodes (radial and circumferential orientations) to analyze the source-related direct and Lorentz fields and to establish that these fields can be understood in terms of simple analytical models. Using an off-line geometry, we create a multifold, 2D image of our trenches as dipping layers, and we also produce a complementary synthetic image through numerical modeling. These images demonstrate that off-line geometry (e.g., crosswell) surveys offer a particularly promising application of the seismoelectric method because they effectively separate the interface response signal from the (generally much stronger) coseismic and source-related fields. ?? 2007 Society of Exploration Geophysicists.

Geophysics↗
USGS2019

Shear-wave seismic reflection studies of unconsolidated sediments in the near surface

We have successfully applied of SH-wave seismic reflection methods to two different near-surface problems targeting unconsolidated sediments. At the former Fort Ord, where the water table is approximately 30m deep, we imaged aeolian and marine aquifer and aquitard stratigraphy to a depth of approximately 80m. We identified reflections from sand/clay and sand/silt interfaces and we mapped these interfaces along our transects. At an aggregate study site in Indiana, where the water table is at a depth of 1to2m, we imaged stratigraphy in alluvial sand and gravel, and observe a strong reflection from the 32-m-deep bedrock surface. In both cases, we exploited the high resolution potential of SH waves, their insensitivity to water content, and the possibility of reducing Love wave contamination by working along a roadway. We accomplished our results using only sledgehammer sources and simple data processing flows.

California, Indiana↗
USGS2015

Ground motion models used in the 2014 U.S. National Seismic Hazard Maps

The National Seismic Hazard Maps (NSHMs) are an important component of seismic design regulations in the United States. This paper compares hazard using the new suite of ground motion models (GMMs) relative to hazard using the suite of GMMs applied in the previous version of the maps. The new source characterization models are used for both cases. A previous paper ( Rezaeian et al. 2014 ) discussed the five NGA-West2 GMMs used for shallow crustal earthquakes in the Western United States (WUS), which are also summarized here. Our focus in this paper is on GMMs for earthquakes in stable continental regions in the Central and Eastern United States (CEUS), as well as subduction interface and deep intraslab earthquakes. We consider building code hazard levels for peak ground acceleration (PGA), 0.2-s, and 1.0-s spectral accelerations (SAs) on uniform firm-rock site conditions. The GMM modifications in the updated version of the maps created changes in hazard within 5% to 20% in WUS; decreases within 5% to 20% in CEUS; changes within 5% to 15% for subduction interface earthquakes; and changes involving decreases of up to 50% and increases of up to 30% for deep intraslab earthquakes for most U.S. sites. These modifications were combined with changes resulting from modifications in the source characterization models to obtain the new hazard maps.

Earthquake Spectra↗
USGS2024

The role of clay minerals in the preservation of organic matter in sediments of Qinghai Lake, NW China

The role of saline lake sediments in preserving organic matter has long been recognized. In order to further understand the preservation mechanisms, the role of clay minerals was studied. Three sediment cores, 25, 57, and 500 cm long, were collected from Qinghai Lake, NW China, and dissected into multiple subsamples. Multiple techniques were employed, including density fractionation, X-ray diffraction, scanning and transmission electron microscopy (SEM and TEM), total organic carbon (TOC) and carbon compound analyses, and surface area determination. The sediments were oxic near the water-sediment interface, but became anoxic at depth. The clay mineral content was as much as 36.8%, consisting mostly of illite, chlorite, and halloysite. The TEM observations revealed that organic matter occurred primarily as organic matter-clay mineral aggregates. The TOC and clay mineral abundances are greatest in the mid-density fraction, with a positive correlation between the TOC and mineral surface area. The TOC of the bulk sediments ranges from 1 to 3% with the non-hydrocarbon fraction being predominant, followed by bitumen, saturated hydrocarbon, aromatic hydrocarbons, and chloroform-soluble bitumen. The bimodal distribution of carbon compounds of the saturated hydrocarbon fraction suggests that organic matter in the sediments was derived from two sources: terrestrial plants and microorganisms/algae. Depthrelated systematic changes in the distribution patterns of the carbon compounds suggest that the oxidizing conditions and microbial abundance near the water-sediment interface promote degradation of labile organic matter, probably in adsorbed form. The reducing conditions and small microbial biomass deeper in the sediments favor preservation of organic matter, because of the less labile nature of organic matter, probably occurring within clay mineral-organic matter aggregates that are inaccessible to microorganisms. These results have important implications for our understanding of mechanisms of organic matter preservation in saline lake sediments.

Qinghai Lake↗
USGS2012

A global earthquake discrimination scheme to optimize ground-motion prediction equation selection

We present a new automatic earthquake discrimination procedure to determine in near-real time the tectonic regime and seismotectonic domain of an earthquake, its most likely source type, and the corresponding ground-motion prediction equation (GMPE) class to be used in the U.S. Geological Survey (USGS) Global ShakeMap system. This method makes use of the Flinn–Engdahl regionalization scheme, seismotectonic information (plate boundaries, global geology, seismicity catalogs, and regional and local studies), and the source parameters available from the USGS National Earthquake Information Center in the minutes following an earthquake to give the best estimation of the setting and mechanism of the event. Depending on the tectonic setting, additional criteria based on hypocentral depth, style of faulting, and regional seismicity may be applied. For subduction zones, these criteria include the use of focal mechanism information and detailed interface models to discriminate among outer-rise, upper-plate, interface, and intraslab seismicity. The scheme is validated against a large database of recent historical earthquakes. Though developed to assess GMPE selection in Global ShakeMap operations, we anticipate a variety of uses for this strategy, from real-time processing systems to any analysis involving tectonic classification of sources from seismic catalogs.

Bulletin of the Seismological Society of America↗
USGS2014

Seismic‐wave attenuation determined from tectonic tremor in multiple subduction zones

Tectonic tremor provides a new source of observations that can be used to constrain the seismic attenuation parameter for ground‐motion prediction and hazard mapping. Traditionally, recorded earthquakes of magnitude ∼3–8 are used to develop ground‐motion prediction equations; however, typical earthquake records may be sparse in areas of high hazard. In this study, we constrain the distance decay of seismic waves using measurements of the amplitude decay of tectonic tremor, which is plentiful in some regions. Tectonic tremor occurs in the frequency band of interest for ground‐motion prediction (i.e., ∼2–8 Hz) and is located on the subducting plate interface, at the lower boundary of where future large earthquakes are expected. We empirically fit the distance decay of peak ground velocity from tremor to determine the attenuation parameter in four subduction zones: Nankai, Japan; Cascadia, United States–Canada; Jalisco, Mexico; and southern Chile. With the large amount of data available from tremor, we show that in the upper plate, the lower crust is less attenuating than the upper crust. We apply the same analysis to intraslab events in Nankai and show the possibility that waves traveling from deeper intraslab events experience more attenuation than those from the shallower tremor due to ray paths that pass through the subducting and highly attenuating oceanic crust. This suggests that high pore‐fluid pressure is present in the tremor source region. These differences imply that the attenuation parameter determined from intraslab earthquakes may underestimate ground motion for future large earthquakes on the plate interface.

Cascadia↗
USGS2016

Toppling analysis of the Echo Cliffs precariously balanced rock

Toppling analysis of a precariously balanced rock (PBR) can provide insight into the nature of ground motion that has not occurred at that location in the past and, by extension, can constrain peak ground motions for use in engineering design. Earlier approaches have targeted 2D models of the rock or modeled the rock–pedestal contact using spring‐damper assemblies that require recalibration for each rock. Here, a method to model PBRs in 3D is presented through a case study of the Echo Cliffs PBR. The 3D model is created from a point cloud of the rock, the pedestal, and their interface, obtained using terrestrial laser scanning. The dynamic response of the model under earthquake excitation is simulated using a rigid‐body dynamics algorithm. The veracity of this approach is demonstrated through comparisons against data from shake‐table experiments. Fragility maps for toppling probability of the Echo Cliffs PBR as a function of various ground‐motion parameters, rock–pedestal interface friction coefficient, and excitation direction are presented. These fragility maps indicate that the toppling probability of this rock is low (less than 0.2) for peak ground acceleration (PGA) and peak ground velocity (PGV) lower than 3 m/s 2 and 0.75 m/s, respectively, suggesting that the ground‐motion intensities at this location from earthquakes on nearby faults have most probably not exceeded the above‐mentioned PGA and PGV during the age of the PBR. Additionally, the fragility maps generated from this methodology can also be directly coupled with existing probabilistic frameworks to obtain direct constraints on unexceeded ground motion at a PBR’s location.

California↗
USGS2021

The impact of 3D finite‐fault information on ground‐motion forecasting for earthquake early warning

We identify aspects of finite‐source parameterization that strongly affect the accuracy of estimated ground motion for earthquake early warning (EEW). EEW systems aim to alert users to impending shaking before it reaches them. The U.S. West Coast EEW system, ShakeAlert, currently uses two algorithms based on seismic data to characterize the earthquake’s location, magnitude, and origin time, treating it as a point or line source. From this information, ShakeAlert calculates shaking intensity and alerts locations where shaking estimates exceed a threshold. Several geodetic EEW algorithms under development would provide 3D finite‐fault information. We investigate conditions under which this information produces sufficiently better intensity estimates to potentially improve alerting. Using scenario crustal and subduction interface sources, we (1) identify the most influential source geometry parameters for an EEW algorithm’s shaking forecast, and (2) assess the intensity alert thresholds and magnitude ranges for which more detailed source characterization affects alert accuracy. We find that alert regions determined using 3D‐source representations of correct magnitude and faulting mechanism are generally more accurate than those obtained using line sources. If a line‐source representation is used and magnitude is calculated from the estimated length, then incorrect length estimates significantly degrade alert region accuracy. In detail, the value of 3D‐source characterization depends on the user’s chosen alert threshold, tectonic regime, and faulting style. For the suite of source models we tested, the error in shaking intensity introduced by incorrect geometry could reach levels comparable to the intrinsic uncertainty in ground‐motion calculations (e.g., 0.5–1.3 modified Mercalli intensity [MMI] units for MMI 4.5) but, especially for crustal sources, was often less. For subduction interface sources, 3D representations substantially improved alert area accuracy compared to line sources, and incorrect geometry parameters were more likely to cause error in calculated shaking intensity that exceeded uncertainties.

Bulletin of the Seismological Society of America↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗
USGS1995

Near real-time monitoring of seismic events and status of portable digital recorders using satellite telemetry

Near real-time monitoring of seismic events and status of portable 16-bit digital recorders has been established for arrays near Parkfield, Mammoth Lakes, and San Francisco, California. This monitoring system provides near real-time seismic event identification (rough location and magnitude) and a cost-effective means to maintain arrays at near 100% operational level. Principal objectives in the design of this system have been portability and low-cost telemetry. The system has been developed to use portable digital seismic recorders ( GEOS —General Earthquake Observation System) and portable data collection platforms ( DCP 's) for the Geostationary Operational Environmental Satellite (GEOS) telemetry system. Data are transferred asynchronously from the GEOS seismic system through a microprocessor-controlled interface every 10 min. The interface stores, determines priority, converts, and synchronously transfers these data to a Sutron Corp. model 8004 DCP for transmission through the GEOS satellite telemetry system. Event parameters include trigger time, peak amplitude, time of peak amplitude, and event duration. Instrument configuration parameters, transmitted at system start-up time and every 24 hr, include recording parameters, trigger parameters, GEOS software version, clock reference, and location parameter. Instrument status includes battery voltage, number of events, and percentage of tape usage. These data are transmitted as appropriate to the U.S. Geological Survey satellite downlink and computers located in Menlo Park, California, where they are processed and displayed.

Bulletin of the Seismological Society of America↗
USGS2007

Human influence on California fire regimes

Periodic wildfire maintains the integrity and species composition of many ecosystems, including the mediterranean-climate shrublands of California. However, human activities alter natural fire regimes, which can lead to cascading ecological effects. Increased human ignitions at the wildland-urban interface (WUI) have recently gained attention, but fire activity and risk are typically estimated using only biophysical variables. Our goal was to determine how humans influence fire in California and to examine whether this influence was linear, by relating contemporary (2000) and historic (1960-2000) fire data to both human and biophysical variables. Data for the human variables included fine-resolution maps of the WUI produced using housing density and land cover data. Interface WUI, where development abuts wildland vegetation, was differentiated from intermix WUI, where development intermingles with wildland vegetation. Additional explanatory variables included distance to WUI, population density, road density, vegetation type, and ecoregion. All data were summarized at the county level and analyzed using bivariate and multiple regression methods. We found highly significant relationships between humans and fire on the contemporary landscape, and our models explained fire frequency (R2 = 0.72) better than area burned (R2 = 0.50). Population density, intermix WUI, and distance to WUI explained the most variability in fire frequency, suggesting that the spatial pattern of development may be an important variable to consider when estimating fire risk. We found nonlinear effects such that fire frequency and area burned were highest at intermediate levels of human activity, but declined beyond certain thresholds. Human activities also explained change in fire frequency and area burned (1960-2000), but our models had greater explanatory power during the years 1960-1980, when there was more dramatic change in fire frequency. Understanding wildfire as a function of the spatial arrangement of ignitions and fuels on the landscape, in addition to nonlinear relationships, will be important to fire managers and conservation planners because fire risk may be related to specific levels of housing density that can be accounted for in land use planning. With more fires occurring in close proximity to human infrastructure, there may also be devastating ecological impacts if development continues to grow farther into wildland vegetation. ?? 2007 by the Ecological Society of America.

Ecological Applications↗
USGS2011

Cougar space use and movements in the wildland-urban landscape of western Washington

The wildland-urban interface lies at the confluence of human-dominated and wild landscapes, creating a number of management and conservation challenges. Because wildlife ecology, behavior, and evolution at this interface are shaped by both natural and human phenomena, this requires greater understanding of how diverse factors affect ecosystem and population processes. We illustrate the challenge of understanding and managing a frequent and often undesired inhabitant of the wildland-urban landscape, the cougar (Puma concolor). In wildland and residential areas of western Washington State, USA, we captured and radiotracked 27 cougars to model space use and understand the role of landscape features in interactions (sightings, encounters, and depredations) between cougars and humans. Resource utilization functions (RUFs) identified cougar use of areas with features that were probably attractive to prey, influential on prey vulnerability, and associated with limited or no residential development. Early-successional forest (+), conifer forest (+), distance to road (-), residential density (-), and elevation (-) were significant positive and negative predictors of use for the population, whereas use of other landscape features was highly variable. Space use and movement rates in wildland and residential areas were similar because cougars used wildland-like forest patches, reserves, and corridors in residential portions of their home range. The population RUF was a good predictor of confirmed cougar interactions, with 72% of confirmed reports occurring in the 50% of the landscape predicted to be medium-high and high cougar use areas. We believe that there is a threshold residential density at which the level of development modifies the habitat but maintains enough wildland characteristics to encourage moderate levels of cougar use and maximize the probability of interaction. Wildlife managers trying to reduce interactions between cougars and people should incorporate information on spatial ecology and landscape characteristics to identify areas with the highest overlap of human and cougar use to focus management, education, and landscape planning. Resource utilization functions provide a proactive tool to guide these activities for improved coexistence with wildlife using both wildland and residential portions of the landscape. ??2011 by the Ecological Society of America.

Ecological Applications↗