Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Oceans”

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 325 records · Page 18Linked to original sources

Stress sensitivity of fault seismicity: A comparison between limited-offset oblique and major strike-slip faults

We present a new three-dimensional inventory of the southern San Francisco Bay area faults and use it to calculate stress applied principally by the 1989 M = 7.1 Loma Prieta earthquake and to compare fault seismicity rates before and after 1989. The major high-angle right-lateral faults exhibit a different response to the stress change than do minor oblique (right-lateral/thrust) faults. Seismicity on oblique-slip faults in the southern Santa Clara Valley thrust belt increased where the faults were undamped. The strong dependence of seismicity change on normal stress change implies a high coefficient of static friction. In contrast, we observe that faults with significant offset (>50–100 km) behave differently; microseismicity on the Hayward fault diminished where right-lateral shear stress was reduced and where it was undamped by the Loma Prieta earthquake. We observe a similar response on the San Andreas fault zone in southern California after the Landers earthquake sequence. Additionally, the offshore San Gregorio fault shows a seismicity rate increase where right-lateral/oblique shear stress was increased by the Loma Prieta earthquake despite also being clamped by it. These responses are consistent with either a low coefficient of static friction or high pore fluid pressures within the fault zones. We can explain the different behavior of the two styles of faults if those with large cumulative offset become impermeable through gouge buildup; coseismically pressurized pore fluids could be trapped and negate imposed normal stress changes, whereas in more limited offset faults, fluids could rapidly escape. The difference in behavior between minor and major faults may explain why frictional failure criteria that apply intermediate coefficients of static friction can be effective in describing the broad distributions of aftershocks that follow large earthquakes, since many of these events occur both inside and outside major fault zones.

Journal of Geophysical Research B: Solid Earth

Paleomagnetism of an east-west transect across the Cascade arc in southern Washington: Implications for regional tectonism

Paleomagnetic data from a transect across the Cascade arc in southern Washington were collected to test models of vertical axis rotations for this part of the Pacific Northwest and to provide information on the tectonic history of the St. Helens seismic zone. The 75 site-mean directions are divided into three groups based on isotopic, fission track, and stratigraphic ages. Group 1 consists of samples from rocks deposited between 39 and 30 Ma and shows a mean clockwise vertical axis rotation of 34° ± 13°, group 2 consists of samples from rocks deposited between 30 and 24 Ma and shows a mean clockwise rotation of 17° ± 11°, and group 3 consists of samples from rocks between 24 and 16 Ma and shows a mean clockwise rotation of 20° ± 12°. Although these three values of rotation are statistically indistinguishable at the 95% confidence level, we interpret them in combination with other data to indicate differential rotation across the St. Helens seismic zone (SHZ). The available paleomagnetic data for Eocene-Oligocene rocks west of the seismic zone show clockwise vertical axis rotations of 30° ± 8° and 35° ± 9° consistent with the value for group 1. The rotational values for groups 2 and 3 and the intrusive suite of Kidd Creek (13 Ma), despite their different ages, have similar values east of the SHZ. Comparing these groups of data east and west of the seismic zone indicates a differential rotation of 10° ± 3° across it. In addition, the rates of rotation in southern Washington are similar to those for rocks in the western Cascades of Oregon and indicate that rotational deformation of Miocene age extends northward into southern Washington and eastward into the axis of the arc.

Journal of Geophysical Research B: Solid Earth

Stress transfer by the 1988-1989 M=5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: Unclamping at the site of peak mainshock slip

We study the stress transferred by the June 27, 1988, M = 5.3 and August 8, 1989, M = 5.4 Lake Elsman earthquakes, the largest events to strike within 15 km of the future Loma Prieta rupture zone during 74 years before the 1989 M = 6.9 Loma Prieta earthquake. We find that the first Lake Elsman event brought the rupture plane of the second event 0.3–1.6 bars (0.03–0.16 MPa) closer to Coulomb failure but that the Lake Elsman events did not bring the future Loma Prieta hypocentral zone closer to failure. Instead, the Lake Elsman earthquakes are calculated to have reduced the normal stress on (or “undamped”) the Loma Prieta rupture surface by 0.5–1.0 bar (0.05–0.10 MPa) at the site where the greatest slip subsequently occurred in the Loma Prieta earthquake. This association between the sites of peak unclamping and slip suggests that the Lake Elsman events did indeed influence the Loma Prieta rupture process. Unclamping the fault would have locally lowered the resistance to sliding. Such an effect could have been enhanced if the lowered normal stress permitted fluid infusion into the undamped part of the fault. Although less well recorded, the M L = 5.0 1964 and M L = 5.3 1967 Corralitos events struck within 10 km of the southwest end of the future Loma Prieta rupture. No similar relationship between the normal stress change and subsequent Loma Prieta slip is observed, although the high-slip patch southwest of the Loma Prieta epicenter corresponds roughly to the site of calculated Coulomb stress increase for a low coefficient of friction. The Lake Elsman-Loma Prieta result is similar to that for the 1987 M = 6.2 Elmore Ranch and M = 6.7 Superstition Hills earthquakes, suggesting that foreshocks might influence the distribution of mainshock slip rather than the site of mainshock nucleation.

