Geology ReportsSearch

SEARCH · Geology Reports

Results for “Earth Imaging Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8Linked to original sources

U.S. Geological Survey spatial data access

The U.S. Geological Survey (USGS) has done a progress review on improving access to its spatial data holdings over the Web. The USGS EROS Data Center has created three major Web-based interfaces to deliver spatial data to the general public; they are Earth Explorer, the Seamless Data Distribution System (SDDS), and the USGS Web Mapping Portal. Lessons were learned in developing these systems, and various resources were needed for their implementation. The USGS serves as a fact-finding agency in the U.S. Government that collects, monitors, analyzes, and provides scientific information about natural resource conditions and issues. To carry out its mission, the USGS has created and managed spatial data since its inception. Originally relying on paper maps, the USGS now uses advanced technology to produce digital representations of the Earth’s features. The spatial products of the USGS include both source and derivative data. Derivative datasets include Digital Orthophoto Quadrangles (DOQ), Digital Elevation Models, Digital Line Graphs, land-cover Digital Raster Graphics, and the seamless National Elevation Dataset. These products, created with automated processes, use aerial photographs, satellite images, or other cartographic information such as scanned paper maps as source data. With Earth Explorer, users can search multiple inventories through metadata queries and can browse satellite and DOQ imagery. They can place orders and make payment through secure credit card transactions. Some USGS spatial data can be accessed with SDDS. The SDDS uses an ArcIMS map service interface to identify the user’s areas of interest and determine the output format; it allows the user to either download the actual spatial data directly for small areas or place orders for larger areas to be delivered on media. The USGS Web Mapping Portal provides views of national and international datasets through an ArcIMS map service interface. In addition, the map portal posts news about new map services available from the USGS, many simultaneously published on the Environmental Systems Research Institute Geography Network. These three information systems use new software tools and expanded hardware to meet the requirements of the users. The systems are designed to handle the required workload and are relatively easy to enhance and maintain. The software tools give users a high level of functionality and help the system conform to industry standards. The hardware and software architecture is designed to handle the large amounts of spatial data and Internet traffic required by the information systems. Last, customer support was needed to answer questions, monitor e-mail, and report customer problems.

Journal of Geospatial Engineering

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth

Overview of the Mars Pathfinder Mission: Launch through landing, surface operations, data sets, and science results

Mars Pathfinder successfully landed at Ares Vallis on July 4, 1997, deployed and navigated a small rover about 100 m clockwise around the lander, and collected data from three science instruments and ten technology experiments. The mission operated for three months and returned 2.3 Gbits of data, including over 16,500 lander and 550 rover images, 16 chemical analyses of rocks and soil, and 8.5 million individual temperature, pressure and wind measurements. Path‐finder is the best known location on Mars, having been clearly identified with respect to other features on the surface by correlating five prominent horizon features and two small craters in lander images with those in high‐resolution orbiter images and in inertial space from two‐way ranging and Doppler tracking. Tracking of the lander has fixed the spin pole of Mars, determined the precession rate since Viking 20 years ago, and indicates a polar moment of inertia, which constrains a central metallic core to be between 1300 and ∼2000 km in radius. Dark rocks appear to be high in silica and geochemically similar to anorogenic andesites; lighter rocks are richer in sulfur and lower in silica, consistent with being coated with various amounts of dust. Rover and lander images show rocks with a variety of morphologies, fabrics and textures, suggesting a variety of rock types are present. Rounded pebbles and cobbles on the surface as well as rounded bumps and pits on some rocks indicate these rocks may be conglomerates (although other explanations are also possible), which almost definitely require liquid water to form and a warmer and wetter past. Air‐borne dust is composed of composite silicate particles with a small fraction of a highly magnetic mineral, interpreted to be most likely maghemite; explanations suggest iron was dissolved from crustal materials during an active hydrologic cycle with maghemite freeze dried onto silicate dust grains. Remote sensing data at a scale of a kilometer or greater and an Earth analog correctly predicted a rocky plain safe for landing and roving with a variety of rocks deposited by catstrophic floods, which are relatively dust free. The surface appears to have changed little since it formed billions of years ago, with the exception that eolian activity may have deflated the surface by ∼3–7 cm, sculpted wind tails, collected sand into dunes, and eroded ventifacts (fluted and grooved rocks). Pathfinder found a dusty lower atmosphere, early morning water ice clouds, and morning near‐surface air temperatures that changed abruptly with time and height. Small scale vortices, interpreted to be dust devils, were observed repeatedly in the afternoon by the meteorology instruments and have been imaged.

