Geology Reports⌕ Search

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 181 records · Page 10Linked to original sources
USGS2014

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

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

Shallow conduit system at Kilauea Volcano, Hawaii, revealed by seismic signals associated with degassing bursts

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in March, 2008 and continuing to the present time is characterized by episodic explosive bursts of gas and ash from a vent within Halemaumau Pit Crater. These bursts are accompanied by seismic signals that are well recorded by a broadband network deployed in the summit caldera. We investigate in detail the dimensions and oscillation modes of the source of a representative burst in the 1−10 s band. An extended source is realized by a set of point sources distributed on a grid surrounding the source centroid, where the centroid position and source geometry are fixed from previous modeling of very-long-period (VLP) data in the 10–50 s band. The source time histories of all point sources are obtained simultaneously through waveform inversion carried out in the frequency domain. Short-scale noisy fluctuations of the source time histories between adjacent sources are suppressed with a smoothing constraint, whose strength is determined through a minimization of the Akaike Bayesian Information Criterion (ABIC). Waveform inversions carried out for homogeneous and heterogeneous velocity structures both image a dominant source component in the form of an east trending dike with dimensions of 2.9 × 2.9 km. The dike extends ∼2 km west and ∼0.9 km east of the VLP centroid and spans the depth range 0.2–3.1 km. The source model for a homogeneous velocity structure suggests the dike is hinged at the source centroid where it bends from a strike E 27°N with northern dip of 85° west of the centroid, to a strike E 7°N with northern dip of 80° east of the centroid. The oscillating behavior of the dike is dominated by simple harmonic modes with frequencies ∼0.2 Hz and ∼0.5 Hz, representing the fundamental mode ν 11 and first degenerate mode ν 12 = ν 21 of the dike. Although not strongly supported by data in the 1–10 s band, a north striking dike segment is required for enhanced compatibility with the model elaborated in the 10–50 s band. This dike provides connectivity between the east trending dike and the new vent within Halemaumau Pit Crater. Waveform inversions with a dual-dike model suggest dimensions of 0.7 × 0.7 km to 2.6 × 2.6 km for this segment. Further elaboration of the complex dike system under Halemaumau does not appear to be feasible with presently available data.

Hawai'i↗
USGS2014

An analysis of the global spatial variability of column-averaged CO 2 from SCIAMACHY and its implications for CO 2 sources and sinks

Satellite observations of carbon dioxide (CO 2 ) are important because of their potential for improving the scientific understanding of global carbon cycle processes and budgets. We present an analysis of the column-averaged dry air mole fractions of CO 2 (denoted XCO 2 ) of the Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY) retrievals, which were derived from a satellite instrument with relatively long-term records (2003–2009) and with measurements sensitive to the near surface. The spatial-temporal distributions of remotely sensed XCO 2 have significant spatial heterogeneity with about 6–8% variations (367–397 ppm) during 2003–2009, challenging the traditional view that the spatial heterogeneity of atmospheric CO 2 is not significant enough (<4%) to have any large effect on terrestrial ecosystem carbon cycles. By comparison with surface measurements from the NOAA Earth System Research Laboratory (ESRL) GLOBALVIEW network, significant statistical relationships between XCO 2 and surface CO 2 were found for major ecosystems, with the exception of tropical forest. In addition, when compared with a simulated terrestrial carbon uptake from the Integrated Biosphere Simulator (IBIS) and the Emissions Database for Global Atmospheric Research (EDGAR) carbon emission inventory, the latitudinal gradient of XCO 2 seasonal amplitude was influenced by the combined effect of terrestrial carbon uptake, carbon emission, and atmospheric transport, suggesting no direct implications for terrestrial carbon sinks. From the investigation of the growth rate of XCO 2 we found that the increase of CO 2 concentration was dominated by temperature in the northern hemisphere (20–90°N) and by precipitation in the southern hemisphere (20–90°S), with the major contribution to global average occurring in the northern hemisphere. These findings indicated that the satellite measurements of atmospheric CO 2 improve not only the estimations of atmospheric inversion, but also the understanding of the terrestrial ecosystem carbon dynamics and its feedback to atmospheric CO 2 .

