SEARCH · USGS Science

Results for “One Earth”

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

USGS2002

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗

USGS2013

Assessment of spectral, misregistration, and spatial uncertainties inherent in the cross-calibration study

Cross-calibration of satellite sensors permits the quantitative comparison of measurements obtained from different Earth Observing (EO) systems. Cross-calibration studies usually use simultaneous or near-simultaneous observations from several spaceborne sensors to develop band-by-band relationships through regression analysis. The investigation described in this paper focuses on evaluation of the uncertainties inherent in the cross-calibration process, including contributions due to different spectral responses, spectral resolution, spectral filter shift, geometric misregistrations, and spatial resolutions. The hyperspectral data from the Environmental Satellite SCanning Imaging Absorption SpectroMeter for Atmospheric CartograpHY and the EO-1 Hyperion, along with the relative spectral responses (RSRs) from the Landsat 7 Enhanced Thematic Mapper (TM) Plus and the Terra Moderate Resolution Imaging Spectroradiometer sensors, were used for the spectral uncertainty study. The data from Landsat 5 TM over five representative land cover types (desert, rangeland, grassland, deciduous forest, and coniferous forest) were used for the geometric misregistrations and spatial-resolution study. The spectral resolution uncertainty was found to be within 0.25%, spectral filter shift within 2.5%, geometric misregistrations within 0.35%, and spatial-resolution effects within 0.1% for the Libya 4 site. The one-sigma uncertainties presented in this paper are uncorrelated, and therefore, the uncertainties can be summed orthogonally. Furthermore, an overall total uncertainty was developed. In general, the results suggested that the spectral uncertainty is more dominant compared to other uncertainties presented in this paper. Therefore, the effect of the sensor RSR differences needs to be quantified and compensated to avoid large uncertainties in cross-calibration results.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2014

Seismicity of the Earth 1900-2013 offshore British Columbia-southeastern Alaska and vicinity

The tectonics of the Pacific margin of North America between Vancouver Island and south-central Alaska are dominated by the northwest motion of the Pacific plate with respect to the North America plate at a velocity of approximately 50 mm/yr. In the south of this mapped region, convergence between the northern extent of the Juan de Fuca plate (also known as the Explorer microplate) and North America plate dominate. North from the Explorer, Pacific, and North America plate triple junction, Pacific:North America motion is accommodated along the ~650-km-long Queen Charlotte fault system. Offshore of Haida Gwaii and to the southwest, the obliquity of the Pacific:North America plate motion vector creates a transpressional regime, and a complex mixture of strike-slip and convergent (underthrusting) tectonics. North of the Haida Gwaii islands, plate motion is roughly parallel to the plate boundary, resulting in almost pure dextral strike-slip motion along the Queen Charlotte fault. To the north, the Queen Charlotte fault splits into multiple structures, continuing offshore of southwestern Alaska as the Fairweather fault, and branching east into the Chatham Strait and Denali faults through the interior of Alaska. The plate boundary north and west of the Fairweather fault ultimately continues as the Alaska-Aleutians subduction zone, where Pacific plate lithosphere subducts beneath the North America plate at the Aleutians Trench. The transition is complex, and involves intraplate structures such as the Transition fault. The Pacific margin offshore British Columbia is one of the most active seismic zones in North America and has hosted a number of large earthquakes historically.

Alaska, British Columbia↗

USGS2007

Integrating laboratory creep compaction data with numerical fault models: A Bayesian framework

[1] We developed a robust Bayesian inversion scheme to plan and analyze laboratory creep compaction experiments. We chose a simple creep law that features the main parameters of interest when trying to identify rate-controlling mechanisms from experimental data. By integrating the chosen creep law or an approximation thereof, one can use all the data, either simultaneously or in overlapping subsets, thus making more complete use of the experiment data and propagating statistical variations in the data through to the final rate constants. Despite the nonlinearity of the problem, with this technique one can retrieve accurate estimates of both the stress exponent and the activation energy, even when the porosity time series data are noisy. Whereas adding observation points and/or experiments reduces the uncertainty on all parameters, enlarging the range of temperature or effective stress significantly reduces the covariance between stress exponent and activation energy. We apply this methodology to hydrothermal creep compaction data on quartz to obtain a quantitative, semiempirical law for fault zone compaction in the interseismic period. Incorporating this law into a simple direct rupture model, we find marginal distributions of the time to failure that are robust with respect to errors in the initial fault zone porosity.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Seismic wave triggering of nonvolcanic tremor, episodic tremor and slip, and earthquakes on Vancouver Island