Journal of Geophysical Research E: Planets

Displacement fields from point cloud data: Application of particle imaging velocimetry to landslide geodesy

Acquiring spatially continuous ground-surface displacement fields from Terrestrial Laser Scanners (TLS) will allow better understanding of the physical processes governing landslide motion at detailed spatial and temporal scales. Problems arise, however, when estimating continuous displacement fields from TLS point-clouds because reflecting points from sequential scans of moving ground are not defined uniquely, thus repeat TLS surveys typically do not track individual reflectors. Here, we implemented the cross-correlation-based Particle Image Velocimetry (PIV) method to derive a surface deformation field using TLS point-cloud data. We estimated associated errors using the shape of the cross-correlation function and tested the method's performance with synthetic displacements applied to a TLS point cloud. We applied the method to the toe of the episodically active Cleveland Corral Landslide in northern California using TLS data acquired in June 2005–January 2007 and January–May 2010. Estimated displacements ranged from decimeters to several meters and they agreed well with independent measurements at better than 9% root mean squared (RMS) error. For each of the time periods, the method provided a smooth, nearly continuous displacement field that coincides with independently mapped boundaries of the slide and permits further kinematic and mechanical inference. For the 2010 data set, for instance, the PIV-derived displacement field identified a diffuse zone of displacement that preceded by over a month the development of a new lateral shear zone. Additionally, the upslope and downslope displacement gradients delineated by the dense PIV field elucidated the non-rigid behavior of the slide.

California

A whole image approach using field measurements for transforming EO1 Hyperion hyperspectral data into canopy reflectance spectra

To maximize the spectral distinctiveness (information) of the canopy reflectance, an atmospheric correction strategy was implemented to provide accurate estimates of the intrinsic reflectance from the Earth Observing 1 (EO1) satellite Hyperion sensor signal. In rendering the canopy reflectance, an estimate of optical depth derived from a measurement of downwelling irradiance was used to drive a radiative transfer simulation of atmospheric scattering and attenuation. During the atmospheric model simulation, the input whole-terrain background reflectance estimate was changed to minimize the differences between the model predicted and the observed canopy reflectance spectra at 34 sites. Lacking appropriate spectrally invariant scene targets, inclusion of the field and predicted comparison maximized the model accuracy and, thereby, the detail and precision in the canopy reflectance necessary to detect low percentage occurrences of invasive plants. After accounting for artifacts surrounding prominent absorption features from about 400nm to 1000nm, the atmospheric adjustment strategy correctly explained 99% of the observed canopy reflectance spectra variance. Separately, model simulation explained an average of 88%??9% of the observed variance in the visible and 98% ?? 1% in the near-infrared wavelengths. In the 34 model simulations, maximum differences between the observed and predicted reflectances were typically less than ?? 1% in the visible; however, maximum reflectance differences higher than ?? 1.6% (<??2.3%) at more than a few wavelengths were observed at three sites. In the near-infrared wavelengths, maximum reflectance differences remained less than ??3% for 68% of the comparisons (??1 standard deviation) and less than ??6% for 95% of the comparisons (??2 standard deviation). Higher reflectance differences in the visible and near-infrared wavelengths were most likely associated with problems in the comparison, not in the model generation. ?? 2005 US Government.

International Journal of Remote Sensing

Overview of spirit microscopic imager results