International Journal of Remote Sensing↗
USGS2009

Rocky Mountain evolution: Tying Continental Dynamics of the Rocky Mountains and Deep Probe seismic experiments with receiver functions

In this study, we have determined the crustal structure using three different receiver function methods using data collected from the northern transect of the Continental Dynamics of the Rocky Mountains (CD-ROM) experiment. The resulting migrated image and crustal thickness determinations confirm and refine prior crustal thickness measurements based on the CD-ROM and Deep Probe experiment data sets. The new results show a very distinct and thick lower crustal layer beneath the Archean Wyoming province. In addition, we are able to show its termination at 42??N latitude, which provides a seismic tie between the CD-ROM and Deep Probe seismic experiments and thus completes a continuous north-south transect extending from New Mexico into Alberta, Canada. This new tie is particularly important because it occurs close to a major tectonic boundary, the Cheyenne belt, between an Archean craton and a Proterozoic terrane. We used two different stacking techniques, based on a similar concept but using two different ways to estimate uncertainties. Furthermore, we used receiver function migration and common conversion point (CCP) stacking techniques. The combined interpretation of all our results shows (1) crustal thinning in southern Wyoming, (2) strong northward crustal thickening beginning in central Wyoming, (3) the presence of an unusually thick and high-velocity lower crust beneath the Wyoming province, and (4) the abrupt termination of this lower crustal layer north of the Cheyenne belt at 42??N latitude. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Mapping cropland extent of Southeast and Northeast Asia using multi-year time-series Landsat 30-m data using Random Forest classifier on Google Earth Engine

Cropland extent maps are useful components for assessing food security. Ideally, such products are a useful addition to countrywide agricultural statistics since they are not politically biased and can be used to calculate cropland area for any spatial unit from an individual farm to various administrative unites (e.g., state, county, district) within and across nations, which in turn can be used to estimate agricultural productivity as well as degree of disturbance on food security from natural disasters and political conflict. However, existing cropland extent maps over large areas (e.g., Country, region, continent, world) are derived from coarse resolution imagery (250 m to 1 km pixels) and have many limitations such as missing fragmented and\or small farms with mixed signatures from different crop types and\or farming practices that can be, confused with other land cover. As a result, the coarse resolution maps have limited useflness in areas where fields are small (<1 ha), such as in Southeast Asia. Furthermore, coarse resolution cropland maps have known uncertainties in both geo-precision of cropland location as well as accuracies of the product. To overcome these limitations, this research was conducted using multi-date, multi-year 30-m Landsat time-series data for 3 years chosen from 2013 to 2016 for all Southeast and Northeast Asian Countries (SNACs), which included 7 refined agro-ecological zones (RAEZ) and 12 countries (Indonesia, Thailand, Myanmar, Vietnam, Malaysia, Philippines, Cambodia, Japan, North Korea, Laos, South Korea, and Brunei). The 30-m (1 pixel = 0.09 ha) data from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) were used in the study. Ten Landsat bands were used in the analysis (blue, green, red, NIR, SWIR1, SWIR2, Thermal, NDVI, NDWI, LSWI) along with additional layers of standard deviation of these 10 bands across 1 year, and global digital elevation model (GDEM)-derived slope and elevation bands. To reduce the impact of clouds, the Landsat imagery was time-composited over four time-periods (Period 1: January- April, Period 2: May-August, and Period 3: September-December) over 3-years. Period 4 was the standard deviation of all 10 bands taken over all images acquired during the 2015 calendar year. These four period composites, totaling 42 band data-cube, were generated for each of the 7 RAEZs. The reference training data (N = 7849) generated for the 7 RAEZ using sub-meter to 5-m very high spatial resolution imagery (VHRI) helped generate the knowledge-base to separate croplands from non-croplands. This knowledge-base was used to code and run a pixel-based random forest (RF) supervised machine learning algorithm on the Google Earth Engine (GEE) cloud computing environment to separate croplands from non-croplands. The resulting cropland extent products were evaluated using an independent reference validation dataset (N = 1750) in each of the 7 RAEZs as well as for the entire SNAC area. For the entire SNAC area, the overall accuracy was 88.1% with a producer’s accuracy of 81.6% (errors of omissions = 18.4%) and user’s accuracy of 76.7% (errors of commissions = 23.3%). For each of the 7 RAEZs overall accuracies varied from 83.2 to 96.4%. Cropland areas calculated for the 12 countries were compared with country areas reported by the United Nations Food and Agriculture Organization and other national cropland statistics resulting in an R 2 value of 0.93. The cropland areas of provinces were compared with the province statistics that showed an R 2 = 0.95 for South Korea and R 2 = 0.94 for Thailand. The cropland products are made available on an interactive viewer at www.croplands.org and for download at National Aeronautics and Space Administration’s (NASA) Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/node/1281 .