[1] We explore the physical conditions that enable triggering of nonvolcanic tremor and earthquakes by considering local seismic activity on Vancouver Island, British Columbia during and immediately after the arrival of large-amplitude seismic waves from 30 teleseismic and 17 regional or local earthquakes. We identify tremor triggered by four of the teleseismic earthquakes. The close temporal and spatial proximity of triggered tremor to ambient tremor and aseismic slip indicates that when a fault is close to or undergoing failure, it is particularly susceptible to triggering of further events. The amplitude of the triggering waves also influences the likelihood of triggering both tremor and earthquakes such that large amplitude waves triggered tremor in the absence of detectable aseismic slip or ambient tremor. Tremor and energy radiated from regional/local earthquakes share the same frequency passband so that tremor cannot be identified during these smaller, more frequent events. We confidently identify triggered local earthquakes following only one teleseism, that with the largest amplitude, and four regional or local events that generated vigorous aftershock sequences in their immediate vicinity. Earthquakes tend to be triggered in regions different from tremor and with high ambient seismicity rates. We also note an interesting possible correlation between large teleseismic events and episodic tremor and slip (ETS) episodes, whereby ETS events that are “late” and have built up more stress than normal are susceptible to triggering by the slight nudge of the shaking from a large, distant event, while ETS events that are “early” or “on time” are not.

Vancouver Island↗

USGS1998

Cretaceous plutonic rocks in the Donner Lake-Cisco Grove area, northern Sierra Nevada, California

The northernmost occurrences of extensive, glaciated exposures of the Sierra Nevada batholith occur in the Donner Lake-Cisco Grove area of the northern Sierra Nevada. The plutonic rocks in this area, which are termed here the Castle Valley plutonic assemblage, crop out over an area of 225 km2 and for the most part are shown as a single undifferentiated mass on previously published geological maps. In the present work, the plutonic assemblage is divided into eight separate intrusive units or lithodemes, two of which each consist of two separate plutons. Compositions are dominantly granodiorite and tonalite, but diorite and granite form small plutons in places. Spectacular examples of comb layering and orbicular texture occur in the diorites. U-Pb zircon ages have been obtained for all but one of the main units and range from ~120 to 114 Ma, indicating that the entire assemblage was emplaced in a narrow time frame in the Early Cretaceous. This is consistent with abundant field evidence that many of the individual phases were intruded penecontemporaneously. The timing of emplacement correlates with onset of major Cretaceous plutonism in the main part of the Sierra Nevada batholith farther south. The emplacement ages also are similar to isotopic ages for gold-quartz mineralization in the Sierran foothills west of the study area, suggesting a direct genetic relationship between the voluminous Early Cretaceous plutonism and hydrothermal gold mineralization.

The Compass: Earth Science Journal of Sigma Gamma ↗

USGS2015

Gravity data from the Sierra Vista Subwatershed, Upper San Pedro Basin, Arizona

Observations of very small changes of Earth’s gravitational field (time-lapse gravity) provide a direct, non-invasive method for measuring changes in aquifer storage change. An existing network of gravity stations in the Sierra Vista Subwatershed was revised in 2014 to better understand the spatial distribution of changes in aquifer storage, especially with relation to ephemeral channel recharge and a groundwater cone of depression associated with pumping in the greater Sierra Vista area. In addition, the network was extended to provide baseline data for possible future enhanced-recharge projects. This report (1) summarizes changes to the Sierra Vista Subwatershed regional time-lapse gravity network with respect to station locations and (2) presents 2014 and 2015 gravity measurements and gravity values at each station. A prior gravity network, established between 2000 and 2005, was revised in 2014 to cover a larger number of stations over a smaller geographic area in order to decrease measurement and interpolation uncertainty. The network currently consists of 59 gravity stations, including 14 absolute-gravity stations. Following above-average rainfall during summer 2014, gravity increased at all but one of the absolute-gravity stations that were observed in both June 2014 and January 2015. This increase in gravity indicates increased groundwater storage in the aquifer and (or) unsaturated zone as a result of rainfall and infiltration.

Arizona↗

