Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 811 records · Page 45Linked to original sources

Search for life on Mars in surface samples: Lessons from the 1999 Marsokhod rover field experiment

The Marsokhod 1999 field experiment in the Mojave Desert included a simulation of a rover-based sample selection mission. As part of this mission, a test was made of strategies and analytical techniques for identifying past or present life in environments expected to be present on Mars. A combination of visual clues from high-resolution images and the detection of an important biomolecule (chlorophyll) with visible/near-infrared (NIR) spectroscopy led to the successful identification of a rock with evidence of cryptoendolithic organisms. The sample was identified in high-resolution images (3 times the resolution of the Imager for Mars Pathfinder camera) on the basis of a green tinge and textural information suggesting the presence of a thin, partially missing exfoliating layer revealing the organisms. The presence of chlorophyll bands in similar samples was observed in visible/NIR spectra of samples in the field and later confirmed in the laboratory using the same spectrometer. Raman spectroscopy in the laboratory, simulating a remote measurement technique, also detected evidence of carotenoids in samples from the same area. Laboratory analysis confirmed that the subsurface layer of the rock is inhabited by a community of coccoid Chroococcidioposis cyanobacteria. The identification of minerals in the field, including carbonates and serpentine, that are associated with aqueous processes was also demonstrated using the visible/NIR spectrometer. Other lessons learned that are applicable to future rover missions include the benefits of web-based programs for target selection and for daily mission planning and the need for involvement of the science team in optimizing image compression schemes based on the retention of visual signature characteristics.

California↗

Nomad Rover field experiment, Atacama Desert, Chile 1. Science results overview

Nomad was deployed for a 45 day traverse in the Atacama Desert, Chile, during the summer of 1997. During this traverse, 1 week was devoted to science experiments. The goal of the science experiments was to test different planetary surface exploration strategies that included (1) a Mars mission simulation, (2) a science on the fly experiment, where the rover was kept moving 75% of the operation time. (The goal of this operation was to determine whether or not successful interpretation of the environment is related to the time spent on a target. The role of mobility in helping the interpretation was also assessed.) (3) a meteorite search using visual and instrumental methods to remotely identify meteorites in extreme environments, and (4) a time-delay experiment with and without using the panospheric camera. The results were as follow: the remote science team positively identified the main characteristics of the test site geological environment. The science on the fly experiment showed that the selection of appropriate targets might be even more critical than the time spent on a study area to reconstruct the history of a site. During the same operation the science team members identified and sampled a rock from a Jurassic outcrop that they proposed to be a fossil. The presence of paleolife indicators in this rock was confirmed later by laboratory analysis. Both visual and instrumental modes demonstrated the feasibility, in at least some conditions, of carrying out a field search for meteorites by using remote-controlled vehicles. Finally, metrics collected from the observation of the science team operations, and the use team members made of mission data, provided critical information on what operation sequences could be automated on board rovers in future planetary surface explorations.

Atacama Desert↗

Observation of the geology and geomorphology of the 1999 Marsokhod test site

The Marsokhod rover returned data from six stations that were used to decipher the geomorphology and geology of a region not previously visited by members of the geomorphology field team. Satellite images and simulated descent images provided information about the regional setting. The landing zone was on an alluvial apron flanking a mountain block to the west and a playa surface to the east. Rover color images, infrared spectra analysis of the mountains, and the apron surface provided insight into the rock composition of the nearby mountains. From the return data the geomorphology team interpreted the region to consist of compressionally deformed, ancient marine sediments and igneous rocks exposed by more recent extensional tectonics. Unconsolidated alluvial materials blanket the lower flanks of the mountains. An ancient shoreline cut into alluvial material marks a high stand of water during a past, wetter climate period. Playa sediments floor a present-day, seasonally, dry lake. Observations made by the rover using panoramic and close-up (hand specimens-scale) image data and color scene data confirmed the presence of boulders, cobbles, and fines of various provinces. Rover traverses to sites identified as geologically distinct, such as fan, channel, shoreline, and playa, provided useful clues to the geologic interpretations. Analysis of local rocks was given context only through comparison with distant geologic features. These results demonstrated the importance of a multifaceted approach to site interpretation through comparison of interpretations derived by differing geologic techniques.

California↗

Results of the Imager for Mars Pathfinder windsock experiment