International Journal of Applied Earth Observation↗
USGS1978

EROS Data Center Landsat digital enhancement techniques and imagery availability

The US Geological Survey's EROS Data Center (EDC) is experimenting with the production of digitally enhanced Landsat imagery. Advanced digital image processing techniques are used to perform geometric and radiometric corrections and to perform contrast and edge enhancements. The enhanced image product is produced from digitally preprocessed Landsat computer compatible tapes (CCTs) on a laser beam film recording system. Landsat CCT data have several geometric distortions which are corrected when NASA produces the standard film products. When producing film images from CCT's, geometric correction of the data is required. The EDC Digital Image Enhancement System (EDIES) compensates for geometric distortions introduced by Earth's rotation, variable line length, non-uniform mirror scan velocity, and detector misregistration. Radiometric anomalies such as bad data lines and striping are common to many Landsat film products and are also in the CCT data. Bad data lines or line segments with more than 150 contiguous bad pixels are corrected by inserting data from the previous line in place of the bad data. Striping, caused by variations in detector gain and offset, is removed with a destriping algorithm applied after digitally enhancing the data. Image enhancement is performed by applying a linear contrast stretch and an edge enhancement algorithm. The linear contrast enhancement algorithm is designed to expand digitally the full range of useful data recorded on the CCT over the range of 256 digital counts. This minimizes the effect of atmospheric scattering and saturates the relative brightness of highly reflecting features such as clouds or snow. It is the intent that no meaningful terrain data are eliminated by the digital processing. The edge enhancement algorithm is designed to enhance boundaries between terrain features that exhibit subtle differences in brightness values along edges of features. After the digital data have been processed, data for each Landsat band are recorded on black-and-white film with a laser beam film recorder (LBR). The LBR corrects for aspect ratio distortions as the digital data are recorded on the recording film over a preselected density range. Positive transparencies of MSS bands 4, 5, and 7 produced by the LBR are used to make color composite transparencies. Color film positives are made photographically from first generation black-and-white products generated on the LBR.

Canadian Journal of Remote Sensing↗
USGS1996

Analysis of fractures intersecting Kahi Puka Well 1 and its relation to the growth of the island of Hawaii

