USGS ScienceSearch

SEARCH · USGS Science

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 37 records · Page 2Linked to original sources

Combinatorial optimization of earthquake spatial distributions under minimum cumulative stress constraints

We determine optimal on‐fault earthquake spatial distributions using a combinatorial method that minimizes the long‐term cumulative stress resolved on the fault. An integer‐programming framework was previously developed to determine the optimal arrangement of a millennia‐scale earthquake sample that minimizes the misfit to a target slip rate determined from geodetic data. The resulting cumulative stress from just slip‐rate optimization, however, can greatly exceed fault strength estimates. Therefore, we add an objective function that minimizes cumulative stress and broad stress constraints to limit the solution space. We find that there is a trade‐off in the two objectives: minimizing the cumulative stress on a fault within fault strength limits concentrates earthquakes in specific areas of the fault and results in excursions from the target slip rate. Both slip‐rate and stress objectives can be combined in either a weighted or lexicographic (hierarchical) method. Using a combination of objectives, we demonstrate that a Gutenberg–Richter sample of earthquakes can be arranged on a constant slip‐rate finite fault with minimal stress and slip‐rate residuals. We apply this method to determine the optimal arrangement of earthquakes on the variable slip‐rate Nankai megathrust over 5000 yr. The sharp decrease in slip rate at the Tokai section of the fault results in surplus cumulative stress under all scenarios. Using stress optimization alone restricts this stress surplus to the northeast end of the fault at the expense of decreasing the slip rate away from the target slip rate at the southwest end of the fault. A combination of both slip‐rate and stress objectives provides an adequate fit to the data, although alternate model formulations for the fault are needed at the Tokai section to explain persistent excess cumulative stress. In general, incorporating stress objectives and constraints into the integer‐programming framework adds an important aspect of fault physics to the resulting earthquake rupture forecasts.

Bulletin of the Seismological Society of America

River meanders

Most river curves have nearly the same value of the ratio of curvature radius to channel width, in the range of 2 to 3. Meanders formed by meltwater on the surface of glaciers, and by the main current of the Gulf Stream, have a relation of meander length to channel width similar to rivers. Because such meanders carry no sediment, the shapes of curves in rivers are evidently determined primarily by the dynamics of flow rather than by relation to debris load. Velocity distributions along river curves provide a generalized picture of flow characteristics. Evidence on flow resistance in curved channels suggests that a basic aspect of meander mechanics may be related to the distribution of energy loss provided by a particular configuration or curvature. No general theory of meanders is as yet satisfactory, however; in fact, present evidence suggests that no single theory will explain the formation and characteristics of all meanders and that few of the physical principles involved have yet been clearly identified.

Geological Society of America Bulletin

Estimating the magnitude and frequency of floods at ungaged locations on streams in Tennessee through the 2013 water year

To improve estimates of the frequency of annual peak flows for ungaged locations on non-urban, unregulated streams in Tennessee, generalized least-squares multiple linear-regression techniques were used to relate annual peak flows from streamgages operated by the U.S. Geological Survey to physical, climatic, and land-use characteristics of their drainage basins. Geospatial data acquired since the previous study in 2003, annual peak-streamflow data through the 2013 water year, and Bulletin 17C methods for frequency analysis of annual peak-streamflow data were used in the study. Generalized least-squares regression equations were developed for four hydrologic areas with distinct hydrologic, geologic, and topographic characteristics. Drainage area was used as an explanatory variable in equations developed for each hydrologic area. In addition to drainage area, a 2-year recurrence-interval climate factor was used for hydrologic area 1, a 10–85 channel slope was used for hydrologic area 2, and percent imperviousness was used for hydrologic area 4. The regression equations can be used to estimate annual exceedance probability streamflows for ungaged locations on non-urban, unregulated streams in Tennessee. The term “unregulated” indicates that streamflow is not appreciably influenced by regulation from reservoirs or other impoundments. Average standard errors of prediction for the regression equations ranged from 44.4 to 51.4 percent for hydrologic area 1; 30.4 to 42.9 percent for hydrologic area 2; 35.7 to 42.6 percent for hydrologic area 3; and 32.5 to 47.4 percent for hydrologic area 4.

Tennessee

The size of earthquakes