USGS2020

Measuring basal force fluctuations of debris flows using seismic recordings and empirical green's functions

We present a novel method for measuring the fluctuating basal normal and shear stresses of debris flows by using along‐channel seismic recordings. Our method couples a simple parameterization of a debris flow as a seismic source with direct measurements of seismic path effects using empirical Green's functions generated with a force hammer. We test this method using two large‐scale (8 and 10 m 3 ) experimental flows at the U.S. Geological Survey debris‐flow flume that were recorded by dozens of three‐component seismic sensors. The seismically derived basal stress fluctuations compare well in amplitude and timing to independent force plate measurements within the valid frequency range (15–50 Hz). We show that although the high‐frequency seismic signals provide band‐limited forcing information, there are systematic relations between the fluctuating stresses and independently measured flow properties, especially mean basal shear stress and flow thickness. However, none of the relationships are simple, and since the flow properties also correlate with one another, we cannot isolate a single factor that relates in a simple way to the fluctuating forces. Nevertheless, our observations, most notably the gradually declining ratio of fluctuating to mean basal stresses during flow passage and the distinctive behavior of the coarse, unsaturated flow front, imply that flow style may be a primary control on the conversion of translational to vibrational kinetic energy. This conversion ultimately controls the radiation of high‐frequency seismic waves. Thus, flow style may provide the key to revealing the nature of the relationship between fluctuating forces and other flow properties.

Journal of Geophysical Research Earth Surface↗

USGS2012

The 1979 Homestead Valley earthquake sequence, California: Control of aftershocks and postseismic deformation

The coseismic slip and geometry of the March 15, 1979, Homestead Valley, California, earthquake sequence are well constrained by precise horizontal and vertical geodetic observations and by data from a dense local seismic network. These observations indicate 0.52±0.10 m of right-lateral slip and 0.17±0.04 m of reverse slip on a buried vertical 6-km-long and 5-km-deep fault and yield a mean static stress drop of 7.2±1.3 MPa. The largest shock had M S = 5.6. Observations of the ground rupture revealed up to 0.1 m of right-lateral slip on two mapped faults that are subparallel to the modeled seismic slip plane. In the 1.9 years since the earthquakes, geodetic network displacements indicate that an additional 60±10 mm of postseismic creep took place. The rate of postseismic shear strain (0.53±0.13 μrad/yr) measured within a 30×30-km network centered on the principal events was anomalously high compared to its preearthquake value and the postseismic rate in the adjacent network. This transient cannot be explained by postseismic slip on the seismic fault but rather indicates that broadscale release of strain followed the earthquake sequence. We have calculated the postearthquake stress field caused by the modeled coseismic slip. We assume that failure is promoted when the sum of the shear stress plus 0.75 times the faultopening stress increases. Most aftershocks concentrate at points where the stresses are enhanced by 0.3 MPa (3 bars) or more; aftershocks are nearly absent where postearthquake stresses decrease by 0.3–0.5 MPa. Isolated off-fault clusters of aftershocks that locate at one fault length from the rupture plane are explainable by this hypothesis. We find that ground rupture and postseismic creep take place where near-surface stresses are calculated to increase within the preexisting fault zones. Two patches that extend 4 km from both ends of the seismic fault exhibited neither aftershocks nor measurable postseismic creep. The sensitivity of aftershocks and ground rupture to changes in stress that are less than 5% of the earthquake stress drop demonstrates that the region around the earthquakes was within a few percent of its failure threshold before the main shocks. The preearthquake stress field and the stress required for failure must also have been nearly uniform.

Journal of Geophysical Research Solid Earth↗

USGS1969

(RLC-14) Geologic map of the Alphonsus region of the moon

This 1:250,000- scale geologic map is one of a series prepared largely from photographs transmitted by Ranger IX (reproduced in a report by Jet Propulsion Lab., 1966). It depicts the geology of the crater Alphonsus (in which Ranger IX impacted) and environs. The primary objective of the mapping was to apply extant lunar mapping techniques used on relatively small scale telescopic photographs to the larger scale Ranger photographs, in preparation for extensive analysis of Lunar Orbiter photography in support of the Apollo program. An additional objective was to shed more light on the formation of Alphonsus and its associated features. The map outlines rock units that are inferred from surface features and characteristics, such as topography and albedo. Each rock unit is assigned an age and arranged in the explanation so that the youngest unit are at the top and the oldest at the bottom. The age classification is based on the early work of Shoemaker and Hackman (1962), and subsequent revisions reported by McCauley (1967) and Wilhelms (1966). Geology of the region has been mapped at a scale of 1:1,000,000 on the basis of Earth-based photographs and observations (Howard and Masursky , 1968). A more detailed, 1;50,000-scale map of part of the Alphonsus floor was prepared by McCauley (1969), who used Ranger IX photographs as his source of geologic information.

