Search USGS⌕ Search

SEARCH · Search USGS

Results for “Physics Bulletin”

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 163 records · Page 9Linked to original sources
USGS1978

Differential use of fresh water environments by wintering waterfowl of coastal Texas

A comparative study of the environmental relationships among 14 species of wintering waterfowl was conducted at the Welder Wildlife Foundation, San Patricia County, near Sinton, Texas during the fall and early winter of 1973. Measurements of 20 environmental factors (social, vegetational, physical, and chemical) were subjected to multivariate statistical methods to determine certain niche characteristics and environmental relationships of waterfowl wintering in the aquatic community.....Each waterfowl species occupied a unique realized niche by responding to distinct combinations of environmental factors identified by principal component analysis. One percent confidence ellipses circumscribing the mean scores plotted for the first and second principal components gave an indication of relative niche width for each species. The waterfowl environments were significantly different interspecifically and water depth at feeding site and % emergent vegetation were most important in the separation. This was shown by subjecting the transformed data to multivariate analysis of variance with an associated step-down procedure. The species were distributed along a community cline extending from shallow water with abundant emergent vegetation to open deep water with little emergent vegetation of any kind. Four waterfowl subgroups were significantly separated along the cline, as indicated by one-way analysis of variance with Duncan?s multiple range test. Clumping of the bird species toward the middle of the available habitat hyperspace was shown in a plot of the principal component scores for the random samples and individual species.....Naturally occurring relationships among waterfowl were clarified using principal comcomponent analysis and related multivariate procedures. These techniques may prove useful in wetland management for particular groups of waterfowl based on habitat preferences.

The Wilson Bulletin↗
USGS2025

Development and evaluation of the remote passive integrated transponder tag reader for customizable monitoring of wildlife

Increasing the quantity and quality of mark-reencounter data can be important when individual identification data are needed to address research or management questions. Physical recapture methods tend to be labor-intensive and therefore expensive, while remote detection methods are not easily applied to all taxa or are proprietary and thus difficult to customize. However, partnerships between wildlife scientists and engineers have the potential to provide innovative solutions that improve data collection while reducing costs. We describe a collaborative effort to improve the collection of mark-reencounter data on an invasive reptile, the brown treesnake ( Boiga irregularis ). In laboratory trials conducted on Guåhan (CHamoru; Guam in English) in 2021, we assessed the feasibility of detecting snakes fitted with passive integrated transponder (PIT) tags using customized remote PIT-tag readers (RePTaR) paired with baited snake traps. We evaluated differences in scanning success as a function of snake traits, location of PIT tags within snakes, and the distance of the PIT tag from the reader. We successfully scanned all individuals, on average 529 times but ranging from 3 to 6,436 times, during 12-hour trials using a designed customizable RePTaR unit. Snake characteristics and marking location explained little of the variability in scanning success, with distance from the reader as the best-supported covariate. By building cross-disciplinary partnerships, we can develop customizable technological solutions and provide additional tools with which to address challenges in wildlife research and management.

Guam↗
USGS2025

New insights into gas-driven phase segregation in andesitic enclaves from Mt. Mazama (Crater Lake), USA

A key process in active magmatic systems is the “recharge” of deep-sourced mafic magma into cooler, more evolved, and crystal-rich shallow reservoirs; recharge may be the cause of, or response to, eruptive activity. Although compositional evidence for recharge has been extensively documented, physical models of recharge are limited, particularly processes that separate exsolving volatiles and melts from rapidly growing crystals. To improve constraints on phase separation behaviors, we re-examine andesitic enclaves in silicic andesite lava flows of Mt. Mazama (Crater Lake), USA, that provided early evidence of gas-driven filter pressing (Bacon, 1986). 2D and 3D imaging shows that enclaves have a sample-spanning crystal framework that is disrupted by melt patches, indicating that initially deformable crystal networks were subject to early phase reorganization. Small enclaves are poorly vesicular and require early gas loss. Large enclaves have porous cores with angular (diktytaxitic) voids that are well-connected in 3D and denser rinds with isolated pores. Large enclave rinds have similar bulk compositions to small enclaves but their less evolved cores require ~ 20% melt removal. In the large enclave, diktytaxitic core textures and gas fingering structures at the core–rind boundary suggest relatively slow late-stage outward gas migration. Both scaling arguments and evidence of outward gas/melt migration require a resistant rind. Rind formation is best explained by differential cooling and demonstrates the importance of thermal gradients for gas-driven filter pressing. A corollary is a limited time scale of recharge, enclave formation, and vesiculation to produce diktytaxitic textures, suggesting that recharge was (near) synchronous with eruption.

