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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35Linked to original sources

USGS2014

Upwelling rebound, ephemeral secondary pycnoclines, and the creation of a near-bottom wave guide over the Monterey Bay continental shelf

Several sequential upwelling events were observed in fall 2012, using measurements from the outer half of the continental shelf in Monterey Bay, during which the infiltration of dense water onto the shelf created a secondary, near-bottom pycnocline. This deep pycnocline existed in concert with the near-surface pycnocline and enabled the propagation of near-bottom, cold, semidiurnal internal tidal bores, as well as energetic, high-frequency, nonlinear internal waves of elevation (IWOE). The IWOE occurred within 20 m of the bottom, had amplitudes of 8–24 m, periods of 6–45 min, and depth-integrated energy fluxes up to 200 W m −1 . Iribarren numbers (<0.03) indicate that these IWOE were nonbreaking in this region of the shelf. These observations further demonstrate how regional upwelling dynamics and the resulting bulk, cross-margin hydrography is a first-order control on the ability of internal waves, at tidal and higher frequencies, to propagate through continental shelf waters.

California↗

USGS2000

Seismovolcanic signals at Deception Island volcano, Antarctica: Wave field analysis and source modeling

The seismovolcanic signals associated with the volcanic activity of Deception Island (Antarctica), recorded during three Antarctic summers (1994–1995, 1995–1996 and 1996–1997), are analyzed using a dense small-aperture (500 m) seismic array. The visual and spectral classification of the seismic events shows the existence of long-period and hybrid isolated seismic events, and of low-frequency, quasi-monochromatic and spasmodic continuous tremors. All spectra have the highest amplitudes in the frequency band between 1 and 4 Hz, while hybrids and spasmodic tremors have also significant amplitudes in the high-frequency band (4–10 Hz). The array analysis indicates that almost all the well-correlated low-frequency signals share similar array parameters (slowness and back azimuth) and have the same source area, close to the array site. The polarization analysis shows that phases at high-frequency are mostly composed of P waves, and those phases dominated by low frequencies can be interpreted as surface waves. No clear shear waves are evidenced. From the energy evaluation, we have found that the reduced displacement values for surface and body waves are confined in a narrow interval. Volcano-tectonic seismicity is located close to the array, at a depth shallower than 1 km. The wave-field properties of the seismovolcanic signals allow us to assume a unique source model, a shallow resonating fluid-filled crack system at a depth of some hundreds of meters. All of the seismic activity is interpreted as the response of a reasonably stable stationary geothermal process. The differences observed in the back azimuth between low and high frequencies are a near-field effect. A few episodes of the degassification process in an open conduit were observed and modeled with a simple organ pipe.

Journal of Geophysical Research: Solid Earth↗

USGS2025

Surface-wave relocation and characterization of the October 2023 tsunamigenic seismic unrest near Sofugan volcano, Izu Islands, Japan

A moderate-magnitude earthquake swarm occurred in the remote Izu Islands region of Japan between October 1 and 8, 2023. The swarm included 151 shallow earthquakes cataloged by the U.S. Geological Survey, which notably included a roughly 2.5-hr episode of 15 successive magnitude (M) < 5.5 earthquakes. Origin times were coincident with regionally recorded tsunami waves, but tsunamigenesis for moderate-magnitude earthquakes is uncommon, indicating that volcanic activity generated the ocean displacements. Leveraging a surface-wave relative relocation approach, we estimate precise epicentroid locations for the remote swarm. Final epicentroids and caldera analogs indicate a three-stage model to explain swarm activity: (a) caldera pressurization due to magma intrusion, (b) depressurization via dike propagation away from the caldera, and (c) eruption corresponding with caldera reactivation either by collapse or additional intrusion.

Sofugan volcano, Izu Islands↗

USGS2002

Interpreting the earthquake source of the Wabash Valley seismic zone (Illinois, Indiana, and Kentucky) from seismic-reflection, gravity, and magnetic-intensity data

