Search USGS⌕ Search

SEARCH · Search USGS

Results for “Applied Physics Letters”

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.

15 recordsLinked to original sources
USGS2013

Focused ion beam and field-emission microscopy of metallic filaments in memory devices based on thin films of an ambipolar organic compound consisting of oxadiazole, carbazole, and fluorene units

We report on the mechanism of operation of organic thin film resistive memory architectures based on an ambipolar compound consisting of oxadiazole, carbazole, and fluorene units. Cross-sections of the devices have been imaged by electron microscopy both before and after applying a voltage. The micrographs reveal the growth of filaments, with diameters of 50 nm–100 nm, on the metal cathode. We suggest that these are formed by the drift of aluminium ions from the anode and are responsible for the observed switching and negative differential resistance phenomena in the memory devices.

Applied Physics Letters↗
USGS2005

Seismic precursory patterns before a cliff collapse and critical point phenomena

We analyse the statistical pattern of seismicity before a 1-2 103 m3 chalk cliff collapse on the Normandie ocean shore, Western France. We show that a power law acceleration of seismicity rate and energy in both 40 Hz-1.5 kHz and 2 Hz-10kHz frequency range, is defined on 3 orders of magnitude, within 2 hours from the collapse time. Simultaneously, the average size of the seismic events increases toward the time to failure. These in situ results are derived from the only station located within one rupture length distance from the rock fall rupture plane. They mimic the "critical point" like behavior recovered from physical and numerical experiments before brittle failures and tertiary creep failures. Our analysis of this first seismic monitoring data of a cliff collapse suggests that the thermodynamic phase transition models for failure may apply for cliff collapse. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2024

Identifying Pareto-efficient eradication strategies for invasive populations

Invasive species are a major cause of biodiversity loss and are notoriously expensive and challenging to manage. We developed a decision-analytic framework for evaluating invasive species removal strategies, given objectives of maximizing eradication probability and minimizing costs. The framework uses an existing estimation model for spatially referenced removal data—one of the most accessible types of invasive species data—to obtain estimates of population growth rate, movement probability, and detection probability. We use these estimates in simulations to identify Pareto-efficient strategies—strategies where increases in eradication probability cannot be obtained without increases in cost—from a set of proposed strategies. We applied the framework post hoc to a successful eradication of veiled chameleons ( Chamaeleo calyptratus ) and identified the potential for substantial improvements in efficiency. Our approach provides managers and policymakers with tools to identify cost-effective strategies for a range of invasive species using only prior knowledge or data from initial physical removals.

Hawaii↗
USGS2012

Energy constraints on kinematic models of oblique faulting: Loma Prieta versus Parkfield-Coalinga

The Loma Prieta earthquake was a combination of right lateral strike-slip and reverse faulting. This oblique motion has been explained by a kinematic model that has the fault accommodate both the translation and convergence required by the Pacific-North American plate motions and a local bend in the San Andreas Fault However, in other areas of California, such as Parkfield-Coalinga, this combination of motions is taken up by decoupled strike-slip and thrust faults. In this paper I develop a simple physical model, based on applying Hamilton's principle to the energy released in earthquakes, that allows us to determine if a region should undergo oblique or decoupled faulting. This determination depends on the dip of available faults and the amount of convergence relative to the amount of translation. These energy constraints correctly predict the different types of faulting at Loma Prieta and Parkfield-Coalinga.

California↗
USGS2017

Characterizing local variability in long‐period horizontal tilt noise

Horizontal seismic data are dominated by atmospherically induced tilt noise at long periods (i.e., 30 s and greater). Tilt noise limits our ability to use horizontal data for sensitive seismological studies such as observing free earth modes. To better understand the local spatial variability of long‐period horizontal noise, we observe horizontal noise during quiet time periods in the Albuquerque Seismological Laboratory (ASL) underground vault using four small‐aperture array configurations. Each array comprises eight Streckeisen STS‐2 broadband seismometers. We analyze the spectral content of the data using power spectral density and magnitude‐squared coherence (γ2‐coherence). Our results show a high degree of spatial variability and frequency dependence in the long‐period horizontal wavefield. The variable nature of long‐period horizontal noise in the ASL vault suggests that it might be highly local in nature and not easily characterized by simple physical models when overall noise levels are low, making it difficult to identify locations in the vault with lower horizontal noise. This variability could be limiting our ability to apply coherence analysis for estimating horizontal sensor self‐noise and could also complicate various indirect methods for removing long‐period horizontal noise (e.g., collocated rotational sensor or microbarograph).

New Mexico↗
USGS2016

