SEARCH · Search USGS

Results for “Developments in Earth Surface Processes”

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.

269 records · Page 15Linked to original sources

USGS2009

EAARL Coastal Topography-Assateague Island National Seashore, 2008: First Surface

These remotely sensed, geographically referenced elevation measurements of lidar-derived first-surface (FS) topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL; the National Park Service (NPS), Northeast Coastal and Barrier Network, Kingston, RI; and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. This project provides highly detailed and accurate datasets of the Assateague Island National Seashore in Maryland and Virginia, acquired March 24-25, 2008. The datasets are made available for use as a management tool to research scientists and natural-resource managers. An innovative airborne lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multi-spectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of lidar data in an interactive or batch mode. Modules for pre-survey flight-line definition, flight-path plotting, lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is used routinely to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of last return elevations.

Data Series↗

USGS2010

EAARL coastal topography-Cape Hatteras National Seashore, North Carolina, post-Nor'Ida, 2009: first surface

These remotely sensed, geographically referenced elevation measurements of lidar-derived first-surface (FS) topography datasets were produced collaboratively by the U.S. Geological Survey (USGS), St. Petersburg Coastal and Marine Science Center, St. Petersburg, FL, and the National Park Service (NPS), Northeast Coastal and Barrier Network, Kingston, RI. This project provides highly detailed and accurate datasets of a portion of the National Park Service Southeast Coast Network's Cape Hatteras National Seashore in North Carolina, acquired post-Nor'Ida (November 2009 nor'easter) on November 27 and 29 and December 1, 2009. The datasets are made available for use as a management tool to research scientists and natural-resource managers. An innovative airborne lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multispectral color-infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine aircraft, but the instrument was deployed on a Pilatus PC-6. A single pilot, a lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of lidar data in an interactive or batch mode. Modules for presurvey flight-line definition, flight-path plotting, lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is used routinely to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of last return elevations. For more information about similar projects, please visit the Decision Support for Coastal Science and Management website.

Data Series↗

USGS2011

Bedform response to flow variability

Laboratory observations and computational results for the response of bedform fields to rapid variations in discharge are compared and discussed. The simple case considered here begins with a relatively low discharge over a flat bed on which bedforms are initiated, followed by a short high-flow period with double the original discharge, during which the morphology of the bedforms adjusts, followed in turn by a relatively long period of the original low discharge. For the grain size and hydraulic conditions selected, the Froude number remains subcritical during the experiment, and sediment moves predominantly as bedload. Observations show rapid development of quasi-two-dimensional bedforms during the initial period of low flow with increasing wavelength and height over the initial low-flow period. When the flow increases, the bedforms rapidly increase in wavelength and height, as expected from other empirical results. When the flow decreases back to the original discharge, the height of the bedforms quickly decreases in response, but the wavelength decreases much more slowly. Computational results using an unsteady two-dimensional flow model coupled to a disequilibrium bedload transport model for the same conditions simulate the formation and initial growth of the bedforms fairly accurately and also predict an increase in dimensions during the high-flow period. However, the computational model predicts a much slower rate of wavelength increase, and also performs less accurately during the final low-flow period, where the wavelength remains essentially constant, rather than decreasing. In addition, the numerical results show less variability in bedform wavelength and height than the measured values; the bedform shape is also somewhat different. Based on observations, these discrepancies may result from the simplified model for sediment particle step lengths used in the computational approach. Experiments show that the particle step length varies spatially and temporally over the bedforms during the evolution process. Assuming a constant value for the step length neglects the role of flow alterations in the bedload sediment-transport process, which appears to result in predicted bedform wavelength changes smaller than those observed. However, observations also suggest that three-dimensional effects play at least some role in the decrease of bedform wavelength, so incorporating better models for particle hop lengths alone may not be sufficient to improve model predictions.

Earth Surface Processes and Landforms↗

USGS2018

Dust on a Hawaiian volcano: A regional model using field measurements to estimate transport and deposition

The western slopes of Hawaii's Mauna Kea volcano are mantled by fine-grained soils, the record of volcanic airfall and eolian deposition. Where exposed, strong winds transport this sediment across West Hawaii, affecting tourism and local communities with decreased air and water quality. Operations on US Army's Ke'amuku Maneuver Area (KMA) have the potential to increase dust flux from these deposits. The USGS established 18 ground monitoring sites and sampling locations surrounding KMA. For over 3 years, each station measured vertical and horizontal dust flux, while co-located anemometers measured wind speed and direction. We used these datasets to develop a parsimonious regional model for dust supply and transport to assess whether KMA is a net dust sink or source. We found that dust transport is most highly correlated with threshold wind speeds of 8 m/s. We used this value as the regional average threshold wind speed for dust entrainment. Using a model that partitions measured horizontal dust flux into inward- and outward-directed components, we estimate that KMA is currently a net dust sink. Geochemical analysis of dust samples illustrates that local organics and carbonate make up 64% of dust mass, the remainder being volcanic silt and fine sand. Measured vertical dust deposition rates of 0.006 mm/yr are similar to 0.004 mm/yr of deposition predicted from taking the divergence of dust across KMA's boundary. These rates are low compared with pre-historic rates of ~0.2–0.3 mm/yr, from radiocarbon dating of buried soils. KMA's soils record persistent deposition over millennia, at rates that imply episodic dust storms. Such events created a soil-mantled landscape in the middle of a largely Pleistocene rocky landscape. A substantial portion of fine-grained soils in other leeward Hawaiian Island landscapes may have formed from similar eolian deposition, and not direct weathering of parent rock. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.

Hawaii↗

USGS2012

Developing a new stream metric for comparing stream function using a bank-floodplain sediment budget: a case study of three Piedmont streams

A bank and floodplain sediment budget was created for three Piedmont streams tributary to the Chesapeake Bay. The watersheds of each stream varied in land use from urban (Difficult Run) to urbanizing (Little Conestoga Creek) to agricultural (Linganore Creek). The purpose of the study was to determine the relation between geomorphic parameters and sediment dynamics and to develop a floodplain trapping metric for comparing streams with variable characteristics. Net site sediment budgets were best explained by gradient at Difficult Run, floodplain width at Little Conestoga Creek, and the relation of channel cross-sectional area to floodplain width at Linganore Creek. A correlation for all streams indicated that net site sediment budget was best explained by relative floodplain width (ratio of channel width to floodplain width). A new geomorphic metric, the floodplain trapping factor, was used to compare sediment budgets between streams with differing suspended sediment yields. Site sediment budgets were normalized by floodplain area and divided by the stream's sediment yield to provide a unitless measure of floodplain sediment trapping. A floodplain trapping factor represents the amount of upland sediment that a particular floodplain site can trap (e.g. a factor of 5 would indicate that a particular floodplain site traps the equivalent of 5 times that area in upland erosional source area). Using this factor we determined that Linganore Creek had the highest gross and net (floodplain deposition minus bank erosion) floodplain trapping factor (107 and 46, respectively) that Difficult Run the lowest gross floodplain trapping factor (29) and Little Conestoga Creek had the lowest net floodplain trapping factor (–14, indicating that study sites were net contributors to the suspended sediment load). The trapping factor is a robust metric for comparing three streams of varied watershed and geomorphic character, it promises to be a useful tool for future stream assessments.

Piedmont;Difficult Run;Little Conestoga Creek;Ling↗

USGS2025

River ice controls permafrost bank erosion across an Arctic delta

Bank erosion in Arctic rivers helps shape channel geometry, mobilizes carbon from permafrost and influences sediment delivery to the Arctic Ocean. On Alaska's Arctic coastal plain, rivers begin flowing during snowmelt in late spring while extensive river ice persists in channels, such that hydraulics are altered and water is kept cool. The effects of river ice on permafrost bank erosion are poorly understood, primarily due to a dearth of field observations and a lack of river ice in existing models. To address this knowledge gap, we developed a numerical model to simulate the melt of substrate interstitial ice and bank collapse along individual permafrost river banks. We parameterize the model with field observations from riverbanks in three different channels on the Canning River delta, which are disparately impacted by river ice during snowmelt. We explore the bank erosion produced without river ice in the model and with modern river ice model scenarios that we drive with different stages and water temperature boundary conditions. We also compare predicted erosion rates to observations from satellite imagery to validate this approach. In the model, banks are idealized as vertical profiles that rise 1–2 m above the river bed and are comprised of silt- to sand-sized sediment with dense roots in the active layer. Underneath, we generalize bank ice content underneath the active layer to represent ice-rich permafrost on the river corridor boundaries. The model predicts that these ice-rich river banks can erode by 2–6 m/yr. Scenarios without ice underpredict erosion in the distributary channels. Scenarios with varying river ice for different deltaic channels produce erosion rates similar to observations. Our results suggest that the prolonged melt of thick river ice in a delta nonlinearly impacts permafrost bank erosion by blocking river discharge to certain branches, heightening stage across the distributary network and locally limiting river water warming. Given expected changes in air temperature and hydrology, future estimates of Arctic river bank erosion could be improved by considering river ice.

Alaska↗

USGS2012

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗

USGS2009

Changes in the chemistry of shallow groundwater related to the 2008 injection of CO2 at the ZERT field site, Bozeman, Montana

Approximately 300 kg/day of food-grade CO 2 was injected through a perforated pipe placed horizontally 2–2.3 m deep during July 9–August 7, 2008 at the MSU-ZERT field test to evaluate atmospheric and near-surface monitoring and detection techniques applicable to the subsurface storage and potential leakage of CO 2 . As part of this multidisciplinary research project, 80 samples of water were collected from 10 shallow monitoring wells (1.5 or 3.0 m deep) installed 1–6 m from the injection pipe, at the southwestern end of the slotted section (zone VI), and from two distant monitoring wells. The samples were collected before, during, and following CO 2 injection. The main objective of study was to investigate changes in the concentrations of major, minor, and trace inorganic and organic compounds during and following CO 2 injection. The ultimate goals were (1) to better understand the potential of groundwater quality impacts related to CO 2 leakage from deep storage operations, (2) to develop geochemical tools that could provide early detection of CO 2 intrusion into underground sources of drinking water (USDW), and (3) to test the predictive capabilities of geochemical codes against field data. Field determinations showed rapid and systematic changes in pH (7.0–5.6), alkalinity (400–1,330 mg/l as HCO 3 ), and electrical conductance (600–1,800 μS/cm) following CO 2 injection in samples collected from the 1.5 m-deep wells. Laboratory results show major increases in the concentrations of Ca (90–240 mg/l), Mg (25–70 mg/l), Fe (5–1,200 ppb), and Mn (5–1,400 ppb) following CO 2 injection. These chemical changes could provide early detection of CO 2 leakage into shallow groundwater from deep storage operations. Dissolution of observed carbonate minerals and desorption-ion exchange resulting from lowered pH values following CO 2 injection are the likely geochemical processes responsible for the observed increases in the concentrations of solutes; concentrations generally decreased temporarily following four significant precipitation events. The DOC values obtained are 5 ± 2 mg/l, and the variations do not correlate with CO 2 injection. CO 2 injection, however, is responsible for detection of BTEX (e.g. benzene, 0–0.8 ppb), mobilization of metals, the lowered pH values, and increases in the concentrations of other solutes in groundwater. The trace metal and BTEX concentrations are all significantly below the maximum contaminant levels (MCLs). Sequential leaching of core samples is being carried out to investigate the source of metals and other solutes.

Montana↗

USGS2014

Evaluation of statistically downscaled GCM output as input for hydrological and stream temperature simulation in the Apalachicola–Chattahoochee–Flint River Basin (1961–99)

The accuracy of statistically downscaled general circulation model (GCM) simulations of daily surface climate for historical conditions (1961–99) and the implications when they are used to drive hydrologic and stream temperature models were assessed for the Apalachicola–Chattahoochee–Flint River basin (ACFB). The ACFB is a 50 000 km 2 basin located in the southeastern United States. Three GCMs were statistically downscaled, using an asynchronous regional regression model (ARRM), to ⅛° grids of daily precipitation and minimum and maximum air temperature. These ARRM-based climate datasets were used as input to the Precipitation-Runoff Modeling System (PRMS), a deterministic, distributed-parameter, physical-process watershed model used to simulate and evaluate the effects of various combinations of climate and land use on watershed response. The ACFB was divided into 258 hydrologic response units (HRUs) in which the components of flow (groundwater, subsurface, and surface) are computed in response to climate, land surface, and subsurface characteristics of the basin. Daily simulations of flow components from PRMS were used with the climate to simulate in-stream water temperatures using the Stream Network Temperature (SNTemp) model, a mechanistic, one-dimensional heat transport model for branched stream networks. The climate, hydrology, and stream temperature for historical conditions were evaluated by comparing model outputs produced from historical climate forcings developed from gridded station data (GSD) versus those produced from the three statistically downscaled GCMs using the ARRM methodology. The PRMS and SNTemp models were forced with the GSD and the outputs produced were treated as “truth.” This allowed for a spatial comparison by HRU of the GSD-based output with ARRM-based output. Distributional similarities between GSD- and ARRM-based model outputs were compared using the two-sample Kolmogorov–Smirnov (KS) test in combination with descriptive metrics such as the mean and variance and an evaluation of rare and sustained events. In general, precipitation and streamflow quantities were negatively biased in the downscaled GCM outputs, and results indicate that the downscaled GCM simulations consistently underestimate the largest precipitation events relative to the GSD. The KS test results indicate that ARRM-based air temperatures are similar to GSD at the daily time step for the majority of the ACFB, with perhaps subweekly averaging for stream temperature. Depending on GCM and spatial location, ARRM-based precipitation and streamflow requires averaging of up to 30 days to become similar to the GSD-based output. Evaluation of the model skill for historical conditions suggests some guidelines for use of future projections; while it seems correct to place greater confidence in evaluation metrics which perform well historically, this does not necessarily mean those metrics will accurately reflect model outputs for future climatic conditions. Results from this study indicate no “best” overall model, but the breadth of analysis can be used to give the product users an indication of the applicability of the results to address their particular problem. Since results for historical conditions indicate that model outputs can have significant biases associated with them, the range in future projections examined in terms of change relative to historical conditions for each individual GCM may be more appropriate.

Florida, Georgia↗

USGS2024

A framework to facilitate development and testing of image-based river velocimetry algorithms

Image-based methods have compelling, demonstrated potential for characterizing flow fields in rivers, but algorithms like particle image velocimetry (PIV) must be further tested and improved to enable more effective use of these techniques. This paper presents a framework designed for this exact purpose: Simulating Hydraulics and Images for Velocimetry Evaluation and Refinement (SHIVER). The approach involves coupling a hydrodynamic model with a synthetic particle generator to advect particles between frames, as dictated by local velocity vectors and thus construct a plausible image sequence specific to the reach of interest. The resulting time series can then be used as input to a velocimetry algorithm to compare image-derived estimates with known (modelled) velocities to perform an exhaustive, spatially distributed accuracy assessment. As an example application of SHIVER, we examined the effects of interrogation area (IA) size, frame rate, flow velocity, and image sequence duration on the performance of a standard PIV algorithm. This analysis indicated that image-derived velocities were generally in close agreement with those from the flow model (root mean square error <10% and mean bias <3%), except when small IAs were coupled with low frame rates. Velocity estimates were most accurate for the lowest modelled discharge ( R 2 =0.97 at baseflow) and became less reliable as the mean flow velocity increased ( R 2 =0.92 for an intermediate discharge and R 2 =0.86 at bankfull). Accuracy was essentially independent of image sequence duration, implying that long occupations might not be necessary. Errors were concentrated along channel margins, where PIV-based velocities tended to be greater than those from the flow model. Small IAs led to underpredictions of velocity, while larger IAs led to overpredictions. SHIVER is highly modular and could be updated to make use of different hydrodynamic models or image simulators. The framework could also facilitate more thorough sensitivity analyses and comparison of various velocimetry algorithms.

Earth Surface Processes and Landforms↗

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗

USGS1995

U.S. Geological Survey Subsidence Interest Group conference, Edwards Air Force Base, Antelope Valley, California, November 18-19, 1992; Abstracts and summary

Land subsidence, the loss of surface elevation as a result of the removal of subsurface support, affects every state in the United States. More than 17,000 mi 2 of land in the United States has been lowered by the various processes that produce land subsidence with annual costs from resulting flooding and structural damage that exceed $125 million. It is estimated that an additional $400 million is spent nationwide in attempts to control subsidence. Common causes of land subsidence include the removal of oil, gas, and water from underground reservoirs; dissolution of limestone aquifers (sinkholes); underground mining activities; drainage of organic soils; and hydrocompaction (the initial wetting of dry soils). Overdrafting of aquifers is the major cause of areally extensive land subsidence, and as ground-water pumping increases, land subsidence also will increase. Land subsidence and its effects on engineering structures have been recognized for centuries, but it was not until this century that the processes that produce land subsidence were identified and understood. In 1928, while working with field data from a test of the Dakota Sandstone aquifer, O.E. Meinzer of the U.S. Geological Survey recognized the compressibility of aquifers. Around the same time, Karl Terzaghi, a soil scientist working at Harvard University, developed the one-dimensional consolidation theory that provided a quantitative means of predicting soil compaction resulting from the drainage of compressible soils. Thus, with the recognition of the compressibility of aquifers (Meinzer), and the development of a quantitative means of predicting soil compaction as a consequence of the reduction of intergranular pore pressure (Terzaghi), the theory of aquifer-system compaction was formed. With the widespread availability of electric power in rural areas, and the advent of the deep turbine pump, ground-water withdrawals increased dramatically throughout the country in the 1940's and 1950's. Along with this unprecedented increase in pumpage, substantial amounts of land subsidence were observed in several areas of the United States, most notably in Arizona, California, and Texas. Beginning in 1955, under the direction of Joseph Poland, the Geological Survey began the "Mechanics of Aquifers Project," which focused largely on the processes that resulted in land subsidence due to the withdrawal of ground water. This research team gained international renown as they advanced the scientific understanding of aquifer mechanics and land-subsidence theory. The results of field studies by members of this research group not only verified the validity of the application of Terzaghi's consolidation theory to compressible aquifers, but they also provided definitions, methods of quantification, and confirmation of the interrelation among hydraulic head declines, aquifer-system compaction, and land subsidence. In addition to conducting pioneering research, this group also formed a "center of expertise," providing a focal point within the Geological Survey for the dissemination of technology and scientific understanding in aquifer mechanics. However, when the "Mechanics of Aquifers Project" was phased out in 1984, the focal point for technology transfer no longer existed. Interest among various state and local agencies in land subsidence has persisted, and the Geological Survey has continued to participate in a broad spectrum of cooperative and Federally funded projects in aquifer mechanics and land subsidence. These projects are designed to identify and monitor areas with the potential for land subsidence, to conduct basic research in the processes that control land subsidence and the development of earth fissures, as well as to develop new quantitative tools to predict aquifer-system deformation. In 1989 an ad hoc "Aquifer Mechanics and Subsidence Interest Group" (referred to herein as the "Subsidence Interest Group") was formed

California↗

USGS2017

Platinum-group elements

The platinum-group elements (PGEs)—platinum, palladium, rhodium, ruthenium, iridium, and osmium—are metals that have similar physical and chemical properties and tend to occur together in nature. PGEs are indispensable to many industrial applications but are mined in only a few places. The availability and accessibility of PGEs could be disrupted by economic, environmental, political, and social events. The United States net import reliance as a percentage of apparent consumption is about 90 percent. PGEs have many industrial applications. They are used in catalytic converters to reduce carbon monoxide, hydrocarbon, and nitrous oxide emissions in automobile exhaust. The chemical industry requires platinum or platinum-rhodium alloys to manufacture nitric oxide, which is the raw material used to manufacture explosives, fertilizers, and nitric acid. In the petrochemical industry, platinum-supported catalysts are needed to refine crude oil and to produce aromatic compounds and high-octane gasoline. Alloys of PGEs are exceptionally hard and durable, making them the best known coating for industrial crucibles used in the manufacture of chemicals and synthetic materials. PGEs are used by the glass manufacturing industry in the production of fiberglass and flat-panel and liquid crystal displays. In the electronics industry, PGEs are used in computer hard disks, hybridized integrated circuits, and multilayer ceramic capacitors. Aside from their industrial applications, PGEs are used in such other fields as health, consumer goods, and finance. Platinum, for example, is used in medical implants, such as pacemakers, and PGEs are used in cancer-fighting drugs. Platinum alloys are an ideal choice for jewelry because of their white color, strength, and resistance to tarnish. Platinum, palladium, and rhodium in the form of coins and bars are also used as investment commodities, and various financial instruments based on the value of these PGEs are traded on major exchanges. PGEs are among the rarest metals; Earth’s upper crust contains only about 0.0005 part per million (ppm) platinum. Today, the average grade of PGEs in ores that are mined primarily for their PGE concentrations varies from 5 to 15 ppm, although the concentration of PGEs in hand-picked ore specimens may range from tens to hundreds of parts per million. More than 100 different minerals have one of the PGEs as an essential component. PGE minerals occur as native metals. They also occur as compounds with other transition metals (copper, iron, mercury, nickel, and silver), post-transition metals (bismuth, lead, and tin), metalloids (antimony, arsenic, and tellurium), and nonmetals (selenium and sulfur). From 1900 to 2011, approximately 14,200 metric tons of PGEs was produced, and roughly 95 percent of that production (13,500 metric tons) took place between 1960 and 2011. The breakdown of production by country shows that, since 1900, about 90 percent of the production came from South Africa and Russia. The secondary supply of platinum, palladium, and rhodium is obtained through the recycling of catalytic converters from end-of-life vehicles, jewelry, and electronic equipment. Recycled platinum, palladium, and rhodium provide a significant proportion of the world’s total supply; these secondary sources are sufficient to close the gap between world mine production and consumption. Exploration and mining companies report resources of about 104,000 metric tons of PGEs (including minor amounts of gold) in mineral deposits around the world that could be developed. For PGEs, almost all the reported production and identified resources are associated with deposits in three geologic features—the Bushveld Complex, which is a layered mafic-to-ultramafic intrusion in South Africa; the Great Dyke, which is a layered mafic-to-ultramafic intrusion in Zimbabwe; and sill-like intrusions associated with flood basalts in the Noril’sk-Talnakh area of Russia. The metallic forms of PGEs are generally considered to be inert. PGEs pose a risk to human health only in cases where individuals are occupationally exposed to synthetic PGE compounds, especially workers in precious-metal refineries. In the natural environment, background PGE concentrations are low in water, sediment, soil, and plants. Anthropogenic sources of PGEs in the environment include catalytic converters used in modern automobiles, platinum-based chemotherapy drugs, and smelter emissions. The abundance of sulfide minerals defines the environmental and geologic characteristics of PGE-enriched magmatic sulfide deposits; those deposits with the highest amount of sulfide minerals could have the highest environmental impact. Acid rock drainage from reef-type and contact-type deposits is unlikely because the ores and their host rocks contain low proportions of sulfide minerals. For some conduit-type orebodies with massive ores, mineral-processing techniques separate and produce concentrates of copper-, iron-, and nickel-bearing sulfide minerals; those with copper and nickel are processed to extract metal, but the iron-sulfide minerals, mainly pyrrhotite, are discarded as waste. This results in waste material with a high acid-generating potential. The most significant primary source of PGEs in the United States is a deposit in the Stillwater Complex, which is a layered igneous intrusion in Montana. Approximately 305 metric tons of platinum and palladium have been mined from the Stillwater Complex deposit since 1986. Exploration and development drilling indicate that another 2,200 metric tons are present. Mining has progressed to depths of 1,800 meters below the surface, but the bottom of the ore deposit has not been reached; geologic estimates suggest that another 1,000 to 6,200 metric tons of PGEs could be present at depth. In the future, PGEs may be mined from deposits found near the base of the Duluth Complex, which is a group of igneous intrusions in Minnesota.

Professional Paper↗

USGS1999

Cobalt-rich ferromanganese crusts in the Pacific

Co-rich Fe-Mn crusts occur throughout the Pacific on seamounts, ridges, and plateaus where currents have kept the rocks swept clean of sediments at least intermittently for millions of years. Crusts precipitate out of cold ambient sea water onto hard-rock substrates forming pavements up to 250 mm thick. Crusts are important as a potential resource for Co, Ni, Pt, Mn, Tl, Te, and other metals, as well as for the paleoclimate signals stored in their stratigraphic layers. Crusts form at water depths of about 400 to 4000 m, with the thickest and most Co-rich crusts occurring at depths of about 800 to 2500 m, which may vary on a regional scale. Gravity processes, sediment cover, submerged and emergent reefs, and currents control the distribution and thickness of crusts on seamounts. Crusts occur on a variety of substrate rocks that generally decrease in the order, breccia, basalt, phosphorite, limestone, hyaloclastite, and mudstone. Because of this wide variety of substrate types, crusts are difficult to distinguish from the substrate using remotely sensed data, such as geophysical measurements, but are generally weaker and lighter-weight than the substrate. Crusts can be distinguished from the substrates, however, by their much higher gamma radiation levels. The mean dry bulk density of crusts is 1.3 g/cm3, the mean porosity is 60%, and the mean surface area is extremely high, 300 m2/g. Crusts generally grow at rates of 1 to 10 mm/Ma. Crust surfaces are botryoidal, which may be modified to a variety of forms by current erosion. In cross-section, crusts are generally layered, with individual layers displaying massive, botryoidal, laminated, columnar, or mottled textures. Characteristic layering is persistent regionally in the Pacific. Crusts are composed of ferruginous vernadite (δ-MnO2) and X-ray amorphous Fe oxyhydroxide, with moderate amounts of carbonate fluorapatite (CFA) in thick crusts and minor amounts of quartz and feldspar in most crusts. Elements most commonly associated with the vernadite phase include Mn, Co, Ni, Cd, and Mo, whereas those most commonly associated with Fe oxyhydroxide are Fe and As. Detrital phases are represented by Si, Al, K, Ti, Cr, Mg, Fe, and Na; the CFA phase by Ca, P, Sr, Y, and CO2; and a residual biogenic phase by Ba, Sr, Ce, Cu, V, Ca, and Mg. Crusts contain Co contents up to about 2.3%, Ni to 1%, and Pt to 3 ppm, with mean Fe/Mn ratios of 0.6 to 1.3. Fe/Mn decreases, whereas Co, Ni, Ti, and Pt increase in central Pacific crusts and Fe/Mn, Si, and Al increase in continental margin crusts and in crusts with proximity to west Pacific volcanic arcs. Vernadite and CFA-related elements decrease, whereas Fe, Cu, and detrital-related elements increase with increasing water depth of crust occurrence. Cobalt, Ce, Tl, and maybe also Ti, Pb, and Pt are strongly concentrated in crusts over other metals because of oxidation reactions. Total rare earth elements (REEs) commonly vary between 0.1% and 0.3% and are derived from sea water along with other hydrogenetic elements, Co, Mn, Ni, etc. Platinum, Rh, Ir, and some Ru in crusts are also derived from sea water, whereas Pd and the remainder of the Ru derive from detrital minerals. The older parts of thick crusts were phosphatized during at least two global phosphogenic events during the Tertiary, which mobilized and redistributed elements in those parts of the crusts. 240Silicon, Fe, Al, Th, Ti, Co, Mn, Pb, and U are commonly depleted, whereas Ni, Cu, Zn, Y, REEs, Sr, and Pt are commonly enriched in phosphatized layers compared to younger nonphosphatized layers. The dominant controls on the concentration of elements in crusts include the concentration of metals in sea water and their ratios, colloid surface charge, types of complexing agents, surface area, and growth rates. Crusts act as closed systems with regard to the isotopic ratios of Be, Nd, Pb, Hf, Os, and U-series, which in part have been used to date crusts and in part used as isotopic tracers of paleoceanographic and paleoclimatic conditions. Those tracers are especially useful in delineating temporal changes in deep-ocean circulation. Research and development on the technology of mining crusts are only in their infancy. Detailed maps of crust deposits and a better understanding of small-scale seamount topography are required to design the most appropriate mining equipment.

Book chapter↗

USGS2012

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗

USGS2016

Users and uses of Landsat 8 satellite imagery—2014 survey results

Executive Summary In 2013, Landsat 8 began adding high quality, global, moderate-resolution imagery to the more than 40-year archive of Landsat imagery. To assess the potential effects of the availability of Landsat 8 imagery on users and their work, the U.S. Geological Survey (USGS) Land Remote Sensing Program (LRS) initiated a survey of Landsat users. The objectives of the survey were to 1. Characterize various Landsat user groups, such as United States (U.S.) and international users and Landsat 8 and non-Landsat 8 users; 2. Identify any differences among user groups in uses and preferences; 3. Measure the importance of and satisfaction with Landsat 8 attributes; 4. Assess the importance to users of the frequency of usable imagery; and 5. Determine any challenges in using Landsat 8. The online survey was sent to 51,617 Landsat users registered with USGS in May 2014. Almost 13,000 people responded to the survey for a response rate of 25 percent (n = 12,966). Current Landsat users (users who had used Landsat in their work in the year prior to the survey) composed 89 percent of the sample (n = 11,549) and past Landsat users composed 11 percent (n = 1,417). The results reported here apply to current Landsat users registered with the USGS Earth Resources Observation and Science (EROS) Center. Users from 161 countries responded to the survey. Of those, 19 percent were citizens or permanent residents of the United States and 81 percent resided in other countries. More than 70 percent of current users had used Landsat 8 in the year prior to the survey. The majority of Landsat 8 users (65 percent) were established users who used Landsat imagery regularly both before and after Landsat 8 imagery became available. The average current Landsat user was male, 36 years old, and highly educated, with 9 years of experience using satellite imagery or geographic information system (GIS) software. Landsat 8 users had, on average, two more years of experience than non-Landsat 8 users. Users were employed predominantly by academic institutions (65 percent), followed by private businesses (13 percent), Federal governments (10 percent), State and local governments (6 percent), and nonprofit organizations (6 percent). Of the Landsat imagery obtained in the past year by current users, on average 31 percent came from a Landsat 8 sensor. An equivalent amount came from the Landsat 7 ETM+ sensor (33 percent); slightly less came from Landsats 4 and 5 TM sensors (27 percent). Much less came from Landsats 1 through 5 MSS sensors (5 percent). Overall, more than a third of users’ work used Landsat imagery (38 percent). Of this work, on average, 37 percent of the work was operational. Landsat 8 users considered a greater proportion of their work operational than non-Landsat 8 users (39 percent compared with 29 percent). Environmental sciences and management were the most commonly selected primary applications (selected by 42 percent of users). Land use/land cover (23 percent) was the second most commonly selected primary application, followed by education (12 percent), agriculture (9 percent), and planning and development (6 percent). Landsat 8 users were asked to rank the importance of certain attributes in determining whether to use Landsat 8 imagery in their work. The archive was ranked most important, followed by cost, spatial resolution, extent of coverage, data quality, and frequency of revisit. Users were asked how satisfied they were with these same attributes as they currently apply to Landsat 8 imagery. On average, users were most satisfied with lack of cost, extent of coverage, data quality, and the archive, but they were satisfied with all attributes. Users were asked how often they needed Landsat imagery to meet various requirements for their primary application. The survey question specifically asked how often users needed usable imagery, which differs from how often they would like the Landsat satellites to acquire an image. Users were asked to identify their needed frequency of usable imagery for the following levels: 1. Threshold level—the minimum frequency of usable imagery needed to be of any value to their primary application. 2. Breakthrough level—the frequency of usable imagery that would result in a significant improvement for their primary application of the imagery. 3. Target level—the frequency of usable imagery that would only provide a limited additional increase in the expected performance for their primary application. To meet the threshold level, three-quarters of users needed usable imagery every 17 days or less frequently. At the breakthrough level, two-thirds of users (64 percent) needed a usable image every 5–16 days. The current constellation of two satellites (Landsat 7 and 8) is capable of meeting the threshold and breakthrough needs of most users at least some of the time, but a single satellite would be highly unlikely to do so. Two-fifths of users (40 percent) felt that usable imagery provided every 4 days or more frequently would meet their target level which the current Landsat constellation cannot provide. Landsat 8 users were significantly more likely than non-Landsat 8 users to need usable imagery more frequently to meet their target levels. Additionally, U.S. Landsat 8 users were significantly more likely than other Landsat users to need usable imagery more frequently in order meet both their breakthrough and target levels. To explore the effect of the availability of Landsat 8 imagery on Landsat imagery use in general, established users (those who had consistently used Landsat imagery both before and after Landsat 8 imagery became available) using Landsat 8 imagery were asked about changes in the amount of Landsat imagery they used. The majority of established users using Landsat 8 imagery (60 percent) reported an average increase of 51 percent in the number of scenes obtained after Landsat 8 imagery became available. Landsat 8 users were asked if they had encountered challenges in using Landsat 8 whereas non-Landsat 8 users were asked if such challenges had played a role in why they were not using Landsat 8 imagery. Although many users did not encounter challenges when using or trying to use Landsat 8 data, slightly less than 30 percent did encounter issues with processing the data to a usable point. The most common issue reported was not being able to create or have access to a surface reflectance corrected product. Other challenges were related to the file sizes of images being too large to download, store, or analyze. There were no statistically significant differences between Landsat 8 and non-Landsat 8 users in terms of challenges encountered when using or trying to use the imagery, which indicates that users were not unduly discouraged by the challenges they may have encountered. When asked about potential consequences of not using Landsat 8, more than half of the non-Landsat 8 users did not report detrimental effects on their work from not using the imagery. Of those who did report detrimental effects, decreased quality of work, decreased scope of work, and increased time spent on work were the most common.

Open-File Report↗