The Imager for Mars Pathfinder (IMP) windsock experiment measured wind speeds at three heights within 1.2 m of the Martian surface during Pathfinder landed operations. These wind data allowed direct measurement of near-surface wind profiles on Mars for the first time, including determination of aerodynamic roughness length and wind friction speeds. Winds were light during periods of windsock imaging, but data from the strongest breezes indicate aerodynamic roughness length of 3 cm at the landing site, with wind friction speeds reaching 1 m/s. Maximum wind friction speeds were about half of the threshold-of-motion friction speeds predicted for loose, fine-grained materials on smooth Martian terrain and about one third of the threshold-of-motion friction speeds predicted for the same size particles over terrain with aerodynamic roughness of 3 cm. Consistent with this, and suggesting that low wind speeds prevailed when the windsock array was not imaged and/or no particles were available for aeolian transport, no wind-related changes to the surface during mission operations have been recognized. The aerodynamic roughness length reported here implies that proposed deflation of fine particles around the landing site, or activation of duneforms seen by IMP and Sojourner, would require wind speeds >28 m/s at the Pathfinder top windsock height (or >31 m/s at the equivalent Viking wind sensor height of 1.6 m) and wind speeds >45 m/s above 10 m. These wind speeds would cause rock abrasion if a supply of durable particles were available for saltation. Previous analyses indicate that the Pathfinder landing site probably is rockier and rougher than many other plains units on Mars, so aerodynamic roughness length elsewhere probably is less than the 3-cm value reported for the Pathfinder site.

Journal of Geophysical Research E: Planets↗

Hydrated salt minerals on Europa's surface from the Galileo near-infrared mapping spectrometer (NIMS) investigation

We reported evidence of heavily hydrated salt minerals present over large areas of Europa's surface from analysis of reflectance spectra returned by the Galileo mission near infrared mapping spectrometer (NIMS) [ McCord et al. , 1997a, b, 1998a, b]. Here we elaborate on this earlier evidence, present spatial distributions of these minerals, examine alternate water-ice interpretations, expand on our hydrated-salts interpretation, consider salt mineral stability on Europa, and discuss the implications. Extensive well-defined areas on Europa show distinct, asymmetric water-related absorption bands in the 1 to 2.5-μm region. Radiative transfer modeling of water ice involving different particle sizes and layers at Europa temperatures does not reproduce the distinctive Europa water bands. However, ice near its melting temperature, such as in terrestrial environments, does have some characteristics of the Europa spectrum. Alternatively, some classes of heavily hydrated minerals do exhibit such water bands. Among plausible materials, heavily hydrated salt minerals, such as magnesium and sodium sulfates, sodium carbonate and their mixtures, are preferred. All Europa spectral features are present in some salt minerals and a very good match to the Europa spectrum can be achieved by mixing several salt spectra. However, no single or mix of salt mineral spectra from the limited library available has so far been found to perfectly match the Europa spectrum in every detail. The material is concentrated at the lineaments and in chaotic terrain, which are technically disrupted areas on the trailing side. Since the spectrum of the material on Europa is nearly the same everywhere so-far studied, the salt or salt-mixture composition may be nearly uniform. This suggests similar sources and processes over at least a near-hemispheric scale. This would suggest that an extensive subsurface ocean containing dissolved salts is the source, and several possible mechanisms for deposit emplacement are considered. The hydrogen bonds associated with hydration of these salts are similar or greater in strength and energy to those in pure water ice. Thus, once on the surface, the salt minerals should be as stable to disruption as water ice at the Europa temperatures, and mechanisms are suggested to enhance the stability of both materials. Spectra obtained of MgSO 4 ·6H 2 O at 77 K show only small differences from room temperature spectra. The main difference is the appearance of the individual absorptions composing the broad, composite water features and associated with the several different H 2 O sites in the salt hydrate molecule. This suggests that the Europa absorption bands are also composites. Thus higher spectral resolution may reveal these diagnostic features in Europa's spectrum. The specific salts present and their relative abundances would be indicators of the chemistry and conditions of an ocean environment, and areas of fresh, heavy concentration of these minerals should make ideal lander mission sampling sites.

Journal of Geophysical Research E: Planets↗

The failure of earthquake failure models

In this study I show that simple heuristic models and numerical calculations suggest that an entire class of commonly invoked models of earthquake failure processes cannot explain triggering of seismicity by transient or "dynamic" stress changes, such as stress changes associated with passing seismic waves. The models of this class have the common feature that the physical property characterizing failure increases at an accelerating rate when a fault is loaded (stressed) at a constant rate. Examples include models that invoke rate state friction or subcritical crack growth, in which the properties characterizing failure are slip or crack length, respectively. Failure occurs when the rate at which these grow accelerates to values exceeding some critical threshold. These accelerating failure models do not predict the finite durations of dynamically triggered earthquake sequences (e.g., at aftershock or remote distances). Some of the failure models belonging to this class have been used to explain static stress triggering of aftershocks. This may imply that the physical processes underlying dynamic triggering differs or that currently applied models of static triggering require modification. If the former is the case, we might appeal to physical mechanisms relying on oscillatory deformations such as compaction of saturated fault gouge leading to pore pressure increase, or cyclic fatigue. However, if dynamic and static triggering mechanisms differ, one still needs to ask why static triggering models that neglect these dynamic mechanisms appear to explain many observations. If the static and dynamic triggering mechanisms are the same, perhaps assumptions about accelerating failure and/or that triggering advances the failure times of a population of inevitable earthquakes are incorrect.

