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,477 records · Page 82Linked to original sources
USGS1998

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↗
USGS1998

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↗
USGS1998

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↗
USGS1998

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↗
USGS1998

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↗
USGS1998

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↗
USGS1999

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↗
USGS1998

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↗
USGS1998

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↗
USGS1998

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↗
USGS1998

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↗
USGS1999

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↗
USGS1998

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↗
USGS1998

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↗
USGS1998

Classification of surface types using SIR-C/X-SAR, Mount Everest Area, Tibet

Imaging radar is a promising tool for mapping snow and ice cover in alpine regions. It combines a high-resolution, day or night, all-weather imaging capability with sensitivity to hydrologic and climatic snow and ice parameters. We use the spaceborne imaging radar-C/X-band synthetic aperture radar (SIR-C/X-SAR) to map snow and glacial ice on the rugged north slope of Mount Everest. From interferometrically derived digital elevation data, we compute the terrain calibration factor and cosine of the local illumination angle. We then process and terrain-correct radar data sets acquired on April 16, 1994. In addition to the spectral data, we include surface slope to improve discrimination among several surface types. These data sets are then used in a decision tree to generate an image classification. This method is successful in identifying and mapping scree/talus, dry snow, dry snow-covered glacier, wet snow-covered glacier, and rock-covered glacier, as corroborated by comparison with existing surface cover maps and other ancillary information. Application of the classification scheme to data acquired on October 7 of the same year yields accurate results for most surface types but underreports the extent of dry snow cover.

Journal of Geophysical Research E: Planets↗
USGS1998

Estimation of Mars radar backscatter from measured surface rock populations

Reanalysis of rock population data at the Mars Viking Lander sites has yielded updated values of rock fractional surface coverage (about 0.16 at both sites, including outcrops) and new estimates of rock burial depths and axial ratios. These data are combined with a finite difference time domain (FDTD) numerical scattering model to estimate diffuse backscatter due to rocks at both the Lander 1 (VL1) and Lander 2 (VL2) sites. We consider single scattering from both surface and subsurface objects of various shapes, ranging from an ideal sphere to an accurate digitized model of a terrestrial rock. The FDTD cross-section calculations explicitly account for the size, shape, composition, orientation, and burial state of the scattering object, the incident wave angle and polarization, and the composition of the surface. We calculate depolarized specific cross sections at 12.6 cm wavelength due to lossless rock-like scatterers of about 0.014 at VL1 and 0.023 at VL2, which are comparable to the measured ranges of 0.019–0.032 and 0.012–0.018, respectively. We also discuss the variation of the diffuse cross section as the local angle of incidence, θ i , changes. Numerical calculations for a limited set of rock shapes indicate a marked difference between the angular backscattering behavior of wavelength-scale surface and subsurface rocks: while subsurface rocks scatter approximately as a cosine power law, surface rocks display a complex variation, often with peak backscattering at high incidence angles (θ i = 70°–75°).

Journal of Geophysical Research E: Planets↗
USGS1999

Preliminary results on photometric properties of materials at the Sagan Memorial Station, Mars

Reflectance measurements of selected rocks and soils over a wide range of illumination geometries obtained by the Imager for Mars Pathfinder (IMP) camera provide constraints on interpretations of the physical and mineralogical nature of geologic materials at the landing site. The data sets consist of (1) three small “photometric spot” subframed scenes, covering phase angles from 20° to 150°; (2) two image strips composed of three subframed images each, located along the antisunrise and antisunset lines (photometric equator), covering phase angles from ∼0° to 155°; and (3) full‐image scenes of the rock “Yogi,” covering phase angles from 48° to 100°. Phase functions extracted from calibrated data exhibit a dominantly backscattering photometric function, consistent with the results from the Viking lander cameras. However, forward scattering behavior does appear at phase angles >140°, particularly for the darker gray rock surfaces. Preliminary efforts using a Hapke scattering model are useful in comparing surface properties of different rock and soil types but are not well constrained, possibly due to the incomplete phase angle availability, uncertainties related to the photometric function of the calibration targets, and/or the competing effects of diffuse and direct lighting. Preliminary interpretations of the derived Hapke parameters suggest that (1) red rocks can be modeled as a mixture of gray rocks with a coating of bright and dark soil or dust, and (2) gray rocks have macroscopically smoother surfaces composed of microscopically homogeneous, clear materials with little internal scattering, which may imply a glass‐like or varnished surface.

Journal of Geophysical Research E: Planets↗
USGS1999

Mars Pathfinder spectral measurements of Phobos and Deimos: Comparison with previous data

The Imager for Mars Pathfinder (IMP) acquired four spectra of parts of the sub‐Mars hemispheres of Phobos and Deimos. The measured region of Phobos is expected to be a mixture of the two spectral units identified on that satellite from Phobos 2 data, and the IMP spectra of Phobos are intermediate to the two units as expected. The derived geometric albedo is consistent with the value for that part of Phobos determined from Viking imagery. The IMP spectrum of Deimos is generally consistent with previous measurements acquired from the ground and from the Hubble Space Telescope (HST), but the signal‐to‐noise ratio is lower than that of the Phobos data. The spectral contrast between the two moons is similar to that seen in HST and Phobos 2 data. Mars Pathfinder measurements therefore substantiate recent results which indicate that Phobos and Deimos are not, as previously believed, analogous to C‐type asteroids. They also provide some indications for an absorption near 700 nm, perhaps like that seen in other low‐albedo asteroids. Both Martian moons are redder than most asteroids, and most closely resemble two analog materials believed to have undergone very dissimilar histories: primitive D‐like asteroids, and highly space‐weathered, mafic‐rich assemblages, such as are present in lunar mare soils.

Journal of Geophysical Research E: Planets↗