Scaling relation between earthquake magnitude and the departure time from P wave similar growth

We introduce a new scaling relation between earthquake magnitude ( M ) and a characteristic of initial P wave displacement. By examining Japanese K-NET data averaged in bins partitioned by M w and hypocentral distance, we demonstrate that the P wave displacement briefly displays similar growth at the onset of rupture and that the departure time ( T dp ), which is defined as the time of departure from similarity of the absolute displacement after applying a band-pass filter, correlates with the final M in a range of 4.5 ≤ M w ≤ 7. The scaling relation between M w and T dp implies that useful information on the final M can be derived while the event is still in progress because T dp occurs before the completion of rupture. We conclude that the scaling relation is important not only for earthquake early warning but also for the source physics of earthquakes.

Geophysical Research Letters↗
USGS2023

Magnitude conversion and earthquake recurrence rate models for the central and eastern United States

Development of Seismic Source Characterization (SSC) models, which is an essential part of Probabilistic Seismic Hazard Analyses (PSHA), can help forecast the temporal and spatial distribution of future damaging earthquakes (𝑀 w ≥ 5) in seismically active regions. Because it is impossible to associate all earthquakes with known faults, seismic source models for PSHA often include sources of diffuse seismicity in which future earthquake scenarios are not localized on mapped faults. These sources of diffuse seismicity are referred to as area source zones, distributed seismicity zones, or just source zones. During the early years of PSHA studies, it was assumed that earthquakes in seismotectonic zones have (1) uniform spatial distribution, (2) Poisson temporal distribution, and (3) exponential magnitude distribution (NRC, 2012). In seismically active regions (e.g., the Western United States), where active faults are readily identified, models of the spatial distribution of earthquakes include both the fault source geometries and the distributed seismicity (background) source zones. Source characterization of active faults is complemented by paleoseismic studies with estimates of earthquake magnitudes, dates of occurrences, and slip rates, which provide important information for PSHA studies. In the Central and Eastern United States (CEUS) very few Quaternary-active faults have the requisite information for use in PSHA (i.e., fault geometry and dimensions, event rates or slip rates, etc.), and we lack knowledge about the causative faults for most observed seismicity in the region. As a result, area source zones are frequently used in site-specific PSHA in the CEUS to represent diffuse seismicity that cannot be associated with faults. However, there are examples of active fault sources in the CEUS, such as the Meers fault, the Cheraw fault, and New Madrid region, where individual faults can be characterized. The source characterization models for background seismicity are based, to a large extent, on an assumption that spatial distribution of historical and recorded seismicity will not change substantially for time periods of interest for PSHA (approximately the next 50-100 years for engineered structures). Furthermore, studies such as those by Kafka (2007, 2009) found a correlation between the locations of small- to moderate-magnitude earthquakes and the locations of large-magnitude earthquakes, indicating that we can, with some level of confidence, use the spatial pattern of smaller earthquakes to forecast the future pattern of damaging earthquakes. Within background seismicity zones, the earthquake rate forecast is developed using spatial smoothing of the small to moderate magnitude events in earthquake catalogs. Different methodologies are used for this purpose and can predict varying distributions of seismicity rates. This in turn affects the results of a seismic hazard analysis. The U.S. Geological Survey (USGS) and Nuclear Regulatory Commission (NRC) use different methods for computing spatially smoothed seismicity rates in the CEUS; the USGS uses kernel-based spatial smoothing methods in developing the National Seismic Hazard Model (NSHM), and the method adopted in the Central and Eastern United States Seismic Source Characterization (CEUS-SSC) project is used when evaluating seismic hazard for nuclear power plant siting. These methods are described and the impact on seismic hazard are evaluated in this Research Information Letter (RIL). Another important input to estimating the rate of distributed seismicity is event magnitudes listed in earthquake catalogs. A substantial source of uncertainty in catalogs is the magnitude assigned to a given earthquake. Numerous different magnitude types exist, with each magnitude type computed in a different way. Therefore, for the sake of consistency, both the CEUS-SSC and the USGS NSHM have attempted to assemble a complete catalog with a uniform magnitude determination. To this end, moment magnitude, 𝑀 w , which is a physics-based measurement, has been adopted as the standard. However, 𝑀 w was not computed routinely until the past few decades. To address this issue, the CEUS-SSC conducted extensive analyses to determine conversion equations from which to take a routinely computed network (e.g., 𝑀 L or 𝑚 bLg ) and convert it into 𝑀 w . Another issue with using 𝑀 w is that it becomes increasingly difficult to compute for earthquakes with 𝑀 less than ~4. This study investigates the effects of moment magnitude estimation and spatial smoothing methods on estimation of the earthquake rate forecast and on seismic hazard. We investigate the validity of the magnitude conversion equations and their associated uncertainties by applying them to a case study for induced earthquakes in southern Kansas and northern Oklahoma, and summarize the use of the decay of the seismic coda to estimate 𝑀 w for small earthquakes (𝑀 w < 4. Furthermore, the study documents a comparison and assessment of background seismicity smoothing methods implemented by the USGS for the NSHM and used by the CEUS-SSC for siting nuclear facilities based on probabilistic seismic hazard estimates from multiple source zones in the CEUS and for multiple sites.

central and eastern United States↗
USGS1999

Dynamic 3D simulations of earthquakes on en echelon faults

One of the mysteries of earthquake mechanics is why earthquakes stop. This process determines the difference between small and devastating ruptures. One possibility is that fault geometry controls earthquake size. We test this hypothesis using a numerical algorithm that simulates spontaneous rupture propagation in a three-dimensional medium and apply our knowledge to two California fault zones. We find that the size difference between the 1934 and 1966 Parkfield, California, earthquakes may be the product of a stepover at the southern end of the 1934 earthquake and show how the 1992 Landers, California, earthquake followed physically reasonable expectations when it jumped across en echelon faults to become a large event. If there are no linking structures, such as transfer faults, then strike-slip earthquakes are unlikely to propagate through stepover s >5 km wide.

Geophysical Research Letters↗
USGS2018

Identifying physics‐based thresholds for rainfall‐induced landsliding

Most regional landslide warning systems utilize empirically derived rainfall thresholds that are difficult to improve without recalibration to additional landslide events. To address this limitation, we explored the use of synthetic rainfall to generate thousands of possible storm patterns and coupled them with a physics‐based hydrology and slope stability model for various antecedent soil saturation scenarios to analyze pore‐water pressure and factor of safety metrics. We used these metrics to generate two‐tiered alert thresholds that can be employed to assess shallow landslide potential for any given combination of storm and antecedent wetness. When applied to the San Francisco Bay region (California, USA), the results are consistent with events that caused widespread landsliding. Our deterministic modeling approach, which accounts for plausible ranges in soil hydraulic and mechanical properties, can inform the development of the next generation of warning systems for rainfall‐induced landsliding.

Geophysical Research Letters↗
USGS2005

The soil physics contributions of Edgar Buckingham

During 1902 to 1906 as a soil physicist at the USDA Bureau of Soils (BOS), Edgar Buckingham originated the concepts of matric potential, soil–water retention curves, specific water capacity, and unsaturated hydraulic conductivity ( K ) as a distinct property of a soil. He applied a formula equivalent to Darcy's law (though without specific mention of Darcy's work) to unsaturated flow. He also contributed significant research on quasi-empirical formulas for K as a function of water content, water flow in capillary crevices and in thin films, and scaling. Buckingham's work on gas flow in soils produced paradigms that are consistent with our current understanding. His work on evaporation elucidated the concept of self-mulching and produced sound and sometimes paradoxical generalizations concerning conditions that favor or retard evaporation. Largely overshadowing those achievements, however, is that he launched a theory, still accepted today, that could predict transient water content as a function of time and space. Recently discovered documents reveal some of the arguments Buckingham had with BOS officials, including the text of a two-paragraph conclusion of his famous 1907 report on soil water, and the official letter documenting rejection of that text. Strained interpersonal relations motivated the departure of Buckingham and other brilliant physicists (N.E. Dorsey, F.H. King, and Lyman Briggs) from the BOS during 1903 to 1906. Given that Buckingham and his BOS colleagues had been rapidly developing the means of quantifying unsaturated flow, these strained relations probably slowed the advancement of unsaturated flow theory.

