SEARCH · USGS Science
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.
An overview is given of examples showing the utility of the Comprehensive Model (CM) approach in the processing of near-surface magnetic anomalies. Effectiveness of CM depends on the spatial and temporal matching of the main and external field variations that it attempts to approximate. Hence, data distribution as well as continuity are particularly important.
With phase inversion, one can estimate subsurface velocities using the phases of first-arriving waves, which are the frequency-domain equivalents of the traveltimes. Phase inversion is modified to make it suitable for processing traveltimes from near-surface refraction surveys. The modifications include parameterizing the model, correcting the observed phases, and selecting the complex frequency. This modified phase inversion is compared to traveltime tomography. For two comparisons using computer-simulated traveltimes, the difference between the estimated and correct models, the residual mean, and the residual standard deviation are smaller for the phase inversion than they are for the traveltime tomography. For a comparison using field data from an S-wave refraction survey, both methods estimate models that are consistent with the known geology. Nonetheless, the phase-inversion model includes small-scale features in the bedrock that are geologically plausible; the residual mean and the residual standard deviation are smaller for the phase inversion than they are for the traveltime tomography.
Invasive species threaten both ecological systems and economies. Successful establishment is governed partly by demographic stochasticity, where random individual-level variation can significantly influence invasion dynamics. A critical factor is the Allee effect, which describes reduced population growth at low densities, potentially serving as a barrier to persistence. Size-dependent growth and survival further shape outcomes based on the size distribution of individuals within the population. To address these complexities, we developed a population-level model integrating an Integral Projection Model (IPM) for the expected size distribution with an Allee effect for recruitment dynamics and a branching process to account for demographic stochasticity. We applied this framework to Silver Carp ( Hypophthalmichthys molitrix ), a species established in the Upper Mississippi River system that poses significant risk of expanding into the Great Lakes. Our analysis compared temporal dynamics and population-size distributions under hypothetical introduction scenarios for large adult and smaller subadult fish. In both cases, initial population size relative to the Allee threshold was a critical determinant of establishment probability, with introductions below the threshold requiring stochastic growth across a critical biomass boundary before self-sustaining dynamics could take hold. Sensitivity analysis further indicated that per-spawner recruitment played a compounding role in invasion risk: higher recruitment rates reduced the critical introduction size for establishment and compressed the time to self-sustaining population size, jointly broadening the window of invasion risk. These findings highlight the framework’s potential for invasion assessments and offer insights into how Allee effects, size structure, and recruitment rates interact to govern establishment outcomes, while underscoring the need for species-specific parameter estimates before quantitative outputs can guide management decisions.
A comprehensive mapping technology was developed utilizing standard image processing and available GIS procedures to automate shoreline identification and mapping from 2 m synthetic aperture radar ( SAR ) HH amplitude data. The development used four NASA Uninhabited Aerial Vehicle SAR ( UAVSAR ) data collections between summer 2009 and 2012 and a fall 2012 collection of wetlands dominantly fronted by vegetated shorelines along the Mississippi River Delta that are beset by severe storms, toxic releases, and relative sea-level rise. In comparison to shorelines interpreted from 0.3 m and 1 m orthophotography, the automated GIS 10 m alongshore sampling found SAR shoreline mapping accuracy to be ±2 m, well within the lower range of reported shoreline mapping accuracies. The high comparability was obtained even though water levels differed between the SAR and photography image pairs and included all shorelines regardless of complexity. The SAR mapping technology is highly repeatable and extendable to other SAR instruments with similar operational functionality.
This chapter explores fire as a physical process, including combustion, fuel characteristics, fuel models, fire weather, ignition sources, mechanisms for fire spread, and fire effects. In wildland fuels, combustion occurs in three phases: preheating, gaseous, and smoldering. Fuel is characterized by physical and chemical properties that affect combustion and fire behavior. Its characteristic classes are defined for a vegetation type and contain data for fuels in up to six strata representing potentially independent combustion environments. Fire weather includes air temperature, atmospheric moisture, atmospheric stability, and clouds and precipitation. Sufficient fuel, conducive weather, and an ignition are necessary ingredients for a fire. In line with this, this chapter investigates how these factors, combined with topography, cause a fire to spread. The chapter also introduces the physical parameters of fire behavior that affect fire severity, spotting, tree scorch height, plant mortality, biomass consumption, and microclimate.
This chapter focuses on the hydrological setting of Mirror Lake and its water budget. It first describes the glacial deposits and bedrock topography in the Mirror Lake area. It then provides an overview of the hydrologic processes associated with Mirror Lake and examines the field and analytical methods used to determine its water budget. It presents results from the hydrologic studies, which are based on monthly and annual water budgets for the calendar years 1981 through 2000.
The Middle – Upper Triassic Shublik Formation is an organic-rich heterogeneous carbonate-siliciclastic-phosphatic unit that generated much of the oil in the Prudhoe Bay field and other hydrocarbon accumulations in northern Alaska. A large dataset, including total organic carbon (TOC), X-ray diffraction (XRD), X-ray fluorescence (XRF) and inductively coupled plasma – mass spectrometry (ICP-MS) measurements, has been built from core and outcrop samples of the Shublik, with a focus on the organic-rich intervals. In addition, two core intervals from the Shublik were analyzed using a hyperspectral imaging system in the visible, near-infrared and shortwave-infrared range. Integration of the hyperspectral results with core descriptions, microfacies interpretations, and analytical data is being used to decipher mudstone depositional and diagenetic processes. Petrographic analysis of Upper Triassic organic-rich intervals within the Shublik suggests that the main microfacies is a laminated bioclastic wackestone/packstone that was episodically disrupted by energetic events of variable intensity. These energetic events produced transitional and sparry calcite bioclastic packstone to grainstone intervals, depending on the depth of sediment column disturbance. By using hyperspectral imaging data from the Ikpikpuk core, individual distribution maps for minerals of interest have been generated and corroborate the microfacies interpretations. These maps also illustrate small-scale vertical changes in mineralogy. The laminated bioclastic wackestone/packstone intervals contain less calcite than the adjacent sparry bioclastic packstone to grainstone intervals. The calcite in these laminated intervals is more iron rich. This interpretation suggests that lower iron concentrations should be expected in the disrupted intervals than in nearby laminated intervals. Textural features are also enhanced in the hyperspectral images relative to visual description of the cores by combining the extraction of the average reflectance in the visible part of the electromagnetic spectrum and the depth of the main carbonate-related feature belonging to calcite. Examples noted in the enhanced imagery include low-angle features, calcite grain-size, and the size, shape and orientation of phosphatic nodules. This enhancement is being used to differentiate laminated from sparry bioclastic packstone to grainstone-rich intervals and provides a more comprehensive assessment of the microfacies than is practical by thin-section analysis.
An experimental computer-generated shaded-relief map of a topographic quadrangle of Nome, Alaska, at 1:1,000,000 scale was created to demonstrate current capabilities of map production. This paper describes the techniques used to produce the map product. The data consisted of forty-eight 1° by 1° blocks of resampled digital elevation model (DEM) data. These data were digitally mosaicked and assigned colors based on intervals of elevation values. The color-coded data set was then used to create a shaded-relief image that was photographically composited with cartographic line information to produce a shaded-relief map. The majority of the processing was completed at the National Mapping Division EROS Data Center in Sioux Falls, South Dakota.
Review of: The recovery process in damaged ecosystems; J. Cairnsed.; 1980; Ann Arbor Science Publisher, Ann Arbor, Michigan.
Habitat use and movements of 25 adult cutthroat trout Oncorhynchus clarkii and 25 adult brook trout Salvelinus fontinalis from fall through winter 2002-2003 were assessed by means of radiotelemetry in a 7-km reach of a Rocky Mountains foothills stream. Temporal dynamics of winter habitat conditions were evaluated by regularly measuring the features of 30 pools and 5 beaver Castor canadensis ponds in the study reach. Groundwater inputs at three locations raised mean daily water temperatures in the stream channel during winter to 0.2-0.6??C and kept at least 250 m of the downstream channel free of ice, but the lack of surface ice further downstream led to the occurrence of frazil ice and anchor ice in pools and unstable habitat conditions for trout. Pools in segments that were not affected by groundwater inputs and beaver ponds tended to be stable and snow accumulated on the surface ice. Pools throughout the study reach tended to become more stable as snow accumulated. Both cutthroat trout and brook trout selected beaver ponds as winter progressed but tended to use lateral scour pools in proportion to their availability. Tagged fish not in beaver ponds selected lateral scour pools that were deeper than average and stable during winter. Movement frequencies by tagged fish decreased from fall through winter, but some individuals of both species moved during winter. Ice processes affected both the habitat use and movement patterns of cutthroat trout and brook trout in this foothills stream.
There is little information on the winter features of salmonid habitats associated with constructed instream structures to provide guidance when planning habitat improvement projects. We assessed winter habitat features for trout of the genera Oncorhynchus and Salvelinus in pools associated with two types of instream structures constructed on a low-gradient reach of a mountain stream in western Wyoming with a mean wetted width of 6.4 m. Pool habitat was affected by temporal variability in ice formations from fall into winter. As surface ice and snow accumulated with the progression of winter, variation in ice formations was less frequent and winter habitat conditions became more stable. However, groundwater inflow that maintained water temperatures at 0.2-0.6??C in a portion of the study reach appeared to contribute to incomplete surface ice cover and variation in ice formations in pools through most of the winter. Hanging dams and anchor ice dams were the primary ice features that affected winter habitat in pools associated with constructed instream structures. Trout were observed in these pools in the fall but tended to abandon pools with variation in ice formations as winter progressed. The potential impacts of groundwater inflow and winter ice processes on trout habitat in pools associated with instream structures should be considered when planning habitat improvement projects. ?? Copyright by the American Fisheries Society 2005.
Yellow perch Perca flavescens support economically important recreational and commercial fisheries in Lake Erie and are intensively managed. Age estimation represents an integral component in the management of Lake Erie yellow perch stocks, as age-structured population models are used to set safe harvest levels on an annual basis. We compared the precision associated with yellow perch (N = 251) age estimates from scales, sagittal otoliths, and anal spine sections and evaluated the time required to process and estimate age from each structure. Three readers of varying experience estimated ages. The precision (mean coefficient of variation) of estimates among readers was 1% for sagittal otoliths, 5-6% for anal spines, and 11-13% for scales. Agreement rates among readers were 94-95% for otoliths, 71-76% for anal spines, and 45-50% for scales. Systematic age estimation differences were evident among scale and anal spine readers; less-experienced readers tended to underestimate ages of yellow perch older than age 4 relative to estimates made by an experienced reader. Mean scale age tended to underestimate ages of age-6 and older fish relative to otolith ages estimated by an experienced reader. Total annual mortality estimates based on scale ages were 20% higher than those based on otolith ages; mortality estimates based on anal spine ages were 4% higher than those based on otolith ages. Otoliths required more removal and preparation time than scales and anal spines, but age estimation time was substantially lower for otoliths than for the other two structures. We suggest the use of otoliths or anal spines for age estimation in yellow perch (regardless of length) from Lake Erie and other systems where precise age estimates are necessary, because age estimation errors resulting from the use of scales could generate incorrect management decisions. ?? Copyright by the American Fisheries Society 2008.
Coastal erosion rates locally exceeding 30 m y-1 have been documented along Alaska's Beaufort Sea coastline, and a number of studies suggest that these erosion rates have accelerated as a result of climate change. However, a lack of direct observational evidence has limited our progress in quantifying the specific processes that connect climate change to coastal erosion rates in the Arctic. In particular, while longer ice-free periods are likely to lead to both warmer surface waters and longer fetch, the relative roles of thermal and mechanical (wave) erosion in driving coastal retreat have not been comprehensively quantified. We focus on a permafrost coastline in the northern National Petroleum Reserve-Alaska (NPR-A), where coastal erosion rates have averaged 10-15 m y-1 over two years of direct monitoring. We take advantage of these extraordinary rates of coastal erosion to observe and quantify coastal erosion directly via time-lapse photography in combination with meteorological observations. Our observations indicate that the erosion of these bluffs is largely thermally driven, but that surface winds play a crucial role in exposing the frozen bluffs to the radiatively warmed seawater that drives melting of interstitial ice. To first order, erosion in this setting can be modeled using formulations developed to describe iceberg deterioration in the open ocean. These simple models provide a conceptual framework for evaluating how climate-induced changes in thermal and wave energy might influence future erosion rates in this setting.
Celastrus orbiculatus is a non-native, invasive liana that was introduced to the United States in the 1860s and has spread rapidly throughout the Northeast Several attributes contribute to the invasiveness of C. orbiculatus, including tolerance to a wide range of light levels and habitat types. We compared soil characteristics in seven sets of adjacent, paired plots, spanning a range of habitats and soil types, with and without C. orbiculatus. The paired plots were similar other than the presence or absence of Celastrus. Plots with C. orbiculatus had significantly higher soil pH, potassium, calcium and magnesium levels. Furthermore, nitrogen mineralization and litter decomposition rates were higher in plots with C. orbiculatus. Phosphorus levels were not significantly different between the paired plots. The results of this study contribute to the growing body of research of the effects of invasive species on ecosystem processes.
Processing of strong-motion data in many cases can be as straightforward as filtering the acceleration time series and integrating to obtain velocity and displacement. To avoid the introduction of spurious low-frequency noise in quantities derived from the filtered accelerations, however, care must be taken to append zero pads of adequate length to the beginning and end of the segment of recorded data. These padded sections of the filtered acceleration need to be retained when deriving velocities, displacements, Fourier spectra, and response spectra. In addition, these padded and filtered sections should also be included in the time series used in the dynamic analysis of structures and soils to ensure compatibility with the filtered accelerations.
Displacement, velocity, and velocity-squared records of P and SH body waves recorded at teleseismic distances are analyzed to determine the rupture characteristics of the Denali fault, Alaska, earthquake of 3 November 2002 ( M W 7.9, M e 8.1 ). Three episodes of rupture can be identified from broadband ( ∼0.1–5.0 Hz ) waveforms. The Denali fault earthquake started as a M W 7.3 thrust event. Subsequent right-lateral strike-slip rupture events with centroid depths of 9 km occurred about 22 and 49 sec later. The teleseismic P waves are dominated by energy at intermediate frequencies (0.1-1 Hz) radiated by the thrust event, while the SH waves are dominated by energy at lower frequencies (0.05-0.2 Hz) radiated by the strike-slip events. The strike-slip events exhibit strong directivity in the teleseismic SH waves. Correcting the recorded P -wave acceleration spectra for the effect of the free surface yields an estimate of 2.8 × 10 15 N m for the energy radiated by the thrust event. Correcting the recorded SH -wave acceleration spectra similarly yields an estimate of 3.3 × 10 16 N m for the energy radiated by the two strike-slip events. The average rupture velocity for the strike-slip rupture process is 1.1 β –1.2 β . The strike-slip events were located 90 and 188 km east of the epicenter. The rupture length over which significant or resolvable energy is radiated is, thus, far shorter than the 340-km fault length over which surface displacements were observed. However, the seismic moment released by these three events, 4 × 10 20 N m, was approximately half the seismic moment determined from very low-frequency analyses of the earthquake. The difference in seismic moment can be reasonably attributed to slip on fault segments that did not radiate significant or coherent seismic energy. These results suggest that very large and great strike-slip earthquakes can generate stress pulses that rapidly produce substantial slip with negligible stress drop and little discernible radiated energy on fault segments distant from the initial point of nucleation. The existence of this energy-deficient rupture mode has important implications for the evaluation of the seismic hazard of very large strike-slip earthquakes.
Displacement waveforms and high-frequency acceleration envelopes from stations at distances of 3-300 km were inverted to determine the source process of the M 7.9 Denali fault earthquake. Fitting the initial portion of the displacement waveforms indicates that the earthquake started with an oblique thrust subevent (subevent # 1) with an east-west-striking, north-dipping nodal plane consistent with the observed surface rupture on the Susitna Glacier fault. Inversion of the remainder of the waveforms (0.02-0.5 Hz) for moment release along the Denali and Totschunda faults shows that rupture proceeded eastward on the Denali fault, with two strike-slip subevents (numbers 2 and 3) centered about 90 and 210 km east of the hypocenter. Subevent 2 was located across from the station at PS 10 (Trans-Alaska Pipeline Pump Station #10) and was very localized in space and time. Subevent 3 extended from 160 to 230 km east of the hypocenter and had the largest moment of the subevents. Based on the timing between subevent 2 and the east end of subevent 3, an average rupture velocity of 3.5 km/sec, close to the shear wave velocity at the average rupture depth, was found. However, the portion of the rupture 130-220 km east of the epicenter appears to have an effective rupture velocity of about 5.0 km/ sec, which is supershear. These two subevents correspond approximately to areas of large surface offsets observed after the earthquake. Using waveforms of the M 6.7 Nenana Mountain earthquake as empirical Green's functions, the high-frequency (1-10 Hz) envelopes of the M 7.9 earthquake were inverted to determine the location of high-frequency energy release along the faults. The initial thrust subevent produced the largest high-frequency energy release per unit fault length. The high-frequency envelopes and acceleration spectra (>0.5 Hz) of the M 7.9 earthquake can be simulated by chaining together rupture zones of the M 6.7 earthquake over distances from 30 to 180 km east of the hypocenter. However, the inversion indicates that there was relatively little high-frequency energy generated along the 60-km portion of the Totschunda fault on the east end of the rupture.
With the implementation of the USGS National Earthquake Information Center Prompt Assessment of Global Earthquakes for Response system (PAGER), rapid determination of earthquake moment magnitude is essential, especially for earthquakes that are felt within the contiguous United States. We report an implementation of moment tensor processing for application to broad, seismically active areas of North America. This effort focuses on the selection of regional crustal velocity models, codification of data quality tests, and the development of procedures for rapid computation of the seismic moment tensor. We systematically apply these techniques to earthquakes with reported magnitude greater than 3.5 in continental North America that are not associated with a tectonic plate boundary. Using the 0.02-0.10 Hz passband, we can usually determine, with few exceptions, moment tensor solutions for earthquakes with M w as small as 3.7. The threshold is significantly influenced by the density of stations, the location of the earthquake relative to the seismic stations and, of course, the signal-to-noise ratio. With the existing permanent broadband stations in North America operated for rapid earthquake response, the seismic moment tensor of most earthquakes that are M w 4 or larger can be routinely computed. As expected the nonuniform spatial pattern of these solutions reflects the seismicity pattern. However, the orientation of the direction of maximum compressive stress and the predominant style of faulting is spatially coherent across large regions of the continent.