Journal of Geophysical Research B: Solid Earth↗

Permeability reduction in granite under hydrothermal conditions

The formation of impermeable fault seals between earthquake events is a feature of many models of earthquake generation, suggesting that earthquake recurrence may depend in part on the rate of permeability reduction of fault zone materials under hydrothermal conditions. In this study, permeability measurements were conducted on intact, fractured, and gouge-bearing Westerly granite at an effective pressure of 50 MPa and at temperatures from 150° to 500°C, simulating conditions in the earthquake-generating portions of fault zones. Pore fluids were cycled back and forth under a 2 MPa pressure differential for periods of up to 40 days. Permeability of the granite decreased with time t , following the exponential relation k = c (10 − rt ). For intact samples run between 250° and 500°C the time constant for permeability decrease r was proportional to temperature and ranged between 0.001 and 0.1 days −1 (i.e., between 0.4 and 40 decades year −1 loss of permeability). Values of r for the lower-temperature experiments differed little from the 250°C runs. In contrast, prefractured samples showed higher rates of permeability decrease at a given temperature. The surfaces of the fractured samples showed evidence of dissolution and mineral growth that increased in abundance with both temperature and time. The experimentally grown mineral assemblages varied with temperature and were consistent with a rock-dominated hydrothermal system. As such mineral deposits progressively seal the fractured samples, their rates of permeability decrease approach the rates for intact rocks at the same temperature. These results place constraints on models of precipitation sealing and suggest that fault rocks may seal at a rate consistent with earthquake recurrence intervals of typical fault zones.

California↗

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗

Strain accumulation and rotation in the Eastern California Shear Zone

Although the Eastern California Shear Zone (ECSZ) (strike ∼N25°W) does not quite coincide with a small circle drawn about the Pacific-North America pole of rotation, trilateration and GPS measurements demonstrate that the motion within the zone corresponds to right-lateral simple shear across a vertical plane (strike N33°W±5°) roughly parallel to the tangent to that local small circle (strike ∼N40°W). If the simple shear is released by slip on faults subparallel to the shear zone, the accumulated rotation is also released, leaving no secular rotation. South of the Garlock fault the principal faults (e.g., Calico-Blackwater fault) strike ∼N40°W, close enough to the strike of the vertical plane across which maximum right-lateral shear accumulates to almost wholly accommodate that accumulation of both strain and rotation by right-lateral slip. North of the Garlock fault dip slip as well as strike slip on the principal faults (strike ∼N20°W) is required to accommodate the simple shear accumulation. In both cases the accumulated rotation is released with the shear strain. The Garlock fault, which transects the ECSZ, is not offset by north-northwest striking faults nor, despite geological evidence for long-term left-lateral slip, does it appear at the present time to be accumulating left-lateral simple shear strain across the fault due to slip at depth. Rather the motion is explained by right-lateral simple shear across the orthogonal ECSZ. Left-lateral slip on the Garlock fault will release the shear strain accumulating there but would augment the accumulating rotation, resulting in a secular clockwise rotation rate ∼80 nrad yr −1 (4.6° Myr −1 ).

California↗

Complex earthquake rupture and local tsunamis

