Search USGS⌕ Search

SEARCH · Search USGS

Results for “Exploration Geophysics”

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 541 records · Page 30Linked to original sources
USGS2021

Coseismic surface displacement in the 2019 ridgecrest earthquakes: Comparison of field measurements and optical image correlation results

A fundamental topic in earthquake studies is understanding the extent to which fault rupture at the surface is localized on primary fault strands as opposed to distributed tens to hundreds of meters away from primary ruptures through off‐fault deformation (OFD) via a combination of discrete secondary faulting and bulk deformation. The 2019 Ridgecrest, CA M w 6.4 and M w 7.1 earthquakes provide an opportunity to explore this problem via comparison of published field‐based and mostly on‐fault offset measurements, with new lateral displacement measurements made over length scales of hundreds of meters across the primary rupture using WorldView satellite images collected before and after the earthquakes. A mean fit to the field observations underestimates net slip (pixel‐correlation results) by an average of 0.4 m along the M w 6.4 rupture (∼41% of net left‐lateral displacement) and 0.7 m along the M w 7.1 rupture (∼65% of net right‐lateral displacement). We attribute these differences to substantial OFD along the lengths of the M w 6.4 (∼59%) and M w 7.1 (∼35%) ruptures. A maximum fit to the field observations provides an improved match to the pixel correlation results (difference of 0.1–0.2 m or ∼76%–98% of net displacement). OFD may in part depend on displacement magnitude, where locations along the rupture with smaller displacements (<2 m) are characterized by a greater percentage of OFD. Furthermore, OFD increases with rupture zone width, especially at the scale of individual rupture strands (hundreds of meters). These findings contribute to a growing understanding of the important role of OFD in accommodating deformation near the surface.

California↗
USGS2017

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗
USGS1969

Application of deep electrical soundings for groundwater exploration in Hawaii

Forty-five resistivity soundings, using Schlumberger and equatorial dipole electrode configurations, were made on the islands of Oahu and Hawaii to determine the applicability of direct current resistivity methods for locating freshwater aquifers in the State of Hawaii. The soundings were made on the northwestern part of the island of Oahu near the town of Waialua and on the island of Hawaii on the "saddle" area near Pohakuloa and Humuula. Interpretation of 32 sounding curves obtained on the island of Oahu indicates that it is possible to correlate five stratigraphic units underlain by a vesicular basalt basement and that the determination of the approximate depth to the freshwater-saline-water interface within the basalt is feasible. Two of these Schlumberger soundings with electrode spacings AB/2 reaching 6000 ft yielded sounding curves of the maximum and minimum types whose terminal branches asymptotically approach a resistivity of about 30 ohm, which is believed to be the true resistivity of basalt saturated with sea water. Near the town of Waialua the aquifer is a coral zone as well as parts of the weathered vesicular basalt basement. On the island of Hawaii, near Pohakuloa, an exploratory well drilled in basalt to a depth of 1001 ft (prior to the resistivity survey) proved to be dry. Interpretation of thirteen deep soundings made with Schlumberger and equatorial arrays suggests that the minimum depth to a conductive layer, which may represent basalt saturated with fresh water, is about 2700 ft below land surface. The groundwater appears to be dike impounded.

Hawaii↗
USGS2010

Mineralogy and chemistry of cobbles at Meridiani Planum, Mars, investigated by the Mars Exploration Rover Opportunity

