SEARCH · USGS Science

Results for “Journal of Maps”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,135 records · Page 63Linked to original sources

USGS2024

Deformation by pressure solution and grain boundary sliding in a retrograde shear zone in southern New England, USA

Alleghanian phyllonites in a shear zone in southern New England were formed by the retrogression and hydration of a high-grade Acadian pelitic schist. The retrogression was locally incomplete, resulting in both heterogeneous mineralogy and mineral compositions, and in many arrested reaction textures. These features, documented by backscattered electron and element map imaging, and by 40 Ar/ 39 Ar dating, provide a record of the processes active as the schist was converted to a phyllonite. On the margins of the shear zone relic biotite and garnet grains survive but are partially replaced by chlorite. Relic high-grade muscovite flakes containing up to 35 % paragonite (pg) also persist but are entrained in swarms of pg-poor muscovite flakes producing folia with an anastomosing network of muscovite grains of variable composition. Single crystals of pg-rich muscovite are truncated by muscovite in successively younger cross-cutting folia that contain decreasing pg contents, falling to <5% pg in the youngest folia. In the core of the shear zone all high-grade minerals are destroyed, but recrystallization continues with earlier greenschist facies folia replaced by lower grade and pg-poorer muscovite in younger folia. Muscovite grains in truncated folia are commonly kinked and selective recrystallization of muscovite in kink bands to pg-poor-compositions demonstrates that strain energy helped drive recrystallization. The boundaries of these new pg-poor muscovite folia are sharp, and truncations indicate that grain boundary sliding (GBS) involved dissolution. GBS also enabled folding by flexural slip along muscovite grain boundaries. Fold amplitudes increased as chlorite dissolved from fold limbs precipitated in fold hinges producing crescent-shaped saddle-reef-like structures. Together these observations of truncation and replacement by dissolution and precipitation demonstrate that pressure solution and GBS facilitated both the formation of the phyllonitic shear zone and slip along it.

Connecticut↗

USGS2018

Methodology for correcting bottomhole temperatures acquired from wireline logging measurements in the onshore U.S. Gulf of Mexico Basin to characterize the thermal regime of total petroleum systems

Characterization of the subsurface thermal regime is critical for understanding many facets of the petroleum system, from thermal maturation of organic-rich source rocks to thermal preservation and non-degradation of hydrocarbon accumulations. On a broad scale, paleo-heatflow has been mapped for the North American continent (Blackwell and Richards, 2004) as well as the contiguous United States (Blackwell et al., 2011). However, in situ reservoir temperature is a fundamental property (Cooper and Jones, 1959) that is difficult to accurately measure in the subsurface (Deming, 1989). Previous work has described the thermal regime of the offshore U.S. Gulf of Mexico Basin (Waples et al., 2004; Forrest et al., 2005; Nagihara and Jones, 2005; Husson et al., 2008); however, due to the lack of an applicable bottomhole temperature (BHT) correction method, virgin rock temperatures of the onshore portion of the basin remains largely uncharacterized in a regional or subregional context. The abundance of BHT measurements offers a useful way to characterize the subsurface thermal environment, provided that they are corrected to reflect the reservoir temperature. This study develops BHT correction methods that are specifically calibrated for the onshore U.S. Gulf of Mexico Basin. These BHT corrections are empirically derived and are based on a newly compiled database of temperatures obtained from BHT wireline measurements and, to a lesser extent, from drill stem test (DST) data. The results of this investigation provide a unified BHT correction methodology for the onshore U.S. Gulf of Mexico Basin as well as provide 12 distinct BHT correction equations for each of the 12 physiographic provinces within the onshore Gulf Coast region. This study also characterizes the geothermal gradient regime across the onshore U.S. Gulf Coast, which ranges from 1.89ºF/100 ft in Sabine Uplift area to 1.39ºF/100 ft in the Southern Louisiana Salt Basin.

Gulf of Mexico Basin↗

USGS2020

Permafrost mapping with electrical resistivity tomography in two wetland systems north of the Tanana River, Interior Alaska