Reprocessing of seismic-reflection data reveals new images of upper- to middle-crustal structures beneath the Wabash Valley seismic zone, located north of the New Madrid seismic zone within the seismically active southern Illinois basin. Four intersecting deep seismic profiles (243 km total) indicate an anomalous, 5–10-km-wide zone of dipping reflections and diffractions below the western flank of the Wabash Valley fault system (WVFS). The zone corresponds in places to gently arched regions of Paleozoic strata. The reflector zone can be interpreted as a result of either (or a combination of) magmatic intrusion or structural deformation. The area encompassing the reflection profiles has experienced several moderate magnitude (3.0 ≤ m bLg ≤ 5.5) earthquakes during the past 50 years, defining the central part of the Wabash Valley seismic zone. The hypocenter of the largest 20th-century earthquake in the central USA midcontinent (9 November 1968, m bLg 5.5) corresponds to the most prominent zone of dipping middle-crustal reflections, just west of the WVFS. Both the focal mechanism (moderately dipping reverse fault) and the expected rupture zone size (∼2.9 km fault length) of this earthquake are consistent with the orientation and size of observed reflectors. Dipping reflector patterns in the Precambrian crust are not collinear with fault surfaces updip in the Paleozoic sedimentary section. This indicates that shallow Paleozoic structures are effectively “decoupled” from deeper, possibly seismogenic structure, which suggests that understanding Paleozoic structure is not the key to understanding the earthquake source. The complex dipping crustal reflectivity beneath the WVFS is typical of Paleozoic continental convergent zones observed elsewhere ( e.g. , Appalachian orogen) and thus may suggest a preserved Proterozoic suture, possibly associated with the distal Grenville orogeny or an older event. Although magnetic intensity, Bouguer gravity, and seismic-reflection data present different means of understanding the deep geology of the area, their integration aids in limiting the number of admissible interpretations. The reflection profiles indicate a variable zone of anomalous crustal structure, including the dipping reflector zone, along a trend of northeast-trending gravity and magnetic highs locally defining the Commerce geophysical lineament (CGL), which is a suspected source of seismic hazard in the central USA midcontinent. Three-dimensional inverse modeling of the residual isostatic gravity anomaly values indicates that the upper part of the dipping reflector zone beneath the CGL lies near an important density boundary in the upper Precambrian crust. The results of our study suggest that the seismogenic source just north of the New Madrid seismic zone consists, in part, of a pre-existing fabric of blind thrusts localized along pre-existing igneous intrusions, locally coincident with the CGL. This suggests that the CGL may be seismogenic in places and thus a potential seismic hazard. The variation in the expression of the CGL using reflection and potential-field data sets is probably partly related to the differing geologic features by which it is expressed, but would be also consistent with its reactivation numerous times under varying stress regimes.

Illinois, Indiana, Kentucky↗

USGS2015

Laboratory electrical resistivity analysis of geologic samples from Fort Irwin, California

Correlating laboratory resistivity measurements with geophysical resistivity models helps constrain these models to the geology and lithology of an area. Throughout the Fort Irwin National Training Center area, 111 samples from both cored boreholes and surface outcrops were collected and processed for laboratory measurements. These samples represent various lithologic types that include plutonic and metamorphic (basement) rocks, lava flows, consolidated sedimentary rocks, and unconsolidated sedimentary deposits that formed in a series of intermountain basins. Basement rocks, lava flows, and some lithified tuffs are generally resistive (≥100 ohm-meters [Ω·m]) when saturated. Saturated unconsolidated samples are moderately conductive to conductive, with resistivities generally less than 100 Ω·m, and many of these samples are less than 50 Ω·m. The unconsolidated samples can further be separated into two broad groups: (1) younger sediments that are moderately conductive, owing to their limited clay content, and (2) older, more conductive sediments with a higher clay content that reflects substantial amounts of originally glassy volcanic ash subsequently altered to clay. The older sediments are believed to be Tertiary. Time-domain electromagnetic (TEM) data were acquired near most of the boreholes, and, on the whole, close agreements between laboratory measurements and resistivity models were found.

California↗

USGS1999

Dark material in the polar layered deposits and dunes on Mars

Viking infrared thermal mapping and bistatic radar data suggest that the bulk density of the north polar erg material is much lower than that of the average Martian surface or of dark dunes at lower latitudes. We have derived a thermal inertia of 245–280 J m −2 s −1/2 K −1 (5.9–6.7×10 −3 cal cm −2 s −1/2 K −1 ) for the Proctor dune field and 25–150 J m −2 s −1/2 K −1 (0.6–3.6×10 −3 cal cm −2 s −1/2 K −1 ) for the north polar erg. The uniqueness of the thermophysical properties of the north polar erg material may be due to a unique polar process that has created them. The visible and near‐infrared spectral reflectance of the erg suggests that the dark material may be composed of basalt or ferrous clays. These data are consistent with the dark material being composed of basaltic ash or filamentary sublimate residue (FSR) particles derived from erosion of the layered deposits. Dark dust may be preferentially concentrated at the surface of the layered deposits by the formation of FSR particles upon sublimation of water ice. Further weathering and erosion of these areas of exposed layered deposits may form the dark, saltating material that is found in both polar regions. Dark FSR particles may saltate for great distances before eventually breaking down into dust grains, re‐mixing with the global dust reservoir, and being recycled into the polar layered deposits via atmospheric suspension.