In contrast to far-field tsunami amplitudes that are fairly well predicted by the seismic moment of subduction zone earthquakes, there exists significant variation in the scaling of local tsunami amplitude with respect to seismic moment. From a global catalog of tsunami runup observations this variability is greatest for the most frequently occurring tsunamigenic subduction zone earthquakes in the magnitude range of 7 < M w < 8.5. Variability in local tsunami runup scaling can be ascribed to tsunami source parameters that are independent of seismic moment: variations in the water depth in the source region, the combination of higher slip and lower shear modulus at shallow depth, and rupture complexity in the form of heterogeneous slip distribution patterns. The focus of this study is on the effect that rupture complexity has on the local tsunami wave field. A wide range of slip distribution patterns are generated using a stochastic, self-affine source model that is consistent with the falloff of far-field seismic displacement spectra at high frequencies. The synthetic slip distributions generated by the stochastic source model are discretized and the vertical displacement fields from point source elastic dislocation expressions are superimposed to compute the coseismic vertical displacement field. For shallow subduction zone earthquakes it is demonstrated that self-affine irregularities of the slip distribution result in significant variations in local tsunami amplitude. The effects of rupture complexity are less pronounced for earthquakes at greater depth or along faults with steep dip angles. For a test region along the Pacific coast of central Mexico, peak nearshore tsunami amplitude is calculated for a large number ( N = 100) of synthetic slip distribution patterns, all with identical seismic moment ( M w = 8.1). Analysis of the results indicates that for earthquakes of a fixed location, geometry, and seismic moment, peak nearshore tsunami amplitude can vary by a factor of 3 or more. These results indicate that there is substantially more variation in the local tsunami wave field derived from the inherent complexity subduction zone earthquakes than predicted by a simple elastic dislocation model. Probabilistic methods that take into account variability in earthquake rupture processes are likely to yield more accurate assessments of tsunami hazards.

Journal of Geophysical Research B: Solid Earth↗

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗

Geochemical evidence for a magmatic CO 2 degassing event at Mammoth Mountain, California, September-December 1997

Recent time series soil CO 2 concentration data from monitoring stations in the vicinity of Mammoth Mountain, California, reveal strong evidence for a magmatic degassing event during the fall of 1997 lasting more than 2 months. Two sensors at Horseshoe Lake first recorded the episode on September 23, 1997, followed 10 days later by a sensor on the north flank of Mammoth Mountain. Direct degassing from shallow intruding magma seems an implausible cause of the degassing event, since the gas released at Horseshoe Lake continued to be cold and barren of other magmatic gases, except for He. We suggest that an increase in compressional strain on the area south of Mammoth Mountain driven by movement of major fault blocks in Long Valley caldera may have triggered an episode of increased degassing by squeezing additional accumulated CO 2 from a shallow gas reservoir to the surface along faults and other structures where it could be detected by the CO 2 monitoring network. Recharge of the gas reservoir by CO 2 emanating from the deep intrusions that probably triggered deep long-period earthquakes may also have contributed to the degassing event. The nature of CO 2 discharge at the soil-air interface is influenced by the porous character of High Sierra soils and by meteorological processes. Solar insolation is the primary source of energy for the Earth atmosphere and plays a significant role in most diurnal processes at the Earth surface. Data from this study suggest that external forcing due largely to local orographic winds influences the fine structure of the recorded CO 2 signals.

California↗

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↗

Dating of major normal fault systems using thermochronology: An example from the Raft River detachment, Basin and Range, western United States

Application of thermochronological techniques to major normal fault systems can resolve the timing of initiation and duration of extension, rates of motion on detachment faults, timing of ductile mylonite formation and passage of rocks through the crystal-plastic to brittle transition, and multiple events of extensional unroofing. Here we determine the above for the top-to-the-east Raft River detachment fault and shear zone by study of spatial gradients in 40 Ar/ 39 Ar and fission track cooling ages of footwall rocks and cooling histories and by comparison of cooling histories with deformation temperatures. Mica 40 Ar/ 39 Ar cooling ages indicate that extension-related cooling began at ∼25–20 Ma, and apatite fission track ages show that motion on the Raft River detachment proceeded until ∼7.4 Ma. Collective cooling curves show acceleration of cooling rates during extension, from 5–10°C/m.y. to rates in excess of 70–100°C/m.y. The apparent slip rate along the Raft River detachment, recorded in spatial gradients of apatite fission track ages, is 7 mm/yr between 13.5 and 7.4 Ma and is interpreted to record the rate of migration of a rolling hinge. Microstructural study of footwall mylonite indicates that deformation conditions were no higher than middle greenschist facies and that deformation occurred during cooling to cataclastic conditions. These data show that the shear zone and detachment fault represent a continuum produced by progressive exhumation and shearing during Miocene extension and preclude the possibility of a Mesozoic age for the ductile shear zone. Moderately rapid cooling in middle Eocene time likely records exhumation resulting from an older, oppositely rooted, extensional shear zone along the west side of the Grouse Creek, Raft River, and Albion Mountains.

Idaho, Utah↗

Strain accumulation across the Eastern California Shear Zone at latitude 36°30'N