Oregon↗
USGS2017

Uncertainty, variability, and earthquake physics in ground‐motion prediction equations

Residuals between ground‐motion data and ground‐motion prediction equations (GMPEs) can be decomposed into terms representing earthquake source, path, and site effects. These terms can be cast in terms of repeatable (epistemic) residuals and the random (aleatory) components. Identifying the repeatable residuals leads to a GMPE with reduced uncertainty for a specific source, site, or path location, which in turn can yield a lower hazard level at small probabilities of exceedance. We illustrate a schematic framework for this residual partitioning with a dataset from the ANZA network, which straddles the central San Jacinto fault in southern California. The dataset consists of more than 3200 1.15≤ M ≤3 earthquakes and their peak ground accelerations (PGAs), recorded at close distances ( R ≤20 km). We construct a small‐magnitude GMPE for these PGA data, incorporating V S 30 site conditions and geometrical spreading. Identification and removal of the repeatable source, path, and site terms yield an overall reduction in the standard deviation from 0.97 (in ln units) to 0.44, for a nonergodic assumption, that is, for a single‐source location, single site, and single path. We give examples of relationships between independent seismological observables and the repeatable terms. We find a correlation between location‐based source terms and stress drops in the San Jacinto fault zone region; an explanation of the site term as a function of kappa, the near‐site attenuation parameter; and a suggestion that the path component can be related directly to elastic structure. These correlations allow the repeatable source location, site, and path terms to be determined a priori using independent geophysical relationships. Those terms could be incorporated into location‐specific GMPEs for more accurate and precise ground‐motion prediction.

Bulletin of the Seismological Society of America↗
USGS1996

The initial cooling of pahoehoe flow lobes

In this paper we describe a new thermal model for the initial cooling of pahoehoe lava flows. The accurate modeling of this initial cooling is important for understanding the formation of the distinctive surface textures on pahoehoe lava flows as well as being the first step in modeling such key pahoehoe emplacement processes as lava flow inflation and lava tube formation. This model is constructed from the physical phenomena observed to control the initial cooling of pahoehoe flows and is not an empirical fit to field data. We find that the only significant processes are (a) heat loss by thermal radiation, (b) heat loss by atmospheric convection, (c) heat transport within the flow by conduction with temperature and porosity-dependent thermal properties, and (d) the release of latent heat during crystallization. The numerical model is better able to reproduce field measurements made in Hawai'i between 1989 and 1993 than other published thermal models. By adjusting one parameter at a time, the effect of each of the input parameters on the cooling rate was determined. We show that: (a) the surfaces of porous flows cool more quickly than the surfaces of dense flows, (b) the surface cooling is very sensitive to the efficiency of atmospheric convective cooling, and (c) changes in the glass forming tendency of the lava may have observable petrographic and thermal signatures. These model results provide a quantitative explanation for the recently observed relationship between the surface cooling rate of pahoehoe lobes and the porosity of those lobes (Jones 1992, 1993). The predicted sensitivity of cooling to atmospheric convection suggests a simple field experiment for verification, and the model provides a tool to begin studies of the dynamic crystallization of real lavas. Future versions of the model can also be made applicable to extraterrestrial, submarine, silicic, and pyroclastic flows.

Bulletin of Volcanology↗
USGS1994

Correlation of basinal carbonate cycles to nearshore parasequences in the Late Cretaceous Greenhorn seaway, Western Interior U.S.A.

