USGS ScienceSearch

SEARCH · USGS Science

Results for “Earth”

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

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

At least 289 records · Page 16Linked to original sources

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

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

International Journal of Applied Earth Observation

Book review: Earth’s magnetism in the age of sail

For many centuries, the source, behavior, and even the essential nature of geomagnetism were enigmatic. Despite this, the effect of geomagnetism was familiar, by imparting a directional preference on the magnetized needle of the compass and providing a useful, if somewhat annoyingly complicated, reference for navigators. Although the compass seems to have first been invented in China, it was the Europeans who made the most systematic early studies of magnetism, who made the first elaborate and practical usage of the compass, and who developed most of the early theories as to the cause of the compass needle’s north-seeking tendency. From the centuries of the Middle Ages, through the late 16th century of the Renaissance, to the 17th century of philosophical enlightenment and the 18th century of discovery, the subject of magnetism and, more specifically, geomagnetism, evolved from a hodgepodge of mystical beliefs into something that we can today recognize as the object of modern scientific pursuit. Those same centuries witnessed the great transoceanic sailing voyages undertaken by European nations for reasons of exploration, territorial claim, religious mission, and mercantile trade. Naturally, the navigator’s compass, and therefore geomagnetism, played an important role in these developments. This romantic intersection of science and history is the subject of Earth’s Magnetism in the Age of Sail , a pleasantly written and scholarly book by A.R.T. Jonkers.

Physics of the Earth and Planetary Interiors

Deep earthquakes beneath the Fiji Basin, SW Pacific: Earth's most intense deep seismicity in stagnant slabs

Previous work has suggested that many of the deep earthquakes beneath the Fiji Basin occur in slab material that has been detached and foundered to the bottom of the transition zone or has been laid down by trench migration in a similar recumbent position. Since nowhere else in the Earth do so many earthquakes occur in slabs stagnated in the transition zone, these earthquakes merit closer study. Accordingly, we have assembled from historical and modern data a comprehensive catalogue of the relocated hypocenters and focal mechanisms of well-located deep events in the geographic area between the bottoms of the main Vanuatu and Tonga Wadati-Benioff zones. Two regions of deep seismogenesis are recognized there: (i) 163 deep shocks have occurred north of 15??S in the Vityaz Group from 1949 to 1996. These seismological observations and the absence of other features characteristic of active subduction suggest that the Vityaz group represents deep failure in a detached slab that has foundered to a horizontal orientation near the bottom of the transition zone. (ii) Another group of nearly 50 'outboard' deep shocks occur between about 450 and 660 km depth, west of the complexly buckled and offset western edge of the Tonga Wadati-Benioff zone. Their geometry is in the form of two or possibly three small-circle arcs that roughly parallel the inferred motion of Tonga trench migration. Earthquakes in the southernmost of these arcs occur in a recumbent high-seismic-wavespeed slab anomaly that connects both to the main inclined Tonga anomaly to the east and a lower mantle anomaly to the west [Van der Hilst, R., 1995. Complex morphology of subducted lithosphere in the mantle beneath the Tonga trench. Nature, Vol. 374, pp. 154-157.]. Both groups show complexity in their focal mechanisms. The major question raised by these observations is the cause of this apparent temporary arrest in the descent of the Tonga slab into the lower mantle. We approach these questions by considering the effects of buoyant metastable peridotite in cold slab material that was detached and rapidly foundered, or was buckled, segmented and laid out in the transition zone.

Physics of the Earth and Planetary Interiors

Human-in-the-Loop segmentation of earth surface imagery