Numerous loose rocks with dimensions of a few centimeters to tens of centimeters and with no obvious physical relationship to outcrop rocks have been observed along the traverse of the Mars Exploration Rover Opportunity. To date, about a dozen of these rocks have been analyzed with Opportunity's contact instruments, providing information about elemental chemistry (Alpha Particle X‐ray Spectrometer), iron mineralogy and oxidation states (Mössbauer Spectrometer) and texture (Microscopic Imager). These “cobbles” appear to be impact related, and three distinct groups can be identified on the basis of chemistry and mineralogy. The first group comprises bright fragments of the sulfate‐rich bedrock that are compositionally and texturally indistinguishable from outcrop rocks. All other cobbles are dark and are divided into two groups, referred to as the “Barberton group” and the “Arkansas group,” after the first specimen of each that was encountered by Opportunity. Barberton group cobbles are interpreted as meteorites with an overall chemistry and mineralogy consistent with a mesosiderite silicate clast composition. Arkansas group cobbles appear to be related to Meridiani outcrop and contain an additional basaltic component. They have brecciated textures, pointing to an impact‐related origin during which local bedrock and basaltic material were mixed.

Journal of Geophysical Research E: Planets↗
USGS1990

Schlumberger soundings near Medicine Lake, California

The use of direct current resistivity soundings to explore the geothermal potential of the Medicine Lake area in northern California proved to be challenging because of high contact resistances and winding roads. Deep Schlumberger soundings were made by expanding current electrode spacings along the winding roads. Measured apparent resistivities were corrected using the geometric factor for the exact array geometry instead of a linear array geometry. For horizontally stratified, laterally homogeneous media, the apparent resistivities measured with a nonlinear Schlumberger array are equal to those measured with a linear Schlumberger array provided that (a) distances from the current electrodes to the center of the array are equal, and (b) the proper geometric factor is used to calculate the apparent resistivity. Corrected sounding data were interpreted using an automatic interpretation method. Forty-two maps of interpreted resistivity were calculated for depths extending from 20 to 1000 m. Computer animation of these 42 maps revealed that (a) certain subtle anomalies migrate laterally with depth and can be traced to their origin, (b) an extensive volume of low-resistivity material underlies the survey area, and (c) the three areas (east of Bullseye Lake, southwest of Glass Mountain, and northwest of Medicine Lake) may be favorable geothermal targets. Six interpreted resistivity maps and three cross-sections illustrate the above findings.

Geophysics↗
USGS2017

A decade of induced slip on the causative fault of the 2015 Mw 4.0 Venus earthquake, northeast Johnson County, Texas

On 7 May 2015, a M w 4.0 earthquake occurred near Venus, northeast Johnson County, Texas, in an area of the Bend Arch-Fort Worth Basin that reports long-term, high-volume wastewater disposal and that has hosted felt earthquakes since 2009. In the weeks following the M w 4.0 earthquake, we deployed a local seismic network and purchased nearby active-source seismic reflection data to capture additional events, characterize the causative fault, and explore potential links between ongoing industry activity and seismicity. Hypocenter relocations of the resulting local earthquake catalog span ~4–6 km depth and indicate a fault striking ~230 ° , dipping to the west, consistent with a nodal plane of the M w 4.0 regional moment tensor. Fault plane solutions indicate normal faulting, with B axes striking parallel to maximum horizontal compressive stress. Seismic reflection data image the reactivated basement fault penetrating the Ordovician disposal layer and Mississippian production layer, but not displacing post-Lower Pennsylvanian units. Template matching at regional seismic stations indicates that low-magnitude earthquakes with similar waveforms began in April 2008, with increasing magnitude over time. Pressure data from five saltwater disposal wells within 5 km of the active fault indicate a disposal formation that is 0.9–4.8 MPa above hydrostatic. We suggest that the injection of 28,000,000 m 3 of wastewater between 2006 and 2015 at these wells led to an increase in subsurface pore fluid pressure that contributed to inducing this long-lived earthquake sequence. The 2015 M w 4.0 event represents the largest event in the continuing evolution of slip on the causative fault.

Texas↗
USGS2021

Geodetic constraints on a 25-year magmatic inflation episode near Three Sisters, central Oregon