Upper Cretaceous limestone-shale couplets developed within the late transgressive stage of the Greenhorn cyclothem may be correlated from carbonate-dominated (basinal) sequences in central Kansas and Colorado westward to clastic cycles in southern Utah. Six such basinal couplets have been traced to corresponding upward-coarsening progradational cycles developed on the western margin of the Western Interior basin. In the central basin in Colorado and Kansas, these sedimentary cycles are represented by limestone-shale and marlstone-shale couplets ∼0.5-1.0 m in thickness. More calcareous parts of these couplets may be correlated westward into condensed, fossiliferous concretion and shell beds in proximal offshore lithofacies of Arizona and Utah. These concretion and shell beds are physically traceable farther landward (westward) into bioturbated, fossil-rich, transgressive lag deposits that bound 10- to 20-m-thick coarsening-upward progradational strand-plain deposits (parasequences) in southwestern Utah. Thus, the progra-dational phase of parasequence deposition correlates with accumulation of clay-rich sediment in the central basin, and the transgressive phase is characterized by reduced terrigenous input and deposition of carbonate-rich sediment. We consider Milankovitch-style orbital forcing of climate and tectonically induced fluctuations in rates of foredeep basin subsidence as possible forcing mechanisms for these basinwide events. Based on the widespread distribution of the limestone-shale couplets, as well as on estimated sedimentation rates and geochronology, it has been widely speculated that these carbonate cycles reflect Milankovitch cycles with periodicities on the order of 20 k.y. to 100 k.y. If so, then stratigraphic data suggest that orbital forcing of climate affected eustasy and/or sediment input and biogenic production in the Western Interior basin. Alternatively, thrusting events in the Sevier orogenic belt may have produced episodic changes in the rates of foredeep basin subsidence and consequent changes in base level, which could have controlled the deposition of the Greenhorn parasequences and carbonate cycles. In either case, correlation of these units demonstrates a consistent basinwide sedimentary response to high-frequency base level or sediment input changes in the Western Interior epicontinental basin.

Arizona, Colorado, Kansas, Utah↗
USGS1960

Beryl-bearing pegmatites in the Ruby Mountains and other areas in Nevada and northwestern Arizona