IMAP↗

USGS1982

Arizona Vegetation Resource Inventory (AVRI) accuracy assessment

A quantitative accuracy assessment was performed for the vegetation classification map produced as part of the Arizona Vegetation Resource Inventory (AVRI) project. This project was a cooperative effort between the Bureau of Land Management (BLM) and the Earth Resources Observation Systems (EROS) Data Center. The objective of the accuracy assessment was to estimate (with a precision of ?10 percent at the 90 percent confidence level) the comission error in each of the eight level II hierarchical vegetation cover types. A stratified two-phase (double) cluster sample was used. Phase I consisted of 160 photointerpreted plots representing clusters of Landsat pixels, and phase II consisted of ground data collection at 80 of the phase I cluster sites. Ground data were used to refine the phase I error estimates by means of a linear regression model. The classified image was stratified by assigning each 15-pixel cluster to the stratum corresponding to the dominant cover type within each cluster. This method is known as stratified plurality sampling. Overall error was estimated to be 36 percent with a standard error of 2 percent. Estimated error for individual vegetation classes ranged from a low of 10 percent ?6 percent for evergreen woodland to 81 percent ?7 percent for cropland and pasture. Total cost of the accuracy assessment was $106,950 for the one-million-hectare study area. The combination of the stratified plurality sampling (SPS) method of sample allocation with double sampling provided the desired estimates within the required precision levels. The overall accuracy results confirmed that highly accurate digital classification of vegetation is difficult to perform in semiarid environments, due largely to the sparse vegetation cover. Nevertheless, these techniques show promise for providing more accurate information than is presently available for many BLM-administered lands.

Open-File Report↗

USGS1953

Radioactive source materials in Los Estados Unidos de Venezuela

This report summarizes the data available on radioactive source materials in Los Estados Unidos de Venezuela accumulated by geologists of the Direccions Tecnica de Geolgia and antecedent agencies prior to June 1951, and the writers from June to November 1951. The investigation comprised preliminary study, field examination, office studies, and the preparation of this report, in which the areas and localities examined are described in detail, the uranium potentialities of Venezuela are summarized, and recommendations are made. Preliminary study was made to select areas and rock types that were known or reported to be radioactive or that geologic experience suggests would be favorable host for uranium deposits, In the office, a study of gamma-ray well logs was started as one means of amassing general radiometric data and of rapidly scanning many of ye rocks in northern Venezuela; gamma-ray logs from about 140 representative wells were examined and their peaks of gamma intensity evaluated; in addition samples were analyzed radiometrically, and petrographically. Radiometic reconnaissance was made in the field during about 3 months of 1951, or about 12 areas, including over 100 localities in the State of Miranda, Carabobo, Yaracuy, Falcon, Lara, Trujillo, Zulia, Merida, Tachira, Bolivar, and Territory Delta Amacuro. During the course of the investigation, both in the filed and office, information was given about geology of uranium deposits, and in techniques used in prospecting and analysis. All studies and this report are designed to supplement and to strengthen the Direccion Tecnica de Geologias's program of investigation of radioactive source in Venezuela now in progress. The uranium potentialities of Los Estados de Venezuela are excellent for large, low-grade deposits of uraniferous phospahtic shales containing from 0.002 to 0.027 percent uranium; fair, for small or moderate-sized, low-grade placer deposits of thorium, rare-earth, and uranium minerals; poor, for high-grade hydrothermal pitchblende deposits; and highly possible for small, medium- to high-grade despots of carnotite-or copper-uranium bearing sandstone. Recommendations for the Venezuelan uranium program include 1) the systematic collection of a mass general radiometric data by examining sample collections, expanding the gamma-ray program, encouraging the use of Geiger counter by field geologists, and by enlisting the aid of the general public; 2) , the examination of specific areas or localities, chosen on the basis of geologic favorability from the results of the amassing of data, or obtained by hints and rumors; 3), the organization of a unit within the Direccion Tecnica de Geologica to direct, collection, and collate metric data. It is emphasized that to be most fruitful the program requires the application of sounds and imaginative geologic theory.

