USGS Science⌕ Search

SEARCH · USGS Science

Results for “Near Surface 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 505 records · Page 28Linked to original sources

Isolated cases of remote dynamic triggering in Canada detected using cataloged earthquakes combined with a matched-filter approach

Here we search for dynamically triggered earthquakes in Canada following global main shocks between 2004 and 2014 with M S &thinsp;>&thinsp;6, depth&thinsp;<&thinsp;100&thinsp;km, and estimated peak ground velocity&thinsp;>&thinsp;0.2&thinsp;cm/s. We use the Natural Resources Canada (NRCan) earthquake catalog to calculate &beta; statistical values in 1&deg;&thinsp;&times;&thinsp;1&deg; bins in 10&thinsp;day windows before and after the main shocks. The statistical analysis suggests that triggering may occur near Vancouver Island, along the border of the Yukon and Northwest Territories, in western Alberta, western Ontario, and the Charlevoix seismic zone. We also search for triggering in Alberta where denser seismic station coverage renders regional earthquake catalogs with lower completeness thresholds. We find remote triggering in Alberta associated with three main shocks using a matched-filter approach on continuous waveform data. The increased number of local earthquakes following the passage of main shock surface waves suggests local faults may be in a critically stressed state.

Alberta↗

Possible deep fault slip preceding the 2004 Parkfield earthquake, inferred from detailed observations of tectonic tremor

Earthquake predictability depends, in part, on the degree to which sudden slip is preceded by slow aseismic slip. Recently, observations of deep tremor have enabled inferences of deep slow slip even when detection by other means is not possible, but these data are limited to certain areas and mostly the last decade. The region near Parkfield, California, provides a unique convergence of several years of high-quality tremor data bracketing a moderate earthquake, the 2004 magnitude 6.0 event. Here, I present detailed observations of tectonic tremor from mid-2001 through 2008 that indicate deep fault slip both before and after the Parkfield earthquake that cannot be detected with surface geodetic instruments. While there is no obvious short-term precursor, I find unidirectional tremor migration accompanied by elevated tremor rates in the 3 months prior to the earthquake, which suggests accelerated creep on the fault ∼16 km beneath the eventual earthquake hypocenter.

Geophysical Research Letters↗

Currents, drag, and sediment transport induced by a tsunami

We report observations of water surface elevation, currents, and suspended sediment concentration (SSC) from a 10-m deep site on the inner shelf in northern Monterey Bay during the arrival of the 2010 Chile tsunami. Velocity profiles were measured from 3.5 m above the bed (mab) to the surface at 2 min intervals, and from 0.1 to 0.7 mab at 1 Hz. SSC was determined from the acoustic backscatter of the near-bed profiler. The initial tsunami waves were directed cross shore and had a period of approximately 16 min. Maximum wave height was 1.1 m, and maximum current speed was 0.36 m/s. During the strongest onrush, near-bed velocities were clearly influenced by friction and a logarithmic boundary layer developed, extending more than 0.3 mab. We estimated friction velocity and bed shear stress from the logarithmic profiles. The logarithmic structure indicates that the flow can be characterized as quasi-steady at these times. At other phases of the tsunami waves, the magnitude of the acceleration term was significant in the near-bed momentum equation, indicating unsteady flow. The maximum tsunami-induced bed shear stress (0.4 N/m 2 ) exceeded the critical shear stress for the medium-grained sand on the seafloor. Cross-shore sediment flux was enhanced by the tsunami. Oscillations of water surface elevation and currents continued for several days. The oscillations were dominated by resonant frequencies, the most energetic of which was the fundamental longitudinal frequency of Monterey Bay. The maximum current speed (hourly-timescale) in 18 months of observations occurred four hours after the tsunami arrived.

California↗

Silica in a Mars analog environment: Ka u Desert, Kilauea Volcano, Hawaii

