Search USGSSearch

SEARCH · Search USGS

Results for “Earth's Future”

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 775 records · Page 43Linked to original sources

Soil biogeochemical responses of a tropical forest to warming and hurricane disturbance

Tropical forests represent <15% of Earths terrestrial surface yet support >50% of the planets species and play a disproportionately large role in determining climate due to the vast amounts of carbon they store and exchange with the atmosphere. Currently, disturbance patterns in tropical ecosystems are changing due to factors such as increased land use pressure and an occurrence of hurricanes. At the same time, these regions are expected to experience unprecedented warming before 2100. Despite the importance of these ecosystems for forecasting the global consequences of multiple stressors, our understanding of how projected changes in climate and disturbance will affect the biogeochemical cycling of tropical forests remains in its infancy. Until now, no studies to our knowledge have evaluated forest recovery following hurricane disturbance within the context of concurrent climatic change. Here, we present soil biogeochemical results from a tropical forest field warming experiment in Puerto Rico where, a year after experimental warming began, Hurricanes Irma and Mara greatly altered the forest, allowing a unique opportunity to explore the interacting effects of hurricane disturbance and warming. We tracked post-hurricane forest recovery for a year without warming to assess legacy effects of prior warming on the disturbance response, and then reinitiated warming treatments to further evaluate interactions between forest recovery and warmer temperatures. The data showed that warming affected multiple aspects of soil biogeochemical cycling even in the first year of treatment, with particularly large positive effects on soil microbial biomass pools (e.g., increases of 54, 43, and 46% relative to the control plots were observed for microbial biomass carbon, nitrogen, and phosphorus, respectively after 6 months of warming). We also observed significant effects of the hurricanes on soil biogeochemical cycling, as well as interactive controls of warming and disturbance. Taken together, our results showed dynamic soil responses that suggest the future of biogeochemical cycling in this tropical wet forest will be strongly shaped by the directional effects of warming and the episodic effects of hurricanes.

Puerto Rico

Reconnaissance study of the major and trace element content of bauxite deposits in the Arkansas bauxite region, Saline and Pulaski Counties, central Arkansas

The Arkansas bauxite district, which comprises about 275 square miles (710 square kilometers) of central Arkansas, produced an order of magnitude more bauxite and alumina than the other bauxite districts in the United States combined. Bauxite was mined in the region continuously from 1898 to 1982. These bauxites are laterite deposits, formed from intensive in-place weathering of the exposed surface of the Granite Mountain pluton, a Late Cretaceous batholith composed mainly of nepheline syenite and lesser amounts of syenite. Nepheline syenite was the aluminum source for the bauxite and clay deposits that blanket the pluton. The early Eocene continental sedimentary rocks that contain and overlie the bauxite deposits indicate that central Arkansas had a warm tropical environment during bauxite formation. Bauxite ores are the principal sources of aluminum. Some of the global bauxite deposits have been found to contain co-occurring metals that have essential applications in modern technologies. For example, bauxite is the largest global source of gallium (Ga), used in semiconductors, which is recovered as a byproduct of processing bauxite to recover alumina. Other critical metal commodities within some bauxites that reportedly have potential for byproduct recovery include niobium (Nb), scandium (Sc), and rare earth elements (REEs). Currently (2021), the United States is wholly dependent on imports for its supplies of bauxite for processing to produce alumina. The United States is also dependent on foreign sources of gallium, niobium, and scandium, as well for most of its domestic requirements of REEs. For these reasons, samples were collected from Arkansas bauxite deposits, associated clays, mill residue wastes (respectively referred to as red muds and black sands), and the parent nepheline syenite to determine their elemental content, with a particular focus on gallium, niobium, scandium, and REEs. Each sample was analyzed for 60 elements; these data and the methods used are published as a U.S. Geological Survey data release. The results indicate that, of the critical metals in bauxites, gallium is a potential byproduct from the central Arkansas bauxite deposits. The highest gallium concentrations occur in the raw bauxite ore, with an average concentration of 76 parts per million (ppm). Gallium partitions with alumina (the product) rather than into mine waste residues. Results indicate an average niobium content of 662 ppm in the Arkansas bauxite ores. Niobium progressively increases in concentration from parent syenite (247 ppm) to clays (315 ppm) and further from bauxite (662 ppm) to processed residues (1,075 ppm). Low concentrations of scandium were found in all samples, averaging 10 ppm or less in the parent rock (syenite), bauxite, clays, and processing residues. Modest concentrations of the light and heavy REEs were found in samples of bauxite ores, bauxitic clays and interbedded clays, syenite, and the residues of ore. The highest REE values were found in processed residues, with average concentrations of 613 ppm total light REEs and 130 ppm total heavy REEs. These concentrations suggest that additional processing to recover REEs is unlikely to be economic in the foreseeable future.

Arkansas