Journal of Geophysical Research B: Solid Earth

Three-dimensional velocity structure of Siletzia and other accreted terranes in the Cascadia forearc of Washington

Eocene mafic crust with high seismic velocities underlies much of the Oregon and Washington forearc and acts as a backstop for accretion of marine sedimentary rocks from the obliquely subducting Juan de Fuca slab. Arc-parallel migration of relatively strong blocks of this terrane, known as Siletzia, focuses upper crustal deformation along block boundaries, which are potential sources of earthquakes. In a three-dimensional velocity model of coastal Washington, we have combined surface geology, well data, and travel times from earthquakes and controlled source seismic experiments to resolve the major boundaries of the Siletz terrane with the adjacent accreted sedimentary prism and volcanic arc. In southern Washington and northern Oregon the Siletz terrane appears to be a thick block (???20 km) that extends west of the coastline and makes a high-angle contact with the offshore accreted sedimentary prism. On its east flank the high-velocity Siletz terrane boundary coincides with an en echelon zone of seismicity in the arc. In northern Washington the western edge of Siletzia makes a lower-angled, fault-bound contact with the accretionary prism. In addition, alternating, east-west trending uplifts and downwarps of the Siletz terrane centered on the antiformal Olympic Mountains may reflect focusing of north-south compression in the northern part of the Siletz terrane. This compressional strain may result from northward transport and clockwise rotation of the Siletz terrane into the relatively fixed Canadian Coast Mountains restraining bend along the coast.

Oregon, Washington

Small scatterers in the lower mantle observed at German broadband arrays

Seismograms of earthquakes from the South Pacific recorded at a German broadband array and network show precursors to PKPdf . These precursors mainly originate from off-path scattering of PKPab or a nearby PKPbc to P (for receiver-side scattering) or from scattering of P to PKPab or PKPbc on the PKPdf path (for source-side scattering). Standard array processing techniques based on plane wave approximations (such as vespagram or frequency-wavenumber analysis) are inadequate for investigating these precursors since scattered waves cannot be approximated as plane waves for arrays and networks larger than 300 × 300 km for short-period waves. We therefore develop a migration method to estimate the location of scatterers in the mantle, at the core-mantle boundary and at the top of the outer core. With our method we are able to find isolated scatterers at the source side and the receiver side, although the depth of the scatterer is not well constrained. However, from looking at the first possible arrival time of precursors at different depth and the region where scattering can take place (scattering volume), we believe that the location of the scatterers is in the lowermost mantle. Since we have detected scatterers in regions where ultralow-velocity zones have been discovered recently, we think that the precursor energy possibly originates from scattering at partial melt at the base of the mantle. Comparing results from broadband and band-pass-filtered data the detection of small-scale structure of the ultralow-velocity zones becomes possible.

Journal of Geophysical Research B: Solid Earth

Formation of natural gas hydrates in marine sediments 1. Conceptual model of gas hydrate growth conditioned by host sediment properties