Trace Elements Investigations↗

USGS2012

Soil-water dynamics and unsaturated storage during snowmelt following wildfire

Many forested watersheds with a substantial fraction of precipitation delivered as snow have the potential for landscape disturbance by wildfire. Little is known about the immediate effects of wildfire on snowmelt and near-surface hydrologic responses, including soil-water storage. Montane systems at the rain-snow transition have soil-water dynamics that are further complicated during the snowmelt period by strong aspect controls on snowmelt and soil thawing. Here we present data from field measurements of snow hydrology and subsurface hydrologic and temperature responses during the first winter and spring after the September 2010 Fourmile Canyon Fire in Colorado, USA. Our observations of soil-water content and soil temperature show sharp contrasts in hydrologic and thermal conditions between north- and south-facing slopes. South-facing burned soils were ∼1–2 °C warmer on average than north-facing burned soils and ∼1.5 °C warmer than south-facing unburned soils, which affected soil thawing during the snowmelt period. Soil-water dynamics also differed by aspect: in response to soil thawing, soil-water content increased approximately one month earlier on south-facing burned slopes than on north-facing burned slopes. While aspect and wildfire affect soil-water dynamics during snowmelt, soil-water storage at the end of the snowmelt period reached the value at field capacity for each plot, suggesting that post-snowmelt unsaturated storage was not substantially influenced by aspect in wildfire-affected areas. Our data and analysis indicate that the amount of snowmelt-driven groundwater recharge may be larger in wildfire-impacted areas, especially on south-facing slopes, because of earlier soil thaw and longer durations of soil-water contents above field capacity in those areas.

Hydrology and Earth System Sciences↗

USGS2024

The National Ocean Biodiversity Strategy

President Biden has been clear that the ocean is central to life on Earth. As he has proclaimed, “the ocean powers millions of jobs; feeds and sustains us; and is a rejuvenating source of inspiration, exploration, and recreation.” The Biden-Harris Administration has worked hard to fulfill the President’s goal to protect and conserve at least 30% of U.S. waters by 2030. The ocean faces increased threats from warming, overfishing, increased acidity, and loss of biodiversity. It is now more important than ever to sustain the many benefits that the ocean, coasts, and Great Lakes provide, including food, a favorable climate, recreation, physical and mental health, and for many, a sense of cultural identity. Ocean life represents an irreplaceable heritage, the foundation of a habitable planet, and a vast trove of resources. Keeping our ocean healthy requires reliable information on the changing status of these living organisms, the drivers of biodiversity change, and options for effectively addressing those drivers. Over 2 million species are estimated to live in the ocean, yet only about 240,000 species have been described by scientists. Most of the ocean’s benefits result from those diverse species interacting with one another and the environment they create. To protect and conserve the ocean, we as a nation need to make better use of existing knowledge and prioritize acquiring new biodiversity knowledge to enable better policy and management decisions. The ability to monitor ocean species and habitats has expanded dramatically over the past decade, with innovations in technology, genomics, taxonomy, big data management and sharing, artificial intelligence, and machine learning. Yet large fractions of the U.S. ocean remain almost unknown. The National Ocean Biodiversity Strategy (strategy) reflects the urgent need to leverage these advances. The goals of this strategy must be guided by the nation’s diverse voices and ways of knowing, in order to maximize effective and equitable stewardship of the ocean’s diverse life and its benefits to people. The strategy is intended as a guiding document for government to advance three overarching goals: ● Goal 1: Drive delivery of ocean biodiversity knowledge at the national scale. Objectives include developing an Implementation Plan for achieving the strategy’s three goals; establishing a coordination mechanism to manage the implementation; and documenting gaps in biodiversity knowledge and the benefits of ocean biodiversity to people and economies. ● Goal 2: Strengthen tools and institutions to deliver ocean biodiversity knowledge. Objectives include establishing a robust information pipeline to support indicators and dynamic maps of ocean biodiversity, from the coasts to the deep sea. This pipeline should include expanded observing systems and comprehensive data management; science and technology solutions to accelerate the availability of biodiversity information; and plans to leverage previous investments to rebuild and expand the nation’s human capital and infrastructure to sustain foundational taxonomy and biodiversity science. ● Goal 3: Protect, conserve, restore, and sustainably use ocean biodiversity. Objectives include expanding the collection, delivery, and use of biodiversity knowledge to inform actions that advance ocean protection, conservation, restoration, and sustainable development. Government should lead in establishing and incentivizing diverse partnerships across scales and sectors to implement those actions and should educate and involve the public to discover and value the nation’s diverse ocean life. Achieving these goals will require commitments across society: new federal and private investments, coordination across sectors to address climate and equity challenges, and engagement of Indigenous Knowledge holders and frontline communities as full partners throughout planning and implementation. The Subcommittee on Ocean Science and Technology (SOST) IWG-Biodiversity will begin developing an Implementation Plan to describe and direct specific actions to implement the strategy. Successful implementation of the strategy will harmonize and expand collection and delivery of timely knowledge on ocean life to all of society. The strategy will also enable evidence-based management and protection of the ocean. Advancing the strategy will build human and institutional capital and partnerships that support both existing mandates and new needs to rebuild and sustain biodiversity, achieve healthy ocean ecosystems, and manage living resources. Implementing the strategy will deliver knowledge for monitoring, modeling, forecasting, and assessments that support food security, public health, and cultural values, and that more effectively protect, conserve, and restore nature.