Segmentation, or the classification of pixels (grid cells) in imagery, is ubiquitously applied in the natural sciences. Manual methods are often prohibitively time-consuming, especially those images consisting of small objects and/or significant spatial heterogeneity of colors or textures. Labeling complicated regions of transition that in Earth surface imagery are represented by collections of mixed-pixels, -textures, and -spectral signatures, can be especially error-prone because it is difficult to reliably unmix, identify and delineate consistently. However, the success of supervised machine learning (ML) approaches is entirely dependent on good label data. We describe a fast, semi-automated, method for interactive segmentation of N-dimensional (x, y, N) images into two-dimensional (x, y) label images. It uses human-in-the-loop ML to achieve consensus between the labeler and a model in an iterative workflow. The technique is reproducible; the sequence of decisions made by human labeler and ML algorithms can be encoded to file, so the entire process can be played back and new outputs generated with alternative decisions and/or algorithms. We illustrate the scientific potential of segmentation of imagery of diverse settings and image types using six case studies from river, estuarine, and open coast environments. These photographic and non-photographic imagery consist of 1- and 3-bands on regular and irregular grids ranging from centimeters to tens of meters. We demonstrate high levels of agreement in label images generated by several labelers on the same imagery, and make suggestions to achieve consensus and measure uncertainty, ideal for widespread application in training supervised ML for image segmentation.

Earth and Space Science

Multi-year tracing of spatial and temporal dynamics of post-fire aeolian sediment transport using rare earth elements provide insights into grassland management

Aeolian sediment transport occurs as a function of, and with feedback to ecosystem changes and disturbances. Many desert grasslands are undergoing rapid changes in vegetation, including the encroachment of woody plants, which alters fire regimes and in turn can change the spatial and temporal patterns of aeolian sediment transport. We investigated aeolian sediment transport and spatial distribution of sediment in the surface soil for 7 years following a prescribed fire using a multiple rare earth element (REE) tracer-based approach in a shrub-encroached desert grassland in the northern Chihuahuan desert. Results indicate that even though the aeolian horizontal sediment mass flux increased approximately three-fold in the first windy season in the burned areas compared to control areas, there were no significant differences after three windy seasons. The soil surface of bare microsites was the major contributor of aeolian sediments in unburned areas (87%), while the shrub microsites contributed the least (<2%) during the observation period. However, after the prescribed fire, the contribution of aeolian sediments from shrub microsites increased considerably (∼40%), indicating post-fire microsite-scale sediment redistribution. The findings of this study, which is the first to use multiple REE tracers for multi-year analysis of the spatial and temporal dynamics of aeolian sediment transport, illustrate how disturbance by prescribed fire can influence aeolian processes and alters dryland soil geomorphology in which distinct soils develop over time at very fine spatial scales of individual plants.

New Mexico

Earth’s mantle composition revealed by mantle plumes

Mantle plumes originate at depths near the core−mantle boundary (~2,800 km). As such, they provide invaluable information about the composition of the deep mantle and insight into convection, crustal formation, and crustal recycling, as well as global heat and volatile budgets. In this Review, we discuss the effectiveness and challenges of using isotopic analyses of plume-generated rocks to infer mantle composition and to constrain geodynamic models. Isotopic analyses of plume-derived ocean island basalts, including radiogenic (Sr, Nd, Pb, Hf, W, noble gas) and stable isotopes (Li, C, O, S, Fe, Tl), permit determination of mantle plume composition, which in turn generate insight into mantle plume origins, dynamics, mantle heterogeneities, early-formed mantle reservoirs, crustal recycling processes, core−mantle interactions and mantle evolution. Nevertheless, the magmatic flux, temperature, tectonic environment and compositions of mantle plumes can vary. Consequently, plumes and their melts are best evaluated along a spectrum that acknowledges their different properties, particularly mantle flux, before making interpretations about the interior of the Earth. To provide insight into specific mantle and plume processes, future work should document correlations across elemental and isotopic data sets on the same sample powder, coordinate targeting sampling strategies, and refine stable isotopic fractionation factors through experiments. Such work will benefit from collaboration across geochemical laboratories, as well as among geochemists, mineral physicists, seismologists and geodynamicists.

Nature Reviews in Earth and Environment

Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model

Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

Hydrology and Earth System Sciences

Mobile laser scanning applied to the earth sciences