The stability of submarine gas hydrates is largely dictated by pressure and temperature, gas composition, and pore water salinity. However, the physical properties and surface chemistry of deep marine sediments may also affect the thermodynamic state, growth kinetics, spatial distributions, and growth forms of clathrates. Our conceptual model presumes that gas hydrate behaves in a way analogous to ice in a freezing soil. Hydrate growth is inhibited within fine-grained sediments by a combination of reduced pore water activity in the vicinity of hydrophilic mineral surfaces, and the excess internal energy of small crystals confined in pores. The excess energy can be thought of as a "capillary pressure" in the hydrate crystal, related to the pore size distribution and the state of stress in the sediment framework. The base of gas hydrate stability in a sequence of fine sediments is predicted by our model to occur at a lower temperature (nearer to the seabed) than would be calculated from bulk thermodynamic equilibrium. Capillary effects or a build up of salt in the system can expand the phase boundary between hydrate and free gas into a divariant field extending over a finite depth range dictated by total methane content and pore-size distribution. Hysteresis between the temperatures of crystallization and dissociation of the clathrate is also predicted. Growth forms commonly observed in hydrate samples recovered from marine sediments (nodules, and lenses in muds; cements in sands) can largely be explained by capillary effects, but kinetics of nucleation and growth are also important. The formation of concentrated gas hydrates in a partially closed system with respect to material transport, or where gas can flush through the system, may lead to water depletion in the host sediment. This "freeze-drying" may be detectable through physical changes to the sediment (low water content and overconsolidation) and/or chemical anomalies in the pore waters and metastable presence of free gas within the normal zone of hydrate stability.

Journal of Geophysical Research B: Solid Earth

Air blasts generated by rockfall impacts: Analysis of the 1996 Happy Isles event in Yosemite National Park

The July 10, 1996, Happy Isles rockfall in Yosemite National Park, California, released 23,000 to 38,000 m 3 of granite in four separate events. The impacts of the first two events which involved a 550-m free fall, generated seismic waves and atmospheric pressure waves (air blasts). We focus on the dynamic behavior of the second air blast that downed over 1000 trees, destroyed a bridge, demolished a snack bar, and caused one fatality and several injuries. Calculated velocities for the air blast from a two-phase, finite difference model are compared to velocities estimated from tree damage. From tornadic studies of tree damage, the air blast is estimated to have traveled <108–120 m/s within 50 m from the impact and decreased to <10–20 m/s within 500 m from the impact. The numerical model simulates the two-dimensional propagation of an air blast through a dusty atmosphere with initial conditions defined by the impact velocity and pressure. The impact velocity (105–107 m/s) is estimated from the Colorado Rockfall Simulation Program that simulates rockfall trajectories. The impact pressure (0.5 MPa) is constrained by the kinetic energy of the impact (10 10 –10 12 J) estimated from the seismic energy generated by the impact. Results from the air blast simulations indicate that the second Happy Isles air blast (weak shock wave) traveled with an initial velocity above the local sound speed. The size and location of the first impact are thought to have injected <50 wt% dust into the atmosphere. This amount of dust lowered the local atmospheric sound speed to ∼220 m/s. The discrepancy between calculated velocity data and field estimated velocity data (∼220 m/s versus ∼110 m/s) is attributed to energy dissipated by the downing of trees and additional entrainment of debris into the atmosphere not included in the calculations.

Journal of Geophysical Research B: Solid Earth

Gravitational stability of three-dimensional stratovolcano edifices

Catastrophic flank collapses have occurred at many stratovolcanoes worldwide. We present a three-dimensional (3-D) slope stability analysis for assessing and quantifying both the locations of minimum edifice stability and the expected volumes of potential failure. Our approach can search the materials underlying a topographic surface, represented as a digital elevation model (DEM), and determine the relative stability of all parts of the edifice. Our 3-D extension of Bishop's [1955] simplified limit-equilibrium analysis incorporates spherical failure surfaces, variable material properties, pore fluid pressures, and earthquake shaking. Although a variety of processes can trigger collapse, we focus here on gravitationally induced instability. Even homogeneous rock properties strongly influence the depth and volume of the least stable potential failure. For large failures in complex topography, patterns of potential instability do not mimic local ground surface slope alone. The May 18, 1980, catastrophic failure of the north flank of Mount St. Helens provides the best documented case history to test our method. Using the undeformed edifice topography of Mount St. Helens in an analysis of dry, static slope stability with homogeneous materials, as might be conducted in a precollapse hazard analysis, our method identified the northwest flank as the least stable region, although the north flank stability was within 5% of the minimum. Using estimates of the conditions that existed 2 days prior to collapse, including deformed topography with a north flank bulge and combined pore pressure and earthquake shaking effects, we obtained good estimates of the actual failure location and volume. Our method can provide estimates of initial failure volume and location to aid in assessing downslope or downstream hazards.

Washington

The crustal thickness of Australia