Journal of Geophysical Research E: Planets↗

USGS2019

Three-dimensional basin and fault structure from a detailed seismic velocity model of Coachella Valley, Southern California

The Coachella Valley in the northern Salton Trough is known to produce destructive earthquakes, making it a high seismic hazard area. Knowledge of the seismic velocity structure and geometry of the sedimentary basins and fault zones is required to improve earthquake hazard estimates in this region. We simultaneously inverted first P wave travel times from the Southern California Seismic Network (39,998 local earthquakes) and explosions (251 land/sea shots) from the 2011 Salton Seismic Imaging Project to obtain a 3-D seismic velocity model. Earthquakes with focal depths ≤10 km were selected to focus on the upper crustal structure. Strong lateral velocity contrasts in the top ~3 km correlate well with the surface geology, including the low-velocity (<5 km/s) sedimentary basin and the high-velocity crystalline basement rocks outside the valley. Sediment thickness is ~4 km in the southeastern valley near the Salton Sea and decreases to <2 km at the northwestern end of the valley. Eastward thickening of sediments toward the San Andreas fault within the valley defines Coachella Valley basin asymmetry. In the Peninsular Ranges, zones of relatively high seismic velocities (~6.4 km/s) between 2 to 4 km depth may be related to Late Cretaceous mylonite rocks or older inherited basement structures. Other high-velocity domains exist in the model down to 9 km depth and help define crustal heterogeneity. We identify a potential fault zone in Lost Horse Valley unassociated with mapped faults in Southern California from the combined interpretation of surface geology, seismicity, and lateral velocity changes in the model.

California↗

USGS2008

Mars Exploration Rover Navigation Camera in-flight calibration

The Navigation Camera (Navcam) instruments on the Mars Exploration Rover (MER) spacecraft provide support for both tactical operations as well as scientific observations where color information is not necessary: large-scale morphology, atmospheric monitoring including cloud observations and dust devil movies, and context imaging for both the thermal emission spectrometer and the in situ instruments on the Instrument Deployment Device. The Navcams are a panchromatic stereoscopic imaging system built using identical charge-coupled device (CCD) detectors and nearly identical electronics boards as the other cameras on the MER spacecraft. Previous calibration efforts were primarily focused on providing a detailed geometric calibration in line with the principal function of the Navcams, to provide data for the MER navigation team. This paper provides a detailed description of a new Navcam calibration pipeline developed to provide an absolute radiometric calibration that we estimate to have an absolute accuracy of 10% and a relative precision of 2.5%. Our calibration pipeline includes steps to model and remove the bias offset, the dark current charge that accumulates in both the active and readout regions of the CCD, and the shutter smear. It also corrects pixel-to-pixel responsivity variations using flat-field images, and converts from raw instrument-corrected digital number values per second to units of radiance (W m-2 nm-1 sr-1), or to radiance factor (I/F). We also describe here the initial results of two applications where radiance-calibrated Navcam data provide unique information for surface photometric and atmospheric aerosol studies. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

USGS2023

Accuracy of finite fault slip estimates in subduction zone regions with topographic Green's functions and seafloor geodesy

Until recently, the lack of seafloor geodetic instrumentation and the use of unrealistically simple, half-space based forward models have resulted in poor resolution of near-trench slip in subduction zone settings. Here, we use a synthetic framework to investigate the impact of topography and geodetic data distribution on coseismic slip estimates in various subduction zone settings. We calculate surface displacements in two synthetic topographic domains that have topography similar to that of Chile and Japan, respectively. We then attempt to image target slip distributions by using a Bayesian approach to solve for slip with two sets of Green's functions—one that accounts for topography and one that does not—and five sets of 50 or more observation points selected from the synthetic surface displacements. Three of these sets of observation points are entirely onland, and two include 5–10 seafloor geodetic sites. We find that the use of topographic Green's functions always improves inferred slip models, and with seafloor geodetic data, it enables an almost perfect recovery of a target slip model, even in the near-trench region. Critically, our results demonstrate that it would be impossible for non-topographic Green's functions to properly recover the true slip distribution, particularly in the near-trench region. We also perform a parameter study with approximately 4,000 slip models estimated using a least-square approach, and find that topographic Green's functions yield significantly more accurate slip models in cases where good data (well distributed and reasonably dense) are available, even in the absence of seafloor geodetic sites.