Lidar (light detection and ranging), a method by which the precise time of flight of emitted pulses of laser energy is measured and converted to distance for reflective targets, has helped scientists make topographic maps of Earth's surface at scales as fine as centimeters. These maps have allowed the discovery and analysis of myriad otherwise unstudied features, such as fault scarps, river channels, and even ancient ruins [Glennie et al., 2013b].

Eos, Transactions, American Geophysical Union

Entrainment of bed material by Earth-surface mass flows: review and reformulation of depth-integrated theory

Earth-surface mass flows such as debris flows, rock avalanches, and dam-break floods can grow greatly in size and destructive potential by entraining bed material they encounter. Increasing use of depth-integrated mass- and momentum-conservation equations to model these erosive flows motivates a review of the underlying theory. Our review indicates that many existing models apply depth-integrated conservation principles incorrectly, leading to spurious inferences about the role of mass and momentum exchanges at flow-bed boundaries. Model discrepancies can be rectified by analyzing conservation of mass and momentum in a two-layer system consisting of a moving upper layer and static lower layer. Our analysis shows that erosion or deposition rates at the interface between layers must in general satisfy three jump conditions. These conditions impose constraints on valid erosion formulas, and they help determine the correct forms of depth-integrated conservation equations. Two of the three jump conditions are closely analogous to Rankine-Hugoniot conditions that describe the behavior of shocks in compressible gasses, and the third jump condition describes shear traction discontinuities that necessarily exist across eroding boundaries. Grain-fluid mixtures commonly behave as compressible materials as they undergo entrainment, because changes in bulk density occur as the mixtures mobilize and merge with an overriding flow. If no bulk density change occurs, then only the shear-traction jump condition applies. Even for this special case, however, accurate formulation of depth-integrated momentum equations requires a clear distinction between boundary shear tractions that exist in the presence or absence of bed erosion.

Reviews of Geophysics

Time causal operational estimation of electric fields induced in the Earth's lithosphere during magnetic storms

In support of projects for monitoring geomagnetic hazards for electric power grids, we develop a simple mathematical formalism, consistent with the time causality of deterministic physics, for estimating electric fields that are induced in the Earth's lithosphere during magnetic storms. For an idealized model of the lithosphere, an infinite half‐space having uniform electrical conductivity properties described by a galvanic tensor, we work in the Laplace‐transformed frequency domain to obtain a transfer function which, when convolved with measured magnetic field time series, gives an estimated electric field time series. Using data collected at the Kakioka, Japan observatory, we optimize lithospheric conductivity parameters by minimizing the discrepancy between model‐estimated electric field variation and that actually measured. With our simple model, we can estimate 87% of the variance in storm time Kakioka electric field data; a more complicated model of lithospheric conductivity would be required to estimate the remaining 13% of the variance. We discuss how our estimation formalism might be implemented for geographically coordinated real‐time monitoring of geoelectric fields.

Geophysical Research Letters

High-latitude dust in the Earth system

Natural dust is often associated with hot, subtropical deserts, but significant dust events have been reported from cold, high latitudes. This review synthesizes current understanding of high-latitude (≥50°N and ≥40°S) dust source geography and dynamics and provides a prospectus for future research on the topic. Although the fundamental processes controlling aeolian dust emissions in high latitudes are essentially the same as in temperate regions, there are additional processes specific to or enhanced in cold regions. These include low temperatures, humidity, strong winds, permafrost and niveo-aeolian processes all of which can affect the efficiency of dust emission and distribution of sediments. Dust deposition at high latitudes can provide nutrients to the marine system, specifically by contributing iron to high-nutrient, low-chlorophyll oceans; it also affects ice albedo and melt rates. There have been no attempts to quantify systematically the expanse, characteristics, or dynamics of high-latitude dust sources. To address this, we identify and compare the main sources and drivers of dust emissions in the Northern (Alaska, Canada, Greenland, and Iceland) and Southern (Antarctica, New Zealand, and Patagonia) Hemispheres. The scarcity of year-round observations and limitations of satellite remote sensing data at high latitudes are discussed. It is estimated that under contemporary conditions high-latitude sources cover >500,000 km2 and contribute at least 80–100 Tg yr−1 of dust to the Earth system (~5% of the global dust budget); both are projected to increase under future climate change scenarios.