Pegmatite occurs widely in Nevada and northwestern Arizona, but little mining has been done for such pegmatite minerals as mica, feldspar, beryl, and lepidolite. Reconnaissance for beryl-bearing pegmatite in Nevada and in part of Mohave County, Ariz., and detailed studies in the Dawley Canyon area, Elko County, Nev., have shown that beryl occurs in at least 11 districts in the region. Muscovite has been prospected or mined in the Ruby and Virgin Mountains, Nev., and in Mohave County, Ariz. Feldspar has been mined in the southern part of the region near Kingman, Ariz., and in Clark County, Nev. The pegmatites in the region range in age from Precambrian to late Mesozoic or Tertiary. Among the pegmatite minerals found or reported in the districts studied are beryl, chrysoberyl, scheelite, wolframite, garnet, tourmaline, fluorite, apatite, sphene, allanite, samarskite, euxenite, gadolinite, monazite, autunite, columbite-tantalite, lepidolite, molybdenite, and pyrite and other sulflde minerals. The principal beryl-bearing pegmatites examined are in the Oreana and Lakeview (Humboldt Canyon) areas, Pershing County; the Dawley Canyon area in the Ruby Mountains, Elko County, Nev.; and on the Hummingbird claims in the Virgin Mountains, Mohave County, Ariz. Beryl has also been reported in the Marietta district, Mineral County; the Sylvania district, Esmeralda County; near Crescent Peak and near Searchlight, Clark County, Nev.; and in the Painted Desert near Hoover Dam, Mohave County, Ariz. Pegmatites are abundant in the Ruby Mountains, chiefly north of the granite stock at Harrison Pass. In the Dawley Canyon area of 2.6 square miles at least 350 pegmatite dikes more than 1 foot thick were mapped, and beryl was found in small quantities in at least 100 of these dikes. Four of these dikes exceed 20 feet in thickness, and 1 is 55 feet thick. A few pegmatites were also examined in the Corral Creek, Gilbert Canyon, and Hankins Canyon areas in the Ruby Mountains. The pegmatite dikes in the Dawley Canyon area intrude granite and metamorphic rocks which consist chiefly of quartzite and schist of probable Early Cambrian age. The granite is of two types: a biotite-muscovite granite that forms the main mass of the stock and albite granite that occurs in the metamorphic rocks near the borders of the stock. The pegmatites were emplaced chiefly along fractures in the granite and along schistosity or bedding planes in the metamorphic rocks. Many of the Dawley Canyon pegmatite dikes are zoned, having several rock units of contrasting mineralogy or grain size formed successively from the walls inward. Aplitic units occur either as zones or in irregular positions in the pegmatite dikes and are a distinctive feature of the Dawley Canyon pegmatites. Some of the aplitic and fine-grained pegmatite units are characterized by thin layers of garnet crystals, forming many parallel bands on outcrop surfaces. The occurrence of aplitic and pegmatitic textures in the same dike presumably indicates abrupt changes in physical-chemical conditions during crystallization, such as changes in viscosity and in content of volatile constituents. Concentrations of 0.1 percent or more beryl, locally more than 1 percent, occur in certain zones in the Dawley Canyon pegmatites. Spectrographic analyses of 23 samples indicate that the BeO content ranges from 0.0017 to 0.003 percent in the albite granite, from ,0.0013 to 0.039 percent in aplitic units in pegmatite, from 0.0005 to 0.10 percent in coarse-grained pegmatite, and from less than 0.0001 to 0.0004 percent in massive quartz veins. The scheelite-beryl deposits at Oreana and in Humboldt Canyon, Pershing County, are rich in beryllium. Twelve samples from the Lakeview (Humboldt Canyon) deposit range from 0.018 to 0.11 percent BeO, but underground crosscuts have failed to intersect similar rock at depth. Beryl locally constitutes as much as 10 percent of the pegmatitic ore at Oreana. The beryl was not recovered during tungsten mining at Oreana and is now in the tailings of the mill at Toulon, Nev. The percentage of beryl is lower than the Oreana ore because of dilution by tailings from other ores milled at Toulon. Beryl has been found in many pegmatite dikes in the Virgin Mountains. Both beryl and chrysoberyl occur in dikes on the Hummingbird claims, north of Virgin Peak, in Mohave County, Ariz. Spectrographic analyses of 5 representative samples of the principal dike on the Hummingbird claims range from 0.055 to 0.11 percent BeO.

Arizona, Nevada↗
USGS2014

Geologic and physiographic controls on bed-material yield, transport, and channel morphology for alluvial and bedrock rivers, western Oregon

The rivers of western Oregon have diverse forms and characteristics, with channel substrates ranging from continuous alluvial gravel to bare bedrock. Analysis of several measurable morphologic attributes of 24 valley reaches on 17 rivers provides a basis for comparing nonalluvial and alluvial channels. Key differences are that alluvial reaches have greater bar area, greater migration rates, and show systematic correlation among variables relating grain size to bed-material transport capacity. We relate these differences between channel types to bed-material transport rates as derived from a coupled regional analysis of empirical sediment yield measurements and physical experiments of clast attrition during transport. This sediment supply analysis shows that overall bed-material transport rates for western Oregon are chiefly controlled by (1) lithology and basin slope, which are the key factors for bed-material supply into the stream network, and (2) lithologic control of bed-material attrition from in-transport abrasion and disintegration. This bed-material comminution strongly affects bed-material transport in the study area, reducing transport rates by 50%–90% along the length of the larger rivers in the study area. A comparison of the bed-material transport estimates with the morphologic analyses shows that alluvial gravel-bed channels have systematic and bounding relations between bed-material transport rate and attributes such as bar area and local transport capacity. By contrast, few such relations are evident for nonalluvial rivers with bedrock or mixed-bed substrates, which are apparently more influenced by local controls on channel geometry and sediment supply. At the scale of western Oregon, the physiographic and lithologic controls on the balance between bed-material supply and transport capacity exert far-reaching influence on the distribution of alluvial and nonalluvial channels and their consequently distinctive morphologies and behaviors—differences germane for understanding river response to tectonics and environmental perturbations, as well as for implementing effective restoration and monitoring strategies.

