Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysics Research”

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 1,405 records · Page 78Linked to original sources

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth↗

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗

Introduction to special section: Stress triggers, stress shadows, and implications for seismic hazard

Many aspects of earthquake mechanics remain an enigma as we enter the closing years of the twentieth century. One potential bright spot is the realization that simple calculations of stress changes may explain some earthquake interactions, just as previous and on going studies of stress changes have begun to explain human-induced seismicity. This paper, which introduces the special section “Stress Triggers, Stress Shadows, and Implications for Seismic Hazard,” reviews many published works and presents a compilation of quantitative earthquake interaction studies from a stress change perspective. This synthesis supplies some clues about certain aspects of earthquake mechanics. It also demonstrates that much work remains before we can understand the complete story of how earthquakes work.

Journal of Geophysical Research B: Solid Earth↗

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Episodic fluid flow in the Nankai accretionary complex: Timescale, geochemistry, flow rates, and fluid budget

Down-hole geochemical anomalies encountered in active accretionary systems can be used to constrain the timing, rates, and localization of fluid flow. Here we combine a coupled flow and solute transport model with a kinetic model for smectite dehydration to better understand and quantify fluid flow in the Nankai accretionary complex offshore of Japan. Compaction of sediments and clay dehydration provide fluid sources which drive the model flow system. We explicitly include the consolidation rate of underthrust sediments in our calculations to evaluate the impact that variations in this unknown quantity have on pressure and chloride distribution. Sensitivity analysis of steady state pressure solutions constrains bulk and flow conduit permeabilities. Steady state simulations with 30% smectite in the incoming sedimentary sequence result in minimum chloride concentrations at site 808 of 550 m M , but measured chlorinity is as low as 447 m M . We simulate the transient effects of hydrofracture or a strain event by assuming an instantaneous permeability increase of 3–4 orders of magnitude along a flow conduit (in this case the décollement), using steady state results as initial conditions. Transient results with an increase in décollement permeability from 10 −16 m 2 to 10 −13 m 2 and 20% smectite reproduce the observed chloride profile at site 808 after 80–160 kyr. Modeled chloride concentrations are highly sensitive to the consolidation rate of underthrust sediments, such that rapid compaction of underthrust material leads to increased freshening. Pressures within the décollement during transient simulations rise rapidly to a significant fraction of lithostatic and remain high for at least 160 kyr, providing a mechanism for maintaining high permeability. Flow rates at the deformation front for transient simulations are in good agreement with direct measurements, but steady state flow rates are 2–3 orders of magnitude smaller than observed. Fluid budget calculations indicate that nearly 71% of the incoming water in the sediments leaves the accretionary wedge via diffuse flow out the seafloor, 0–5% escapes by focused flow along the décollement, and roughly 1% is subducted.

Journal of Geophysical Research B: Solid Earth↗

Seismotectonics of the Loma Prieta, California, region determined from three-dimensional Vp, Vp/Vs, and seismicity

Three-dimensional V p and V p / V s velocity models for the Loma Prieta region were developed from the inversion of local travel time data (21,925 P arrivals and 1,116 S arrivals) from earthquakes, refraction shots, and blasts recorded on 1700 stations from the Northern California Seismic Network and numerous portable seismograph deployments. The velocity and density models and microearthquake hypocenters reveal a complex structure that includes a San Andreas fault extending to the base of the seismogenic layer. A body with high V p extends the length of the rupture and fills the 5 km wide volume between the Loma Prieta mainshock rupture and the San Andreas and Sargent faults. We suggest that this body controls both the pattern of background seismicity on the San Andreas and Sargent faults and the extent of rupture during the mainshock, thus explaining how the background seismicity outlined the along-strike and depth extent of the mainshock rupture on a different fault plane 5 km away. New aftershock focal mechanisms, based on three-dimensional ray tracing through the velocity model, support a heterogeneous postseismic stress field and can not resolve a uniform fault normal compression. The subvertical (or steeply dipping) San Andreas fault and the fault surfaces that ruptured in the 1989 Loma Prieta earthquake are both parts of the San Andreas fault zone and this section of the fault zone does not have a single type of characteristic event.