This paper provides an overview of Mars Exploration Rover Spirit Microscopic Imager (MI) operations and the calibration, processing, and analysis of MI data. The focus of this overview is on the last five Earth years (2005–2010) of Spirit 's mission in Gusev crater, supplementing the previous overview of the first 450 sols of the Spirit MI investigation. Updates to radiometric calibration using in‐flight data and improvements in high‐level processing are summarized. Released data products are described, and a table of MI observations, including target/feature names and associated data sets, is appended. The MI observed natural and disturbed exposures of rocks and soils as well as magnets and other rover hardware. These hand‐lens‐scale observations have provided key constraints on interpretations of the formation and geologic history of features, rocks, and soils examined by Spirit . MI images complement observations by other Spirit instruments, and together show that impact and volcanic processes have dominated the origin and evolution of the rocks in Gusev crater, with aqueous activity indicated by the presence of silica‐rich rocks and sulfate‐rich soils. The textures of some of the silica‐rich rocks are similar to terrestrial hot spring deposits, and observations of subsurface cemented layers indicate recent aqueous mobilization of sulfates in places. Wind action has recently modified soils and abraded many of the rocks imaged by the MI, as observed at other Mars landing sites.

Journal of Geophysical Research E: Planets

Resolution of fault slip along the 470-km-long rupture of the great 1906 San Francisco earthquake and its implications

Data from all available triangulation networks affected by the 1906 earthquake have been combined to assess the trade-off between slip resolution and its uncertainty and to construct a conservative image of coseismic slip along the rupture. Because of varying network aperture and station density, slip resolution is very uneven. Although slip is determined within uncertainties of ??1.0 m along 60% of the fault, constraints are poor on the remaining, mostly offshore portions of the rupture. Slip decreases from maxima of 8.6 and 7.5 m at Shelter Cove and Tomales Bay to 4.5 m near Mount Tamalpais and 2.7 m at Loma Prieta. The geodetically derived slip distribution is in poor agreement with estimates based on analysis of S wave seismograms, probably because these waves register only 20-30% of the total seismic moment obtained from longer-period surface waves. Consideration of a range of fault geometries for 1906 slip near Loma Prieta indicates right-lateral motions lie between 2.3 and 3.1 m. These values are considerably greater than the 1.5 m of measured surface slip on which several assessments of high earthquake hazard for this fault segment were based. This factor, along with the absence of 1989 slippage where 1906 surface slip was used to make the forecasts, casts doubt on some claims of success in predicting the 1989 M = 6.9 Loma Prieta earthquake.

Journal of Geophysical Research B: Solid Earth

Characterization of intra-annual reflectance properties of land cover classes in southeastern South Dakota using Landsat TM and ETM+ data

Landsat-7 and Landsat-5 have orbits that are offset from each other by 8 days. During the time that the sensors on both satellites are operational, there is an opportunity for conducting analyses that incorporate multiple intra-annual high spatial resolution data sets for characterizing the Earth's land surface. In the current study, nine Landsat thematic mapper (TM) and enhanced thematic mapper plus (ETM+) data sets, covering the same path and row on different dates, were acquired during a 1-year time interval for a region in southeastern South Dakota and analyzed. Scenes were normalized using pseudoinvariant objects, and digital data from a series of test sites were extracted from the imagery and converted to surface reflectance. Sunphotometer data acquired on site were used to atmospherically correct the data. Ground observations that were made throughout the growing season by a large group of volunteers were used to help interpret spectroradiometric patterns and trends. Normalized images were found to be very effective in portraying the seasonal patterns of reflectance change that occurred throughout the region. Many of the radiometric patterns related to plant growth and development, but some also related to different background properties. The different kinds of land cover in the region were spectrally and radiometrically characterized and were found to have different seasonal patterns of reflectance. The degree to which the land cover classes could be separated spectrally and radiometrically, however, depended on the time of year during which the data sets were acquired, and no single data set appeared to be adequate for separating all types of land cover. This has practical implications for classification studies because known patterns of seasonal reflectance properties for the different types of land cover within a region will facilitate selection of the most appropriate data sets for producing land cover classifications.

Canadian Journal of Remote Sensing

Performance and effects of land cover type on synthetic surface reflectance data and NDVI estimates for assessment and monitoring of semi-arid rangeland