We investigate the crustal structure of the Australian continent using the temporary broadband stations of the Skippy and Kimba projects and permanent broadband stations. We isolate near-receiver information, in the form of crustal P-to-S conversions, using the receiver function technique. Stacked receiver functions are inverted for S velocity structure using a Genetic Algorithm approach to Receiver Function Inversion (GARFI). From the resulting velocity models we are able to determine the Moho depth and to classify the width of the crust-mantle transition for 65 broadband stations. Using these results and 51 independent estimates of crustal thickness from refraction and reflection profiles, we present a new, improved, map of Moho depth for the Australian continent. The thinnest crust (25 km) occurs in the Archean Yilgarn Craton in Western Australia; the thickest crust (61 km) occurs in Proterozoic central Australia. The average crustal thickness is 38.8 km (standard deviation 6.2 km). Interpolation error estimates are made using kriging and fall into the range 2.5–7.0 km. We find generally good agreement between the depth to the seismologically defined Moho and xenolith-derived estimates of crustal thickness beneath northeastern Australia. However, beneath the Lachlan Fold Belt the estimates are not in agreement, and it is possible that the two techniques are mapping differing parts of a broad Moho transition zone. The Archean cratons of Western Australia appear to have remained largely stable since cratonization, reflected in only slight variation of Moho depth. The largely Proterozoic center of Australia shows relatively thicker crust overall as well as major Moho offsets. We see evidence of the margin of the contact between the Precambrian craton and the Tasman Orogen, referred to as the Tasman Line.

Journal of Geophysical Research B: Solid Earth

Style and age of late Oligocene-early Miocene deformation in the southern Stillwater Range, west central Nevada: Paleomagnetism, geochronology, and field relations

Paleomagnetic and geochronologic data combined with geologic mapping tightly restrict the timing and character of a late Oligocene to early Miocene episode of large magnitude extension in the southern Stillwater Range and adjacent regions of west central Nevada. The southern Stillwater Range was the site of an Oligocene to early Miocene volcanic center comprising (1) 28.3 to 24.3 Ma intracaldera ash flow tuffs, lava flows, and subjacent plutons associated with three calderas, (2) 24.8 to 20.7 Ma postcaldera silicic dikes and domes, and (3) unconformably overlying 15.3 to 13.0 Ma dacite to basalt lava flows, plugs, and dikes. The caldera-related tuffs, lava flows, and plutons were tilted 60°-70° either west or east during the initial period of Cenozoic deformation that accommodated over 100% extension. Directions of remanent magnetization obtained from these extrusive and intrusive, caldera-related rocks are strongly deflected from an expected Miocene direction in senses appropriate for their tilt. A mean direction for these rocks after tilt correction, however, suggests that they were also affected by a moderate (33.4° ± 11.8°) component of counterclockwise vertical axis rotation. Paleomagnetic data indicate that the episode of large tilting occurred during emplacement of 24.8 to 20.7 Ma postcaldera dikes and domes. In detail, an apparent decrease in rotation with decreasing age of individual, isotopically dated bodies of the postcaldera group indicates that most tilting occurred between 24.4 and 24.2 Ma. The onset of tilting immediately following after the final caldera eruptions suggests that the magmatism and deformation were linked. Deformation was not driven by magma buoyancy, however, because tilting equally affected the caldera systems of different ages, including their plutonic roots. It is more likely that regional extension was focused in the southern Stillwater Range due to magmatic warming and reduction of tensile strength of the brittle crust. Faults that accommodated deformation in the southern Stillwater Range initially dipped steeply and cut deeply to expose more than 9 km of crustal section. The exposed crustal sections are probably rotated blocks above an unexposed basal detachment that lay near the early Miocene brittle-ductile transition.

Journal of Geophysical Research B: Solid Earth

Crustal deformation associated with glacial fluctuations in the eastern Chugach Mountains, Alaska