Journal of Geophysical Research: Solid Earth↗

USGS2006

Gravity, magnetic, and physical property data in the Smoke Creek Desert area, northwest Nevada

The Smoke Creek Desert, located approximately 100 km (60 mi) north of Reno near the California-Nevada border, is a large basin situated along the northernmost parts of the Walker Lane Belt (Stewart, 1988), a physiographic province defined by northwest-striking topographic features and strike-slip faulting. Because geologic framework studies play an important role in understanding the hydrology of the Smoke Creek Desert, a geologic and geophysical effort was begun to help determine basin geometry, infer structural features, and estimate depth to Pre-Cenozoic rocks, or basement. In May and June of 2004, and June of 2005, the U.S. Geological Survey (USGS) collected 587 new gravity stations, more than 160 line-kilometers (100 line-miles) of truck-towed magnetometer data, and 111 rock property samples in the Smoke Creek Desert and vicinity in northwest Nevada, as part of an effort to characterize its hydrogeologic framework. In the Smoke Creek Desert area, gravity highs occur over rocks of the Skedaddle Mountains, Fox Range, Granite Range, and over portions of Tertiary volcanic rocks in the Buffalo Hills. These gravity highs likely reflect basement rocks, either exposed at the surface or buried at shallow depths. The southern Smoke Creek Desert corresponds to a 25-mGal isostatic gravity low, which corresponds with a basin depth of approximately 2 km. Magnetic highs are likely due to granitic, andesitic, and metavolcanic rocks, whereas magnetic lows are probably associated with less magnetic gneiss and metasedimentary rocks in the region. Three distinctive patterns of magnetic anomalies occur throughout the Smoke Creek Desert and Squaw Creek Valley, likely reflecting three different geological and structural settings.

Open-File Report↗

USGS1984

Evaluation of the ground-water resources of coastal Georgia; preliminary report on the data available as of July 1983

A compilation of ground-water data that have been collected for nearly 100 years in the coastal area of Georgia as part of cooperative activities between the U.S. Geological Survey and other agencies is presented in this report. The compilation of pertinent data indicates that information is available for use in the evaluation of the ground-water resources of the 13 counties of coastal Georgia. Included in this report is a fairly complete discussion of previous and ongoing investigations and monitoring networks, and an extensive list of references. Maps at 1:24,000, 1:100,000; and 1:1000,000 scales contain well locations and identifers for all wells in the Ground Water Site Inventory (GWSI) data base of the National Water Data Storage and Retrieval System (WATSTORE). Tabular summaries of selected site information from GWSI, including well identifiers and names , latitude-longitude location, depth of well, altitude of land surface, and use of water are presented. Water-use data from the National Water Use Data System, and water use for irrigation from the University of Georgia, Department of Agriculture survey , are tabulated. Also included are pertinent information on geophysical surveys and data obtained, and proposed project activities, particularly test-monitor well drilling.

Georgia↗

USGS2018

Characterizing the Los Angeles Aqueduct crossing of the San Andreas Fault for improved earthquake resilience

The five-mile-long Elizabeth Tunnel, which crosses the San Andreas fault (SAF) zone near Lake Hughes, California, is part of the Los Angeles Aqueduct (LAA) that delivers water from Owens Valley to the City of Los Angeles. Geologic characterization of the Elizabeth Tunnel alignment is focused on developing a better understanding of fault displacement hazards at the SAF crossing to support design of both short- and long-term strategies to increase the earthquake resilience of the City’s water supply system. The results of this study define a fault zone that is wider and more complex at the surface than at tunnel depth. A 750-ft transect of deep, angled core borings, located ~ 100 ft west of, and parallel to the tunnel alignment, penetrated a wide zone (>550 ft) of fault-damaged granitic and gneissic rocks containing a few thick fault zones that can be correlated confidently between borings as well as a multitude of very thin faults and shears throughout the rock mass. Fault rocks are predominantly defined as coherent cataclasites and ultracataclasites, as opposed to incoherent clay gouge, even at this relatively shallow sampling depth. Integrating geomorphic, geologic, petrographic, and geophysical data defines a steeply south-dipping SAF that splays upward and increases in complexity toward the surface, and provides the basis for engineering design considerations at this critical SAF crossing.