Federal land management agencies provide stewardship over much of the rangelands in the arid andsemi-arid western United States, but they often lack data of the proper spatiotemporal resolution andextent needed to assess range conditions and monitor trends. Recent advances in the blending of com-plementary, remotely sensed data could provide public lands managers with the needed information.We applied the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) to five Landsat TMand concurrent Terra MODIS scenes, and used pixel-based regression and difference image analyses toevaluate the quality of synthetic reflectance and NDVI products associated with semi-arid rangeland. Pre-dicted red reflectance data consistently demonstrated higher accuracy, less bias, and stronger correlationwith observed data than did analogous near-infrared (NIR) data. The accuracy of both bands tended todecline as the lag between base and prediction dates increased; however, mean absolute errors (MAE)were typically &le;10%. The quality of area-wide NDVI estimates was less consistent than either spectra lband, although the MAE of estimates predicted using early season base pairs were &le;10% throughout the growing season. Correlation between known and predicted NDVI values and agreement with the 1:1regression line tended to decline as the prediction lag increased. Further analyses of NDVI predictions,based on a 22 June base pair and stratified by land cover/land use (LCLU), revealed accurate estimates through the growing season; however, inter-class performance varied. This work demonstrates the successful application of the STARFM algorithm to semi-arid rangeland; however, we encourage evaluation of STARFM&rsquo;s performance on a per product basis, stratified by LCLU, with attention given to the influence of base pair selection and the impact of the time lag.

Idaho, Utah, Wyoming

An improved procedure for detection and enumeration of walrus signatures in airborne thermal imagery

In recent years, application of remote sensing to marine mammal surveys has been a promising area of investigation for wildlife managers and researchers. In April 2006, the United States and Russia conducted an aerial survey of Pacific walrus (Odobenus rosmarus divergens) using thermal infrared sensors to detect groups of animals resting on pack ice in the Bering Sea. The goal of this survey was to estimate the size of the Pacific walrus population. An initial analysis of the U.S. data using previously-established methods resulted in lower detectability of walrus groups in the imagery and higher variability in calibration models than was expected based on pilot studies. This paper describes an improved procedure for detection and enumeration of walrus groups in airborne thermal imagery. Thermal images were first subdivided into smaller 200 x 200 pixel "tiles." We calculated three statistics to represent characteristics of walrus signatures from the temperature histogram for each the. Tiles that exhibited one or more of these characteristics were examined further to determine if walrus signatures were present. We used cluster analysis on tiles that contained walrus signatures to determine which pixels belonged to each group. We then calculated a thermal index value for each walrus group in the imagery and used generalized linear models to estimate detection functions (the probability of a group having a positive index value) and calibration functions (the size of a group as a function of its index value) based on counts from matched digital aerial photographs. The new method described here improved our ability to detect walrus groups at both 2 m and 4 m spatial resolution. In addition, the resulting calibration models have lower variance than the original method. We anticipate that the use of this new procedure will greatly improve the quality of the population estimate derived from these data. This procedure may also have broader applicability to thermal infrared surveys of other wildlife species. Published by Elsevier B.V.

International Journal of Applied Earth Observation

Tomographic image of a seismically active volcano: Mammoth Mountain, California