California, Nevada, Oregon, Washington↗
USGS2020

Repeatable source, path, and site effects from the 2019 Ridgecrest M7.1 earthquake sequence

We use a large instrumental dataset from the 2019 Ridgecrest earthquake sequence ( Rekoske et al. , 2019 , 2020 ) to examine repeatable source‐, path‐, and site‐specific ground motions. A mixed‐effects analysis is used to partition total residuals relative to the Boore et al. (2014 ; hereafter, BSSA14) ground‐motion model. We calculate the Arias intensity stress drop for the earthquakes and find strong correlation with our event terms, indicating that they are consistent with source processes. We look for physically meaningful trends in the partitioned residuals and test the ability of BSSA14 to capture the behavior we observe in the data. We find that BSSA14 is a good match to the median observations for M > 4 "> M > 4 M>4 ⁠ . However, we find bias for individual events, especially those with small magnitude and hypocentral depth ≥ 7    km "> depth ≥ 7 km depth≥7 km ⁠ , for which peak ground acceleration is underpredicted by a factor of 2.5. Although the site amplification term captures the median site response when all sites are considered together, it does not capture variations at individual stations across a range of site conditions. We find strong basin amplification in the Los Angeles, Ventura, and San Gabriel basins. We find weak amplification in the San Bernardino basin, which is contrary to simulation‐based findings showing a channeling effect from an event with a north–south azimuth. This and an additional set of ground motions from earthquakes southwest of Los Angeles suggest that there is an azimuth‐dependent southern California basin response related to the orientation of regional structures when ground motion from waves traveling south–north are compared with those in the east–west direction. These findings exhibit the power of large, spatially dense ground‐motion datasets and make clear that nonergodic models are a way to reduce bias and uncertainty in ground‐motion estimation for applications like the U.S. Geological Survey National Seismic Hazard Model and the ShakeAlert earthquake early warning System.

California↗
USGS2021

Solar-system-wide significance of Mars polar science

1. The North Polar Layered Deposits contain thousands of ice layers that record accumulation and climate history for at least several million years, making the most accessible and most complete climate record aside from the Earth’s – and the only one to record the impact of large obliquity shifts. 2. Mars Polar Science is a diverse and integrated system spanning much of the planet, above and below the surface: geographic poles, mid- and high-latitudes, numerous active volatile processes that interact with the near-surface, and the lower and middle atmosphere circulations, even influencing the upper atmosphere and atmospheric escape. 3. Mars polar processes serves as an analogue and scientific link to numerous solar system bodies. Mars has two volatile cycles that influence surface processes, atmospheric processes, and climate cycles. Mars, with many data sets and observations, serves as a proxy for outer planets and moons that undergo the same physics but with other volatile species. Additionally, technology designed for Mars polar exploration benefits exploration elsewhere. 4. Mars Polar and Atmospheric Science results are critical for future human exploration and in situ resources and the study of planetary habitability.

Bulletin of the AAS↗
USGS2020

A new technique to calculate earthquake stress transfer and to forecast aftershocks

Coseismic stress changes have been the primary physical principle used to explain aftershocks and triggered earthquakes. However, this method does not adequately forecast earthquake rates and diverse rupture populations when subjected to formal testing. We show that earthquake forecasts can be impaired by assumptions made in physics-based models, such as the existence of hypothetical optimal faults, and regional scale invariability of the stress field. We compare calculations made under these assumptions along with different realizations of a new conceptual triggering model that features a complete assay of all possible ruptures. In this concept, there always exists a set of theoretical planes that has positive failure stress conditions under a combination of background and coseismic static stress change. In the Earth, all of these theoretical planes may not exist, and if they do, they may not be ready to fail. Thus the actual aftershock plane may not correspond to the plane with the maximum stress change value. This is consistent with observations that mainshocks commonly activate faults with exotic orientations and rakes. Our testing ground is the M=7.2, 2010 El Mayor-Cucapah earthquake sequence that activated multiple diverse fault populations across the USA-Mexico border in California and Baja California. We carry out a retrospective test involving 748 M≥3.0 triggered earthquakes that occurred during a 3-yr period after the mainshock. We find that a probabilistic expression of possible aftershock planes constrained by pre-mainshock rupture patterns is strongly favoured (89% of aftershocks consistent with static stress triggering) versus an optimal fault implementation (35% consistent). Results show that coseismic stress change magnitudes do not necessarily control earthquake triggering, instead we find that the summed background stress and coseismic stress change promotes diverse ruptures. Our model can thus explain earthquake triggering in regions where optimal plane mapping shows coseismic stress reduction.