Airborne Visible/Near-Infrared Imaging Spectrometer (AVIRIS) data acquired over the Ka u Desert are atmospherically corrected to ground reflectance and used to identify the mineralogic components of relatively young basaltic materials, including 250-700 and 200-400 year old lava flows, 1971 and 1974 flows, ash deposits, and solfatara incrustations. To provide context, a geologic surface units map is constructed, verified with field observations, and supported by laboratory analyses. AVIRIS spectral endmembers are identified in the visible (0.4 to 1.2 ??m) and short wave infrared (2.0 to 2.5 ??m) wavelength ranges. Nearly all the spectral variability is controlled by the presence of ferrous and ferric iron in such minerals as pyroxene, olivine, hematite, goethite, and poorly crystalline iron oxides or glass. A broad, nearly ubiquitous absorption feature centered at 2.25 ??m is attributed to opaline (amorphous, hydrated) silica and is found to correlate spatially with mapped geologic surface units. Laboratory analyses show the silica to be consistently present as a deposited phase, including incrustations downwind from solfatara vents, cementing agent for ash duricrusts, and thin coatings on the youngest lava flow surfaces. A second, Ti-rich upper coating on young flows also influences spectral behavior. This study demonstrates that secondary silica is mobile in the Ka u Desert on a variety of time scales and spatial domains. The investigation from remote, field, and laboratory perspectives also mimics exploration of Mars using orbital and landed missions, with important implications for spectral characterization of coated basalts and formation of opaline silica in arid, acidic alteration environments. Copyright 2010 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Geologic implications of topographic, gravity, and aeromagnetic data in the northern Yukon-Koyukuk province and its borderlands, Alaska

The northern Yukon-Koyukuk province is characterized by low elevation and high Bouguer gravity and aeromagnetic anomalies in contrast to the adjacent Brooks Range and Ruby geanticline. Using newly compiled digital topographic, gravity, and aeromagnetic maps, I have divided the province into three geophysical domains. The Koyukuk domain, which is nearly equivalent to the Koyukuk lithotectonic terrane, is a horseshoe-shaped area, open to the south, of low topography, high gravity, and high-amplitude magnetic anomalies caused by an intraoceanic magmatic arc. The Angayucham and Kanuti domains are geophysical subdivisions of the Angayucham lithotectonic terrane that occur along the northern and southeastern margins of the Yukon-Koyukuk province, where oceanic rocks have been thrust over continental rocks of the Brooks Range and Ruby geanticline. Basalt of the Angayucham domain causes strong gravity highs and weak magnetic highs. The Kanuti domain is distinguished from the Angayucham domain by intense magnetic highs caused by cumulus mafic and ultramafic plutonic rocks, abundant ultramafic mantle tectonites, and magnetic syenite and monzonite. Long-wavelength, low-intensity magnetic highs and undulating gravity anomalies indicate an undulating basement surface of varied lithology beneath the Kobuk-Koyukuk and Lower Yukon basins. Modeling of gravity and magnetic anomalies shows that oceanic rocks of the Angayucham and Kanuti domains dip inward beneath the Kobuk-Koyukuk basin. The modeling supports, but does not prove, the hypothesis that the crust of the Kobuk-Koyukuk basin is 32–35 km thick, consisting of a tectonically thickened section of Cretaceous volcanic and sedimentary rocks and older oceanic crust. Plutons of the Brooks Range and the southern Ruby geanticline are nonmagnetic, ilmenite series, S-type granites that cause magnetic lows. Plutons of the northern Ruby geanticline are variable in their magnetic properties and cause both highs and lows. Plutons of both the eastern and western Yukon-Koyukuk province are variable in their magnetic expression but commonly cause magnetic lows in contrast to andesite.

Journal of Geophysical Research Solid Earth↗

Geophysical instrumentation near Parkfield

The geophysical instrumentation operated by the U.S Geological Survey and others near Parkfield is designed to monitor ongoing tectonic processes that generate earthquakes and to record the strong shaking that results from larger shocks and its effects. this discussion focuses on the former objectives; the latter is discussed in the next section "Ground Shaking and Engineering Studies on the Parkfield Section of the San Andreas Fault Zone." Because scientists expect the anticipated earthquake to resemble the historic Parkfield earthquakes, and in particular that in 1966, the data from the 1966 shock were used to site instruments for optimun benefit before, during, and after the next shock. the primary feature used for siting was the "1966 rupture zone," which is shown as the orange fault traces on the maps in this section. This zone defines the extent of surface tectonic cracks in 1966 and includes the source areas for fore shocks to the 1934 and 1966 earthquakes (north end of the zone) as well as for apparent precursory fault creep in 1966 (near center of the zone). Scientists believe that if precursors to the next shock are observed, they most likely will be near the 1966 rupture zone.