As part of the Hawaii Scientific Drilling Project, Kahi Puka Well 1 penetrated about 275 m of Mauna Loa basalts overlying a sequence of Mauna Kea flow units as it was drilled and cored to a total depth of 1053 m below land surface. A borehole televiewer (BHTV) was run in most of the well in successive stages prior to casing in order to obtain magnetically oriented acoustic images of the borehole wall. A total of 283 individual fractures were identified from this log and characterized in terms of strike and dip. These data are divided into three vertical sections based upon age and volcanic source, and lower hemisphere stereographic plots identify two predominant, subparallel fracture subsets common to each section. Assuming that most of the steeply dipping fractures observed in the BHTV log are tensile features generated within basalt flows during deposition and cooling, this fracture information can be combined with models of the evolution of the island of Hawaii to investigate the depositional history of these Mauna Loa and Mauna Kea basalts over the past 400 kyr. The directions of high-angle fractures appear to be generally parallel to topography or to the coastline at the time of deposition, as is supported by surface mapping of modern flows. Consequently, an overall counterclockwise rotation of about 75?? in the strike of these fractures from the bottom to the top of the well represents a systematic change in depositional slope direction over time. We attribute the observed rotation in the orientations of the two predominant fracture subsets over the past 400 kyr to changes in the configurations of volcanic sources during shield building and to the structural interference of adjacent volcanoes that produces shifts in topographic patterns.

Journal of Geophysical Research B: Solid Earth↗
USGS1993

Oblique synoptic images, produced from digital data, display strong evidence of a "new" caldera in southwestern Guatemala

The synoptic view of broad regions of the Earth's surface as displayed in Landsat and other satellite images has greatly aided in the recognition of calderas, ignimbrite plateaus and other geologic landforms. Remote-sensing images that include visual representation of depth are an even more powerful tool for geologic interpretation of landscapes, but their use has been largely restricted to the exploration of planets other than Earth. By combining Landsat images with digitized topography, we have generated regional oblique views that display compelling evidence for a previously undocumented late-Cenozoic caldera within the active volcanic zone of southwestern Guatemala. This "new" caldera, herein called Xela, is a depression about 30 km wide and 400-600 m deep, which includes the Quezaltenango basin. The caldera depression is breached only by a single river canyon. The caldera outline is broadly circular, but a locally scalloped form suggests the occurrence of multiple caldera-collapse events, or local slumping of steep caldera walls, or both. Within its northern part, Xela caldera contains a toreva block, about 500 m high and 2 km long, that may be incompletely foundered pre-caldera bedrock. Xela contains several post-caldera volcanoes, some of which are active. A Bouguer gravity low, tens of milligals in amplitude, is approximately co-located with the proposed caldera. The oblique images also display an extensive plateau that dips about 2?? away from the north margin of Xela caldera. We interpret this landform to be underlain by pyroclastic outflow from Xela and nearby Atitla??n calderas. Field mapping by others has documented a voluminous rhyolitic pumiceous fallout deposit immediately east of Xela caldera. We speculate that Xela caldera was the source of this deposit. If so, the age of at least part of the caldera is between about 84 ka and 126 ka, the ages of deposits that stratigraphically bracket this fallout. Most of the floor of Xela caldera is covered with Los Chocoyos pyroclastics, 84-ka deposits erupted from Atitla??n caldera. Oblique images produced from digital data are unique tools that can greatly facilitate initial geologic interpretation of morphologically young volcanic (and other) terrains where field access is limited, especially because conventional visual representations commonly lack depth perspective and may cover only part of the region of interest. ?? 1993.

Journal of Volcanology and Geothermal Research↗
USGS2022

River bathymetry retrieval from Landsat-9 images based on neural networks and comparison to SuperDove and Sentinel-2

The Landsat mission has kept an eye on our planet, including water bodies, for 50 years. With the launch of Landsat-9 and its onboard Operational Land Imager 2 (OLI-2) in September 2021, more subtle variations in brightness (14-bit dynamic range) can be captured than previous sensors in the Landsat series (e.g., 12-bit Landsat-8). The enhanced radiometric resolution of OLI-2 appeals to the aquatic remote sensing community because the instrument might be capable of resolving smaller differences in water-leaving radiance. This study evaluates the potential to map river bathymetry from Landsat-9 imagery. We employ a neural network (NN)-based regression model for bathymetry retrieval and compare the results with optimal band ratio analysis (OBRA). The effect of Landsat-9 pan-sharpening on depth retrieval is also examined. In addition, we perform an intersensor comparison with Sentinel-2 and newly available 8-band SuperDoves from the PlanetScope constellation. Depth retrieval results from the Colorado and Potomac Rivers imply that Landsat-9 provided more accurate bathymetry across a range of depths up to 20 m, particularly when pan-sharpened. Downsampling the SuperDove data improved bathymetry retrieval due to enhanced signal-to-noise ratio, most notably in deep waters (maximum detectable depth increased from ∼15 to ∼20 m). Similarly, the enhanced spectral resolution of 8-band SuperDoves improved depth retrieval relative to 4-band Doves. The NN-based model outperformed OBRA by incorporating more spectral information.