The motion of a linear array of monuments extending across the Eastern California Shear Zone (ECSZ) has been measured from 1994 to 1999 with the Global Positioning System. The linear array is oriented N54&deg;E, perpendicular to the tangent to the local small circle drawn about the Pacific-North America pole of rotation, and the observed motion across the ECSZ is approximated by differential rotation about that pole. The observations suggest uniform deformation within the ECSZ (strike N23&deg;W) (26 nstrain yr &minus;1 extension normal to the zone and 39 nstrain yr &minus;1 simple right-lateral shear across it) with no significant deformation in the two blocks (the Sierra Nevada mountains and southern Nevada) on either side. The deformation may be imposed by right-lateral slip at depth on the individual major fault systems within the zone if the slip rates are: Death Valley-Furnace Creek fault 3.2&plusmn;0.9 mm yr &minus;1 , Hunter Mountain-Panamint Valley fault 3.3&plusmn;1.6 mm yr &minus;1 , and Owens Valley fault 6.9&plusmn;1.6 mm yr &minus;1 . However, this estimate of the slip rate on the Owens Valley fault is 3 times greater than the geologic estimate.

Journal of Geophysical Research B: Solid Earth↗

On the modified Mercalli intensities and magnitudes of the 1811-1812 New Madrid earthquakes

We reexamine original felt reports from the 1811–1812 New Madrid earthquakes and determine revised isoseismal maps for the three principal mainshocks. In many cases we interpret lower values than those assigned by earlier studies. In some cases the revisions result from an interpretation of original felt reports with an appreciation for site response issues. Additionally, earlier studies had assigned modified Mercalli intensity (MMI) values of V-VII to a substantial number of reports that we conclude do not describe damage commensurate with intensities this high. We investigate several approaches to contouring the MMI values using both analytical and subjective methods. For the first mainshock on 02∶15 LT December 16, 1811, our preferred contouring yields M ω7.2–7.3 using the area-moment regressions of Johnston [1996]. For the 08∶00 LT on January 23, 1812, and 03∶45 LT on February 7, 1812, mainshocks, we obtain M ω7.0 and M ω7.4–7.5, respectively. Our magnitude for the February mainshock is consistent with the established geometry of the Reelfoot fault, which all evidence suggests to have been the causative structure for this event. We note that the inference of lower magnitudes for the New Madrid events implies that site response plays a significant role in controlling seismic hazard at alluvial sites in the central and eastern United States. We also note that our results suggest that thrusting may have been the dominant mechanism of faulting associated with the 1811–1812 sequence.

Arkansas, Illinois, Kentucky, Missouri, Tennessee↗

Pore fluid pressure, apparent friction, and Coulomb failure

Many recent studies of stress-triggered seismicity rely on a fault failure model with a single free parameter, the apparent coefficient of friction, presumed to be a material constant with possible values 0 ≤ μ′ ≤ 1. These studies may present a misleading view of fault strength and the role of pore fluid pressure in earthquake failure. The parameter μ′ is intended to incorporate the effects of both friction and pore pressure, but is a material constant only if changes in pore fluid pressure induced by changes in stress are proportional to the normal stress change across the potential failure plane. Although specific models of fault zones permit such a relation, neither is it known that fault zones within the Earth behave this way, nor is this behavior expected in all cases. In contrast, for an isotropic homogeneous poroelastic model the pore pressure changes are proportional to changes in mean stress, μ′ is not a material constant, and −∞ ≤ μ′ ≤ +∞. Analysis of the change in Coulomb failure stress for tectonically loaded reverse and strike-slip faults shows considerable differences between these two pore pressure models, suggesting that such models might be distinguished from one another using observations of triggered seismicity (e.g., aftershocks). We conclude that using the constant apparent friction model exclusively in studies of Coulomb failure stress is unwise and could lead to significant errors in estimated stress change and seismic hazard.

Journal of Geophysical Research B: Solid Earth↗

Density structure of the lithosphere in the southwestern United States and its tectonic significance

We calculate a density model of the lithosphere of the southwestern United States through an integrated analysis of gravity, seismic refraction, drill hole, and geological data. Deviations from the average upper mantle density are as much as ?? 3%. A comparison with tomographic images of seismic velocities indicates that a substantial part (>50%) of these density variations is due to changes in composition rather than temperature. Pronounced mass deficits are found in the upper mantle under the Basin and Range Province and the northern part of the California Coast Ranges and adjacent ocean. The density structure of the northern and central/southern Sierra Nevada is remarkably different. The central/southern part is anomalous and is characterized by a relatively light crust underlain by a higher-density upper mantle that may be associated with a cold, stalled subducted plate. High densities are also determined within the uppermost mantle beneath the central Transverse Ranges and adjoining continental slope. The average density of the crystalline crust under the Great Valley and western Sierra Nevada is estimated to be up to 200 kg m~3 higher than the regional average, consistent with tectonic models for the obduction of oceanic crust and uppermost mantle in this region.

Journal of Geophysical Research B: Solid Earth↗