Journal of Geophysical Research B: Solid Earth↗

Deformation across the rupture zone of the 1964 Alaska earthquake, 1993–1997

A linear array of 15 geodetic monuments was installed in 1993 across the rupture zone of the 1964 Alaska earthquake ( M w = 9.2). The array extends from Middleton Island (at the edge of the continental shelf and 80 km from the Alaska‐Aleutian trench) to north of Palmer, Alaska (380 km from the trench), in the approximate direction of Pacific‐North American plate convergence (N15.5°W). The array was surveyed in June 1993, May 1995, and June 1997. The changes between surveys are a measure of the deformation of the continental margin across the subduction zone in southern Alaska. Measured relative to the interior of the North American plate, the horizontal velocities on the outer plate margin are parallel to the direction of plate convergence (N15.5°W ) and reach a maximum (58 mm yr −1 ) about 150 km from the trench. Beyond about 300 km from the trench the observed horizontal velocities are small. A narrow (halfwidth 50 km) zone of significant uplift (10 mm yr −1 maximum) is observed about 300 km from the trench, coinciding roughly with the locus of maximum coseismic subsidence associated with the 1964 Alaska earthquake. Although the deformation is roughly described by the conventional model of deformation at a subduction zone (deformation due to virtual back slip on the main thrust zone at the 55 mm yr −1 plate convergence rate), a better fit is given with a 65 mm yr −1 virtual back (normal) slip rate. This higher rate is attributed to continued postseismic relaxation. The model does not explain the relatively high uplift rate and low N15.5°W velocity observed at Middleton Island. That anomalous motion is attributed to continued thrusting on postulated upward trending splays from the subduction zone beneath the island.

Alaska↗

Paleomagnetism of the Miocene intrusive suite of Kidd Creek: Timing of deformation in the Cascade arc, southern Washington

Paleomagnetic study of the intrusive suite of Kidd Creek in the southern Washington Cascades (23 sites in dikes and sills) was undertaken to help determine if these rocks are comagmatic and whether they postdate regional folding of the volcanic arc. Fission track and 40 Ar- 39 Ar age determinations indicate an age of ∼12.7 Ma (middle Miocene) for these rocks. The similarity of normal-polarity characteristic directions for most samples corroborate the available geochemical data indicating that these rocks are most likely comagmatic. Reversed-polarity directions for samples from four sites, however, show that emplacement of Kidd Creek intrusions spanned at least one reversal of the geomagnetic field. The paleomagnetic directions for the dikes and sills fail a fold test at the 99% confidence level indicating that the Kidd Creek rocks postdate regional folding. The mean in situ direction also indicates that the Kidd Creek and older rocks have been rotated 22°±6° clockwise about a vertical or near-vertical axis from the expected Miocene direction. Compression and regional folding of the Cascade arc in southern Washington therefore had ended by ∼12 Ma prior to the onset of deformation resulting in rotation of these rocks.

Journal of Geophysical Research B: Solid Earth↗

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗

Coseismic slip resolution along a plate boundary megathrust: the Nankai Trough, southwest Japan

Geodetic survey measurements are used to estimate the coseismic slip distribution in the 1944 Tonankai (M w =8.1) and 1946 Nankaido (M w =8.3) earthquakes and to assess quantitatively the degree to which this slip is resolved on the plate boundary megathrust. Data used include 798 angle changes from triangulation surveys, 328 leveling section differences, and 5 coseismic tidal gage offsets. Many of the nominally coseismic triangulation data span ∼50 years, nearly half the earthquake cycle, and correction for interseismic deformation using post-1950 observations is applied. Microseismicity is used to define the configuration of the plate boundary interface and approximate it with a continuous, multisegment fault model. Because the onshore geodetic data have very limited resolving power for offshore fault segments, offshore coseismic slip was constrained by Satctke's [1993] estimation based on tsunami data. The majority of the coseismic slip occurs between 15 and 25 km depth. Although resolution declines toward the trench axis, it is sufficiently good to define two distinct high-slip regions, one off southeastern Shikoku Island (11 m maximum) and the other offshore of Kii Peninsula (3 m maximum). The slip magnitude off southeastern Shikoku, coupled with the plate convergence rate, would imply an recurrence interval of about 270 years, much-longer than the average repeat time of ∼120 years for historical great earthquakes on the Nankai Trough. However, the maximum coseismic slip is sensitive to the assumed fault geometry, and slippage on trough-parallel splay faults could significantly decrease the maximum slip to about 6 m.