High-resolution tomographic P wave, S wave, and V P / V S velocity structure models are derived for Mammoth Mountain, California, using phase data from the Northern California Seismic Network and a temporary deployment of broadband seismometers. An anomalous volume (5.1 &times; 10 9 to 5.9 &times; 10 10 m 3 ) of low P and low S wave velocities is imaged beneath Mammoth Mountain, extending from near the surface to a depth of &sim;2 km below sea level. We infer that the reduction in seismic wave velocities is due to the presence of CO 2 distributed in oblate spheroid pores with mean aspect ratio &alpha; = 1.6 &times; 10 &minus;3 to 7.9 &times; 10 &minus;3 (crack-like pores) and mean gas volume fraction ϕ = 8.1 &times; 10 &minus;4 to 3.4 &times; 10 &minus;3 . The pore density parameter &kappa; = 3 ϕ /(4&pi; &alpha; ) = na 3 =0.11, where n is the number of pores per cubic meter and a is the mean pore equatorial radius. The total mass of CO 2 is estimated to be 4.6 &times; 10 9 to 1.9 &times; 10 11 kg. The local geological structure indicates that the CO 2 contained in the pores is delivered to the surface through fractures controlled by faults and remnant foliation of the bedrock beneath Mammoth Mountain. The total volume of CO 2 contained in the reservoir suggests that given an emission rate of 500 tons day &minus;1 , the reservoir could supply the emission of CO 2 for &sim;25&ndash;1040 years before depletion. Continued supply of CO 2 from an underlying magmatic system would significantly prolong the existence of the reservoir.

California

Reconstructing the velocity and deformation of a rapid landslide using multiview video

Noncontact measurements of spatially varied ground surface deformation during landslide motion can provide important constraints on landslide mechanics. Here, we present and test a new method for extracting measurements of rapid landslide surface displacement and velocity (accelerations of approximately 1 m/s 2 ) using sequences of stereo images obtained from a pair of inexpensive, stationary 4K video cameras with nominal frame rates of 29.97 Hz. The method combines elements of Structure from Motion with those of optical flow to extract data on 3‐D evolution of the ground surface during slope failure. We apply the method to an experiment at the U.S. Geological Survey debris‐flow flume in which a high‐speed, liquefying landslide was triggered by gradually adding water to a 6‐m 3 prism of loosely packed sediment on a 31° slope. Strip‐scanning lidar measurements made during the experiment corroborate our video‐based measurements, but the latter encompassed the entire landslide surface and were much lower in cost. Our video‐based measurements enabled computation of depth‐integrated landslide dilation/contraction rates. The range of computed rates was within the ranges inferred from independent measurements of evolving pore water pressures and reasonable estimates of the hydraulic permeability of the sediment. Dilation and contraction rates play a crucial role in landslide mechanics. The dilation and contraction we observe contradict the incompressible flow assumption used in many studies that have employed noncontact methods to infer landslide properties.

Oregon

Galileo observations of Post-Imbrium lunar craters during the first Earth-Moon flyby

Copernican‐age craters are among the most conspicuous features seen on the far side and western limb of the Moon in the Galileo multispectral images acquired in December 1990. Among the new morphologic observations of far‐side craters are bright rays, continuous ejecta deposits, and dark rings associated with probable impact‐melt veneers. These observations suggest that the mapped age assignments of several large far‐side craters (Ohm, Robertson, and possibly Lowell and Lenz) need revision. New crater size‐frequency measurements on Lunar Orbiter images suggest the following age reassignments: Hausen (170 km diameter), Pythagoras (120 km), and Bullialdus (61 km) from Eratosthenian to Upper Imbrian, and Carpenter (60 km) and Harpalus (39 km) from Copernican to Eratosthenian. Colors and albedos of craters (away from impact‐melt veneers) are correlated with their geologic emplacement ages as determined from counts of superposed craters; these age‐color relations are used to estimate the emplacement age (time since impact event) for other Copernican‐age craters. These age‐color relations indicate a probable Copernican age for 27 far‐side or western limb craters larger than 10 km diameter that were not previously mapped as Copernican. The apparent deficiency of Copernican craters on the far side compared with the near side in published geologic maps is not present in our data. Age‐color trends differ between mare and highland regions and between the interiors and continuous ejecta of the craters. Similar trends are established for color and albedo versus soil‐maturity indices for the returned lunar samples, with distinct trends for mare and highland soils. However, the mare versus highland offsets are reversed in the two comparisons. These relations can be explained by variations in regolith thicknesses and rates of mixing with relatively fresh, crystalline ejecta. Therefore, the soil‐maturity trends represent longer geologic time periods in regions with thinner regoliths, such as the maria.

Journal of Geophysical Research E: Planets

Volcano dome dynamics at Mount St. Helens: Deformation and intermittent subsidence monitored by seismicity and camera imagery pixel offsets