Detecting sea-level hazards: Simple regression-based methods for calculating the acceleration of sea level

This report documents the development of statistical tools used to quantify the hazard presented by the response of sea-level elevation to natural or anthropogenic changes in climate and ocean circulation. A hazard is a physical process (or processes) that, when combined with vulnerability (or susceptibility to the hazard), results in risk. This study presents the development and comparison of new and existing sea-level analysis methods, exploration of the strengths and weaknesses of the methods using synthetic time series, and when appropriate, synthesis of the application of the method to observed sea-level time series. These reports are intended to enhance material presented in peer-reviewed journal articles where it is not always possible to provide the level of detail that might be necessary to fully support or recreate published results. The purpose of this report is to document and compare three simple methodologies that have previously been used to provide estimates with associated errors of the acceleration of sea-level elevation. These techniques have been used by coastal scientists and planners in assessing coastal risk over a wide range of spatial and temporal scales. Because relative sea-level (SL) elevation time series contain energetic fluctuations at many time scales, extracting what can be relatively small rate and acceleration signals (along with estimates of the error) from much larger &ldquo;noise&rdquo; has proven to be both difficult and controversial. Acceleration is a preferred measure of SL response to recent changes in the Earth&rsquo;s climate because over time scales of 100 years or less slow vertical land motions (such as glacial isostatic adjustment) contribute only to the linear signal and not to acceleration, thus reducing the complexity of the analysis. Hence acceleration is useful if the goal of a study is to characterize and quantify the hazard associated with the changing relative elevation of water with respect to land on decadal time scales. Although in some cases it may be necessary to determine the cause of relative sea level rise, as a first step, it is important to accurately estimate the magnitude of the threat. Most researchers agree that global sea level (GSL) rose persistently through much of the 20th century at about 1.5&ndash;2.0 millimeters per year (mm/yr). There is far less agreement about whether the rate of sea-level rise (SLR) is increasing (that is, an acceleration in SL). Recent studies, and most of their predecessors, use tide gage data to quantify SL acceleration, A SL (t). In the current study, three techniques were used to calculate acceleration from tide gage data, and of those examined, it was determined that the two techniques based on sliding a regression window through the time series are more robust compared to the technique that fits a single quadratic form to the entire time series, particularly if there is temporal variation in the magnitude of the acceleration. The single-fit quadratic regression method has been the most commonly used technique in determining acceleration in tide gage data. The inability of the single-fit method to account for time-varying acceleration may explain some of the inconsistent findings between investigators. Properly quantifying A SL (t) from field measurements is of particular importance in evaluating numerical models of past, present, and future SLR resulting from anticipated climate change.

Open-File Report

Land area change in coastal Louisiana from 1932 to 2010

Coastal Louisiana wetlands make up the seventh largest delta on Earth, contain about 37 percent of the estuarine herbaceous marshes in the conterminous United States, and support the largest commercial fishery in the lower 48 States. These wetlands are in peril because Louisiana currently undergoes about 90 percent of the total coastal wetland loss in the continental United States. Documenting and understanding the occurrence and rates of wetland loss are necessary for effective planning, protection, and restoration activities. The analyses of landscape change presented in this report use historical surveys, aerial data, and satellite data to track landscape changes. Summary data are presented for 1932-2010; trend data are presented for 1985-2010. These later data were calculated separately because of concerns over the comparability of the 1932 and 1956 datasets (which are based on survey and aerial data, respectively) with the later datasets (which are all based on satellite imagery). These analyses show that coastal Louisiana has undergone a net change in land area of about -1,883 square miles (mi 2 ) from 1932 to 2010. This net change in land area amounts to a decrease of about 25 percent of the 1932 land area. Persistent losses account for 95 percent of this land area decrease; the remainder are areas that have converted to water but have not yet exhibited the persistence necessary to be classified as "loss." Trend analyses from 1985 to 2010 show a wetland loss rate of 16.57 mi 2 per year. If this loss were to occur at a constant rate, it would equate to Louisiana losing an area the size of one football field per hour. The use of 17 datasets plus the application of consistent change criteria in this study provide opportunities to better understand the timing and causal mechanisms of wetland loss that are critical for forecasting landscape changes in the future.

Louisiana

Executive summary. In Second State of the Carbon Cycle Report (SOCCR2): A Sustained Assessment Report