Surface-based 2D electrical resistivity tomography (ERT) surveys were used to characterize permafrost distribution at wetland sites on the alluvial plain north of the Tanana River, 20 km southwest of Fairbanks, Alaska, in June and September 2014. The sites were part of an ecologically-sensitive research area characterizing biogeochemical response of this region to warming and permafrost thaw, and the site contained landscape features characteristic of interior Alaska, including thermokarst bog, forested permafrost plateau, and a rich fen. The results show how vegetation reflects shallow (0–10 m depth) permafrost distribution. Additionally, we saw shallow (0–3 m depth) low resistivity areas in forested permafrost plateau potentially indicating the presence of increased unfrozen water content as a precursor to ground instability and thaw. Time-lapse study from June to September suggested a depth of seasonal influence extending several meters below the active layer, potentially as a result of changes in unfrozen water content. A comparison of several electrode geometries (dipole-dipole, extended dipole-dipole, Wenner-Schlumberger) showed that for depths of interest to our study (0–10 m) results were similar, but data acquisition time with dipole-dipole was the shortest, making it our preferred geometry. The results show the utility of ERT surveys to characterize permafrost distribution at these sites, and how vegetation reflects shallow permafrost distribution. These results are valuable information for ecologically sensitive areas where ground-truthing can cause excessive disturbance. ERT data can be used to characterize the exact subsurface geometry of permafrost such that over time an understanding of changing permafrost conditions can be made in great detail. Characterizing the depth of thaw and thermal influence from the surface in these areas also provides important information as an indication of the depth to which carbon storage and microbially-mediated carbon processing may be affected.

Alaska↗

USGS2006

Distribution of recent volcanism and the morphology of seamounts and ridges in the GLIMPSE study area: Implications for the lithospheric cracking hypothesis for the origin of intraplate, non-hot spot volcanic chains

Lithospheric cracking by remotely applied stresses or thermoelastic stresses has been suggested to be the mechanism responsible for the formation of intraplate volcanic ridges in the Pacific that clearly do not form above fixed hot spots. As part of the Gravity Lineations Intraplate Melting Petrology and Seismic Expedition (GLIMPSE) project designed to investigate the origin of these features, we have mapped two volcanic chains that are actively forming to the west of the East Pacific Rise using multibeam echo sounding and side-scan sonar. Side-scan sonar reveals the distribution of rough seafloor corresponding to recent, unsedimented lava flows. In the Hotu Matua volcanic complex, recent flows and volcanic edifices are distributed over a region 450 km long and up to 65 km wide, with an apparent, irregular age progression from older flows in the west to younger in the east. The 550-km-long Southern Cross Seamount/Sojourn Ridge/Brown Ridge chain appears to have been recently active only at its eastern end near the East Pacific Rise. A third region of recent flows is found 120 km north of Southern Cross Seamount in seafloor approximately 9 Myr old. No indication of lithospheric extension in the form of faulting or graben formation paralleling the trend of the volcanic chains is found in the vicinity of recent flows or anywhere else in the study area. Thermoelastic cracking could be a factor in the formation of a few small, very narrow volcanic ridges, but most of the volcanic activity is broadly distributed in wide swaths with no indication of formation along narrow cracks. The Sojourn and Brown chains appear to begin as distributed zones of small seamounts that later develop into segmented ridges, perhaps under the influence of membrane stresses from self-loading. We suggest that the linear volcanic chains are created by moving melting anomalies in the asthenosphere and that lithospheric cracking plays at most a secondary role. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Morphology and distribution of seamounts surrounding Easter Island