How we should measure the size of an earthquake has been historically a very important, as well as a very difficult, seismological problem. For example, figure 1 shows the loss of life caused by earthquakes in recent times and clearly demonstrates that 1976 was the worst year for earthquake casualties in the 20th century. However, the damage caused by an earthquake is due not only to its physical size but also to other factors such as where and when it occurs; thus, figure 1 is not necessarily an accurate measure of the "size" of earthquakes in 1976. the point is that the physical process underlying an earthquake is highly complex; we therefore cannot express every detail of an earthquake by a simple straightforward parameter. Indeed, it would be very convenient if we could find a single number that represents the overall physical size of an earthquake. This was in fact the concept behind the Richter magnitude scale introduced in 1935.

Earthquake Information Bulletin (USGS)

Strain-estimated ground motions associated with recent earthquakes in California

Peak ground velocity (PGV) is a commonly used parameter in earthquake ground‐motion models (GMMs) and hazard analyses, because it is closely related to structural damage and felt ground shaking, and is typically measured on broadband seismometers. Here, we demonstrate that strainmeters, which directly measure in situ strain in the bulk rock, can easily be related to ground velocity by a factor of bulk shear‐wave velocity and, thus, can be used to measure strain‐estimated PGV. We demonstrate the parity of velocity to strain utilizing data from borehole strainmeters deployed along the plate boundaries of the west coast of the United States for nine recent M "> M M 4.4–7.1 earthquakes in California, including the largest two events of the July 2019 Ridgecrest earthquake sequence. PGVs derived from maximum horizontal shear strains fall within the range of seismic‐estimated values recorded at the same distances. We compare the strain‐estimated data with GMMs based on seismic PGVs and find consistency in residual polarity (positive vs. negative; the sign of the difference between observed and modeled data) for certain earthquake–station paths, where some paths indicate an overestimation and others indicate an underestimation of strain‐derived PGVs, as compared with the GMMs. We surmise that this may be indicative of over or underestimation of shear‐wave velocity along those paths, as compared with the average velocity used to derive PGV from strain measurements, or indicative of repeatable site and path effects that are not accounted for in our analyses. This direct comparison of strain with velocity can highlight physical path effects, as well as improve the density and capability of ground‐motion recordings.

California

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

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

California

Seismology in South America; an interview with Alberto Giesecke

Dr. Alberto A. Giesecke is head of the Instituto Geofisico del Peru, in Lima, Peru, and Director of Centro Regional de Sismologia para America del Sur (CERESIS). The center is dedicated to the coordination and promotion of earthquake hazard mitigation. Dr. Giesecke was President of the National Research Council of Peru and currently is a member of the Board of Directors of the National Institute for Industrial Technological Research and Standards and of the National Institute for Research and Training in Telecommunications. He presided over the Organizing Committee for the General Assemblies of the International Association for Seismology and Physics of the Earth's interor and the International Union for Radio Science held in Lima, Peru, in 1973 and 1975, respectively.

Earthquake Information Bulletin (USGS)

Remotely triggered seismicity on the United States west coast following the Mw 7.9 Denali fault earthquake

The M w 7.9 Denali fault earthquake in central Alaska of 3 November 2002 triggered earthquakes across western North America at epicentral distances of up to at least 3660 km. We describe the spatial and temporal development of triggered activity in California and the Pacific Northwest, focusing on Mount Rainier, the Geysers geothermal field, the Long Valley caldera, and the Coso geothermal field. The onset of triggered seismicity at each of these areas began during the Love and Raleigh waves of the M w 7.9 wave train, which had dominant periods of 15 to 40 sec, indicating that earthquakes were triggered locally by dynamic stress changes due to low-frequency surface wave arrivals. Swarms during the wave train continued for ∼4 min (Mount Rainier) to ∼40 min (the Geysers) after the surface wave arrivals and were characterized by spasmodic bursts of small ( M ≤ 2.5) earthquakes. Dynamic stresses within the surface wave train at the time of the first triggered earthquakes ranged from 0.01 MPa (Coso) to 0.09 MPa (Mount Rainier). In addition to the swarms that began during the surface wave arrivals, Long Valley caldera and Mount Rainier experienced unusually large seismic swarms hours to days after the Denali fault earthquake. These swarms seem to represent a delayed response to the Denali fault earthquake. The occurrence of spatially and temporally distinct swarms of triggered seismicity at the same site suggests that earthquakes may be triggered by more than one physical process.

North America

Mineralogy and petrology of the currant creek magnesite deposits and associated rocks of Nevada