Journal of Selected Topics in Applied Earth Observ↗
USGS2002

Sensitive high resolution ion microprobe (SHRIMP) detrital zircon geochronology provides new evidence for a hidden neoproterozoic foreland basin to the Grenville Orogen in the eastern Midwest, U.S.A

A sensitive high resolution ion microprobe (SHRIMP) was used in combination with backscattered electron (BSE) and cathodoluminescence (CL) images to determine the age of detrital zircons from sandstones in the Neoproterozoic Middle Run Formation of the eastern Midwest, United States. Eleven samples from seven drill cores of the upper part of the Middle Run Formation contain detrital zircons ranging in age from 1030 to 1982 Ma (84 analyses), with six distinctive modes at 1.96, 1.63, 1.47, 1.34, 1.15, and 1.08 Ga. This indicates that most, but not all, of the zircon at the top of the Middle Run Formation was derived from the Grenville Orogen. The youngest concordant detrital zircon yields a maximum age of 1048 ?? 22 Ma for the Middle Run Formation, indicating that the formation is younger than ca. 1026 Ma minus the added extra time needed for later uplift, denudation, thrusting, erosion, and transport to southwestern Ohio. Thus, as judged by proximity, composition, thickness, and geochronology, it is a North American equivalent to other Neoproterozoic Grenvillian-derived basins, such as the Torridon Group of Scotland and the Palmeiral Formation of South America. An alternate possibility, although much less likely in our opinion, is that it could be much younger, any time between 1048 ?? 22 Ma and the deposition of the Middle Cambrian Mount Simon Sandstone at about 510 Ma, and still virtually almost all derived from rocks of the Grenville Orogen.

Canadian Journal of Earth Sciences↗
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↗
USGS1996

Seismic tomography and deformation modeling of the junction of the San Andreas and Calaveras faults

Local earthquake P traveltime data is inverted to obtain a three-dimensional tomographic image of the region centered on the junction of the San Andreas and Calaveras faults. The resulting velocity model is then used to relocate more than 17,000 earthquakes and to produce a model of fault structure in the region. These faults serve as the basis for modeling the topography using elastic dislocation methods. The region is of interest because active faults join, it marks the transition zone from creeping to locked fault behavior on the San Andreas fault, it exhibits young topography, and it has a good spatial distribution of seismicity. The tomographic data set is extensive, consisting of 1445 events, 96 stations, and nearly 95,000 travel time readings. Tomographic images are resolvable to depths of 12 km and show significant velocity contrasts across the San Andreas and Calaveras faults, a low-velocity zone associated with the creeping section of the San Andreas fault, and shallow low-velocity sediments in the southern Santa Clara valley and northern Salinas valley. Relocated earthquakes only occur where v p > 5 km/s and indicate that portions of the San Andreas and Calaveras faults are non vertical, although we cannot completely exclude the possibility that all or part of this results from ray tracing problems. The new dips are more consistent with geological observations that dipping faults intersect the surface where surface traces have been mapped. The topographic modeling predicts extensive subsidence in regions characterized by shallow low-velocity material, presumably the result of recent sedimentation. Some details of the topography at the junction of the San Andreas and Calaveras faults are not consistent with the modeling results, suggesting that the current position of this “triple junction” has changed with time. The model also predicts those parts of the fault subject to contraction or extension perpendicular to the fault strike and hence the sense of any dip-slip component. In each locality the relative vertical motion across the fault is consistent with the fault dips found with the new hypocentral locations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Active high-resolution seismic tomography of compressional wave velocity and attenuation structure at Medicine Lake Volcano, northern California Cascade Range