Journal of Geophysical Research B: Solid Earth↗

Greenland Sea Odden sea ice feature: Intra-annual and interannual variability

The “Odden” is a large sea ice feature that forms in the east Greenland Sea that may protrude eastward to 5°E from the main sea ice pack (at about 8°W) between 73° and 77°N. It generally forms at the beginning of the winter season and can cover 300,000 km 2 . Throughout the winter the outer edge of the Odden may advance and retreat by several hundred kilometers on timescales of a few days to weeks. Satellite passive microwave observations from 1978 through 1995 provide a continuous record of the spatial and temporal variations of this extremely dynamic phenomenon. Aircraft synthetic aperture radar, satellite passive microwave, and ship observations in the Odden show that the Odden consists of new ice types, rather than older ice types advected eastward from the main pack. The 17-year record shows both strong interannual and intra-annual variations in Odden extent and temporal behavior. For example, in 1983 the Odden was weak, in 1984 the Odden did not occur, and in 1985 the Odden returned late in the season. An analysis of the ice area and extent time series derived from the satellite passive microwave observations along with meteorological data from the International Arctic Buoy Program (IABP) determined the meteorological forcing associated with Odden growth, maintenance, and decay. The key meteorological parameters that are related to the rapid ice formation and decay associated with the Odden are, in order of importance, air temperature, wind speed, and wind direction. Oceanographic parameters must play an important role in controlling Odden formation, but it is not yet possible to quantify this role because of a lack of long-term oceanographic observations.

Journal of Geophysical Research C: Oceans↗

Winter and early spring CO2 efflux from tundra communities of northern Alaska

Carbon dioxide concentrations through snow were measured in different arctic tundra communities on the North Slope of Alaska during winter and early spring of 1996. Subnivean CO 2 concentrations were always higher than atmospheric CO 2 . A steady state diffusion model was used to generate conservative estimates of CO 2 flux to the atmosphere. The magnitude of CO 2 efflux differed with tundra community type, and rates of carbon release increased from March to May. Winter CO 2 efflux was highest in riparian and snow bed communities and lowest in dry heath, upland tussock, and wet sedge communities. Snow generally accrues earlier in winter and is deeper in riparian and snow bed communities compared with other tundra communities, which are typically windswept and do not accumulate much snow during the winter. These results support the hypothesis that early and deep snow accumulation may insulate microbial populations from very cold temperatures, allowing sites with earlier snow cover to sustain higher levels of activity throughout winter compared to communities that have later developing snow cover. Extrapolating our estimates of CO 2 efflux to the entire snow-covered season indicates that total carbon flux during winter in the Arctic is 13–109 kg CO 2 -C ha −1 , depending on the vegetation community type. Wintertime CO 2 flux is a potentially important, yet largely overlooked, part of the annual carbon cycle of tundra, and carbon release during winter should be accounted for in estimates of annual carbon balance in arctic ecosystems.

Journal of Geophysical Research D: Atmospheres↗

Analysis of simulated advanced spaceborne thermal emission and reflection (ASTER) radiometer data of the Iron Hill, Colorado, study area for mapping lithologies