The surface deformation field measured at volcanic domes provides insights into the effects of magmatic processes, gravity- and gas-driven processes, and the development and distribution of internal dome structures. Here we study short-term dome deformation associated with earthquakes at Mount St. Helens, recorded by a permanent optical camera and seismic monitoring network. We use Digital Image Correlation (DIC) to compute the displacement field between successive images and compare the results to the occurrence and characteristics of seismic events during a 6 week period of dome growth in 2006. The results reveal that dome growth at Mount St. Helens was repeatedly interrupted by short-term meter-scale downward displacements at the dome surface, which were associated in time with low-frequency, large-magnitude seismic events followed by a tremor-like signal. The tremor was only recorded by the seismic stations closest to the dome. We find a correlation between the magnitudes of the camera-derived displacements and the spectral amplitudes of the associated tremor. We use the DIC results from two cameras and a high-resolution topographic model to derive full 3-D displacement maps, which reveals internal dome structures and the effect of the seismic activity on daily surface velocities. We postulate that the tremor is recording the gravity-driven response of the upper dome due to mechanical collapse or depressurization and fault-controlled slumping. Our results highlight the different scales and structural expressions during growth and disintegration of lava domes and the relationships between seismic and deformation signals.

Washington

Seismic images of a tectonic subdivision of the Greenville Orogen beneath lakes Ontario and Erie

New seismic data from marine air-gun and Vibroseis profiles in Lake Ontario and Lake Erie provide images of subhorizontal Phanerozoic sediments underlain by a remarkable series of easterly dipping reflections that extends from the crystalline basement to the lower crust. These reflections are interpreted as structural features of crustal-scale subdivisions within the Grenville Orogen. Broadly deformed, imbricated, and overlapping thrust sheets within the western Central Metasedimentary Belt are succeeded to the west by a complex zone of easterly dipping, apparent thrust faults that are interpreted as a southwest subsurface extension of the boundary zone between the Central Metasedimentary Belt and the Central Gneiss Belt. The interpreted Central Metasedimentary Belt boundary zone has a characteristic magnetic anomaly that provides a link from the adjacent ends of lakes Ontario and Erie to structures exposed 150&ensp;km to the north. Less reflective, west-dipping events are interpreted as structures within the eastern Central Gneiss Belt. The seismic interpretation augments current tectonic models that suggest the exposed ductile structures formed at depth as a result of crustal shortening along northwest-verging thrust faults. Relatively shallow reflections across the boundary region suggest local, Late Proterozoic extensional troughs containing post-Grenville sediments, preserved possibly as a result of pre-Paleozoic reactivation of basement structures.

Lake Erie, Lake Ontario

The global short-period wavefield modelled with a Monte Carlo seismic phonon method

At high frequencies (∼1 Hz), much of the seismic energy arriving at teleseismic distances is not found in the main phases (e.g. P , PP , S , etc.) but is contained in the extended coda that follows these arrivals. This coda results from scattering off small-scale velocity and density perturbations within the crust and mantle and contains valuable information regarding the depth dependence and strength of this heterogeneity as well as the relative importance of intrinsic versus scattering attenuation. Most analyses of seismic coda to date have concentrated on S -wave coda generated from lithospheric scattering for events recorded at local and regional distances. Here, we examine the globally averaged vertical-component, 1-Hz wavefield (>10° range) for earthquakes recorded in the IRIS FARM archive from 1990 to 1999. We apply an envelope-function stacking technique to image the average time–distance behavior of the wavefield for both shallow (≤50 km) and deep (≥500 km) earthquakes. Unlike regional records, our images are dominated by P and P coda owing to the large effect of attenuation on PP and S at high frequencies. Modelling our results is complicated by the need to include a variety of ray paths, the likely contributions of multiple scattering and the possible importance of P -to- S and S -to- P scattering. We adopt a stochastic, particle-based approach in which millions of seismic phonons are randomly sprayed from the source and tracked through the Earth. Each phonon represents an energy packet that travels along the appropriate ray path until it is affected by a discontinuity or a scatterer. Discontinuities are modelled by treating the energy normalized reflection and transmission coefficients as probabilities. Scattering probabilities and scattering angles are computed in a similar fashion, assuming random velocity and density perturbations characterized by an exponential autocorrelation function. Intrinsic attenuation is included by reducing the energy contained in each particle as an appropriate function of traveltime. We find that most scattering occurs in the lithosphere and upper mantle, as previous results have indicated, but that some lower-mantle scattering is likely also required. A model with 3 to 4 per cent rms velocity heterogeneity at 4-km scale length in the upper mantle and 0.5 per cent rms velocity heterogeneity at 8-km scale length in the lower mantle (with intrinsic attenuation of Q α = 450 above 200 km depth and Q α = 2500 below 200 km) provides a reasonable fit to both the shallow- and deep-earthquake observations, although many trade-offs exist between the scale length, depth extent and strength of the heterogeneity.