We investigate the morphology and distribution of a seamount population on a section of seafloor influenced by both superfast seafloor spreading and hotspot volcanism. The population under investigation is part of a broad chain of seamounts extending eastward from the East Pacific Rise, near Easter Island. In order to define the morphological variability of the seamounts, basal shape, cross-sectional area, volume, flatness, and flank slope are plotted against height for 383 seamounts with heights greater than 200 m, based on bathymetry data collected by GLORI-B and SeaBeam 2000, during three cruises onboard the R/V Melville in the spring of 1993. Nearly complete swath mapping coverage of the seamounts is available for the analysis of size and shape distribution. We quantitatively describe the seamount population of this active region, in which seamounts cover ???27% of the seafloor, and account for ???4.2% of the total crustal volume. Over 50% of the total volume (61,000 km3) of seamounts used in this study is made up by the 14 largest seamounts, and the remaining volume is made up by the 369 smaller seamounts (>200 m in height). Our analysis indicates there are at least two seamount populations in the Easter Island-Salas y Gomez Island (25??-29??S, 113??-104??W) study area. One population of seamounts is composed of short seamounts (<1200 m in height) with variable flatness from pointy cones to flattened domes (flatness from 0.01 to 0.57) and predominantly steep flanks (slopes from 5?? to 32??). A second population is of massive (>1200 m), shield-like, pointy cones (flatness < 0.2) and gentle slopes (from 5?? to 15??). An exponential maximum likelihood distribution is fit to the binned raw frequency of height and gives a characteristic height of the seamount population of 308??12m and an expected number of seamounts per 1000 km2 of 2.7??0.15. Many seamounts that have different slope and flatness relationships with height are formed next to each other. We speculate that the larger volcanoes (>???1200 m) originate exclusively from a hotspot source, but only a portion of the smaller volcanoes (<???1200 m) are formed from a hotspot source. The remainder would be presumably formed by a normal mantle or mixed source.

Journal of Geophysical Research B: Solid Earth↗

USGS2008

Estimating watershed level nonagricultural pesticide use from golf courses using geospatial methods

Limited information exists on pesticide use for nonagricultural purposes, making it difficult to estimate pesticide loadings from nonagricultural sources to surface water and to conduct environmental risk assessments. A method was developed to estimate the amount of pesticide use on recreational turf grasses, specifically golf course turf grasses, for watersheds located throughout the conterminous United States (U.S.). The approach estimates pesticide use: (1) based on the area of recreational turf grasses (used as a surrogate for turf associated with golf courses) within the watershed, which was derived from maps of land cover, and (2) from data on the location and average treatable area of golf courses. The area of golf course turf grasses determined from these two methods was used to calculate the percentage of each watershed planted in golf course turf grass (percent crop area, or PCA). Turf-grass PCAs derived from the two methods were used with recommended application rates provided on pesticide labels to estimate total pesticide use on recreational turf within 1,606 watersheds associated with surface-water sources of drinking water. These pesticide use estimates made from label rates and PCAs were compared to use estimates from industry sales data on the amount of each pesticide sold for use within the watershed. The PCAs derived from the land-cover data had an average value of 0.4% of a watershed with minimum of 0.01% and a maximum of 9.8%, whereas the PCA values that are based on the number of golf courses in a watershed had an average of 0.3% of a watershed with a minimum of <0.01% and a maximum of 14.2%. Both the land-cover method and the number of golf courses method produced similar PCA distributions, suggesting that either technique may be used to provide a PCA estimate for recreational turf. The average and maximum PCAs generally correlated to watershed size, with the highest PCAs estimated for small watersheds. Using watershed specific PCAs, combined with label rates, resulted in greater than two orders of magnitude over-estimation of the pesticide use compared to estimates from sales data. ?? 2008 American Water Resources Association.

Journal of the American Water Resources Associatio↗

USGS2006

Volcanic history and 40 Ar/ 39 Ar and 14 C geochronology of Terceira Island, Azores, Portugal

Seven new 40 Ar/ 39 Ar and 23 new radiocarbon ages of eruptive units, in support of new geologic mapping, improve the known chronology of Middle to Late Pleistocene and Holocene volcanic activity on the island of Terceira, Azores and define an east-to-west progression in stratovolcano growth. The argon ages indicate that Cinco Picos Volcano, the oldest on Terceira, completed its main subaerial cone building activity by about 370–380 ka. Collapse of the upper part of the stratovolcanic edifice to form a 7 × 9 km caldera occurred some time after 370 ka. Postcaldera eruptions of basalt from cinder cones on and near the caldera floor and trachytic pyroclastic flow and pumice fall deposits from younger volcanoes west of Cinco Picos have refilled much of the caldera. The southern portion of Guilherme Moniz Volcano, in the central part of the island, began erupting prior to 270 ka and produced trachyte domes, flows, and minor pyroclastic deposits until at least 111 ka. The northern part of Guilherme Moniz Caldera is less well exposed than the southern part, but reflects a similar age range. The northwest portion of the caldera was formed sometime after 44 ka. Several well-studied ignimbrites that blanket much of the island likely erupted from Guilherme Moniz Volcano. The Pico Alto Volcanic Center, a tightly spaced cluster of trachyte domes and short flows, is a younger part of Guilherme Moniz Volcano. Stratigraphic studies and our new radiocarbon ages suggest that most of the Pico Alto eruptions occurred during the period from about 9000 to 1000 years BP. Santa Barbara Volcano is the youngest stratovolcano on Terceira, began erupting prior to 29 ka, and has been active historically.