Central to life on Earth, carbon is essential to the molecular makeup of all living things and plays a key role in regulating global climate. To understand carbon’s role in these processes, researchers measure and evaluate carbon stocks and fluxes. A stock is the quantity of carbon contained in a pool or reservoir in the Earth system (e.g., carbon in forest trees), and a flux is the direction and rate of carbon’s transfer between pools (e.g., the movement of carbon from the atmosphere into forest trees during photosynthesis). This document, the Second State of the Carbon Cycle Report (SOCCR2), examines the patterns of carbon stocks and fluxes—collectively called the “carbon cycle.” Emphasis is given to these patterns in specific sectors (e.g., agriculture and energy) and ecosystems (e.g., forests and coastal waters) and to the response of the carbon cycle to human activity. The purpose of SOCCR2 is to assess the current state of the North American carbon cycle and to present recent advances in understanding the factors that influence it. Concentrating on North America—Canada, the United States, and Mexico—the report describes carbon cycling for air, land, inland waters (streams, rivers, lakes, and reservoirs), and coastal waters (see Figure ES.1, p. 23). The questions framing the publication A U.S. Carbon Cycle Science Plan (Michalak et al., 2011) inspired development of three slightly modified questions that guide SOCCR2’s content and focus on North America in a global context: How have natural processes and human actions affected the global carbon cycle on land, in the atmosphere, in the ocean and other aquatic systems, and at ecosystem interfaces (e.g., coastal, wetland, and urban-rural)? How have socioeconomic trends affected atmospheric levels of the primary carbon-containing gases, carbon dioxide (CO2) and methane (CH4)? How have species, ecosystems, natural resources, and human systems been impacted by increasing greenhouse gas (GHG) concentrations, associated changes in climate, and carbon management decisions and practices? SOCCR2 synthesizes the most recent understanding of carbon cycling in North America, assessing new carbon cycle findings and information, the state of knowledge regarding core methods used to study the carbon cycle, and future research needed to best inform carbon management and policy options. Focusing on scientific developments in the decade since the First State of the Carbon Cycle Report (SOCCR1; CCSP 2007), SOCCR2 summarizes the past, current, and projected state of carbon sources, sinks, and natural processes, as well as contributions by human activities. In addition to CO2 and CH4, the report sometimes discusses nitrous oxide (N2O), a GHG associated with activities and processes that affect fluxes of carbon gases.1 SOCCR2 also describes improvements in analysis tools; developments in decision support; and new insights into ecosystem carbon cycling, human causes of changes in the carbon cycle, and social science perspectives on carbon. Since publication of SOCCR1, coordinated research from agencies in the three North American countries has enabled innovative observational, analytical, and modeling capabilities to further advance understanding of the North American carbon cycle (see Appendix D: Carbon Measurement Approaches and Accounting Frameworks, p. 834). Some of the report’s main conclusions, based on the Key Findings of each chapter, are highlighted in Box ES.1, Main Findings of SOCCR2, p. 24.

Report

Alkalic-type epithermal gold deposit model