Bulletin of the Seismological Society of America↗
USGS2024

Developing a photography-based harvest survey to estimate age and subspecies composition of midcontinent sandhill cranes

Midcontinent sandhill cranes ( Antigone canadensis ) are managed as a single population, but hunting regulations are structured so harvest is targeted towards the more numerous lesser sandhill cranes ( A. c. canadensis ). However, research indicates that greater sandhill cranes ( A. c. tabida ) have been disproportionally exposed to harvest at a rate exceeding their proportion within the midcontinent population. In addition, harvest has increased 22% per year in the U.S. Central Flyway states. The midcontinent population appears to be growing in recent years, but variability in annual abundance estimates has increased substantially. With limited resources and harvest management uncertainty increasing, we developed methods for a citizen science, photography-based harvest survey to estimate age and subspecies composition of harvested midcontinent sandhill cranes. To develop survey methods, we collected physical parts from 284 sandhill cranes in North Dakota in 2019 and 2020. We manually measured the culmen and tarsus using calipers, and digitally measured these parts using photographs and computer software. All digitally derived measurements were 2.5% to 5.9% larger than manual measurements; therefore, we developed linear models that adjusted digital measurements, facilitating subspecies prediction using an existing morphometric-based technique. In 2021, we requested an equal number of hunters to participate using 2 data collection methods to test if hunters could reliably take photographs suitable for digital measurement. Collection method 1 involved photographing the head and leg simultaneously, and Collection method 2 involved photographing the head only. Hunters submitted a total of 239 photographs. Only 80 of these photographs were submitted using Collection method 1, and 72% were suitable for digital measurement. Conversely, hunters submitted twice as many photographs using Collection method 2, and 88% of these photographs were deemed suitable. Although obtaining the tarsus measurement slightly improved subspecies predictability, Collection method 2 increased participation and usable data. We believe our results could be used to develop an operational survey of age and subspecies composition of midcontinent sandhill cranes, wherein a sample of crane hunters throughout the midcontinent population range would be asked to electronically submit a photograph of the head of each bird they harvest. A time series of age and subspecies composition of this population would provide managers with valuable information and improve harvest management at minimal additional cost and burden, compared to a traditional parts collection survey administered by mail.

Wildlife Society Bulletin↗
USGS2017

Alternative rupture-scaling relationships for subduction interface and other offshore environments

Alternative fault-rupture-scaling relationships are developed for Mw 7.1– 9.5 subduction interface earthquakes using a new database of consistently derived finitefault rupture models from teleseismic inversion. Scaling relationships are derived for rupture area, rupture length, rupture width, maximum slip, and average slip. These relationships apply width saturation for large-magnitude interface earthquakes (approximately Mw >8:6) for which the physical characteristics of subduction zones limit the depth extent of seismogenic rupture, and consequently, the down-dip limit of strong ground motion generation. On average, the down-dip rupture width for interface earthquakes saturates near 200 km (196 km on average). Accordingly, the reinterpretation of rupture-area scaling for subduction interface earthquakes through the use of a bilinear scaling model suggests that rupture asperity area is less well correlated with magnitude for earthquakes Mw >8:6. Consequently, the size of great-magnitude earthquakes appears to be more strongly controlled by the average slip across asperities. The sensitivity of the interface scaling relationships is evaluated against geographic region (or subduction zone) and average dip along the rupture interface to assess the need for correction factors. Although regional perturbations in fault-rupture scaling could be identified, statistical significance analyses suggest there is little rationale for implementing regional correction factors based on the limited number of interface rupture models available for each region. Fault-rupture-scaling relationships are also developed for intraslab (within the subducting slab), extensional outer-rise and offshore strike-slip environments. For these environments, the rupture width and area scaling properties yield smaller dimensions than interface ruptures for the corresponding magnitude. However, average and maximum slip metrics yield larger values than interface events. These observations reflect both the narrower fault widths and higher stress drops in these faulting environments. Although expressing significantly different rupture-scaling properties from earthquakes in subduction environments, the characteristics of offshore strike-slip earthquake ruptures compare similarly to commonly used rupture-scaling relationships for onshore strike-slip earthquakes.