Azores↗

USGS2001

Seismic mapping of shallow fault zones in the San Gabriel Mountains from the Los Angeles Region Seismic Experiment, southern California

During the Los Angeles Region Seismic Experiment (LARSE), a reflection/refraction survey was conducted along a profile (line 1) extending from Seal Beach, California, northeastward to the Mojave Desert and crossing the Los Angeles and San Gabriel Valley basins and San Gabriel Mountains. In most shot gathers from the southern and central San Gabriel Mountains, clear secondary arrivals are seen that merge, or appear to merge, with first arrivals at three locations, including the location of the Vincent thrust fault, an exposed late Mesozoic/early Cenozoic megathrust. These secondary arrivals are interpretable as reflections in the shallow crust (<5 km depth) from a concave-upward interface that projects to the surface in the north near the Vincent thrust fault, is offset in its central part at the San Gabriel fault (an old branch of the San Andreas fault), and terminates in the south at 1 to 2 km depth at the southern mountain front. The velocity structure above and below this interface strongly suggests it is the Vincent thrust fault: intermediate velocities (6.2 km/s), consistent with mylonites overlying the Vincent thrust fault, are observed above it; lower velocities (5.8 km/s), consistent with the Pelona Schist underlying the Vincent thrust fault, are observed below it. Problems arise, however, in attempting to match this reflector to the exposed Vincent thrust fault, which is seen in outcrops east of line 1. The Vincent thrust fault is shallower than the reflector in most places. An unmapped structure (steep fault, monocline, or thrust fault) is required between line 1 and the outcrops that either drops the Vincent thrust fault down to the depths of the reflector or repeats the Vincent thrust fault beneath line 1 in the footwall of another thrust fault. An alternative interpretation of the reflector is a deep greenstone horizon within the Pelona Schist, although this alternative is not favored by the velocity structure. Copyright 2001 by the American Geophysical Union.

California↗

USGS2000

Distribution of surficial sediment in Long Island Sound and adjacent waters: Texture and total organic carbon

The surficial sediment distribution within Long Island Sound has been mapped and described using bottom samples, photography, and sidescan sonar, combined with information from the geologic literature. The distributions of sediment type and total organic carbon (TOC) reveal several broad trends that are largely related to the sea-floor geology, the bathymetry, and the effects of modern tidal- and wind-driven currents. Sediment types are most heterogeneous in bathymetrically complex and shallow nearshore areas; the heterogeneity diminishes and the texture fines with decreasing bottom-current energy. Lag deposits of gravel and gravelly sand dominate the surficial sediment texture in areas where bottom currents are the strongest (such as where tidal flow is constricted) and where glacial till crops out at the sea floor. Sand is the dominant sediment type in areas characterized by active sediment transport and in shallow areas affected by fine-grained winnowing. Silty sand and sand-silt-clay mark transitions within the basin from higher- to lower-energy environments, suggesting a diminished hydraulic ability to sort and transport sediment. Clayey silt and silty clay are the dominant sediment types accumulating in the central and western basins and in other areas characterized by long-term depositional environments. The amount of TOC in the sediments of Long Island Sound varies inversely with sediment grain size. Concentrations average more than 1.9% (dry weight) in clayey silt, but are less than 0.4% in sand. Generally, values for TOC increase both toward the west in the Sound and from the shallow margins to the deeper parts of the basin floor. Our data also suggest that TOC concentrations can vary seasonally.

Long Island Sound↗

USGS2026

Aligning legacy NLCD land cover maps based on Landsat Collection 1 to Collection 2