California↗

Reconnaissance stratigraphic studies in the Susitna basin, Alaska, during the 2014 field season

The Susitna basin is a poorly-understood Cenozoic successor basin immediately north of Cook Inlet in south-central Alaska (Kirschner, 1994). The basin is bounded by the Castle Mountain fault and Cook Inlet basin on the south, the Talkeetna Mountains on the east, the Alaska Range on the north, and the Alaska–Aleutian Range on the west (fig. 2-1). The Cenozoic fill of the basin includes coal-bearing nonmarine rocks that are partly correlative with Paleogene strata in the Matanuska Valley and Paleogene and Neogene formations in Cook Inlet (Stanley and others, 2013, 2014). Mesozoic sedimentary rocks are present in widely-scattered uplifts in and around the margins of the basin; these rocks differ significantly from Mesozoic rocks in the forearc basin to the south. Mesozoic strata in the Susitna region were likely part of a remnant ocean basin that preceded the nonmarine Cenozoic basin (Trop and Ridgway, 2007). The presence of coal-bearing strata similar to units that are proven source rocks for microbial gas in Cook Inlet (Claypool and others, 1980) suggests the possibility of a similar system in the Susitna basin (Decker and others, 2012). In 2011 the Alaska Division of Geological & Geophysical Surveys (DGGS) and Alaska Division of Oil and Gas, in collaboration with the U.S. Geological Survey, initiated a study of the gas potential of the Susitna basin (Gillis and others, 2013). This report presents a preliminary summary of the results from 14 days of helicopter-supported field work completed in the basin in August 2014. The goals of this work were to continue the reconnaissance stratigraphic work begun in 2011 aimed at understanding reservoir and seal potential of Tertiary strata, characterize the gas source potential of coals, and examine Mesozoic strata for source and reservoir potential

Alaska↗

Peak ground velocity spatial variability revealed by dense seismic array in southern California

Understanding and modeling variability of ground motion is essential for building accurate and precise ground motion prediction equations, which can net site‐specific characterization and reduced hazard levels. Here, we explore the spatial variability in peak ground velocity (PGV) at Sage Brush Flats along the San Jacinto Fault in southern California. We use data from a dense array (0.6 x 0.6 km 2 , 1,108 geophones, station spacings 10‐30 m) deployed in 2014 for ~1‐month. These data offer an opportunity to study small scale variability in this region. We examine 38 2≤M L ≤4.2 earthquakes within 200 km of the array. Fault strands and a small basin impact the ground motions, producing PGV variations up to 22% of the mean and a 40% reduction in P and S wave near‐surface velocities. We find along‐fault rupture directivity, source, and path effects can increase PGVs by 167%. Surface PGV measurements exceed the co‐located borehole station (depth 148 m) PGV by factors of 3‐10, confirming the impact on PGV from near surface fault structures, basins, topography, and amplifications from soft sediments. Consistently we find high PGVs within the basin structure. A pair of ~co‐located M L 2.6 events produce repeatable PGV values with similar spatial patterns. The average corner frequencies of these two events are 11‐16 Hz and viable measurements of stress drop can differ by 6.45MPa. Within this small array, the PGV values are variable implying spatial extrapolation of PGV to regions of known faults and basins, even across a small area, should be done with caution.

California↗

Named landforms of the World: A geomorphological and physiographic compilation