We determine compressional wave velocity and attenuation structures for the upper crust beneath Medicine Lake volcano in northeast California using a high-resolution active source seismic tomography method. Medicine Lake volcano is a basalt through rhyolite shield volcano of the Cascade Range, lying east of the range axis. The Pg wave from eight explosive sources which has traveled upward through the target volume to a dense array of 140 seismographs provides 1- to 2-km resolution in the upper 5 to 7 km of the crust beneath the volcano. The experiment tests the hypothesis that Cascade Range volcanoes of this type are underlain only by small silicic magma chambers. We image a low-velocity low- Q region not larger than a few tens of cubic kilometers in volume beneath the summit caldera, supporting the hypothesis. A shallower high-velocity high-density feature, previously known to be present, is imaged for the first time in full plan view; it is east-west elongate, paralleling a topographic lineament between Medicine Lake volcano and Mount Shasta. This lineament is interpreted to be the result of an old crustal weakness now affecting the emplacement of magma, both on direct ascent from the lower crust and mantle and in migration from the shallow silicic chamber to summit vents. Differences between this high-velocity feature and the equivalent feature at Newbeny volcano, a volcano in central Oregon resembling Medicine Lake volcano, may partly explain the scarcity of surface hydrothermal features at Medicine Lake volcano. A major low-velocity low- Q feature beneath the southeast flank of the volcano, in an area with no Holocene vents, is interpreted as tephra, flows, and sediments from the volcano deeply ponded on the downthrown side of the Gillem fault, a normal fault mapped at the surface north of the volcano. A high- Q normal-velocity feature beneath the north rim of the summit caldera may be a small, possibly hot, subsolidus intrusion. A high-velocity low- Q region beneath the eastern caldera may be an area of boiling water between the magma chamber and the ponded east flank material. These structural data are useful both for understanding Cascade Range volcanism and for geothermal development in progress on the volcano.

Journal of Geophysical Research Solid Earth↗
USGS2003

Imaging spectroscopy: Earth and planetary remote sensing with the USGS Tetracorder and expert systems

Imaging spectroscopy is a tool that can be used to spectrally identify and spatially map materials based on their specific chemical bonds. Spectroscopic analysis requires significantly more sophistication than has been employed in conventional broadband remote sensing analysis. We describe a new system that is effective at material identification and mapping: a set of algorithms within an expert system decision‐making framework that we call Tetracorder. The expertise in the system has been derived from scientific knowledge of spectral identification. The expert system rules are implemented in a decision tree where multiple algorithms are applied to spectral analysis, additional expert rules and algorithms can be applied based on initial results, and more decisions are made until spectral analysis is complete. Because certain spectral features are indicative of specific chemical bonds in materials, the system can accurately identify and map those materials. In this paper we describe the framework of the decision making process used for spectral identification, describe specific spectral feature analysis algorithms, and give examples of what analyses and types of maps are possible with imaging spectroscopy data. We also present the expert system rules that describe which diagnostic spectral features are used in the decision making process for a set of spectra of minerals and other common materials. We demonstrate the applications of Tetracorder to identify and map surface minerals, to detect sources of acid rock drainage, and to map vegetation species, ice, melting snow, water, and water pollution, all with one set of expert system rules. Mineral mapping can aid in geologic mapping and fault detection and can provide a better understanding of weathering, mineralization, hydrothermal alteration, and other geologic processes. Environmental site assessment, such as mapping source areas of acid mine drainage, has resulted in the acceleration of site cleanup, saving millions of dollars and years in cleanup time. Imaging spectroscopy data and Tetracorder analysis can be used to study both terrestrial and planetary science problems. Imaging spectroscopy can be used to probe planetary systems, including their atmospheres, oceans, and land surfaces.