The changes of the solid Earth in south central Alaska in response to two major glacial fluctuations on different temporal and spatial scales have been estimated and we evaluated their influence on the stress state and ongoing tectonic deformation of the region. During the recent (1993–1995) Bering Glacier surge, a large transfer of ice from the Bagley Ice Field to the Bering Glacier terminus region occurred. We estimated the elastic displacement of the solid Earth due to ice mass redistribution from Global Positioning System (GPS) measurements at sites near the surging glacier. We can account for these displacements by transfer of an ice volume of about 14 km 3 from the surge reservoir area to the terminus region. We examined the background seismicity (M L ≥ 2.5) before, during, and after the surge. We found that the occurrence of small earthquakes (M L ≤ 4.0) in the surge reservoir region increased during the surge time interval possibly in response to a decrease in ice mass. This suggests that a small decrease in the vertical stress, σ 3 , could be enough to modulate the occurrence of small, shallow earthquakes in this dominantly thrust fault setting. During this century the southern Alaska coastal glaciers have been undergoing an overall decrease in volume. Based on our compilation of changes in the extent and thickness of the coastal glaciers between the Malaspina and Bering, we calculated surface displacements due to the Earth's viscoelastic response to annual thinning and to the cumulative retreat over the last 100 years. The uplift of the region due to an average annual thinning rate of 1–6 m/yr in the ablation region is 1–12 mm/yr. For our reference model with a viscosity of 5×10 19 Pa s for depths between ≈ 40 and 200 km the total viscoelastic response due to the retreat over the last century may be as much as a couple of meters within the coastal ablation zone near Icy Bay. The maximum decrease in σv between 0 and 10 km was ≈ 1.0 MPa, which is significant in relation to the stress drops in recent earthquakes (≈ 2 to 10 MPa) but small in relation to the estimated tectonic stress magnitude. Therefore the occurrence of an earthquake such as the St. Elias (1979, M S = 7.2) may have been advanced in time; however, most of the ongoing stress accumulation would be primarily due to tectonic forces.

Journal of Geophysical Research B: Solid Earth

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth

Precipitation structure in the Sierra Nevada of California during winter

Influences of upper air characteristics along the coast of California upon winter-time (November-April) precipitation in the Sierra Nevada are investigated. Precipitation events in the Sierra Nevada region occur mostly during wintertime, irrespective of station location (leeside or windside) and elevation. Most precipitation episodes in the region are associated with moist southwesterly winds (coming from the southwest direction) and also tend to occur when the 700-mbar temperature at the upwind direction is close to −2°C. This favored wind direction and temperature signify the importance of both moisture transport and orographic lifting in augmenting precipitation in the region. By utilizing the observed dependency of the precipitation upon the upper air conditions, a linear model is formulated to quantify the precipitation observed at different sites as a function of moisture transport. The skill of the model increases with timescale of aggregation, reaching more than 50% variance explained at an aggregation period of 5–7 days. This indicates that upstream air moisture transport can be used to estimate the precipitation totals in the Sierra Nevada region.

Journal of Geophysical Research D: Atmospheres

Chamber measurement of surface-atmosphere trace gas exchange: Numerical evaluation of dependence on soil interfacial layer, and source/sink products

We employed a three-dimensional finite difference gas diffusion model to simulate the performance of chambers used to measure surface-atmosphere tace gas exchange. We found that systematic errors often result from conventional chamber design and deployment protocols, as well as key assumptions behind the estimation of trace gas exchange rates from observed concentration data. Specifically, our simulationshowed that (1) when a chamber significantly alters atmospheric mixing processes operating near the soil surface, it also nearly instantaneously enhances or suppresses the postdeployment gas exchange rate, (2) any change resulting in greater soil gas diffusivity, or greater partitioning of the diffusing gas to solid or liquid soil fractions, increases the potential for chamber-induced measurement error, and (3) all such errors are independent of the magnitude, kinetics, and/or distribution of trace gas sources, but greater for trace gas sinks with the same initial absolute flux. Finally, and most importantly, we found that our results apply to steady state as well as non-steady-state chambers, because the slow rate of gas diffusion in soil inhibits recovery of the former from their initial non-steady-state condition. Over a range of representative conditions, the error in steady state chamber estimates of the trace gas flux varied from -30 to +32%, while estimates computed by linear regression from non-steadystate chamber concentrations were 2 to 31% too small. Although such errors are relatively small in comparison to the temporal and spatial variability characteristic of trace gas exchange, they bias the summary statistics for each experiment as well as larger scale trace gas flux estimates based on them.

Journal of Geophysical Research D: Atmospheres

Retention of an atmosphere on early Mars