Prior to the current era of digital geomorphological mapping, global and regional-scale land surface characterization was advanced by qualitative interpretations that relied on human visualization aided by disciplinary knowledge of geophysical processes combined with extensive field study. In the early twentieth century, Fenneman proposed to devise systematic physiographic divisions of the United States and in 1916 produced what is still regarded as an authoritative map of these divisions. His physiographic regions were developed to provide context when describing land surface characteristics of smaller areas using well-known regional characteristics and descriptors. In 1968, geographer Richard E. Murphy published a large-format map of the “Landforms of the World” to fill a gap in the suite of standard classroom maps. In 1990, the British geomorphologist E. M. Bridges published World Geomorphology , providing the first global treatment and description of divisions, provinces, and sections—the same hierarchical land partitioning concepts that Fenneman used decades earlier. In the twenty-first century, geographic information systems (GIS) technologies are nearly ubiquitous, yet neither Murphy’s nor Bridges’s work existed as GIS data. To further illuminate their pioneering work, we (1) recompiled Murphy’s landforms as a spatial combination of modern existing data layers, and (2) used the recompiled Murphy’s landforms as a basis for the boundaries of the divisions, provinces, and sections described by Bridges. Our aggregation yields a new resource, Named Landforms of the World, version 2.0, which provides a reference-level, basemap-quality data layer that can significantly facilitate mapping, assessing, and understanding Earth surface features.

Annals of the AAG↗

Ground water beneath coastal bays of the Delmarva Peninsula: Ages and nutrients

To complement a large-scale geophysical investigation of occurrence and discharge of fresh water beneath Delaware, Maryland, and Virginia (Delmarva) coastal bays, we measured (1) salinity and nutrient concentrations in ground water samples from several offshore coring sites and (2) a suite of chemical and isotopic parameters, including age tracers, in ground water samples from a Delaware site. Samples were collected in a variety of Holocene and Plio-Pleistocene sediments in nearshore and offshore areas of the bays. Ground waters that were significantly fresher than overlying waters were found in plumes up to at least 15 m thick extending to more than 500 m offshore in some areas. Steep salinity and nutrient gradients occur within a few meters of the sediment surface in most locations studied. The zone of transition from deeper fresher waters to shallower brackish waters is generally thin near shore, but thickens and becomes more gradual offshore. Ground water ages at the Delaware site were mostly < 50 yr in both fresh waters and brackish waters up to 22 m below the bay bottom. Water chemistry and age data indicate that fresh water plumes beneath the estuary are active extensions of the surficial aquifer carrying nitrate from recharge areas on land, whereas brackish ground water surrounding the fresh water plumes is recharged beneath the estuary and contains ammonium and phosphate released by diagenesis of shallow estuarine sediments. Denitrification affects some of the fresh water nitrate before it mixes with brackish ground water or discharges to surface water.

Delaware, Maryland, Virginia↗

Geology and geophysics of the southern Raft River Valley geothermal area, Idaho, USA

The Raft River valley, near the boundary of the Snake River plain with the Basin and Range province, is a north-trending late Cenozoic downwarp bounded by faults on the west, south, and east. Pleistocene alluvium and Miocene-Pliocene tuffaceous sediments, conglomerate, and felsic volcanic rocks aggregate 2 km in thickness. Large gravity, magnetic, and total field resistivity highs probably indicate a buried igneous mass that is too old to serve as a heat source. Differing seismic velocities relate to known or inferred structures and to a suspected shallow zone of warm water. Resistivity anomalies reflect differences of both composition and degree of alteration of Cenozoic rocks. Resistivity soundings show a 2 to 5 ohm·m unit with a thickness of 1 km beneath a large part of the valley, and the unit may indicate partly hot water and partly clayey sediments. Observed self-potential anomalies are believed to indicate zones where warm water rises toward the surface. Boiling wells at Bridge, Idaho are near the intersection of north-northeast normal faults which have moved as recently as the late (?) Pleistocene, and an east-northeast structure, probably a right-lateral fault. Deep circulation of ground water in this region of relatively high heat flow and upwelling along faults is the probable cause of the thermal anomaly.

Idaho↗

Geology of the Upheaval Dome impact structure, southeast Utah