Geophysical Journal International

Landsat TM and ETM+ thermal band calibration

Landsat-5 has been imaging the Earth since March 1984, and Landsat-7 was added to the series of Landsat instruments in April 1999. The Landsat Project Science Office and the Landsat-7 Image Assessment System have been monitoring the on-board calibration of Landsat-7 since launch. Additionally, two separate university teams have been evaluating the on-board thermal calibration of Landsat-7 through ground-based measurements since launch. Although not monitored as closely over its lifetime, a new effort is currently being made to validate the calibration of Landsat-5. Two university teams are beginning to collect ground truth under Landsat-5, along with using other vicarious calibration methods to go back into the archive to validate the history of the calibration of Landsat-5. This paper considers the calibration efforts for the thermal band, band 6, of both the Landsat-5 and Landsat-7 instruments. Though stable since launch, Landsat-7 had an initial calibration error of about 3 K, and changes were made to correct for this beginning 1 October 2000 for data processed with the National Landsat Archive Production System (NLAPS) and beginning 20 December 2000 for data processed with the Landsat Product Generation System (LPGS). Recent results from Landsat-5 vicarious calibration efforts show an offset of –0.7 K over the lifetime of the instrument. This suggests that historical calibration efforts may have been detecting errors in processing systems rather than changes in the instrument. A correction to the Landsat-5 processing has not yet been implemented but will be in the near future.

Canadian Journal of Remote Sensing

Radar and photoclinometric studies of wrinkle ridges on Mars

Earth-based radar altimetry and image derived photoclinometric profiles were analyzed to examine both the long- and short-wavelength topography associated with wrinkle ridges on Mars. Photoclinometrically derived elevation data across wrinkle ridges were evaluated to determine the sensitivity of profiles to two empirical photoclinometric parameters, the horizontal digital number (HDN) and the scattered light value (SLV). The photoclinometric profiles are extremely sensitive to small variations in HDN. The sense of slope of a profile can be completely reversed over a range in HDN of as little as ±1. Comparably small variations in the SLV have relatively minor effects on the photoclinometrically derived elevations. The existence of elevation offsets from one side of the ridge to the other, reported in previous photoclinometric studies of martian wrinkles, were not confirmed through photoclinometry. In addition, no evidence of elevation offsets were found in Earth-based radar altimetry profiles across wrinkle ridges. In order to more accurately model wrinkle ridge topography, we controlled photoclinometrically derived elevations with long-wavelength topography obtained from the radar altimetry. The results of this study do not support kinematic models for the origin of planetary wrinkle ridges that involve deeply rooted thrust faults which separate crustal blocks at different elevations. A kinematic model involving buckling of shallow crustal layers into concentric folds that close, leading to the development of thrust faults, is consistent with wrinkle ridge morphology and terrestrial analogs. Recent geophysical studies of terrestrial analogs and the influence of shallow subsurface structures, particularly buried craters, on the localization of many wrinkle ridges on Mars suggest that thrust faults associated with the ridges are confined to the ridged plains material and do not extend into the lithosphere.

Journal of Geophysical Research E: Planets