Bulletin of the Seismological Society of America↗
USGS2023

Crustal block-controlled contrasts in deformation, uplift, and exhumation in the Santa Cruz Mountains, California, USA, imaged through apatite (U-Th)/He thermochronology and 3-D geological modeling

Deformation along strike-slip plate margins often accumulates within structurally partitioned and rheologically heterogeneous crustal blocks within the plate boundary. In these cases, contrasts in the physical properties and state of juxtaposed crustal blocks may play an important role in accommodation of deformation. Near the San Francisco Bay Area, California, USA, the Pacific−North American plate-bounding San Andreas fault bisects the Santa Cruz Mountains (SCM), which host numerous distinct, fault-bounded lithotectonic blocks that surround the San Andreas fault zone. In the SCM, a restraining bend in the San Andreas fault (the SCM bend) caused recent uplift of the mountain range since ca. 4 Ma. To understand how rheologic heterogeneity within a complex fault zone might influence deformation, we quantified plausible bounds on deformation and uplift across two adjacent SCM lithotectonic blocks on the Pacific Plate whose stratigraphic and tectonic histories differ. This was accomplished by combining 31 new apatite (U-Th)/He ages with existing thermochronological datasets to constrain exhumation of these two blocks. Additionally, surface exposures of the latest Miocene to late Pliocene Purisima Formation interpreted in 18 structural cross sections spanning the SCM allowed construction and restoration of Pliocene deformation in a three-dimensional geologic model. We found that rock uplift and deformation concentrated within individual Pacific Plate lithotectonic blocks in the SCM. Since 4 Ma, maximum principal strain computed for the more deformed block adjacent to the fault exceeded that computed for the less deformed block by at least 375%, and cumulative uplift has been more spatially extensive and higher in magnitude. We attribute the difference in uplift and deformation between the two blocks primarily to contrasts in lithotectonic structure, which resulted from diverging geologic histories along the evolving plate boundary.

California↗
USGS1985

General P, type-I S, and type-II S waves in anelastic solids; inhomogeneous wave fields in low-loss solids

The physical characteristics for general plane-wave radiation fields in an arbitrary linear viscoelastic solid are derived. Expressions for the characteristics of inhomogeneous wave fields, derived in terms of those for homogeneous fields, are utilized to specify the characteristics and a set of reference curves for general P and S wave fields in arbitrary viscoelastic solids as a function of wave inhomogeneity and intrinsic material absorption. The expressions show that an increase in inhomogeneity of the wave fields causes the velocity to decrease, the fractional-energy loss ( Q −1 ) to increase, the deviation of maximum energy flow with respect to phase propagation to increase, and the elliptical particle motions for P and type-I S waves to approach circularity. Q −1 for inhomogeneous type-I S waves is shown to be greater than that for type-II S waves, with the deviation first increasing then decreasing with inhomogeneity. The mean energy densities (kinetic, potential, and total), the mean rate of energy dissipation, the mean energy flux, and Q −1 for inhomogeneous waves are shown to be greater than corresponding characteristics for homogeneous waves, with the deviations increasing as the inhomogeneity is increased for waves of fixed maximum displacement amplitude. For inhomogeneous wave fields in low-loss solids, only the tilt of the particle motion ellipse for P and type-I S waves is independent to first order of the degree of inhomogeneity. Quantitative estimates for the characteristics of inhomogeneous plane body waves in layered low-loss solids are derived and guidelines established for estimating the effect of inhomogeneity on seismic body waves and a Rayleigh-type surface wave in low-loss media.