The advanced spaceborne thermal emission and reflection (ASTER) radiometer was designed to record reflected energy in nine channels with 15 or 30 m resolution, including stereoscopic images, and emitted energy in five channels with 90 m resolution from the NASA Earth Observing System AMI platform. A simulated ASTER data set was produced for the Iron Hill, Colorado, study area by resampling calibrated, registered airborne visible/infrared imaging spectrometer (AVIRIS) data, and thermal infrared multispectral scanner (TIMS) data to the appropriate spatial and spectral parameters. A digital elevation model was obtained to simulate ASTER-derived topographic data. The main lithologic units in the area are granitic rocks and felsite into which a carbonatite stock and associated alkalic igneous rocks were intruded; these rocks are locally covered by Jurassic sandstone, Tertiary rhyolitic tuff, and colluvial deposits. Several methods were evaluated for mapping the main lithologic units, including the unsupervised classification and spectral curve-matching techniques. In the five thermalinfrared (TIR) channels, comparison of the results of linear spectral unmixing and unsupervised classification with published geologic maps showed that the main lithologic units were mapped, but large areas with moderate to dense tree cover were not mapped in the TIR data. Compared to TIMS data, simulated ASTER data permitted slightly less discrimination in the mafic alkalic rock series, and carbonatite was not mapped in the TIMS nor in the simulated ASTER TIR data. In the nine visible and near-infrared channels, unsupervised classification did not yield useful results, but both the spectral linear unmixing and the matched filter techniques produced useful results, including mapping calcitic and dolomitic carbonatite exposures, travertine in hot spring deposits, kaolinite in argillized sandstone and tuff, and muscovite in sericitized granite and felsite, as well as commonly occurring illite/muscovite. However, the distinction made in AVIRIS data between calcite and dolomite was not consistently feasible in the simulated ASTER data. Comparison of the lithologie information produced by spectral analysis of the simulated ASTER data to a photogeologic interpretation of a simulated ASTER color image illustrates the high potential of spectral analysis of ASTER data to geologic interpretation.

Journal of Geophysical Research D: Atmospheres↗

Non-water-ice constituents in the surface material of the icy Galilean satellites from the Galileo near-infrared mapping spectrometer investigation

We present evidence for several non‐ice constituents in the surface material of the icy Galilean satellites, using the reflectance spectra returned by the Galileo near infrared mapping spectrometer (NIMS) experiment. Five new absorption features are described at 3.4, 3.88, 4.05, 4.25, and 4.57 μm for Callisto and Ganymede, and some seem to exist for Europa as well. The four absorption bands strong enough to be mapped on Callisto and Ganymede are each spatially distributed in different ways, indicating different materials are responsible for each absorption. The spatial distributions are correlated at the local level in complex ways with surface features and in some cases show global patterns. Suggested candidate spectrally active groups, perhaps within larger molecules, producing the five absorptions include C‐H, S‐H, SO 2 , CO 2 , and C≡N. Organic material like tholins are candidates for the 4.57‐ and 3.4‐μm features. We suggest, based on spectroscopic evidence, that CO 2 is present as a form which does not allow rotational modes and that SO 2 is present neither as a frost nor a free gas. The CO 2 , SO 2 , and perhaps cyanogen (4.57 μm) may be present as very small collections of molecules within the crystal structure, perhaps following models for radiation damage and/or for comet and interstellar grain formation at low temperatures. Some of the dark material on these surfaces may be created by radiation damage of the CO 2 and other carbon‐bearing species and the formation of graphite. These spectra suggest a complex chemistry within the surface materials and an important role for non‐ice materials in the evolution of the satellite surfaces.

Journal of Geophysical Research E: Planets↗

Photogrammetric analysis of horizon panoramas: The Pathfinder landing site in Viking orbiter images

Tiepoint measurements, block adjustment techniques, and sunrise/sunset pictures were used to obtain precise pointing data with respect to north for a set of 33 IMP horizon images. Azimuth angles for five prominent topographic features seen at the horizon were measured and correlated with locations of these features in Viking orbiter images. Based on this analysis, the Pathfinder line/sample coordinates in two raw Viking images were determined with approximate errors of 1 pixel, or 40 m. Identification of the Pathfinder location in orbit imagery yields geological context for surface studies of the landing site. Furthermore, the precise determination of coordinates in images together with the known planet-fixed coordinates of the lander make the Pathfinder landing site the most important anchor point in current control point networks of Mars.

Journal of Geophysical Research E: Planets↗

Erosional valleys in the Thaumasia region of Mars: Hydrothermal and seismic origins