Crustal inflation near the Three Sisters volcanic center documented since the mid-1990s has persisted for more than two decades. We update past analyses of the event through 2020 by simultaneously inverting InSAR interferograms, GPS time series, and leveling data for time-dependent volcanic deformation source parameters. We explore several source models to estimate how the deformation rate varied through time and to identify parameters that can reproduce measured deformation. Our preferred model is a Mogi source 4.1 km below sea level (5.9 km below the surface) about 5 km west of the summit of South Sister. Inflation started in late 1995 or 1996; the rate increased rapidly during 1998–1999, and peaked in late 1999, resulting in maximum surface uplift of about 30 cm by mid-2020. Since 2000, the inflation rate generally declined exponentially with a time constant of about 6 years. Two source inflation scenarios fit the data equally well. In the first, the crust surrounding the source is elastic and the net source-volume increase, which we attribute to persistent magma input, has been about 49 × 10 6 m 3 . The second scenario adds a viscoelastic shell surrounding the Mogi source. In that case, an injection of about 21 × 10 6 m 3 of magma prior to 2000, followed by continuing relaxation of the viscoelastic shell, can account for most of the observed surface deformation. In both scenarios, modeling reveals quasiperiodic increases in the inflation rate (pulses) with a recurrence interval of 3–4 years, both before and after 2000.

Oregon↗
USGS2009

Responses of stream nitrate and dissolved organic carbon loadings to hydrological forcing and climate change in an upland forest of the northeast USA

[1] In coming decades, higher annual temperatures, increased growing season length, and increased dormant season precipitation are expected across the northeastern United States in response to anthropogenic forcing of global climate. We synthesized long-term stream hydrochemical data from the Sleepers River Research Watershed in Vermont, United States, to explore the relationship of catchment wetness to stream nitrate and DOC loadings. We modeled changes in growing season length and precipitation patterns to simulate future climate scenarios and to assess how stream nutrient loadings respond to climate change. Model results for the 2070&ndash;2099 time period suggest that stream nutrient loadings during both the dormant and growing seasons will respond to climate change. During a warmer climate, growing season stream fluxes (runoff +20%, nitrate +57%, and DOC +58%) increase as more precipitation (+28%) and quick flow (+39%) occur during a longer growing season (+43 days). During the dormant season, stream water and nutrient loadings decrease. Net annual stream runoff (+8%) and DOC loading (+9%) increases are commensurate with the magnitude of the average increase of net annual precipitation (+7%). Net annual stream water and DOC loadings are primarily affected by increased dormant season precipitation. In contrast, decreased annual loading of stream nitrate (&minus;2%) reflects a larger effect of growing season controls on stream nitrate and the effects of lengthened growing seasons in a warmer climate. Our findings suggest that leaching of nitrate and DOC from catchment soils will be affected by anthropogenic climate forcing, thereby affecting the timing and magnitude of annual stream loadings in the northeastern United States.

Vermont↗
USGS2017

Delayed seismicity rate changes controlled by static stress transfer

On 15 June 2010, a M w 5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity.

California↗
USGS2012

Archive of digital boomer subbottom data collected during USGS cruise 05FGS01 offshore east-central Florida, July 17-29, 2005