Report↗

USGS2024

Ecological dissimilarity matters more than geographical distance when predicting land surface indicators using machine learning

Supervised training techniques, such as those used in machine learning, use generally large sets of in situ data to train models that can, in turn, be used to make predictions (or prediction maps) about the Earth’s surface in times or places where no in situ data exist. The purpose of the present study is to investigate, using a very large set of in situ data from across the western United States (U.S.), the conditions under which training data from a different geographic region where predictions are desired may be substituted. To do this, we train models using in situ data from level IV ecoregions and test how well these models predict surface conditions in different ecoregions. We characterize the difference between the possible pairs of ecoregion in terms of geographical (centroid-to-centroid) distance and “ecological dissimilarity.” Ecological dissimilarity between pairs of ecoregions is defined in two ways: 1) as the Euclidean distance in multivariate space defined by in situ indicators designed for monitoring purposes and 2) in terms of the difference in temporal behavior from model- and remote sensing-derived datasets. Although, overall, prediction error increases with geographical distance between training and testing ecoregions, our results indicate that ecological dissimilarity can be used to predict the error expected from a model trained with data from one ecoregion when applied in a different ecoregion.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS1971

Aeromagnetic study of the midcontinent gravity high of central United States

A composite map of detailed aeromagnetic surveys over the midcontinent gravity high provides coverage of the 600-mi-long buried belt of mafic rocks of the Keweenawan Series from their outcrop localities in Minnesota and Wisconsin through Iowa and Nebraska. A map of the subsurface extent of the mafic rocks, based on the intricate magnetic patterns, shows that the rocks form a long, semicontinuous block, averaging 40 mi wide and consisting mainly of a sequence of layered flows. This sequence is probably fault-bounded and has been tilted up along the margins, where the linearity of the anomalies indicates steeper dips. The associated clastic rocks, indicated by a smoother magnetic pattern, occur in basins along both sides of the mafic belt and in grabens and a series of axial basins on the upper surface of the block. The well-defined outliers of flows marginal to the main block and the truncation of some of the outermost flow units along a diagonal boundary striking at an angle to them suggest that the present boundaries of the block are postdepositional structural features. The basins and the edges of the block appear to have controlled later, largely vertical movement in the overlying Paleozoic and younger sedimentary cover. Calculated models based on coincident magnetic and detailed gravity profiles along typical cross sections of the midcontinent gravity high show that the block of mafic rocks is steep-sided and as much as several miles thick. The free-air gravity anomaly, which consists of a large positive maximum flanked by minima, averages very close to zero, indicating that this major crustal feature is regionally compensated, although locally each of its components shows a large departure from equilibrium. Remanent magnetization is a primary factor in the interpretation of the magnetic data. Magnetic property studies of Keweenawan mafic rocks in the Lake Superior region show that remanent magnetization may be five times the magnetization induced by the present Earth's field and differs from it radically in direction. This magnetization was acquired before the flows were tilted into their present positions. A computed magnetic profile shows that a trough of flows with such a magnetization and inward-dipping limbs can account for the observed persistent lows along the western edge of the block, the relatively low magnetic values along the axis of the block, and the large positive anomaly along the eastern side of the block. Flows as much as 1 mi thick near the base of the sequence have a remanent magnetization with a nearly opposite polarity. This reverse polarity has been measured on both sides of Lake Superior and is probably also present farther south, particularly in Iowa where the outer units of the block in an area north of Des Moines give rise to a prominent magnetic low. The axis of this long belt of Keweenawan mafic rocks cuts discordantly through the prevailing east-west-trending fabric of the older Precambrian terrane from southern Kansas to Lake Superior. This belt has several major left-lateral offsets, one of which produces a complete hiatus in the vicinity of the 40th parallel where an east-west transcontinental rift or fracture zone has been proposed. The axial basins of clastic rocks are outlined by linear magnetic anomalies and show a concordant relation to the structure of the mafic flows. These basins are oriented at an angle to the main axis, suggesting that the entire feature originated as a major rift composed of a series of short, linear, en echelon segments with offsets similar to the transform faults characterizing the present mid-ocean rift system. This midcontinent rift may well have been part of a Keweenawan global rift system with initial offsets consisting of transform faults along pre-existing fractures, but apparently it never fully developed laterally into an ocean basin, and the upwelling mafic material was localized along a relatively narrow belt.