The transition from Landsat Collection 1 to Collection 2 introduced significant improvements in radiometric and geometric accuracy. However, the improvements cause location misalignment between the existing Landsat-derived land cover products and the new collection. The legacy National Land Cover Database (NLCD) has been used as a cornerstone land cover source for a variety of research. Therefore, a method aligning the legacy NLCD product to Collection 2 is required to ensure its continuity and consistency of service. We developed a strategy to not only align legacy NLCD to match new Collection 2 geometric locations but also improve land cover labeling in the region that was affected by the geometric shifts. The method identifies boundary pixels of homogeneous land cover patches as potential problem areas that are likely impacted by geometric shifts and generates candidate labels from 3 × 3 window with the target pixel at the center and segmentation-derived majority label. Standard phenology patterns of each candidate land cover type are established based on the random samples except boundary pixels within a 1000-pixels × 1000-pixels processing window region. The phenological distance to each standard land cover type pattern is calculated through a penalty dynamic time warping (DTW) method for each target pixel in the boundary region. Finally, the method determines the most suitable label based on the phenological distance from the candidate labels. Both visual and accuracy assessment results demonstrate that the alignment preserves the overall land cover patterns in the original legacy NLCD product while reducing the spatial discrepancies between the Landsat Collection 2 and land cover. In addition, it enhances the accuracy of land cover labeling of boundary pixels. The overall accuracy (OA) was increased by 7% in the land cover boundary regions after alignment. The quality and confusion matrix comparison between the alignment results and the original legacy NLCD confirm the reliability of the method. Our alignment method has the potential to serve as a framework for aligning other Landsat-derived land cover products to future collections.

conterminous United States↗

USGS2010

Irrigated areas of India derived using MODIS 500 m time series for the years 2001-2003

The overarching goal of this research was to develop methods and protocols for mapping irrigated areas using a Moderate Resolution Imaging Spectroradiometer (MODIS) 500 m time series, to generate irrigated area statistics, and to compare these with ground- and census-based statistics. The primary mega-file data-cube (MFDC), comparable to a hyper-spectral data cube, used in this study consisted of 952 bands of data in a single file that were derived from MODIS 500 m, 7-band reflectance data acquired every 8-days during 2001-2003. The methods consisted of (a) segmenting the 952-band MFDC based not only on elevation-precipitation-temperature zones but on major and minor irrigated command area boundaries obtained from India's Central Board of Irrigation and Power (CBIP), (b) developing a large ideal spectral data bank (ISDB) of irrigated areas for India, (c) adopting quantitative spectral matching techniques (SMTs) such as the spectral correlation similarity (SCS) R 2 -value, (d) establishing a comprehensive set of protocols for class identification and labeling, and (e) comparing the results with the National Census data of India and field-plot data gathered during this project for determining accuracies, uncertainties and errors. The study produced irrigated area maps and statistics of India at the national and the subnational (e.g., state, district) levels based on MODIS data from 2001-2003. The Total Area Available for Irrigation (TAAI) and Annualized Irrigated Areas (AIAs) were 113 and 147 million hectares (MHa), respectively. The TAAI does not consider the intensity of irrigation, and its nearest equivalent is the net irrigated areas in the Indian National Statistics. The AIA considers intensity of irrigation and is the equivalent of "irrigated potential utilized (IPU)" reported by India's Ministry of Water Resources (MoWR). The field-plot data collected during this project showed that the accuracy of TAAI classes was 88% with a 12% error of omission and 32% of error of commission. Comparisons between the AIA and IPU produced an R 2 -value of 0.84. However, AIA was consistently higher than IPU. The causes for differences were both in traditional approaches and remote sensing. The causes of uncertainties unique to traditional approaches were (a) inadequate accounting of minor irrigation (groundwater, small reservoirs and tanks), (b) unwillingness to share irrigated area statistics by the individual Indian states because of their stakes, (c) absence of comprehensive statistical analyses of reported data, and (d) subjectivity involved in observation-based data collection process. The causes of uncertainties unique to remote sensing approaches were (a) irrigated area fraction estimate and related sub-pixel area computations and (b) resolution of the imagery. The causes of uncertainties common in both traditional and remote sensing approaches were definitions and methodological issues. ?? 2009 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS).

ISPRS Journal of Photogrammetry and Remote Sensing↗

USGS2016

Optimizing selection of training and auxiliary data for operational land cover classification for the LCMAP initiative