In July of 2005, the U.S. Geological Survey (USGS), in cooperation with the Florida Geological Survey (FGS), conducted a geophysical survey of the Atlantic Ocean offshore of Florida's east coast from Flagler Beach to Daytona Beach. This report serves as an archive of unprocessed digital boomer subbottom data, trackline maps, navigation files, Geographic Information System (GIS) files, Field Activity Collection System (FACS) logs and formal Federal Geographic Data Committee (FGDC) metadata. Filtered and gained (showing a relative increase in signal amplitude) digital images of the seismic profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. The USGS Saint Petersburg Coastal and Marine Science Center (SPCMSC) assigns a unique identifier to each cruise or field activity. For example, 05FGS01 tells us the data were collected in 2005 for cooperative work with the FGS and the data were collected during the first field activity for that project in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity ID. The boomer subbottom processing system consists of an acoustic energy source that is made up of capacitors charged to a high voltage and discharged through a transducer in the water. The transducer is towed on a sled floating on the water surface and when discharged emits a short acoustic pulse, or shot, which propagates through the water column and shallow stratrigraphy below. The acoustic energy is reflected at density boundaries (such as the seafloor or sediment layers beneath the seafloor), detected by the receiver (a hydrophone streamer), and recorded by a PC-based seismic acquisition system. This process is repeated at timed intervals (for example, 0.5 s) and recorded for specific intervals of time (for example, 100 ms). In this way, a two-dimensional (2-D) vertical image of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the acquisition geometry. Refer to table 1 for a summary of acquisition parameters and table 2 for trackline statistics. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG Y format (Barry and others, 1975), except an ASCII format is used for the first 3,200 bytes of the card image header instead of the standard EBCDIC format. For a detailed description about the recorded trace headers, refer to the SEG Y Format page. The SEG Y files may be downloaded and processed with commercial or public domain software such as Seismic Unix (Cohen and Stockwell, 2005). See the How To Download SEG Y Data page for download instructions. The printable profiles provided here are GIF images that were processed and gained using SU software; refer to the Software page for links to example SU processing scripts. The processed SEG Y data were also exported to Chesapeake Technology, Inc. (CTI) SonarWeb software to produce a geospatially interactive version of the profile that allows the user to obtain a geographic location and depth from the profile for a given cursor position; this information is displayed in the status bar of the browser. Please note that clicking on the profile image switches it to "Expanded View" (a compressed image of the entire line) and cursor tracking is not available in this mode.

Florida↗
USGS2011

Extreme volcanism on Io: Latest insights at the end of Galileo era

Galileo has now completed 7 years exploring Jupiter. The spacecraft obtained breathtaking views of the four major satellites, and studied Jupiter's clouds and atmospheric composition, rings, small satellites, and magnetic field. It had five successful close flybys and many distant observations of Io. Scientists already knew from Voyager and Earth‐based astronomy that Io is by far the most volcanically active object in the solar system. Galileo has given us stunning color panoramas of Io's surface and unprecedented close views of erupting volcanoes (Figure 1) and the largest active flows observed anywhere. Among recent discoveries about Io, perhaps most astonishing since Voyager, is that some lavas possess emission temperatures greater than any lavas erupted on Earth today and possibly since the start of Earth's geologic history. The Io science community has identified three alternative interpretations of Io's hottest lavas: (1) ultramafic material similar to komatiite; (2) superheated lava; or (3) an ultra‐refractory substance deficient in silica and rich in Ca‐Al oxides.

Eos, Transactions, American Geophysical Union↗
USGS2020

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗
USGS2021

Warmer winters increase the biomass of phytoplankton in a large floodplain river

Winters are changing rapidly across the globe but the implications for aquatic productivity and food webs are not well understood. In addition, the degree to which winter dynamics in aquatic systems respond to large-scale climate versus ecosystem-level factors is unclear but important for understanding and managing potential changes. We used a unique winter data set from the Upper Mississippi River System to explore spatial and temporal patterns in phytoplankton biomass (chlorophyll a , CHL) and associated environmental covariates across 25 years and ∼1,500 river km. To assess the role of regional climate versus site-specific drivers of winter CHL, we evaluated whether there were coherent long-term CHL dynamics from north to south and across lotic-lentic areas. We then estimated the degree to which these patterns were associated with climate variability (i.e., the Multivariate El Nino-Southern Oscillation Index), winter severity (freezing degree days), river discharge, or site-specific environmental variables (ice depth, snow depth, and nutrient concentrations). We found that winter CHL was typically highest in ice-free reaches and backwater lakes, occasionally exceeding summer values. We did not find highly synchronous CHL dynamics across the basin, but instead show that temporal trends were independent among river reaches and lotic-lentic areas of the river. Moreover, after accounting for these spatial dynamics, we found that CHL was most responsive to winter air temperature, being consistently higher in years with warmer winters across the basin. These results indicate that although productivity dynamics are highly dynamic within large river ecosystems, changes in the duration and severity of winter may uniformly increase wintertime productivity.