Two vastly different phenomena, impact and salt diapirism, have been proposed for the origin of Upheaval Dome, a spectacular scenic feature in southeast Utah. Detailed geologic mapping and seismic refraction data indicate that the dome originated by collapse of a transient cavity formed by impact. Evidence is as follows: (1) sedimentary strata in the center of the structure are pervasively imbricated by top-toward-the-center thrust faulting and are complexly folded as well; (2) top-toward-the-center normal faults are found at the perimeter of the structure; (3) clastic dikes are widespread; (4) the top of the underlying salt horizon is at least 500 m below the surface at the center of the dome, and there are no exposures of salt or associated rocks of the Paradox Formation in the dome to support the possibility that a salt diapir has ascended through it; and (5) planar microstructures in quartz grains, fantailed fracture surfaces (shatter surfaces), and rare shatter cones are present near the center of the structure. We show that the dome formed mainly by centerward motion of rock units along listric faults. Outcrop-scale folding and upturning of beds, especially common in the center, are largely a consequence of this motion. We have also detected some centerward motion of fault-bounded wedges resulting from displacements on subhorizontal faults that conjoin and die out within horizontal bedding near the perimeter of the structure. The observed deformation corresponds to the central uplift and the encircling ring structural depression seen in complex impact craters.

Utah↗

Source model for the Mw 6.7, 23 October 2002, Nenana Mountain earthquake (Alaska) from InSAR

The 23 October 2002 Nenana Mountain Earthquake (M w ∼ 6.7) occurred on the Denali Fault (Alaska), to the west of the M w ∼ 7.9 Denali Earthquake that ruptured the same fault 11 days later. We used 6 interferograms, constructed using radar images from the Canadian Radarsat-1 and European ERS-2 satellites, to determine the coseismic surface deformation and a source model. Data were acquired on ascending and descending satellite passes, with incidence angles between 23 and 45 degrees, and time intervals of 72 days or less. Modeling the event as dislocations in an elastic half space suggests that there was nearly 0.9 m of right-lateral strike-slip motion at depth, on a near-vertical fault, and that the maximum slip in the top 4 km of crust was less than 0.2 m. The Nenana Mountain Earthquake increased the Coulomb stress at the future hypocenter of the 3 November 2002, Denali Earthquake by 30–60 kPa.

Alaska↗

Structural controls on ground-water conditions and estimated aquifer properties near Bill Williams Mountain, Williams, Arizona

As of 1999, surface water collected and stored in reservoirs is the sole source of municipal water for the city of Williams. During 1996 and 1999, reservoirs reached historically low levels. Understanding the ground-water flow system is critical to managing the ground-water resources in this part of the Coconino Plateau. The nearly 1,000-meter-deep regional aquifer in the Redwall and Muav Limestones, however, makes studying or utilizing the resource difficult. Near-vertical faults and complex geologic structures control the ground-water flow system on the southwest side of the Kaibab Uplift near Williams, Arizona. To address the hydrogeologic complexities in the study area, a suite of techniques, which included aeromagnetic, gravity, square-array resistivity, and audiomagnetotelluric surveys, were applied as part of a regional study near Bill Williams Mountain. Existing well data and interpreted geophysical data were compiled and used to estimate depths to the water table and to prepare a potentiometric map. Geologic characteristics, such as secondary porosity, coefficient of anisotropy, and fracture-strike direction, were calculated at several sites to examine how these characteristics change with depth. The 14-kilometer-wide, seismically active northwestward-trending Cataract Creek and the northeastward-trending Mesa Butte Fault systems intersect near Bill Williams Mountain. Several north-south-trending faults may provide additional block faulting north and west of Bill Williams Mountain. Because of the extensive block faulting and regional folding, the volcanic and sedimentary rocks are tilted toward one or more of these faults. These faults provide near-vertical flow paths to the regional water table. The nearly radial fractures allow water that reaches the regional aquifer to move away from the Bill Williams Mountain area. Depth to the regional aquifer is highly variable and depends on location and local structures. On the basis of interpreted audiomagnetotelluric and square-array resistivity sounding curves and limited well data, depths to water may range from 450 to 1,300 meters.

Arizona↗

Timescales and processes of methane hydrate formation and breakdown, with application to geologic systems