California↗

USGS2008

Light-toned salty soils and co-existing Si-rich species discovered by the Mars Exploration Rover Spirit in Columbia Hills

Light-toned soils were exposed, through serendipitous excavations by Spirit Rover wheels, at eight locations in the Columbia Hills. Their occurrences were grouped into four types on the basis of geomorphic settings. At three major exposures, the light-toned soils are hydrous and sulfate-rich. The spatial distributions of distinct types of salty soils vary substantially: with centimeter-scaled heterogeneities at Paso Robles, Dead Sea, Shredded, and Champagne-Penny, a well-mixed nature for light-toned soils occurring near and at the summit of Husband Hill, and relatively homogeneous distributions in the two layers at the Tyrone site. Aeolian, fumarolic, and hydrothermal fluid processes are suggested to be responsible for the deposition, transportation, and accumulation of these light-toned soils. In addition, a change in Pancam spectra of Tyrone yellowish soils was observed after being exposed to current Martian surface conditions for 175 sols. This change is interpreted to be caused by the dehydration of ferric sulfates on the basis of laboratory simulations and suggests a relative humidity gradient beneath the surface. Si-rich nodules and soils were observed near the major exposures of S-rich soils. They possess a characteristic feature in Pancam visible near-infrared (Vis-NIR) spectra that may be diagnostic of hydrated species, and this spectral feature can be used to search for additional Si-rich species. The exposures of hydrated salty soils within various geomorphic settings imply the potential existence of hydrous minerals in similar settings over a much wider area. Hydrous sulfates represent one of the candidates that may contribute the high level of water equivalent hydrogen in equatorial regions detected by the Neutron Spectrometer on Mars Odyssey.

Journal of Geophysical Research↗

USGS2012

Coseismic displacements: 1992 Landers, California, earthquake

The rupture associated with the 1992 Landers earthquake ( M s = 7.5) extended along the west edge of a 60 × 50 km trilateration network that had been surveyed two months before the earthquake and was resurveyed a month after it. The coseismic displacement field within the network was calculated from the changes observed in the trilateration network supplemented by the displacements measured at four GPS stations and other less timely trilateration measurements. Stations near the rupture were displaced by as much as 3 m. The observed changes are reproduced reasonably well by dislocation modeling in which the Landers rupture is represented by 8 vertical fault segments, each extending from the surface to 9 km depth. Two additional segments represent the Big Bear and Eureka Peak faults. The average modeled right-lateral slip on the Landers rupture is about 5 m. However, the residuals for the best-fit slip distribution are about four times greater than expected from measurement error. Given that uncertainty, a broad range of slip models furnish acceptable fits to the data. The large residuals are attributed to inelastic deformation (principally sympathetic slip on numerous unmodeled faults) within the fault blocks assumed to be elastic in the model.

Geophysical Research Letters↗

USGS1995

Finite-fault analysis of the 1979 March 14 Petatlan, Mexico, earthquake using teleseismic P waveforms

Vertical, teleseismic P waves recorded for the 1979 March 14 Petatlan, Mexico, earthquake were used to derive the distribution of coseismic slip using a linear finite-fault inversion scheme that solves for the amount of slip in each of a series of consecutive time windows. Data recorded by six stations of the Global Digital Seismograph Network were inverted in addition to digitized analogue long-period recordings available from nine Worldwide Standardized Seismograph Network stations. The digital data include four broad-band and short-period velocity waveforms reconstructed from the short- and long-period components. The time-window approach allows for a variable rise time on the fault and accounts for the source multiplicity evident in the recorded P waveforms. Synthetic tests conducted using the inversion method on the limited data set, however, reveal that the data are insufficient to identify the exact dislocation duration on the fault. The method is thus implemented by prescribing the fault rise time using five consecutive 1 s time windows. The coseismic slip inferred from the P waves shows a small 70 cm peak near the earthquake hypocentre and a large zone of dislocation (1.2 m maximum) further south-east. The slip pattern covers depths from 3 to 25 km and is located south-east of other recent large interplate ruptures on the Michoacan segment of the Mexican subduction zone. This result indicates that the 1979 Petatlan earthquake broke an independent, adjacent portion of the Cocos-North America plate boundary. The seismic moment of 1.5 × 10 27 dyn cm inferred from the P waves is approximately one-half the long-period moment estimated by other investigators from the observed surface waves. Although the discrepancy is within the uncertainty of the seismic-moment estimates, it may suggest the presence of a component of slow interplate motion that did not radiate significant P -wave energy.

