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,567 records · Page 87Linked to original sources
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

Time functions of deep earthquakes from broadband and short-period stacks

To constrain dynamic source properties of deep earthquakes, we have systematically constructed broadband time functions of deep earthquakes by stacking and scaling teleseismic P waves from U.S. National Seismic Network, TERRAscope, and Berkeley Digital Seismic Network broadband stations. We examined 42 earthquakes with depths from 100 to 660 km that occurred between July 1, 1992 and July 31, 1995. To directly compare time functions, or to group them by size, depth, or region, it is essential to scale them to remove the effect of moment, which varies by more than 3 orders of magnitude for these events. For each event we also computed short-period stacks of P waves recorded by west coast regional arrays. The comparison of broadband with short-period stacks yields a considerable advantage, enabling more reliable measurement of event duration. A more accurate estimate of the duration better constrains the scaling procedure to remove the effect of moment, producing scaled time functions with both correct timing and amplitude. We find only subtle differences in the broadband time-function shape with moment, indicating successful scaling and minimal effects of attenuation at the periods considered here. The average shape of the envelopes of the short-period stacks is very similar to the average broadband time function. The main variations seen with depth are (1) a mild decrease in duration with increasing depth, (2) greater asymmetry in the time functions of intermediate events compared to deep ones, and (3) unexpected complexity and late moment release for events between 350 and 550 km, with seven of the eight events in that depth interval displaying markedly more complicated time functions with more moment release late in the rupture than most events above or below. The first two results are broadly consistent with our previous studies, while the third is reported here for the first time. The greater complexity between 350 and 550 km suggests greater heterogeneity in the failure process in that depth range.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Evolution of the Gorda Escarpment, San Andreas fault and Mendocino triple junction from multichannel seismic data collected across the northern Vizcaino block, offshore northern California

The Gorda Escarpment is a north facing scarp immediately south of the Mendocino transform fault (the Gorda/Juan de Fuca-Pacific plate boundary) between 126°W and the Mendocino triple junction. It elevates the seafloor at the northern edge of the Vizcaino block, part of the Pacific plate, ∼1.5 km above the seafloor of the Gorda/Juan de Fuca plate to the north. Stratigraphy interpreted from multichannel seismic data across and close to the Gorda Escarpment suggests that the escarpment is a relatively recent pop-up feature caused by north-south compression across the plate boundary. Close to 126°W, the Vizcaino block acoustic basement shallows and is overlain by sediments that thin north toward the Gorda Escarpment. These sediments are tilted south and truncated at the seafloor. By contrast, in a localized region at the eastern end of the Gorda Escarpment, close to the Mendocino triple junction, the top of acoustic basement dips north and is overlain by a 2-km-thick wedge of pre-11 Ma sedimentary rocks that thickens north, toward the Gorda Escarpment. This wedge of sediments is restricted to the northeast corner of the Vizcaino block. Unless the wedge of sediments was a preexisting feature on the Vizcaino block before it was transferred from the North American to the Pacific plate, the strong spatial correlation between the sedimentary wedge and the triple junction suggests the entire Vizcaino block, with the San Andreas at its eastern boundary, has been part of the Pacific plate since significantly before 11 Ma.

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

Synthetic aperture radar interferometry coherence analysis over Katmai volcano group, Alaska

The feasibility of measuring volcanic deformation or monitoring deformation of active volcanoes using space-borne synthetic aperture radar (SAR) interferometry depends on the ability to maintain phase coherence over appropriate time intervals. Using ERS 1 C band (λ=5.66 cm) SAR imagery, we studied the seasonal and temporal changes of the interferometric SAR coherence for fresh lava, weathered lava, tephra with weak water reworking, tephra with strong water reworking, and fluvial deposits representing the range of typical volcanic surface materials in the Katmai volcano group, Alaska. For interferograms based on two passes with 35 days separation taken during the same summer season, we found that coherence increases after early June, reaches a peak between the middle of July and the middle of September, and finally decreases until the middle of November when coherence is completely lost for all five sites. Fresh lava has the highest coherence, followed by either weathered lava or fluvial deposits. These surfaces maintain relatively high levels of coherence for periods up to the length of the summer season. Coherence degrades more rapidly with time for surfaces covered with tephra. For images taken in different summers, only the lavas maintained coherence well enough to provide useful interferometric images, but we found only a small reduction in coherence after the first year for surfaces with lava. Measurement of volcanic deformation is possible using summer images spaced a few years apart, as long as the surface is dominated by lavas. Our studies suggest that in order to make volcanic monitoring feasible along the Aleutian arc or other regions with similar climatic conditions, observation intervals of the satellite with C band SAR should be at least every month from July through September, every week during the late spring/early summer or late fall, and every 2–3 days during the winter.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Mapping the rupture process of moderate earthquakes by inverting accelerograms

We present a waveform inversion method that uses recordings of small events as Green's functions to map the rupture growth of moderate earthquakes. The method fits P and S waveforms from many stations simultaneously in an iterative procedure to estimate the subevent rupture time and amplitude relative to the Green's function event. We invert the accelerograms written by two moderate Parkfield earthquakes using smaller events as Green's functions. The first earthquake ( M = 4.6) occurred on November 14, 1993, at a depth of 11 km under Middle Mountain, in the assumed preparation zone for the next Parkfield main shock. The second earthquake ( M = 4.7) occurred on December 20, 1994, some 6 km to the southeast, at a depth of 9 km on a section of the San Andreas fault with no previous microseismicity and little inferred coseismic slip in the 1966 Parkfield earthquake. The inversion results are strikingly different for the two events. The average stress release in the 1993 event was 50 bars, distributed over a geometrically complex area of 0.9 km 2 . The average stress release in the 1994 event was only 6 bars, distributed over a roughly elliptical area of 20 km 2 . The ruptures of both events appear to grow spasmodically into relatively complex shapes: the inversion only constrains the ruptures to grow more slowly than the S wave velocity but does not use smoothness constraints.

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

Seismic structure of the southern Cascadia subduction zone and accretionary prism north of the Mendocino triple junction

Four multichannel-seismic reflection profiles, collected as part of the Mendocino triple junction seismic experiment, image the toe of the southern Cascadia accretionary prism. Today, 250–600 m of sediment is subducting with the Gorda plate, and 1500–3200 m is accreting to the northern California margin. Faults imaged west and east of the deformation front show mixed structural vergence. A north-south trending, 20 km long portion of the central margin is landward vergent for the outer 6–8 km of the toe of the prism. This region of landward vergence exhibits no frontal thrust, is unusually steep and narrow, and is likely caused by a seaward-dipping backstop close to the deformation front. The lack of margin-wide preferred seaward vergence and wedge-taper analysis suggests the prism has low basal shear stress. The three southern lines image wedge-shaped fragments of oceanic crust 1.1–7.3 km in width and 250–700 m thick near the deformation front. These wedges suggest shortening and thickening of the upper oceanic crust. Discontinuities in the seafloor west of the prism provide evidence for mass wasting in the form of slump blocks and debris fans. The southernmost profile extends 75 km west of the prism imaging numerous faults that offset both the Gorda basin oceanic crust and overlying sediments. These high-angle faults, bounding basement highs, are interpreted as strike-slip faults reactivating structures originally formed at the spreading ridge. Northeast or northwest trending strike-slip faults within the basin are consistent with published focal mechanism solutions and are likely caused by north-south Gorda-Pacific plate convergence.

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↗