Gas hydrate is an ice-like form of water and low molecular weight gas stable at temperatures of roughly -10ºC to 25ºC and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one-sixth of Earth’s methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean-atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean-atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth’s near-surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some fields studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

Journal of Geophysical Research-Solid Earth↗

The chromograph, a new analytical tool for laboratory and field use

Water supplies in the Greater Atlanta Region, Georgia, are obtained mainly from surface-water sources, but there is concern that these sources may be unable to meet rising demands. This study indicates that large quantities of good-quality ground water are obtainable from selected sites in the crystalline rocks of the area. Hydrogeologic investigations of 1,051 wells that produce 20 to nearly 500 gallons per minute revealed that large supplies can be developed where favorable structural, stratigraphic, and topographic features result in localized increases in bedrock permeability. These features are described in detail and methods are presented for using them in selecting sites for locating high-producing wells. Most of the site-selection methods are applicable to the north half of the report area, which has rectangular and trellis drainage systems, but their use is restricted in the south half where a dendritic drainage system predominates. Borehole geophysical logs and core drilling revealed that some of the highest production wells derive water from 1- to 8-inch wide horizontal fractures. These fractures occur in a variety of topographic settings and are believed to be stress-relief fractures. (USGS)

Circular↗

Shear-wave velocity model from Rayleigh wave group velocities centered on the Sacramento/San Joaquin Delta

Rayleigh wave group velocities obtained from ambient noise tomography are inverted for an upper crustal model of the Central Valley, California, centered on the Sacramento/San Joaquin Delta. Two methods were tried; the first uses SURF96, a least-squares routine. It provides a good fit to the data, but convergence is dependent on the starting model. The second uses a genetic algorithm, whose starting model is random. This method was tried at several nodes in the model and compared to the output from SURF96. The genetic code is run five times and the variance of the output of all five models can be used to obtain an estimate of error. SURF96 produces a more regular solution mostly because it is typically run with a smoothing constraint. Models from the genetic code are generally consistent with the SURF96 code sometimes producing lower velocities at depth. The full model, calculated using SURF96, employed a 2-pass strategy, which used a variable damping scheme in the first pass. The resulting model shows low velocities near the surface in the Central Valley with a broad asymmetrical sedimentary basin located close to the western edge of the Central Valley near 122°W longitude. At shallow depths the Rio Vista Basin is found nestled between the Pittsburgh/Kirby Hills and Midland faults, but a significant basin also seems to exist to the west of the Kirby Hills fault. There are other possible correlations between fast and slow velocities in the Central Valley and geologic features such as the Stockton Arch, oil or gas producing regions and the fault-controlled western boundary of the Central Valley.

California↗

Premonitory acoustic emissions and stick-slip in natural and smooth-faulted Westerly granite

A stick-slip event was induced in a cylindrical sample of Westerly granite containing a preexisting natural fault by loading at constant confining pressure of 150 MPa. Continuously recorded acoustic emission (AE) data and computer tomography (CT)-generated images of the fault plane were combined to provide a detailed examination of microscale processes operating on the fault. The dynamic stick-slip event, considered to be a laboratory analog of an earthquake, generated an ultrasonic signal that was recorded as a large-amplitude AE event. First arrivals of this event were inverted to determine the nucleation site of slip, which is associated with a geometric asperity on the fault surface. CT images and AE locations suggest that a variety of asperities existed in the sample because of the intersection of branch or splay faults with the main fault. This experiment is compared with a stick-slip experiment on a sample prepared with a smooth, artificial saw-cut fault surface. Nearly a thousand times more AE were observed for the natural fault, which has a higher friction coefficient (0.78 compared to 0.53) and larger shear stress drop (140 compared to 68 MPa). However at the measured resolution, the ultrasonic signal emitted during slip initiation does not vary significantly between the two experiments, suggesting a similar dynamic rupture process. We propose that the natural faulted sample under triaxial compression provides a good laboratory analogue for a field-scale fault system in terms of the presence of asperities, fault surface heterogeneity, and interaction of branching faults. ?? 2009.

Journal of Geophysical Research B: Solid Earth↗