Reviews of Geophysics

Quantifying postfire aeolian sediment transport using rare earth element tracers

Grasslands, which provide fundamental ecosystem services in many arid and semiarid regions of the world, are undergoing rapid increases in fire activity and are highly susceptible to postfire-accelerated soil erosion by wind. A quantitative assessment of physical processes that integrates fire-wind erosion feedbacks is therefore needed relative to vegetation change, soil biogeochemical cycling, air quality, and landscape evolution. We investigated the applicability of a novel tracer technique—the use of multiple rare earth elements (REE)—to quantify soil transport by wind and to identify sources and sinks of wind-blown sediments in both burned and unburned shrub-grass transition zone in the Chihuahuan Desert, NM, USA. Results indicate that the horizontal mass flux of wind-borne sediment increased approximately threefold following the fire. The REE tracer analysis of wind-borne sediments shows that the source of the horizontal mass flux in the unburned site was derived from bare microsites (88.5%), while in the burned site it was primarily sourced from shrub (42.3%) and bare (39.1%) microsites. Vegetated microsites which were predominantly sinks of aeolian sediments in the unburned areas became sediment sources following the fire. The burned areas showed a spatial homogenization of sediment tracers, highlighting a potential negative feedback on landscape heterogeneity induced by shrub encroachment into grasslands. Though fires are known to increase aeolian sediment transport, accompanying changes in the sources and sinks of wind-borne sediments may influence biogeochemical cycling and land degradation dynamics. Furthermore, our experiment demonstrated that REEs can be used as reliable tracers for field-scale aeolian studies.

Journal of Geophysical Research: Biogeosciences

Down to Earth with an electric hazard from space

In reaching across traditional disciplinary boundaries, solid-Earth geophysicists and space physicists are forging new collaborations to map magnetic-storm hazards for electric-power grids. Future progress in evaluation storm time geoelectric hazards will come primarily through monitoring, surveys, and modeling of related data.

Space Weather

Invasive buffelgrass detection using high-resolution satellite and UAV imagery on Google Earth Engine

Methods to detect and monitor the spread of invasive grasses are critical to avoid ecosystem transformations and large economic costs. The rapid spread of non‐native buffelgrass (Pennisetum ciliare) has intensified fire risk and is replacing fire intolerant native vegetation in the Sonoran Desert of the southwestern US. Coarse‐resolution satellite imagery has had limited success in detecting small patches of buffelgrass, whereas ground‐based and aerial survey methods are often cost prohibitive. To improve detection, we trained 2 m resolution DigitalGlobe WorldView‐2 satellite imagery with 12 cm resolution unmanned aerial vehicle (UAV) imagery and classified buffelgrass on Google Earth Engine, a cloud computing platform, using Random Forest (RF) models in Saguaro National Park, Arizona, USA. Our classification models had an average overall accuracy of 93% and producer's accuracies of 94–96% for buffelgrass, although user's accuracies were low. We detected a 2.92 km 2 area of buffelgrass in the eastern Rincon Mountain District (1.07% of the total area) and a 0.46 km 2 area (0.46% of the total area) in the western Tucson Mountain District of Saguaro National Park. Buffelgrass cover was significantly greater in the Sonoran Paloverde‐Mixed Cacti Desert Scrub vegetation type, on poorly developed Entisols and Inceptisol soils and on south‐facing topographic aspects compared to other areas. Our results demonstrate that high‐resolution imagery improve on previous attempts to detect and classify buffelgrass and indicate potential areas where the invasive grass might spread. The methods demonstrated in this study could be employed by land managers as a low‐cost strategy to identify priority areas for control efforts and continued monitoring.

Arizona

Analysis of low levels of rare earths by radiochemical neutron activation analysis