Geophysical Journal International↗

USGS2006

Phosphate oxygen isotope ratios as a tracer for sources and cycling of phosphate in North San Francisco Bay, California

[1] A seasonal analysis assesing variations in the oxygen isotopic composition of dissolved inorganic phosphate (DIP) was conducted in the San Francisco Bay estuarine system, California. Isotopic fractionation of oxygen in DIP (exchange of oxygen between phosphate and environmental water) at surface water temperatures occurs only as a result of enzyme‐mediated, biological reactions. Accordingly, if phospate demand is low relative to input and phosphate is not heavily cycled in the ecosystem, the oxygen isotopic composition of DIP (δ 18 O p ) will reflect the isotopic composition of the source of phosphate to the system. Such is the case for the North San Francisco Bay, an anthropogenically impacted estuary with high surface water phosphate concentrations. Variability in the δ 18 O p in the bay is primarily controlled by mixing of water masses with different δ 18 O p signatures. The δ 18 O p values range from 11.4‰ at the Sacramento River to 20.1‰ at the Golden Gate. Deviations from the two‐component mixing model for the North Bay reflect additional, local sources of phosphate to the estuary that vary seasonally. Most notably, deviations from the mixing model occur at the confluence of a major river into the bay during periods of high river discharge and near wastewater treatment outlets. These data suggest that δ 18 O p can be an effective tool for identifying P point sources and understanding phosphate dynamics in estuarine systems.

Journal of Geophysical Research G: Biogeosciences↗

USGS2021

LiDAR and paleoseismology solve earthquake mystery in the Pacific Northwest, USA

One of the largest historical earthquakes in the U.S. Pacific Northwest occurred on December 15, 1872 near the south end of Lake Chelan. Lack of recognized surface deformation suggested that the earthquake occurred on a blind, perhaps deep, fault. New LiDAR data revealed a NW-side-up scarp along the north side of Spencer Canyon near Entiat, Washington. Landslides triggered during the earthquake impounded small ponds in Spencer Canyon; the larger of the two landslides obliterated a portion of the scarp. Tree-ring counts show that the oldest trees on each landslide are 130 and 128 years old, and lend credence to the idea that the earthquake triggered the landslides. Trenches across the scarp exposed a NW-dipping thrust fault offsetting young soils and Mesozoic bedrock. Radiocarbon and tree ring data shows that the last fault movement was between 1856 and 1873 CE, and was most likely during the 1872 CE earthquake.

Idaho, Montana, Oregon, Washington↗

USGS1997

Creep of water ices at planetary conditions: A compilation

Many constitutive laws for the flow of ice have been published since the advent of the Voyager explorations of the outer solar system. Conflicting data have occasionally come from different laboratories, and refinement of experimental techniques has led to the publication of laws that supersede earlier ones. In addition, there are unpublished data from ongoing research that also amend the constitutive laws. Here we compile the most current laboratory-derived flow laws for water ice phases I, II, III, V, and VI, and ice I mixtures with hard particulates. The rheology of interest is mainly that of steady state, and the conditions reviewed are the pressures and temperatures applicable to the surfaces and interiors of icy moons of the outer solar system. Advances in grain-size-dependent creep in ices I and II as well as in phase transformations and metastability under differential stress are also included in this compilation. At laboratory strain rates the several ice polymorphs are rheologically distinct in terms of their stress, temperature, and pressure dependencies but, with the exception of ice III, have fairly similar strengths. Hard particulates strengthen ice I significantly only at high particulate volume fractions. Ice III has the potential for significantly affecting mantle dynamics because it is much weaker than the other polymorphs and its region of stability, which may extend metastably well into what is nominally the ice II field, is located near likely geotherms of large icy moons.

Journal of Geophysical Research E: Planets↗