The U.S. Geological Survey’s Land Change Monitoring, Assessment, and Projection (LCMAP) initiative is a new end-to-end capability to continuously track and characterize changes in land cover, use, and condition to better support research and applications relevant to resource management and environmental change. Among the LCMAP product suite are annual land cover maps that will be available to the public. This paper describes an approach to optimize the selection of training and auxiliary data for deriving the thematic land cover maps based on all available clear observations from Landsats 4–8. Training data were selected from map products of the U.S. Geological Survey’s Land Cover Trends project. The Random Forest classifier was applied for different classification scenarios based on the Continuous Change Detection and Classification (CCDC) algorithm. We found that extracting training data proportionally to the occurrence of land cover classes was superior to an equal distribution of training data per class, and suggest using a total of 20,000 training pixels to classify an area about the size of a Landsat scene. The problem of unbalanced training data was alleviated by extracting a minimum of 600 training pixels and a maximum of 8000 training pixels per class. We additionally explored removing outliers contained within the training data based on their spectral and spatial criteria, but observed no significant improvement in classification results. We also tested the importance of different types of auxiliary data that were available for the conterminous United States, including: (a) five variables used by the National Land Cover Database, (b) three variables from the cloud screening ‘‘Function of mask” (Fmask) statistics, and (c) two variables from the change detection results of CCDC. We found that auxiliary variables such as a Digital Elevation Model and its derivatives (aspect, position index, and slope), potential wetland index, water probability, snow probability, and cloud probability improved the accuracy of land cover classification. Compared to the original strategy of the CCDC algorithm (500 pixels per class), the use of the optimal strategy improved the classification accuracies substantially (15-percentage point increase in overall accuracy and 4-percentage point increase in minimum accuracy).

ISPRS Journal of Photogrammetry and Remote Sensing↗

USGS1996

Identifying coseismic subsidence in tidal-wetland stratigraphic sequences at the Cascadia subduction zone of western North America

Tidal-wetland stratigraphy reveals that great plate boundary earthquakes have caused hundreds of kilometers of coast to subside at the Cascadia subduction zone. However, determining earthquake recurrence intervals and mapping the coastal extent of past great earthquake ruptures in this region are complicated by the effects of many sedimentologic, hydrographic, and oceanographic processes that occur on the coasts of tectonically passive as well as active continental margins. Tidal-wetland stratigraphy at many Cascadia estuaries differs little from that at similar sites on passive-margin coasts where stratigraphic sequences form through nonseismic processes unrelated to coseismic land level changes. Methods developed through study of similar stratigraphic sequences in Europe provide a framework for investigating the Cascadia estuarine record. Five kinds of criteria must be evaluated when inferring regional coastal subsidence due to great plate boundary earthquakes: the suddenness and amount of submergence, the lateral extent of submerged tidal-wetland soils, the coincidence of submergence with tsunami deposits, and the degree of synchroneity of submergence events at widely spaced sites. Evaluation of such criteria at the Cascadia subduction zone indicates regional coastal subsidence during at least two great earthquakes. Evidence for a coseismic origin remains equivocal, however, for the many peat-mud contacts in Cascadia stratigraphic sequences that lack (1) contrasts in lithology or fossils indicative of more than half a meter of submergence, (2) well-studied tsunami deposits, or (3) precise ages needed for regional correlation. Paleoecologic studies of fossil assemblages are particularly important in estimating the size of sudden sea level changes recorded by abrupt peat-mud contacts and in helping to distinguish erosional and gradually formed contacts from coseismic contacts. Reconstruction of a history of great earthquakes for the Cascadia subduction zone will require rigorous application of the above criteria and many detailed investigations.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Developing framework to constrain the geometry of the seismic rupture plane on subduction interfaces a priori - A probabilistic approach