A procedure for the radiochemical neutron-activation analysis for the rare earth elements (REE) involves the separation of the REE as a group by rapid ion-exchange methods and determination of yields by reactivation or by energy dispersive X-ray fluorescence (EDXRF) spectrometry. The U. S. Geological Survey (USGS) standard rocks, BCR-1 and AGV-1, were analyzed to determine the precision and accuracy of the method. We found that the precision was ??5-10% on the basis of replicate analysis and that, in general the accuracy was within ??5% of accepted values for most REE. Data for USGS standard rocks BIR-1 (Icelandic basalt) and DNC-1 (North Carolina diabase) are also presented. ?? 1985 Akade??miai Kiado??.

Journal of Radioanalytical and Nuclear Chemistry A

Lava flows in mare imbrium: An evaluation of anomalously low earth-based radar reflectivity

The lunar maria reflect two to five times less Earth-based radar power than the highlands, the spectrally blue maria surfaces returning the lowest power levels. This effect of weakening signal return has been attributed to increased signal absorption related to the electrical and magnetic characteristics of the mineral ilmenite (FeTiO 3 ). The surface of Mare Imbrium contains some of the most distinct red-blue colorimetric boundaries and depolarized 70 cm wavelength reflectivity variations on the near side of the Moon. The weakest levels of both 3.8 cm and 70 cm reflectivity within Imbrium are confined to regional mare surfaces of the blue spectral type that can be recognized as stratigraphically unique flow surfaces. Frequency distributions of the 70 cm polarized and depolarized radar return power for five mare surfaces within the basin indicate that signal absorption, and probably the ilmenite content, increases generally from the beginning of the Imbrian Period to the end of the Eratosthenian Period with slight reversal between the end of the Imbrian and beginning of the Eratosthenian. TiO 2 calibrated radar reflectivity curves can be utilized for lunar maria geochemical mapping in the same manner as the TiO 2 calibrated spectral reflectivity curves of Charette et al. (1974). The long wavelength radar data may be a sensitive indicator of mare chemical variations as it is unaffected by the normal surface rock clutter that includes ray materials from large impact craters.

The Moon

Modeling the mid-Piacenzian warm climate using the water isotope-enabled Community Earth System Model (iCESM1.2-ITPCAS)

The mid-Piacenzian Warm Period (MPWP, ~ 3.264–3.025 Ma) is the most recent example of a persistently warmer climate in equilibrium with atmospheric CO 2 concentrations similar to today. Towards studying patterns and dynamics of a warming climate the MPWP is often compared to today. Following the Pliocene Model Intercomparison Project, Phase 2 (PlioMIP2) protocol we prepare a water isotope-enabled Community Earth System Model (iCESM1.2) simulation that is warmer and wetter than the PlioMIP2 multi-model ensemble (MME). While our simulation resembles PlioMIP2 MME in many aspects we find added insights. (1) Considerable warmth at high latitudes exceeds previous simulations. Polar amplification (PA) is comparable to proxies, enabled by iCESM1.2’s high climate sensitivity and a distinct method of ocean initialization. (2) Major driver of warmth is the downward component of clear-sky surface long-wave radiation (Δ T rlds_clearsky ). (3) In iCESM1.2 modulated dominance of dynamic (δDY) processes causes different low-latitude (~ 30 S°–10°N) precipitation response than the PlioMIP2 MME, where thermodynamic processes (δTH) dominate. (4) Modulated local condensation leads to lower δ18O p across tropical Indian Ocean and surrounding Asian-African-Australian monsoon regions. (5) We find contrasting changes in tropical atmospheric circulations (Hadley and Walker cells). Anomalous regional meridional (zonal) circulation, forced by changes in tropical-subtropical (tropical) diabatic processes, presents a more comprehensive perspective than explaining weakened and expanded Hadley circulation (strengthened and westward-shifted Walker circulation) via static stability. (6) Enhanced Atlantic meridional overturning circulation owes to a closed Bering Strait.

Climate Dynamics

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $3.2 per tonne of coal and $186 billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$4.4 per tonne of ash and $1.2 billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

Alabama, Florida, Louisiana, Mississippi, Texas