This report summarizes the primary characteristics of alkalic-type epithermal gold (Au) deposits and provides an updated descriptive model. These deposits, primarily of Mesozoic to Neogene age, are among the largest epithermal gold deposits in the world. Considered a subset of low-sulfidation epithermal deposits, they are spatially and genetically linked to small stocks or clusters of intrusions containing high alkali-element contents. Deposits occur as disseminations, breccia-fillings, and veins and may be spatially and genetically related to skarns and low-grade porphyry copper (Cu) or molybdenum (Mo) systems. Gold commonly occurs as native gold, precious metal tellurides, and as sub-micron gold in arsenian pyrite. Quartz, carbonate, fluorite, adularia, and vanadian muscovite/roscoelite are the most common gangue minerals. Alkalic-type gold deposits form in a variety of geological settings including continent-arc collision zones and back-arc or post-subduction rifts that are invariably characterized by a transition from convergent to extensional or transpressive tectonics. The geochemical compositions of alkaline igneous rocks spatially linked with these deposits span the alkaline-subalkaline transition. Their alkali enrichment may be masked by potassic alteration, but the unaltered or least altered rocks (1) have chondrite normalized patterns that are commonly light rare earth element (LREE) enriched, (2) are heavy rare earth element (HREE) depleted, and (3) have high large ion lithophile contents and variable enrichment of high-field strength elements. Radiogenic isotopes suggest a mantle derivation for the alkalic magmas but allow crustal contamination. Oxygen and hydrogen isotope compositions show that the fluids responsible for deposit formation are dominantly magmatic, although meteoric or other external fluids (seawater, evolved groundwater) also contributed to the ore-forming fluids responsible for these deposits. Carbon and sulfur isotope compositions in vein-hosted carbonates and sulfide gangue minerals, respectively, coincide with magmatic values, although a sedimentary source of carbon and sulfur is evident in several deposits. Deep-seated structures are critical for the upwelling of hydrous alkalic magmas and for focusing magmatic-hydrothermal fluids to the site of precious metal deposition. The source of gold, silver (Ag), tellurium (Te), vanadium (V), and fluorine (F) was probably the alkalic igneous rocks themselves, and the coexistence of native gold, gold tellurides, and roscoelite in several deposits is primarily a function of similar physicochemical conditions during deposition (for example, overlapping pH and oxygen fugacity ( f O2). Potential environmental impacts related to the mining and processing of alkalic-type epithermal gold deposits include acid mine drainage with high levels of metals, especially zinc (Zn), copper, lead (Pb), and arsenic. However, because alkalic-type gold deposits typically contain carbonates, which contribute calcium and magnesium ions that increase water hardness, aquatic life may be afforded some protection. Impacts vary widely as a function of host rocks, climate, topography, and mining methods. Geologic mapping to (1) highlight the distribution of potassic alteration; (2) define fault density and orientation of structures; (3) determine the distribution of alkaline rocks and hydrothermal breccias; and (4) identify uniquely colored gangue minerals, such as fluorite and roscoelite, will be critical to exploration and future discoveries. Geophysical techniques that identify potassium (K) anomalies (for example, radiometric and spectroscopic surveys), as well as magnetic, resistivity, aeromagnetic, and gravity surveys, may help locate zones of high-permeability that control advecting hydrothermal fluids. Geochemical surveys that include analyses for Au, Ag, barium, Te, K, F, V, Mo, and mercury, which are key elements in these deposits, should be undertaken along with the measurement of other pathfinder elements such as arsenic, bismuth, Cu, iron, nickel, Pb, antimony, selenium, and Zn.

Scientific Investigations Report

Earth's magnetic field complex: U.S. National activities during the Decade of Geopotential Field Research

The US geomagnetism community is supported by NASA, NOAA, USGS, NSF, DOD, and US universities. During the Decade of Geopotential Field Research, inaugurated in 1999 with the launch of the Danish satellite Ørsted on a US rocket, the US community has been involved in satellite mission development and analysis, instrument development, model development, and in the discovery and understanding of new processes with satellite magnetic signatures. The ESA Swarm mission has been a primary focus of the US community, with three US scientists on Swarm's Mission Advisory Group. Swarm will measure, for the first time, the E-W gradient of the magnetic field. One of us (T. Sabaka) is involved with the development of a Comprehensive inversion scheme as part of the SMART consortium. This effort is an outgrowth of the Comprehensive Model [1]. Swarm will also provide valuable observations for ionospheric specification and forecast. The geomagnetism group at NOAA (S. Maus, P. Alken and C. Manoj) has developed algorithms to estimate the strength of the eastward electric field (EEF). As the driver of the equatorial plasma fountain, the EEF is an important space weather parameter. ESA is considering the implementation of the EEF as a dedicated inversion chain in the Level-2 Facility. In 2006, NASA launched a minisatellite magnetometer constellation mission (ST-5) to test technologies and software. The ST-5 constellation featured the first along-track gradient measurements. NASA has also initiated efforts to study geomagnetism mission concepts after Swarm. One of the ideas under consideration is the systematic measurement of radial field gradients. Instrument development, and geomagnetic observatories, are also an integral part of the US effort. The past decade has seen significant advances in the development of a self-calibrating vector helium magnetometer, and in the automation of the US observatory network. Working in coordination with Intermagnet, the USGS Geomagnetism Program has made operational 1-second data acquisition at 13 of its magnetic observatories. The Program is also developing a realtime 1-minute and 1-hour Dst service. Within the past decade, US scientists have been leaders in the development of models that describe the global geomagnetic environment, including comprehensive models (the CM series), maps of the lithospheric field from satellite (MF-series), near surface maps of the lithospheric field (WDMAM-series), models of the thickness of the magnetic crust, the IGRF and World Magnetic Model series, ionospheric models such as the EEJM1, JVDM1, and the IRI, and data assimilation-based models (MoSST-series) that predict the future state of the geomagneic field.

Conference Paper

On the influence of biomass burning on the seasonal CO2 signal as observed at monitoring stations

We investigated the role of biomass burning in simulating the seasonal signal in both prognostic and diagnostic analyses. The prognostic analysis involved the High-Resolution Biosphere Model, a prognostic terrestrial biosphere model, and the coupled vegetation fire module, which together produce a prognostic data set of biomass burning. The diagnostic analysis involved the Simple Diagnostic Biosphere Model (SDBM) and the Hao and Liu [1994] diagnostic data set of biomass burning, which have been scaled to global 2 and 4 Pg C yr −1 , respectively. The monthly carbon exchange fields between the atmosphere and the biosphere with a spatial resolution of 0.5° × 0.5°, the seasonal atmosphere-ocean exchange fields, and the emissions from fossil fuels have been coupled to the three-dimensional atmospheric transport model TM2. We have chosen eight monitoring stations of the National Oceanic and Atmospheric Administration network to compare the predicted seasonal atmospheric CO 2 signals with those deduced from atmosphere-biosphere carbon exchange fluxes without any contribution from biomass burning. The prognostic analysis and the diagnostic analysis with global burning emissions of 4 Pg C yr −1 agree with respect to the change in the amplitude of the seasonal CO 2 concentration introduced through biomass burning. We find that the seasonal CO 2 signal at stations in higher northern latitudes (north of 30°N) is marginally influenced by biomass burning. For stations in tropical regions an increase in the CO 2 amplitude of more than 1 ppmv (up to 50% with respect to the observed trough to peak amplitude) has been calculated. Biomass burning at stations farther south accounts for an increase in the CO 2 amplitude of up to 59% (0.6 ppmv). A change in the phase of the seasonal CO 2 signal at tropical and southern stations has been shown to be strongly influenced by the onset of biomass burning in southern tropical Africa and America. Comparing simulated and observed seasonal CO 2 signals, we find higher discrepancies at southern tropical stations if biomass burning emissions are included. This is caused by the additional increase in the amplitude in the prognostic analysis and a phase shift in a diagnostic analysis. In contrast, at the northern tropical stations biomass burning tends to improve the estimates of the seasonal CO 2 signal in the prognostic analysis because of strengthening of the amplitude. Since the SDBM predicts the seasonal CO 2 signal reasonably well for the northern hemisphere tropical stations, no general improvement of the fit occurs if biomass burning emissions are considered.

Global Biogeochemical Cycles

The Cottage Lake aeromagnetic lineament: A possible onshore extension of the southern Whidbey Island fault, Washington

The northwest-striking southern Whidbey Island fault zone (SWIF) was mapped previously using borehole data and potential-field anomalies on Whidbey Island and marine seismic surveys beneath surrounding waterways. Abrupt subsidence at a coastal marsh on south-central Whidbey Island suggests that the SWIF experienced a MW 6.5 to 7.0 earthquake about 3000 years ago. Southeast of Whidbey Island, a hypothesized southeastward projection of the SWIF would make landfall between the cities of Seattle and Everett. As part of systematic, ongoing studies by the U.S. Geological Survey, University of Washington, and other earth science organizations to evaluate potentially active faults and other earth hazards throughout the Puget Lowland, we test this hypothesis using aeromagnetic, lidar, and borehole data. Linear, northwest-striking magnetic anomalies traversing the mainland region project southeastward toward the communities of Woodinville and Maltby, Washington. All of these magnetic anomalies are low in amplitude and best illuminated in residual magnetic fields. The most prominent of the residual magnetic anomalies extends at least 16 km, lies approximately on strike with the SWIF on Whidbey Island, and passes near Crystal and Cottage Lakes, about 27 km southeast of downtown Everett. In places, this magnetic anomaly is associated with topographic lineaments, but spectral analysis indicates that the source of the anomaly extends to depths greater than 2 km and cannot be explained entirely by topographic effects. The Alderwood #1 oil exploration well located on strike with the Cottage Lake aeromagnetic lineament shows evidence of deformation over a total depth range of 3000 m; some beds within this interval exhibit intense fracturing and shearing, although deformation within the well can only be constrained as post-early Oligocene and pre-Pleistocene. Boreholes acquired as part of a wastewater tunnel project show evidence of soil disturbance at locations where some topographic and aeromagnetic lineaments cross the tunnel alignment. Some of the disturbance is likely tectonic in origin, although other explanations are possible. Some of the soil disturbance demonstrably predates the 15-13 ka Fraser glaciation of the Puget Lowland; other samples have inconclusive ages and may be younger. Subtle scarps in Pleistocene surfaces are visible on high-resolution lidar topography at a number of locations along the Cottage Lake aeromagnetic lineament. Collectively, the scarps are parallel to the trend of the aeromagnetic lineament and extend a total distance of 18 km. In the field, scarps exhibit 1 to 5 m of north-side-up offset. The scarps provide targets for future paleoseismic trenching studies to test the hypothesis that they have a tectonic origin.

Washington

Rhenium

Rhenium is one of the rarest elements in Earth’s continental crust; its estimated average crustal abundance is less than 1 part per billion. Rhenium is a metal that has an extremely high melting point and a heat-stable crystalline structure. More than 80 percent of the rhenium consumed in the world is used in high-temperature superalloys, especially those used to make turbine blades for jet aircraft engines. Rhenium’s other major application is in platinum-rhenium catalysts used in petroleum refining. Rhenium rarely occurs as a native element or as its own sulfide mineral; most rhenium is present as a substitute for molybdenum in molybdenite. Annual world mine production of rhenium is about 50 metric tons. Nearly all primary rhenium production (that is, rhenium produced by mining rather than through recycling) is as a byproduct of copper mining, and about 80 percent of the rhenium obtained through mining is recovered from the flue dust produced during the roasting of molybdenite concentrates from porphyry copper deposits. Molybdenite in porphyry copper deposits can contain hundreds to several thousand grams per metric ton of rhenium, although the estimated rhenium grades of these deposits range from less than 0.1 gram per metric ton to about 0.6 gram per metric ton. Continental-arc porphyry copper-(molybdenum-gold) deposits supply most of the world’s rhenium production and have large inferred rhenium resources. Porphyry copper mines in Chile account for about 55 percent of the world’s mine production of rhenium; rhenium is also recovered from porphyry copper deposits in the United States, Armenia, Kazakhstan, Mexico, Peru, Russia, and Uzbekistan. Sediment-hosted strata-bound copper deposits in Kazakhstan (of the sandstone type) and in Poland (of the reduced-facies, or Kupferschiefer, type) account for most other rhenium produced by mining. These types of deposits also have large amounts of identified rhenium resources. The future supply of rhenium is likely to depend largely on the capacity of the specialized processing facilities needed to recover rhenium from molybdenite concentrates. The environmental consequences of rhenium recovery are closely linked to the consequences of mining large porphyry copper and strata-bound copper deposits; no additional environmental impact from recovery of rhenium from these deposits has been identified. No information is available regarding the potential toxic effects of rhenium on humans, partly because of the low natural abundance of rhenium.

Professional Paper

Quartz-pebble-conglomerate gold deposits

Quartz-pebble-conglomerate gold deposits represent the largest repository of gold on Earth, largely due to the deposits of the Witwatersrand Basin, which account for nearly 40 percent of the total gold produced throughout Earth’s history. This deposit type has had a controversial history in regards to genetic models. However, most researchers conclude that they are paleoplacer deposits that have been modified by metamorphism and hydrothermal fluid flow subsequent to initial sedimentation. The deposits are found exclusively within fault-bounded depositional basins. The periphery of these basins commonly consists of granite-greenstone terranes, classic hosts for lode gold that source the detrital material infilling the basin. The gold reefs are typically located along unconformities or, less commonly, at the top of sedimentary beds. Large quartz pebbles and heavy-mineral concentrates are found associated with the gold. Deposits that formed prior to the Great Oxidation Event (circa 2.4 giga-annum [Ga]) contain pyrite, whereas younger deposits contain iron oxides. Uranium minerals and hydrocarbons are also notable features of some deposits. Much of the gold in these types of deposits forms crystalline features that are the product of local remobilization. However, some gold grains preserve textures that are undoubtedly of detrital origin. Other heavy minerals, such as pyrite, contain growth banding that is truncated along broken margins, which indicates that they were transported into place as opposed to forming by in situ growth in a hydrothermal setting. The ore tailings associated with these deposits commonly contain uranium-rich minerals and sulfides. Oxidation of the sulfides releases sulfuric acid and mobilizes various metals into the environment. The neutralizing potential of the tailings is minimal, since carbonate minerals are rare. The continuity of the tabular ore bodies, such as those of the Witwatersrand Basin, has allowed these mines to be the deepest in the world. The extreme depths create engineering complications and safety issues for the miners, such as rock bursts as a result of pressure release. The richness of these deposits makes them a desirable exploration target. However, the likelihood of future discoveries is minimal. Small deposits found in the United States include those found at Nemo in the Black Hills of South Dakota and Deep Lake in the Sierra Madre of Wyoming.

Scientific Investigations Report

Evaluation of Copernicus DEM and comparison to the DEM used for Landsat collection-2 processing

Having highly accurate and reliable Digital Elevation Models (DEMs) of the Earth’s surface is critical to orthorectify Landsat imagery. Without such accuracy, pixel locations reported in the data are difficult to assure as accurate, especially in more mountainous landscapes, where the orthorectification process is the most challenging. To this end, the Landsat Calibration and Validation Team (Cal/Val) compared the Copernicus DEM (CopDEM) to the DEM that is currently used in Collection-2 processing (called “Collection-2 DEM”). NGS ground-surveyed and lidar-based ICESat-2 points were used, and the CopDEM shows improvement to be less than 1 m globally, except in Asia where the accuracy and resolution of the DEM were greater for the CopDEM compared to the Collection-2 DEM. Along with slightly improved accuracy, the CopDEM showed more consistent results globally due to its virtually seamless source and consistent creation methods throughout the dataset. While CopDEM is virtually seamless, having greater than 99% of their data coming from a single source (Tandem-X), there are significantly more voids in the higher elevations which were mostly filled with SRTM derivatives. The accuracy of the CopDEM fill imagery was also compared to the Collection-2 DEM and the results were very similar, showing that the choice of fill imagery used by CopDEM was appropriate. A qualitative assessment using terrain-corrected products processed with different DEMs and viewing them as anaglyphs to evaluate the DEMs proved useful for assessing orbital path co-registration. While the superiority of the CopDEM was not shown to be definitive by the qualitative method for many of the regions assessed, the CopDEM showed a clear advantage in Northern Russia, where the Collection-2 DEM uses some of the oldest and least accurate datasets in the compilation of the Collection-2 DEM. This paper presents results from the comparison study, along with the justification for proceeding with using the Copernicus DEM in future Landsat processing. As of this writing, the Copernicus DEM is planned to be used in Collection-3 processing, which is anticipated to be released no earlier than 2025.

Remote Sensing

Observing coastal wetland transitions using national land cover products

Over the coming century, climate change and sea-level rise are predicted to cause widespread change to coastal wetlands. Estuarine vegetated wetlands can adapt to sea-level rise through both vertical development (i.e., biophysical feedbacks and sedimentation) and upslope/horizontal migration. Quantifying changes to estuarine vegetated wetlands over time can help to inform current and future decisions regarding land management and resource stewardship. In this study, we show how coastal land cover maps readily available in the US can be used to assess and understand estuarine vegetated wetland changes. This assessment involves two steps: (1) identifying the net gain/loss of estuarine vegetated wetlands and (2) determining which land cover types contribute to the net gain/loss. From this information, we developed estuarine vegetated wetland change scenarios that evaluate whether estuarine vegetated wetland gain kept up with loss and whether the contribution was from: (1) estuarine vegetated wetland migration or tidal restoration; (2) land building (i.e., development); or (3) both. We assessed changes from 1996 to 2016 for: (1) the conterminous US; (2) each major US coastline; and (3) focal estuaries with the most change per coast. We found that the change scenario (1, 2, or 3) varied across coastlines. Moving forward, national coastal land cover programs can be informed by utilizing methodologies that leverage contemporary information for delineating the estuarine zone from upslope/adjacent wetlands. We highlight approaches that could be used to address this challenge and provide complementary information related to wetland condition changes.

Progress in Physical Geography: Earth and Environm

Forecast experiment: do temporal and spatial b value variations along the Calaveras fault portend M ≥ 4.0 earthquakes?

The power law distribution of earthquake magnitudes and frequencies is a fundamental scaling relationship used for forecasting. However, can its slope (b value) be used on individual faults as a stress indicator? Some have concluded that b values drop just before large shocks. Others suggested that temporally stable low b value zones identify future large-earthquake locations. This study assesses the frequency of b value anomalies portending M ≥ 4.0 shocks versus how often they do not. I investigated M ≥ 4.0 Calaveras fault earthquakes because there have been 25 over the 37-year duration of the instrumental catalog on the most active southern half of the fault. With that relatively large sample, I conducted retrospective time and space earthquake forecasts. I calculated temporal b value changes in 5-km-radius cylindrical volumes of crust that were significant at 90% confidence, but these changes were poor forecasters of M ≥ 4.0 earthquakes. M ≥ 4.0 events were as likely to happen at times of high b values as they were at low ones. However, I could not rule out a hypothesis that spatial b value anomalies portend M ≥ 4.0 events; of 20 M ≥ 4 shocks that could be studied, 6 to 8 (depending on calculation method) occurred where b values were significantly less than the spatial mean, 1 to 2 happened above the mean, and 10 to 13 occurred within 90% confidence intervals of the mean and were thus inconclusive. Thus spatial b value variation might be a useful forecast tool, but resolution is poor, even on seismically active faults.

California

A Bayesian network to predict coastal vulnerability to sea level rise

Sea level rise during the 21st century will have a wide range of effects on coastal environments, human development, and infrastructure in coastal areas. The broad range of complex factors influencing coastal systems contributes to large uncertainties in predicting long-term sea level rise impacts. Here we explore and demonstrate the capabilities of a Bayesian network (BN) to predict long-term shoreline change associated with sea level rise and make quantitative assessments of prediction uncertainty. A BN is used to define relationships between driving forces, geologic constraints, and coastal response for the U.S. Atlantic coast that include observations of local rates of relative sea level rise, wave height, tide range, geomorphic classification, coastal slope, and shoreline change rate. The BN is used to make probabilistic predictions of shoreline retreat in response to different future sea level rise rates. Results demonstrate that the probability of shoreline retreat increases with higher rates of sea level rise. Where more specific information is included, the probability of shoreline change increases in a number of cases, indicating more confident predictions. A hindcast evaluation of the BN indicates that the network correctly predicts 71% of the cases. Evaluation of the results using Brier skill and log likelihood ratio scores indicates that the network provides shoreline change predictions that are better than the prior probability. Shoreline change outcomes indicating stability (-1 < rate < 1 m/yr) or erosion (rate < -1 m/yr) tend to occur for two sets of input scenarios. Stable shoreline change rates occur mainly for low rates of relative sea level rise and occur in low-vulnerability geomorphic settings. Rates indicating erosion result for cases where the rate of relative sea level rise is high and moderate-to-high vulnerability geomorphic settings occur. In contrast, accretion (rate > 1 m/yr) was not well predicted. We find that BNs can assimilate important factors contributing to coastal change in response to sea level rise and can make quantitative, probabilistic predictions that can be applied to coastal management decisions. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface

Characterizing climate-change impacts on the 1.5-yr flood flow in selected basins across the United States: a probabilistic approach

The U.S. Geological Survey Precipitation-Runoff Modeling System (PRMS) model was applied to basins in 14 different hydroclimatic regions to determine the sensitivity and variability of the freshwater resources of the United States in the face of current climate-change projections. Rather than attempting to choose a most likely scenario from the results of the Intergovernmental Panel on Climate Change, an ensemble of climate simulations from five models under three emissions scenarios each was used to drive the basin models. Climate-change scenarios were generated for PRMS by modifying historical precipitation and temperature inputs; mean monthly climate change was derived by calculating changes in mean climates from current to various future decades in the ensemble of climate projections. Empirical orthogonal functions (EOFs) were fitted to the PRMS model output driven by the ensemble of climate projections and provided a basis for randomly (but representatively) generating realizations of hydrologic response to future climates. For each realization, the 1.5-yr flood was calculated to represent a flow important for sediment transport and channel geomorphology. The empirical probability density function (pdf) of the 1.5-yr flood was estimated using the results across the realizations for each basin. Of the 14 basins studied, 9 showed clear temporal shifts in the pdfs of the 1.5-yr flood projected into the twenty-first century. In the western United States, where the annual peak discharges are heavily influenced by snowmelt, three basins show at least a 10% increase in the 1.5-yr flood in the twenty-first century; the remaining two basins demonstrate increases in the 1.5-yr flood, but the temporal shifts in the pdfs and the percent changes are not as distinct. Four basins in the eastern Rockies/central United States show at least a 10% decrease in the 1.5-yr flood; the remaining two basins demonstrate decreases in the 1.5-yr flood, but the temporal shifts in the pdfs and the percent changes are not as distinct. Two basins in the eastern United States show at least a 10% decrease in the 1.5-yr flood; the remaining basin shows little or no change in the 1.5-yr flood.

Earth Interactions

Forecast model for moderate earthquakes near Parkfield, California

Earthquake instability models have possible application to earthquake forecasting because the models simulate both preseismic and coseismic changes of fault slip and ground deformation. In the forecast procedure proposed here, repeated measurements of preseismic fault slip and ground deformation constrain the values of model parameters. The early part of the model simulation corresponds to the available field data, and the subsequent part constitutes an estimate of future faulting and ground deformation. In particular, the time, location, and size of unstable faulting are estimates of the pending earthquake parameters. The forecast accuracy depends on the model realism and parameter resolution. The forecast procedure is applied to fault creep and trilateration data measured near Parkfield, California, where at least five magnitude 5.5 to 6 earthquakes have occurred regularly since 1881, the last in 1966. The quasi-static model consists of a flat vertical plane embedded in an elastic half space. Spacially variable fault slip of strike-slip sense is driven by an increasing regional shear stress but is impeded by a relatively strong patch of brittle, strain-softening fault. The field data are consistent with these approximate values of patch parameters: radius of 3 km, patch center 5 km deep and 8 km southeast of the 1966 epicenter, and maximum brittle strength of 26 bars. Fluctuations in the available field data prevent estimating the earthquake time with any more precision than use of the 21±8 year recurrence interval. However, the model may later give a more precise estimate of the earthquake time if the fault slip rate near the inferred patch increases before the earthquake, as predicted by the model.

Journal of Geophysical Research Solid Earth

Ground-water research in the U.S.A.

Ground-water reservoirs and the overlying unsaturated zone-collectively, the "subsurface"-have an enormous capacity to supply water to wells and useful plants, to store water to meet future needs for the same purposes, and, under suitable precautions, to accept wastes. This capacity can be exploited on a maximum scale, however, only on the basis of information one or more orders of magnitude greater than that available at present on the distribution, recoverability, and replenishability of subsurface water. Because usable water must be made available, and waste water must be disposed of, at costs of only a cent or a few cents per cubic meter, there is a critical need for research to devise methods of accomplishing these water-management tasks at reasonable cost. Among the chief target areas for research in subsurface hydrology are permeability distribution, including vertical permeability; prediction of the departure of the storage coefficient from the theoretically "instantaneous" property assumed in flow equations; theory of unsaturated flow based on fundamental soil characteristics that can be measured practicably; geochemical relations including the effects of injecting water of one composition into zones occupied by waters of different composition, generation of acid mine water, occurrence of saline water, and salt-fresh-water relations in coastal and other areas; prediction of the fate of wastes injected underground; geophysical techniques both surface and subsurface to extend, at low cost, information obtained by other means; and practical techniques of artificial recharge, especially through wells. ?? 1967.

Earth-Science Reviews