Minnesota, Missouri, Illinois, Iowa, Wisconsin↗
USGS1998

Two-way coupling between Vesuvius eruptions and southern Apennine earthquakes, Italy, by elastic stress transfer

During the past 1000 years, eruptions of Vesuvius have often been accompanied by large earthquakes in the Apennines 50–60 km to the northeast. Statistical investigations had shown that earthquakes often preceded eruptions, typically by less than a decade, but did not provide a physical explanation for the correlation. Here, we explore elastic stress interaction between earthquakes and eruptions under the hypothesis that small stress changes can promote events when the Apennine normal faults and the Vesuvius magma body are close to failure. We show that earthquakes can promote eruptions by compressing the magma body at depth and opening suitably oriented near-surface conduits. Voiding the magma body in turns brings these same normal faults closer to Coulomb failure, promoting earthquakes. Such a coupling is strongest if the magma reservoir is a dike oriented normal to the regional extension axis, parallel to the Apennines, and the near-surface conduits and fissures are oriented normal to the Apennines. This preferred orientation suggests that the eruptions issuing from such fissures should be most closely linked in time to Apennine earthquakes. Large Apennine earthquakes since 1400 are calculated to have transferred more stress to Vesuvius than all but the largest eruptions have transferred to Apennine faults, which may explain why earthquakes more commonly lead than follow eruptions. A two-way coupling may thus link earthquakes and Vesuvius eruptions along a 100-km-long set of faults. We test the statistical significance of the earthquake-eruption correlation in the two-way coupling zone, and find a correlation significant at the 95% confidence level.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Alternative ways of using field-based estimates to calibrate ecosystem models and their implications for carbon cycle studies

Model-data fusion is a process in which field observations are used to constrain model parameters. How observations are used to constrain parameters has a direct impact on the carbon cycle dynamics simulated by ecosystem models. In this study, we present an evaluation of several options for the use of observations in modeling regional carbon dynamics and explore the implications of those options. We calibrated the Terrestrial Ecosystem Model on a hierarchy of three vegetation classification levels for the Alaskan boreal forest: species level, plant-functional-type level (PFT level), and biome level, and we examined the differences in simulated carbon dynamics. Species-specific field-based estimates were directly used to parameterize the model for species-level simulations, while weighted averages based on species percent cover were used to generate estimates for PFT- and biome-level model parameterization. We found that calibrated key ecosystem process parameters differed substantially among species and overlapped for species that are categorized into different PFTs. Our analysis of parameter sets suggests that the PFT-level parameterizations primarily reflected the dominant species and that functional information of some species were lost from the PFT-level parameterizations. The biome-level parameterization was primarily representative of the needleleaf PFT and lost information on broadleaf species or PFT function. Our results indicate that PFT-level simulations may be potentially representative of the performance of species-level simulations while biome-level simulations may result in biased estimates. Improved theoretical and empirical justifications for grouping species into PFTs or biomes are needed to adequately represent the dynamics of ecosystem functioning and structure.

Alaska↗
USGS2013

A fluid-driven earthquake swarm on the margin of the Yellowstone caldera