Soil Science Society of America Journal↗
USGS2023

Fracture-mesh faulting in the swarm-like 2020 Maacama sequence revealed by high-precision earthquake detection, location, and focal mechanisms

In August of 2020, an earthquake sequence initiated within the Maacama fault zone in northern California, raising questions about its relationship with the larger-scale fault. To investigate the faulting geometry and its implications for physical processes driving seismicity, we applied an integrated, multi-faceted seismic analysis including waveform-correlation-based event detection, relative relocation, and high-precision focal mechanisms. To determine mechanisms for a large population of very small earthquakes (predominantly M < 1), we combined correlation-derived polarity analysis with a recently developed technique to determine S/P amplitude ratios from single seismic components. Finally, we applied an iterative stress inversion to distinguish between the likely fault and auxiliary planes. This analysis reveals that although the sequence initiated on a right-lateral fault, plausibly a strand of the Maacama Fault, it also activated numerous en echelon left-lateral conjugate faults. Together, these interlocking faults form a fracture mesh, consistent with structures associated elsewhere with fluid-induced faulting and earthquake swarms.

California↗
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↗
USGS2014

A hydrogeomorphic river network model predicts where and why hyporheic exchange is important in large basins

Hyporheic exchange has been hypothesized to have basin-scale consequences; however, predictions throughout river networks are limited by available geomorphic and hydrogeologic data and by models that can analyze and aggregate hyporheic exchange flows across large spatial scales. We developed a parsimonious but physically based model of hyporheic flow for application in large river basins: Networks with EXchange and Subsurface Storage (NEXSS). We applied NEXSS across a broad range of geomorphic diversity in river reaches and synthetic river networks. NEXSS demonstrates that vertical exchange beneath submerged bed forms rather than lateral exchange through meanders dominates hyporheic fluxes and turnover rates along river corridors. Per kilometer, low-order streams have a biogeochemical potential at least 2 orders of magnitude larger than higher-order streams. However, when biogeochemical potential is examined per average length of each stream order, low- and high-order streams were often found to be comparable. As a result, the hyporheic zone's intrinsic potential for biogeochemical transformations is comparable across different stream orders, but the greater river miles and larger total streambed area of lower order streams result in the highest cumulative impact from low-order streams. Lateral exchange through meander banks may be important in some cases but generally only in large rivers.