A key step in many earthquake source inversions requires knowledge of the geometry of the fault surface on which the earthquake occurred. Our knowledge of this surface is often uncertain, however, and as a result fault geometry misinterpretation can map into significant error in the final temporal and spatial slip patterns of these inversions. Relying solely on an initial hypocentre and CMT mechanism can be problematic when establishing rupture characteristics needed for rapid tsunami and ground shaking estimates. Here, we attempt to improve the quality of fast finite-fault inversion results by combining several independent and complementary data sets to more accurately constrain the geometry of the seismic rupture plane of subducting slabs. Unlike previous analyses aimed at defining the general form of the plate interface, we require mechanisms and locations of the seismicity considered in our inversions to be consistent with their occurrence on the plate interface, by limiting events to those with well-constrained depths and with CMT solutions indicative of shallow-dip thrust faulting. We construct probability density functions about each location based on formal assumptions of their depth uncertainty and use these constraints to solve for the ‘most-likely’ fault plane. Examples are shown for the trench in the source region of the M w 8.6 Southern Sumatra earthquake of March 2005, and for the Northern Chile Trench in the source region of the November 2007 Antofagasta earthquake. We also show examples using only the historic catalogues in regions without recent great earthquakes, such as the Japan and Kamchatka Trenches. In most cases, this method produces a fault plane that is more consistent with all of the data available than is the plane implied by the initial hypocentre and CMT mechanism. Using the aggregated data sets, we have developed an algorithm to rapidly determine more accurate initial fault plane geometries for source inversions of future earthquakes.

Geophysical Journal International↗

USGS2013

Geochemistry, petrography, and zircon U-Pb geochronology of Paleozoic metaigneous rocks in the Mount Veta area of east-central Alaska: implications for the evolution of the westernmost part of the Yukon-Tanana terrane

We report the results of new mapping, whole-rock major, minor, and trace-element geochemistry, and petrography for metaigneous rocks from the Mount Veta area in the westernmost part of the allochthonous Yukon–Tanana terrane (YTT) in east-central Alaska. These rocks include tonalitic mylonite gneiss and mafic metaigneous rocks from the Chicken metamorphic complex and the Nasina and Fortymile River assemblages. Whole-rock trace-element data from the tonalitic gneiss, whose igneous protolith was dated by SHRIMP U–Pb zircon geochronology at 332.6 ± 5.6 Ma, indicate derivation from tholeiitic arc basalt. Whole-rock analyses of the mafic rocks suggest that greenschist-facies rocks from the Chicken metamorphic complex, a mafic metavolcanic rock from the Nasina assemblage, and an amphibolite from the Fortymile River assemblage formed as island-arc tholeiite in a back-arc setting; another Nasina assemblage greenschist has MORB geochemical characteristics, and another mafic metaigneous rock from the Fortymile River assemblage has geochemical characteristics of calc-alkaline basalt. Our geochemical results imply derivation in an arc and back-arc spreading region within the allochthonous YTT crustal fragment, as previously proposed for correlative units in other parts of the terrane. We also describe the petrography and geochemistry of a newly discovered tectonic lens of Alpine-type metaharzburgite. The metaharzburgite is interpreted to be a sliver of lithospheric mantle from beneath the Seventymile ocean basin or from sub-continental mantle lithosphere of the allochthonous YTT or the western margin of Laurentia that was tectonically emplaced within crustal rocks during closure of the Seventymile ocean basin and subsequently displaced and fragmented by faults.

Alaska, British Columbia, Yukon↗

USGS2017

Surface morphology of caldera-forming eruption deposits revealed by lidar mapping of Crater Lake National Park, Oregon- Implications for emplacement and surface modification

Large explosive eruptions of silicic magma can produce widespread pumice fall, extensive ignimbrite sheets, and collapse calderas. The surfaces of voluminous ignimbrites are rarely preserved or documented because most terrestrial examples are heavily vegetated, or severely modified by post-depositional processes. Much research addresses the internal sedimentary characteristics, flow processes, and depositional mechanisms of ignimbrites, however, surface features of ignimbrites are less well documented and understood, except for comparatively small-volume deposits of historical eruptions. The ~7,700 calendar year B.P. climactic eruption of Mount Mazama, USA vented ~50 km3 of magma, deposited first as rhyodacite pumice fall and then as a zoned rhyodacite-to-andesite ignimbrite as Crater Lake caldera collapsed. Lidar collected during summer 2010 reveals the remarkably well-preserved surface of the Mazama ignimbrite and related deposits surrounding Crater Lake caldera in unprecedented detail despite forest cover. The ±1 m lateral and ±4 cm vertical resolution lidar allows surface morphologies to be classified. Surface morphologies are created by internal depositional processes and can point to the processes at work when pyroclastic flows come to rest. We describe nine surface features including furrow-ridge sets and wedge-shaped mounds in pumice fall eroded by high-energy pyroclastic surges, flow- parallel ridges that record the passage of multiple pyroclastic flows, perched benches of marginal deposits stranded by more-mobile pyroclastic-flow cores, hummocks of dense clasts interpreted as lag deposit, transverse ridges that mark the compression and imbrication of flows as they came to rest, scarps indicating ignimbrite remobilization, fields of pit craters caused by phreatic explosions, fractures and cracks caused by extensional processes resulting from ignimbrite volume loss, and stream channels eroded in the newly formed surface. The nine morphologies presented here illustrate a dynamic depositional environment that varied spatially and with time during the eruption, and show that multiple processes modified the ignimbrite after deposition, both during and after the eruption.