Cryptocrystalline magnesite occurs as relatively small dense, white masses in a Tertiary volcanic tuff, here named the Currant tuff, in White Pine and Nye counties, 29 miles southwest of Ely, Nevada. The deposits are small, but some of the magnesite is of very high quality with almost no iron or aluminum. However, much of the magnesite is mixed with dolomite and calcite, and a serpentinelike mineral occurs in one of the deposits in sufficient abundance to be of possible use in ceramics. In scattered areas the tuff has been altered by solutions rich in magnesium and bicarbonic acid to dolomite, magnesite, magnesium silicate, and calcite. Mineral and chemical characteristics are the same throughout the entire rock assemblage of flows and tuffs, suggesting that all the rocks have been derived from the same parent magma. The tuff formation, which ranges in thickness from a few feet to over 400 feet, occurs between two groups of volcanic flows here referred to as the lower volcanics and the upper volcanics. The lower volcanics is composed mainly of flows of hypersthene dacite, but in the southwestern part of the area flows of a hypersthene andesite of balsaltic habit occur between the dacite and the overlying Currant tuff. The upper volcanics which overlie the Currant tuff consist of porphyritic quartz latite and an overlying massive latite crystal tuff. The Currant tuff containing the magnesite has both unmodified and reworked tuff members. The formation of hydrothermal dolomite, magnesite, and deweylite is discussed in detail, and the physical chemistry of solutions containing calcium, magnesium, and carbon dioxide is applied in an explanation of the origin of magnesium-rich solutions.

Nevada

Ground motion simulation for the 23 August 2011, Mineral, Virginia earthquake using physics-based and stochastic broadband methods

Three broadband simulation methods are used to generate synthetic ground motions for the 2011 Mineral, Virginia, earthquake and compare with observed motions. The methods include a physics‐based model by Hartzell et al. (1999, 2005), a stochastic source‐based model by Boore (2009), and a stochastic site‐based model by Rezaeian and Der Kiureghian (2010, 2012). The ground‐motion dataset consists of 40 stations within 600 km of the epicenter. Several metrics are used to validate the simulations: (1) overall bias of response spectra and Fourier spectra (from 0.1 to 10 Hz); (2) spatial distribution of residuals for GMRotI50 peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration (PSA) at various periods; (3) comparison with ground‐motion prediction equations (GMPEs) for the eastern United States. Our results show that (1) the physics‐based model provides satisfactory overall bias from 0.1 to 10 Hz and produces more realistic synthetic waveforms; (2) the stochastic site‐based model also yields more realistic synthetic waveforms and performs superiorly for frequencies greater than about 1 Hz; (3) the stochastic source‐based model has larger bias at lower frequencies (<0.5&thinsp;&thinsp;Hz) and cannot reproduce the varying frequency content in the time domain. The spatial distribution of GMRotI50 residuals shows that there is no obvious pattern with distance in the simulation bias, but there is some azimuthal variability. The comparison between synthetics and GMPEs shows similar fall‐off with distance for all three models, comparable PGA and PSA amplitudes for the physics‐based and stochastic site‐based models, and systematic lower amplitudes for the stochastic source‐based model at lower frequencies (<0.5&thinsp;&thinsp;Hz).

Bulletin of the Seismological Society of America

Seafloor observations eliminate a landslide as the source of the 1918 Puerto Rico Tsunami

The 11 October 1918 devastating tsunami in northwest Puerto Rico had been used as an example for earthquake‐induced landslide tsunami hazard. Three pieces of evidence pointed to a landslide as the origin of the tsunami: the discovery of a large submarine landslide scar from bathymetry data collected by shipboard high‐resolution multibeam sonar, reported breaks of submarine cable within the scar, and the fit of tsunami models to flooding observations. Newly processed seafloor imagery collected by remotely operated vehicle (ROV) show, however, pervasive Fe–Mn crust (patina) on the landslide walls and floor, indicating that the landslide scar is at least several hundred years old. C 14 "> 14 C dates of sediment covering the landslide floor verify this interpretation. Although we have not searched the region systematically for an alternative tsunami source, we propose a possible source—a two‐segment normal‐fault rupture along the eastern wall of Mona rift. The proposed fault location matches the published normal faults with steep bathymetry and is close to the International Seismological Center–Global Earthquake Model catalog locations of the 1918 mainshock and aftershocks. The ROV observations further show fresh vertical slickensides and rock exposure along the proposed fault trace. Hydrodynamic models from an M w "> M w 7.2 earthquake rupture along the eastern wall of the rift faithfully reproduce the reported tsunami amplitudes, polarities, and arrival times. Our analysis emphasizes the value of close‐up observations and physical samples to augment remote sensing data in natural hazard studies.

Puerto Rico