Analysis of erosional valleys, geologic materials and features, and topography through time in the Thaumasia region of Mars using co-registered digital spatial data sets reveals significant associations that relate to valley origin. Valleys tend to originate (1) on Noachian to Early Hesperian (stages 1 and 2) large volcanoes, (2) within 50–100 km of stages 1 and 2 rift systems, and (3) within 100 km of Noachian (stage 1) impact craters >50 km in diameter. These geologic preferences explain observations of higher valley-source densities (VSDs) in areas of higher elevations and regional slopes (>1°) because the volcanoes, rifts, and craters form high, steep topography or occur in terrain of high relief. Other stage 1 and stage 2 high, steep terrains, however, do not show high VSDs. The tendency for valleys to concentrate near geologic features and the overall low drainage densities in Thaumasia compared to terrestrial surfaces rule out widespread precipitation as a major factor in valley formation (as is proposed in wann, wet climate scenarios) except perhaps during the Early Noachian, for which much of the geologic record has been obliterated. Instead, volcanoes and rifts may indicate the presence of shallow crustal intrusions that could lead to local hydrothermal circulation, melting of ground ice and snow, and groundwater sapping. However, impact-crater melt would provide a heat source at the surface that might drive away water, fonning valleys in the process. Post-stage 1 craters mostly have low nearby VSDs, which, for valleys incised in older rocks, suggests burial by ejecta and, for younger valleys, may indicate desiccation of near-surface water and deepening of the cryosphere. Later Hesperian and Amazonian (stages 3 and 4) valleys originate within 100–200 km of three young, large impact craters and near rifts systems at Warrego Valles and the southern part of Coprates rise. These valleys likely developed when the cryosphere was a couple kilometers or more thick, inhibiting valley development by hydrothermal circulation. However, eruption of groundwater may have occurred from impact-induced fracturing and lateral and perhaps minor upward transport of water due to seismic pumping. The two smaller craters formed along the plateau margin where the highest potential hydraulic head would occur in aquifers beneath the plateau. In the case of the larger crater (Lowell, 200 km in diameter), potential aquifers would likely be at depths of kilometers below the cryosphere. Seismic energy generated by the Lowell impactor would have been much greater, pumping both groundwater and perhaps fluidized slurry to the surface from beneath the cryosphere to form the young valleys and flow deposit. Along the margin of Thaumasia, tectonic pressurization of groundwater also may have contributed to valley formation. Dissection of rim materials of the Argyre impact may relate to tectonic activity and the unconsolidated state of basin ejecta.

Journal of Geophysical Research E: Planets↗

Rocky 7 prototype Mars rover field geology experiments: 1. Lavic Lake and Sunshine volcanic field, California

Experiments with the Rocky 7 rover were performed in the Mojave Desert to better understand how to conduct rover-based, long-distance (kilometers) geological traverses on Mars. The rover was equipped with stereo imaging systems for remote sensing science and hazard avoidance and 57 Fe Mössbauer and nuclear magnetic resonance spectrometers for in situ determination of mineralogy of unprepared rock and soil surfaces. Laboratory data were also obtained using the spectrometers and an X ray diffraction (XRD)/XRF instrument for unprepared samples collected from the rover sites. Simulated orbital and descent image data assembled for the test sites were found to be critical for assessing the geologic setting, formulating hypotheses to be tested with rover observations, planning traverses, locating the rover, and providing a regional context for interpretation of rover-based observations. Analyses of remote sensing and in situ observations acquired by the rover confirmed inferences made from orbital and simulated descent images that the Sunshine Volcanic Field is composed of basalt flows. Rover data confirmed the idea that Lavic Lake is a recharge playa and that an alluvial fan composed of sediments with felsic compositions has prograded onto the playa. Rover-based discoveries include the inference that the basalt flows are mantled with aeolian sediment and covered with a dense pavement of varnished basalt cobbles. Results demonstrate that the combination of rover remote sensing and in situ analytical observations will significantly increase our understanding of Mars and provide key connecting links between orbital and descent data and analyses of returned samples.

Journal of Geophysical Research E: Planets↗