The presence of valley networks and indications of high erosion rates in ancient terrains on Mars suggest that Mars was warm and wet during heavy bombardment. Various processes that could occur on early Mars were integrated into a self-consistent model to determine what circumstances might lead to warm temperatures during and at the end of heavy bombardment. Included were weathering and burial of CO 2 as carbonates, impact erosion, sputtering, and recycling of CO 2 back into the atmosphere by burial and heating. The models suggest that despite losses from the atmosphere by weathering and impact erosion, Mars could retain a 0.5 to 1 bar atmosphere at the end of heavy bombardment partly because weathering temporarily sequesters CO 2 in the ground and protects it from impact erosion while the impact rate is declining and impact erosion is becoming less effective. Because of the low output of the early Sun, surface temperatures can be above freezing only for a very efficient greenhouse, such as that suggested by Forget and Pierrehumbert [1997]. With weak greenhouse models, temperatures are below freezing throughout heavy bombardment, and such a large amount of CO 2 is left in the atmosphere at the end of heavy bombardment that it is difficult to eliminate subsequently to arrive at the present surface inventory. With strong greenhouse models, temperatures are well above freezing during heavy bombardment and drop to close to freezing at the end of heavy bombardment, at which time the atmosphere contains 0.5 to 1 bar of CO 2 . This can be largely eliminated subsequently by sputtering and low-temperature weathering. Such a model is consistent with the change in erosion rate and the declining rate of valley formation at the end of heavy bombardment. Conditions that favor warm temperatures at the end of heavy bombardment are an efficient greenhouse, low weathering rates, low impact erosion rates, and a smaller fraction of heat lost by conduction as opposed to transport of lava to the surface.

Journal of Geophysical Research E: Planets

Mineralogic and compositional properties of Martian soil and dust: results from Mars Pathfinder

Mars Pathfinder obtained multispectral, elemental, magnetic, and physical measurements of soil and dust at the Sagan Memorial Station during the course of its 83 sol mission. We describe initial results from these measurements, concentrating on multispectral and elemental data, and use these data, along with previous Viking, SNC meteorite, and telescopic results, to help constrain the origin and evolution of Martian soil and dust. We find that soils and dust can be divided into at least eight distinct spectral units, based on parameterization of Imager for Mars Pathfinder (IMP) 400 to 1000 nm multispectral images. The most distinctive spectral parameters for soils and dust are the reflectivity in the red, the red/blue reflectivity ratio, the near-IR spectral slope, and the strength of the 800 to 1000 nm absorption feature. Most of the Pathfinder spectra are consistent with the presence of poorly crystalline or nanophase ferric oxide(s), sometimes mixed with small but varying degrees of well-crystalline ferric and ferrous phases. Darker soil units appear to be coarser-grained, compacted, and/or mixed with a larger amount of dark ferrous materials relative to bright soils. Nanophase goethite, akaganeite, schwertmannite, and maghemite are leading candidates for the origin of the absorption centered near 900 nm in IMP spectra. The ferrous component in the soil cannot be well-constrained based on IMP data. Alpha proton X-ray spectrometer (APXS) measurements of six soil units show little variability within the landing site and show remarkable overall similarity to the average Viking-derived soil elemental composition. Differences exist between Viking and Pathfinder soils, however, including significantly higher S and Cl abundances and lower Si abundances in Viking soils and the lack of a correlation between Ti and Fe in Pathfinder soils. No significant linear correlations were observed between IMP spectral properties and APXS elemental chemistry. Attempts at constraining the mineralogy of soils and dust using normative calculations involving mixtures of smectites and silicate and oxide minerals did not yield physically acceptable solutions. We attempted to use the Pathfinder results to constrain a number of putative soil and dust formation scenarios, including palagonitization and acid-fog weathering. While the Pathfinder soils cannot be chemically linked to the Pathfinder rocks by palagonitization, this study and McSween et al. [1999] suggest that palagonitic alteration of a Martian basaltic rock, plus mixture with a minor component of locally derived andesitic rock fragments, could be consistent with the observed soil APXS and IMP properties.

Journal of Geophysical Research E: Planets

Compositional analyses of small lunar pyroclastic deposits using Clementine multispectral data