Journal of Geophysical Research E: Planets↗
USGS2017

Spectral matching techniques (SMTs) and automated cropland classification algorithms (ACCAs) for mapping croplands of Australia using MODIS 250-m time-series (2000–2015) data

Mapping croplands, including fallow areas, are an important measure to determine the quantity of food that is produced, where they are produced, and when they are produced (e.g. seasonality). Furthermore, croplands are known as water guzzlers by consuming anywhere between 70% and 90% of all human water use globally. Given these facts and the increase in global population to nearly 10 billion by the year 2050, the need for routine, rapid, and automated cropland mapping year-after-year and/or season-after-season is of great importance. The overarching goal of this study was to generate standard and routine cropland products, year-after-year, over very large areas through the use of two novel methods: (a) quantitative spectral matching techniques (QSMTs) applied at continental level and (b) rule-based Automated Cropland Classification Algorithm (ACCA) with the ability to hind-cast, now-cast, and future-cast. Australia was chosen for the study given its extensive croplands, rich history of agriculture, and yet nonexistent routine yearly generated cropland products using multi-temporal remote sensing. This research produced three distinct cropland products using Moderate Resolution Imaging Spectroradiometer (MODIS) 250-m normalized difference vegetation index 16-day composite time-series data for 16 years: 2000 through 2015. The products consisted of: (1) cropland extent/areas versus cropland fallow areas, (2) irrigated versus rainfed croplands, and (3) cropping intensities: single, double, and continuous cropping. An accurate reference cropland product (RCP) for the year 2014 (RCP2014) produced using QSMT was used as a knowledge base to train and develop the ACCA algorithm that was then applied to the MODIS time-series data for the years 2000–2015. A comparison between the ACCA-derived cropland products (ACPs) for the year 2014 (ACP2014) versus RCP2014 provided an overall agreement of 89.4% (kappa = 0.814) with six classes: (a) producer’s accuracies varying between 72% and 90% and (b) user’s accuracies varying between 79% and 90%. ACPs for the individual years 2000–2013 and 2015 (ACP2000–ACP2013, ACP2015) showed very strong similarities with several other studies. The extent and vigor of the Australian croplands versus cropland fallows were accurately captured by the ACCA algorithm for the years 2000–2015, thus highlighting the value of the study in food security analysis. The ACCA algorithm and the cropland products are released through http://croplands.org/app/map and http://geography.wr.usgs.gov/science/croplands/algorithms/australia_250m.html

International Journal of Digital Earth↗
USGS2002

High-resolution image of Calaveras Fault seismicity

By measuring relative earthquake arrival times using waveform cross correlation and locating earthquakes using the double difference technique, we are able to reduce hypocentral errors by 1 to 2 orders of magnitude over routine locations for nearly 8000 events along a 35-km section of the Calaveras Fault. This represents ∼92% of all seismicity since 1984 and includes the rupture zone of the M 6.2 1984 Morgan Hill, California, earthquake. The relocated seismicity forms highly organized structures that were previously obscured by location errors. There are abundant repeating earthquake sequences as well as linear clusters of earthquakes. Large voids in seismicity appear with dimensions of kilometers that have been aseismic over the 30-year time interval, suggesting that these portions of the fault are either locked or creeping. The area of greatest slip in the Morgan Hill main shock coincides with the most prominent of these voids, suggesting that this part of the fault may be locked between large earthquakes. We find that the Calaveras Fault at depth is extremely thin, with an average upper bound on fault zone width of 75 m. Given the location error, however, this width is not resolvably different from zero. The relocations reveal active secondary faults, which we use to solve for the stress field in the immediate vicinity of the Calaveras Fault. We find that the maximum compressive stress is at a high angle, only 13° from the fault normal, supporting previous interpretations that this fault is weak.

California↗
USGS2000

Fracture process zone in granite

In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 μm/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ∼9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ∼2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J.

Journal of Geophysical Research B: Solid Earth↗