Oregon↗

USGS2008

Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States

Spatial climate data sets of 1971–2000 mean monthly precipitation and minimum and maximum temperature were developed for the conterminous United States. These 30‐arcsec (∼800‐m) grids are the official spatial climate data sets of the U.S. Department of Agriculture. The PRISM (Parameter‐elevation Relationships on Independent Slopes Model) interpolation method was used to develop data sets that reflected, as closely as possible, the current state of knowledge of spatial climate patterns in the United States. PRISM calculates a climate–elevation regression for each digital elevation model (DEM) grid cell, and stations entering the regression are assigned weights based primarily on the physiographic similarity of the station to the grid cell. Factors considered are location, elevation, coastal proximity, topographic facet orientation, vertical atmospheric layer, topographic position, and orographic effectiveness of the terrain. Surface stations used in the analysis numbered nearly 13 000 for precipitation and 10 000 for temperature. Station data were spatially quality controlled, and short‐period‐of‐record averages adjusted to better reflect the 1971–2000 period. PRISM interpolation uncertainties were estimated with cross‐validation (C‐V) mean absolute error (MAE) and the 70% prediction interval of the climate–elevation regression function. The two measures were not well correlated at the point level, but were similar when averaged over large regions. The PRISM data set was compared with the WorldClim and Daymet spatial climate data sets. The comparison demonstrated that using a relatively dense station data set and the physiographically sensitive PRISM interpolation process resulted in substantially improved climate grids over those of WorldClim and Daymet. The improvement varied, however, depending on the complexity of the region. Mountainous and coastal areas of the western United States, characterized by sparse data coverage, large elevation gradients, rain shadows, inversions, cold air drainage, and coastal effects, showed the greatest improvement. The PRISM data set benefited from a peer review procedure that incorporated local knowledge and data into the development process.

International Journal of Climatology↗

USGS1996

Detection and monitoring of H2O and CO2 ice clouds on Mars

We have developed an observational scheme for the detection and discrimination of Mars atmospheric H 2 O and CO 2 clouds using ground-based instruments in the near infrared. We report the results of our cloud detection and characterization study using Mars near IR images obtained during the 1990 and 1993 oppositions. We focused on specific wavelengths that have the potential, based on previous laboratory studies of H 2 O and CO 2 ices, of yielding the greatest degree of cloud detectability and compositional discriminability. We have detected and mapped absorption features at some of these wavelengths in both the northern and southern polar regions of Mars. Compositional information on the nature of these absorption features was derived from comparisons with laboratory ice spectra and with a simplified radiative transfer model of a CO 2 ice cloud overlying a bright surface. Our results indicate that both H 2 O and CO 2 ices can be detected and distinguished in the polar hood clouds. The region near 3.00 μm is most useful for the detection of water ice clouds because there is a strong H 2 O ice absorption at this wavelength but only a weak CO 2 ice band. The region near 3.33 μm is most useful for the detection of CO 2 ice clouds because there is a strong, relatively narrow CO 2 ice band at this wavelength but only broad “continuum” H 2 O ice absorption. Weaker features near 2.30 μm could arise from CO 2 ice at coarse grain sizes, or surface/dust minerals. Narrow features near 2.00 μm, which could potentially be very diagnostic of CO 2 ice clouds, suffer from contamination by Mars atmospheric CO 2 absorptions and are difficult to interpret because of the rather poor knowledge of surface elevation at high latitudes. These results indicate that future ground-based, Earth-orbital, and spacecraft studies over a more extended span of the seasonal cycle should yield substantial information on the style and timing of volatile transport on Mars, as well as a more detailed understanding of the role of CO 2 condensation in the polar heat budget.

Journal of Geophysical Research E: Planets↗