Over the past several decades, the Yellowstone caldera has experienced frequent earthquake swarms and repeated cycles of uplift and subsidence, reflecting dynamic volcanic and tectonic processes. Here, we examine the detailed spatial-temporal evolution of the 2010 Madison Plateau swarm, which occurred near the northwest boundary of the Yellowstone caldera. To fully explore the evolution of the swarm, we integrated procedures for seismic waveform-based earthquake detection with precise double-difference relative relocation. Using cross-correlation of continuous seismic data and waveform templates constructed from cataloged events, we detected and precisely located 8710 earthquakes during the three-week swarm, nearly four times the number of events included in the standard catalog. This high-resolution analysis reveals distinct migration of earthquake activity over the course of the swarm. The swarm initiated abruptly on January 17, 2010 at about 10 km depth and expanded dramatically outward (both shallower and deeper) over time, primarily along a NNW-striking, ~55º ENE-dipping structure. To explain these characteristics, we hypothesize that the swarm was triggered by the rupture of a zone of confined high-pressure aqueous fluids into a pre-existing crustal fault system, prompting release of accumulated stress. The high-pressure fluid injection may have been accommodated by hybrid shear and dilatational failure, as is commonly observed in exhumed hydrothermally affected fault zones. This process has likely occurred repeatedly in Yellowstone as aqueous fluids exsolved from magma migrate into the brittle crust, and it may be a key element in the observed cycles of caldera uplift and subsidence.

Wyoming↗
USGS2011

Analysis of dispersion and attenuation of surface waves in poroelastic media in the exploration-seismic frequency band

We analyse dispersion and attenuation of surface waves at free surfaces of possible vacuum/poroelastic media: permeable-‘open pore’, impermeable-‘closed pore’ and partially permeable boundaries, which have not been previously reported in detail by researchers, under different surface-permeable, viscous-damping, elastic and fluid-flowing conditions. Our discussion is focused on their characteristics in the exploration-seismic frequency band (a few through 200 Hz) for near-surface applications. We find two surface-wave modes exist, R 1 waves for all conditions, and R 2 waves for closed-pore and partially permeable conditions. For R 1 waves, velocities disperse most under partially permeable conditions and least under the open-pore condition. High-coupling damping coefficients move the main dispersion frequency range to high frequencies. There is an f 1 frequency dependence as a constant- Q model for attenuation at high frequencies. R 1 waves for the open pore are most sensitive to elastic modulus variation, but least sensitive to tortuosities variation. R 1 waves for partially permeable surface radiate as non-physical waves (Im( k ) < 0) at low frequencies. For R 2 waves, velocities are slightly lower than the bulk slow P 2 waves. At low frequencies, both velocity and attenuation are diffusive of f 1/2 frequency dependence, as P 2 waves. It is found that for partially permeable surfaces, the attenuation displays - f 1 frequency dependence as frequency increasing. High surface permeability, low-coupling damping coefficients, low Poisson′s ratios, and low tortuosities increase the slope of the - f 1 dependence. When the attenuation coefficients reach 0, R 2 waves for partially permeable surface begin to radiate as non-physical waves.

Geophysical Journal International↗
USGS2014

Preslip and cascade processes initiating laboratory stick slip

Recent modeling studies have explored whether earthquakes begin with a large aseismic nucleation process or initiate dynamically from the rapid growth of a smaller instability in a “cascade-up” process. To explore such a case in the laboratory, we study the initiation of dynamic rupture (stick slip) of a smooth saw-cut fault in a 76mm diameter cylindrical granite laboratory sample at 40–120MPa confining pressure. We use a high dynamic range recording system to directly compare the seismic waves radiated during the stick-slip event to those radiated from tiny (M _6) discrete seismic events, commonly known as acoustic emissions (AEs), that occur in the seconds prior to each large stick slip. The seismic moments, focal mechanisms, locations, and timing of the AEs all contribute to our understanding of their mechanics and provide us with information about the stick-slip nucleation process. In a sequence of 10 stick slips, the first few microseconds of the signals recorded from stick-slip instabilities are nearly indistinguishable from those of premonitory AEs. In this sense, it appears that each stick slip begins as an AE event that rapidly (~20 μs) grows about 2 orders of magnitude in linear dimension and ruptures the entire 150mm length of the simulated fault. We also measure accelerating fault slip in the final seconds before stick slip. We estimate that this slip is at least 98% aseismic and that it both weakens the fault and produces AEs that will eventually cascade-up to initiate the larger dynamic rupture.

Journal of Geophysical Research B: Solid Earth↗