Bulletin of the Seismological Society of America↗
USGS2013

Compaction and gas loss in welded pyroclastic deposits as revealed by porosity, permeability, and electrical conductivity measurements of the Shevlin Park Tuff

Pyroclastic flows produced by large volcanic eruptions commonly densify after emplacement. Processes of gas escape, compaction, and welding in pyroclastic-flow deposits are controlled by the physical and thermal properties of constituent material. Through measurements of matrix porosity, permeability, and electrical conductivity, we provide a framework for understanding the evolution of pore structure during these processes. Using data from the Shevlin Park Tuff in central Oregon, United States, and from the literature, we find that over a porosity range of 0%–70%, matrix permeability varies by almost 10 orders of magnitude (from 10 –20 to 10 –11 m 2 ), with over three orders of magnitude variation at any given porosity. Part of the variation at a given porosity is due to permeability anisotropy, where oriented core samples indicate higher permeabilities parallel to foliation (horizontally) than perpendicular to foliation (vertically). This suggests that pore space is flattened during compaction, creating anisotropic crack-like networks, a geometry that is supported by electrical conductivity measurements. We find that the power law equation: k 1 = 1.3 × 10 –21 × ϕ 5.2 provides the best approximation of dominant horizontal gas loss, where k 1 = permeability, and ϕ = porosity. Application of Kozeny-Carman fluid-flow approximations suggests that permeability in the Shevlin Park Tuff is controlled by crack- or disk-like pore apertures with minimum widths of 0.3 and 7.5 μm. We find that matrix permeability limits compaction over short times, but deformation is then controlled by competition among cooling, compaction, water resorption, and permeable gas escape. These competing processes control the potential for development of overpressure (and secondary explosions) and the degree of welding in the deposit, processes that are applicable to viscous densification of volcanic deposits in general. Further, the general relationships among porosity, permeability, and pore geometry are relevant for flow of any fluid through an ignimbritic host.

Oregon↗
USGS2025

An empirical Green’s function approach for isolating directivity effects in earthquake ground-motion amplitudes

In this study, we apply an empirical Green’s function (eGf) method within a ground‐motion modeling framework to mitigate trade‐offs between source, path, and site effects. Many physical processes contribute to spatial variations in observed ground motions, including earthquake radiation pattern, directivity, variable path attenuation, and site effects. Current nonergodic ground‐motion models use spatially varying coefficients for path and site effects, but they do not address trade‐offs with complex earthquake source effects. To quantify the influence of directivity on ground‐motion amplitudes, we use records from multiple smaller earthquakes with epicenters near that of a larger event. We use these small magnitude events as eGfs and estimate repeatable path and site effects at individual stations, assuming that the average adjustments are not controlled by directivity. We adjust residuals from the larger earthquake using the eGf terms, isolating effects related to the rupture. This method clearly enhances the observed broadband directivity observed in the 2022 M 5.1 and 2007 M 5.4 Alum Rock earthquake ground motions, reinforcing the conclusion that their ruptures were unilateral. For the 2004 M 6.0 Parkfield earthquake, we find a bilateral rupture model better fits the data because variations in rupture velocity, slip rate, and slip distribution seem to have a stronger effect on the ground motions than rupture direction alone. Applying eGf adjustments reduces the standard deviation of the rupture models over the three earthquakes by 32% on average and by up to 57% for the 2022 Alum Rock earthquake, confirming we have effectively removed repeatable effects related to the wave propagation path and site response. We propose a novel measure of the frequency‐dependent directivity amplification strength as the reduction in ground‐motion residual variability gained by fitting a directivity model; for the three earthquakes considered, this parameter varies between 25% and 75%, indicating that directivity can strongly influence ground motions and should be considered in ground‐motion modeling.

California↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