Geophysical Research Letters↗
USGS2012

Casual instrument corrections for short-period and broadband seismometers

Of all the filters applied to recordings of seismic waves, which include source, path, and site effects, the one we know most precisely is the instrument filter. Therefore, it behooves seismologists to accurately remove the effect of the instrument from raw seismograms. Applying instrument corrections allows analysis of the seismogram in terms of physical units (e.g., displacement or particle velocity of the Earth’s surface) instead of the output of the instrument (e.g., digital counts). The instrument correction can be considered the most fundamental processing step in seismology since it relates the raw data to an observable quantity of interest to seismologists. Complicating matters is the fact that, in practice, the term “instrument correction” refers to more than simply the seismometer. The instrument correction compensates for the complete recording system including the seismometer, telemetry, digitizer, and any anti‐alias filters. Knowledge of all these components is necessary to perform an accurate instrument correction. The subject of instrument corrections has been covered extensively in the literature (Seidl, 1980; Scherbaum, 1996). However, the prospect of applying instrument corrections still evokes angst among many seismologists—the authors of this paper included. There may be several reasons for this. For instance, the seminal paper by Seidl (1980) exists in a journal that is not currently available in electronic format and cannot be accessed online. Also, a standard method for applying instrument corrections involves the programs TRANSFER and EVALRESP in the Seismic Analysis Code (SAC) package (Goldstein et al. , 2003). The exact mathematical methods implemented in these codes are not thoroughly described in the documentation accompanying SAC.

Seismological Research Letters↗
USGS2011

Small-scale sediment transport patterns and bedform morphodynamics: New insights from high resolution multibeam bathymetry

New multibeam echosounder and processing technologies yield sub-meter-scale bathymetric resolution, revealing striking details of bedform morphology that are shaped by complex boundary-layer flow dynamics at a range of spatial and temporal scales. An inertially aided post processed kinematic (IAPPK) technique generates a smoothed best estimate trajectory (SBET) solution to tie the vessel motion-related effects of each sounding directly to the ellipsoid, significantly reducing artifacts commonly found in multibeam data, increasing point density, and sharpening seafloor features. The new technique was applied to a large bedform field in 20–30 m water depths in central San Francisco Bay, California (USA), revealing bedforms that suggest boundary-layer flow deflection by the crests where 12-m-wavelength, 0.2-m-amplitude bedforms are superimposed on 60-m-wavelength, 1-m-amplitude bedforms, with crests that often were strongly oblique (approaching 90°) to the larger features on the lee side, and near-parallel on the stoss side. During one survey in April 2008, superimposed bedform crests were continuous between the crests of the larger features, indicating that flow detachment in the lee of the larger bedforms is not always a dominant process. Assessment of bedform crest peakedness, asymmetry, and small-scale bedform evolution between surveys indicates the impact of different flow regimes on the entire bedform field. This paper presents unique fine-scale imagery of compound and superimposed bedforms, which is used to (1) assess the physical forcing and evolution of a bedform field in San Francisco Bay, and (2) in conjunction with numerical modeling, gain a better fundamental understanding of boundary-layer flow dynamics that result in the observed superimposed bedform orientation.

California↗