Geological Society of America Bulletin↗

USGS1990

A numerical study of some potential sources of error in side-by-side seismometer evaluations

This report presents the results of a series of computer simulations of potential errors in test data, which might be obtained when conducting side-by-side comparisons of seismometers. These results can be used as guides in estimating potential sources and magnitudes of errors one might expect when analyzing real test data. First, the derivation of a direct method for calculating the noise levels of two sensors in a side-by-side evaluation is repeated and extended slightly herein. This bulk of this derivation was presented previously (see Holcomb 1989); it is repeated here for easy reference. This method is applied to the analysis of a simulated test of two sensors in a side-by-side test in which the outputs of both sensors consist of white noise spectra with known signal-tonoise ratios (SNR's). This report extends this analysis to high SNR's to determine the limitations of the direct method for calculating the noise levels at signal-to-noise levels which are much higher than presented previously (see Holcomb 1989). Next, the method is used to analyze a simulated test of two sensors in a side-by-side test in which the outputs of both sensors consist of bandshaped noise spectra with known signal-tonoise ratios. This is a much more realistic representation of real world data because the earth's background spectrum is certainly not flat. Finally, the results of the analysis of simulated white and bandshaped side-by-side test data are used to assist in interpreting the analysis of the effects of simulated azimuthal misalignment in side-by-side sensor evaluations. A thorough understanding of azimuthal misalignment errors is important because of the physical impossibility of perfectly aligning two sensors in a real world situation. The analysis herein indicates that alignment errors place lower limits on the levels of system noise which can be resolved in a side-by-side measurement It also indicates that alignment errors are the source of the fact that real data noise spectra tend to follow the earth's background spectra in shape.

Open-File Report↗

USGS1971

40Ar/39Ar technique of KAr dating: A comparison with the conventional technique

K-Ar ages have been determined by the 40 Ar/ 39 Ar total fusion technique on 19 terrestrial samples whose conventional K-Ar ages range from 3.4 my to nearly 1700 my. Sample materials included biotite, muscovite, sanidine, adularia, plagioclase, hornblende, actinolite, alunite, dacite, and basalt. For 18 samples there are no significant differences at the 95% confidence level between the K Ar ages obtained by these two techniques; for one sample the difference is 4.3% and is statistically significant. For the neutron doses used in these experiments (≈4 × 10 18 nvt) it appears that corrections for interfering Ca- and K-derived Ar isotopes can be made without significant loss of precision for samples with K/Ca > 1 as young as about 5 × 10 5 yr, and for samples with K/Ca < 1 as young as about 10 7 yr. For younger samples the combination of large atmospheric Ar corrections and large corrections for Ca- and K-derived Ar may make the precision of the 40 Ar/ 39 Ar technique less than that of the conventional technique unless the irradiation parameters are adjusted to minimize these corrections.

Earth and Planetary Science Letters↗