Clementine ultraviolet-visible (UVVIS) data are used to examine the compositions of 18 pyroclastic deposits (15 small, three large) at 13 sites on the Moon. Compositional variations among pyroclastic deposits largely result from differing amounts of new basaltic (or juvenile) material and reworked local material entrained in their ejecta upon eruption. Characterization of pyroclastic deposit compositions allows us to understand the mechanisms of lunar explosive volcanism. Evidence for compositional differences between small pyroclastic deposits at a single site is observed at Atlas crater. At all sites, compositional variation among the small pyroclastic deposits is consistent with earlier classification based on Earth-based spectra: three compositional groups can be observed, and the trend of increasing mafic absorption band strength from Group 1 to Group 2 to Group 3 is noted. As redefined here, Group 1 deposits include those of Alphonsus West, Alphonsus Southeast, Alphonsus Northeast 2, Atlas South, Crüger, Franklin, Grimaldi, Lavoisier, Oppenheimer, Orientale, and Riccioli. Group 1 deposits resemble lunar highlands, with weak mafic bands and relatively high UV/VIS ratios. Group 2 deposits include those of Alphonsus Northeast 1, Atlas North, Eastern Frigoris East and West, and Aristarchus Plateau; Group 2 deposits are similar to mature lunar maria, with moderate mafic band depths and intermediate UV/VIS ratios. The single Group 3 deposit, J. Herschel, has a relatively strong mafic band and a low UV/VIS ratio, and olivine is a likely juvenile component. Two of the deposits in these groups, Orientale and Aristarchus, are large pyroclastic deposits. The third large pyroclastic deposit, Apollo 17/Taurus Littrow, has a very weak mafic band and a high UV/VIS ratio and it does not belong to any of the compositional groups for small pyroclastic deposits. The observed compositional variations indicate that highland and mare materials are also present in many large and small pyroclastic deposits, and they suggest that volcanic glasses or spheres may not be dominant juvenile components in all large pyroclastic deposits.

Journal of Geophysical Research E: Planets

Detection of crystalline hematite mineralization on Mars by the Thermal Emission Spectrometer: evidence for near-surface water

The Thermal Emission Spectrometer (TES) instrument on the Mars Global Surveyor (MGS) mission has discovered a remarkable accumulation of crystalline hematite (α-Fe 2 O 3 ) that covers an area with very sharp boundaries approximately 350 by 350–750 km in size centered near 2°S latitude between 0° and 5°W longitude (Sinus Meridiani). Crystalline hematite is uniquely identified by the presence of fundamental vibrational absorption features centered near 300, 450, and >525 cm −1 and by the absence of silicate fundamentals in the 1000 cm −1 region. Spectral features resulting from atmospheric CO 2 , dust, and water ice were removed using a radiative transfer model. The spectral properties unique to Sinus Meridiani were emphasized by removing the average spectrum of the surrounding region. The depth and shape of the hematite fundamental bands show that the hematite is crystalline and relatively coarse grained (>5–10 μm). Diameters up to and greater than hundreds of micrometers are permitted within the instrumental noise and natural variability of hematite spectra. Hematite particles <5–10 μm in diameter (as either unpacked or hard-packed powders) fail to match the TES spectra. The spectrally derived areal abundance of hematite varies with particle size from ∼10% (>30 μm diameter) to 40–60% (10 μm diameter). The hematite in Sinus Meridiani is thus distinct from the fine-grained (diameter <5–10 μm), red, crystalline hematite considered, on the basis of visible, near-IR data, to be a minor spectral component in Martian bright regions like Olympus-Amazonis. Sinus Meridiani hematite is closely associated with a smooth, layered, friable surface that is interpreted to be sedimentary in origin. This material may be the uppermost surface in the region, indicating that it might be a late stage sedimentary unit or a layered portion of the heavily cratered plains units. We consider five possible mechanisms for the formation of coarse-grained, crystalline hematite. These processes fall into two classes depending on whether they require a significant amount of near-surface water: the first is chemical precipitation that includes origin by (1) precipitation from standing, oxygenated, Fe-rich water (oxide iron formations), (2) precipitation from Fe-rich hydrothermal fluids, (3) low-temperature dissolution and precipitation through mobile ground water leaching, and (4) formation of surface coatings, and the second is thermal oxidation of magnetite-rich lavas. Weathering and alteration processes, which produce nanophase and red hematite, are not consistent with the coarse, crystalline hematite observed in Sinus Meridiani. We prefer chemical precipitation models and favor precipitation from Fe-rich water on the basis of the probable association with sedimentary materials, large geographic size, distance from a regional heat source, and lack of evidence for extensive groundwater processes elsewhere on Mars. The TES results thus provide mineralogic evidence for probable large-scale water interactions. The Sinus Meridiani region may be an ideal candidate for future landed missions searching for biotic and prebiotic environments, and the physical characteristics of this site satisfy all of the engineering requirements for the missions currently planned.

Journal of Geophysical Research E: Planets