Search USGS⌕ Search

SEARCH · Search USGS

Results for “One 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 757 records · Page 42Linked to original sources

Database for the Geologic Map of Upper Eocene to Holocene Volcanic and Related Rocks of the Cascade Range, Oregon

Since 1979, Earth scientists of the Geothermal Research Program of the U.S. Geological Survey have carried out multidisciplinary research in the Cascade Range. The goal of this research is to understand the geology, tectonics, and hydrology of the Cascades in order to characterize and quantify geothermal resource potential. A major goal of the program is compilation of a comprehensive geologic map of the entire Cascade Range that incorporates modern field studies and that has a unified and internally consistent explanation. This map is one of three in a series that shows Cascade Range geology by fitting published and unpublished mapping into a province-wide scheme of rock units distinguished by composition and age; map sheets of the Cascade Range in Washington (Smith, 1993) and California will complete the series. The complete series forms a guide to exploration and evaluation of the geothermal resources of the Cascade Range and will be useful for studies of volcano hazards, volcanology, and tectonics. This digital release contains all the information used to produce the geologic map published as U.S. Geological Survey Geologic Investigations Series I-2569 (Sherrod and Smith, 2000). The main component of this digital release is a geologic map database prepared using ArcInfo GIS. This release also contains files to view or print the geologic map and accompanying descriptive pamphlet from I-2569.

Data Series↗

Experimental investigations regarding the use of sand as an inhibitor of air convection in deep seismic boreholes

Tilt has been the nemesis of horizontal long period seismology since its inception. Modern horizontal long period seismometers with their long natural periods are incredibly sensitive to tilt. They can sense tilts smaller than 10 -11 radians. To most readers, this is just a very very small number, so we will begin with an example, which should help to illustrate just how small 10 -11 radians is. Suppose we have an absolutely rigid rod which is approximately 4170 kilometers long; this just happens to be the Rand McNally map scaled crow flight distance between Los Angeles and Boston. Tilting this rod 10 -11 radians corresponds to raising one end of the rod 0.0000417 meters. Alas, this is just another very very small number! However, this corresponds to slipping a little less than one third a sheet of ordinary copying paper under one end of this perfectly rigid rod. To clarify, we mean, take a sheet of paper just like the paper this report is printed on and split it a little less than one third in the thickness direction, then put it under the end of the 4170 kilometer long rod! This will tilt the rod 10 -11 radians. Real world seismometers are nowhere near the length of this rod. A KS-54000 is about two meters long. Tilting a rod only two meters long 10~n radians corresponds to moving one end of this rod a mere 0.00000000002 meters or 0.02 millimicrons. As one of the authors old math teachers used to say, "That's PDS" (PDS = Pretty Damn Small). Unfortunately, the long period seismologist does not have the luxury of ignoring PDS numbers when it suits him as the mathematician frequently does. He must live in the real world in which tilts this small create severe contamination of long period seismic data. At periods longer than 20 seconds, tilt noise contaminates the long period data from all instruments installed on or near the earth's surface. Many years of experimentation revealed that installing the sensors at depth in deep mines drastically reduced the level of tilt noise in long period data. However, low levels of tilt noise persisted even at great depth; this noise was caused by air convection in the vault in which the sensors were installed. Over the years, methods were developed to control the air motion with mechanical barriers (boxes) around the sensors and by stratifying (creating a situation in which the air temperature increases with height) the air in the vault near the seismometer. These methods decreased tilt noise in deep mines to very low levels. However, deep mines, that are economically and environmentally suitable and accessible to seismology, are not plentiful and are not evenly distributed over the earth's surface. Therefore, the borehole deployable Teledyne Geotech KS-36000 and later the KS-54000 sensor systems were developed to fulfill the need for instruments that could be installed at depth wherever high quality long period data was desired. Early in the development program, it became evident to the Teledyne Geotech personnel that air convection within the borehole was going to be a significant problem in KS deployments. Experimental and theoretical investigations conducted by Teledyne Geotech (see Douze and Sherwin, 1975, and Sherwin and Cook, 1976) produced a list of recommended installation procedures for reducing the effects of air convection. These procedures consisted of wrapping the sensor in a relatively thin layer of foam insulation, filling the free space volume in the vicinity of the centralizer-bail assembly with foam insulation, and the installation of styrofoam hole plugs immediately above the cable strain relief assembly at the top of the sensor package and at the top of the borehole. This technology has performed quite satisfactorily for over 20 years but evidence of tilt noise in the system output has persisted throughout the KS deployment program (the evidence was that the horizontal components were usually noisier than the vertical components) even in deep boreholes. Some deep borehole sites have been plagued by quite high levels of horizontal noise. Therefore, there has been a definite need for a new technique for controlling low level tilt noise in deep boreholes and the use of sand has been under consideration for several years. Figure 1 contains conceptual illustrations of both the conventional holelock installed KS sensor system and the same sensor installed in sand. This figure demonstrates the major differences between the two installation methods. The curved arrows in the borehole on the left in the figure denote possible air convection cells which are believed to be the source of tilt noise in some of the conventional installations. This air motion is eliminated in a sand installation by filling most of the free air volume surrounding the seismometer with sand as shown in the right hand portion of the figure. The sand actually performs two functions; it prevents air motion and provides a remarkably ridgid clamping of the seismometer in the borehole. This report presents the results of quantitative experimental investigations into the effectiveness of controlling low level air convection in seismic borehole installations with sand. The main body of the experimental effort consisted of installing two KS-540001 sensor systems in closely spaced shallow boreholes, allowing the sensors to reach equilibrium operation, and then pouring sand into both boreholes to observe any changes caused by pouring sand into the holes. The hypothesis of the experiment was that the sand would fill up the entire free air volume between the sensor package and the borehole walls thereby preventing movement of the air in the vicinity of the sensor package. The validity of this hypothesis had been qualitatively proven by earlier experiments at ASL and by the sand installations at the IRIS/ASL stations ANMO in 1995 and COLA in 1996. This experiment documents the degree of improved noise levels to be expected if KS instruments are installed in sand instead of in the conventional manner.

Open-File Report↗

Quantity and quality of seepage from two earthen basins used to store livestock waste in southern Minnesota during the first year of operation, 1997-98

Numerous earthen basins have been constructed in Minnesota for storage of livestock waste. Typically, these basins are excavated pits with partially above-grade, earth-walled embankments and compacted clay liners. Some have drain tile installed around them to prevent shallow ground and soil water to discharge into the basins. Environmental concerns associated with the waste include contamination of ground water by nitrogen compounds and pathogens. The U.S. Geological Survey, in cooperation with the MFC A (Minnesota Pollution Control Agency), studied the quantity and quality of seepage from two earthen basins used to store livestock waste in southern Minnesota during the first year of operation. One basin (site A), located at a small dairy farm, holds a manuresilage mixture, milkhouse wastewater, and local runoff; the other basin (site B), located at a large hog farm, holds a manure-water mixture from a nearby gestation barn. Monitoring systems were installed below compacted clay liners in portions of the sidewalls and bottoms of the basins to determine the quantity and quality of the seepage. Total seepage flow from the site A basin ranged from about 900 to 2,400 gal/d (gallons per day) except during April 1998 when the flow increased to about 4,200 gal/d. Seepage flow in areal units, which closely correlated with flow in gal/d, generally ranged from about 0.07 to 0.18 in./d (inches per day), which exceeded the recommended maximum design rate of 0.018 in./d established by the MPCA. Seepage flow commonly was greater through the sidewalls than through the bottom. Seepage from the site A basin (based on 11 samples each from the bottom and sidewall) had chloride concentrations of 220-350 mg/L (milligrams per liter); ammonium-N (nitrogen) concentrations of 2.40 mg/L or less (except for one concentration of 18.4 mg/L); nitrate-N concentrations of 5.24 mg/L or less; and organic-N concentrations of 6.97 mg/L or less. Ground water would be enriched in chloride and diluted in inorganic-N from mixing with basin seepage. Fecal Coliform bacteria, although abundant in the basin wastewater, were present in very small amounts in the seepage. Total seepage flow from the site B basin generally ranged from 400 to 2,200 gal/d except during 1-month and 3-month periods when the flow ranged from about 3,800 to 6,200 gal/d. Seepage flow in areal units generally ranged from about 0.025 to 0.15 in./d, and, as at the site A basin, exceeded the MPCA recommended maximum design rate of 0.018 in./d. Seepage flow in areal units generally correlated with the flow in gal/d except through the sidewalls when the basin was unfilled. Except during the first three months of the study, seepage flow was greater through the sidewalls than through the bottom. Seepage from the site B basin (based on 10 samples each from the bottom and sidewall) had chloride concentrations of 11 to 100 mg/L; ammonium-N concentrations of 2.58 mg/L or less; nitrate-N concentrations of 25.7 mg/L or less (except for one concentration of 146 mg/L); and organic-N concentrations of 0.92 mg/L or less. Nitrate-N concentrations in the seepage exceeded the U.S. Environmental Protection Agency (1996). MCL (maximum contaminant level) of 10 mg/L in 17 of 22 samples. Background ground-water quality, however, indicated that nitrate-N concentrations were greater than the MCL prior to operation of the basin. Fecal Coliform bacteria, as at the site A basin, were abundant in the basin wastewater, but not in the seepage.

Minnesota↗

Landsat 9 cross calibration under-fly of Landsat 8: Planning, and execution

During the early post-launch phase of the Landsat 9 mission, the Landsat 8 and 9 mission teams conducted a successful under-fly of Landsat 8 by Landsat 9, allowing for the near-simultaneous data collection of common Earth targets by the on-board sensors for cross-calibration. This effort, coordinated by the Landsat Calibration and Validation team, required contributions from various entities across National Aeronautics and Space Administration and U.S. Geological Survey such as Flight Dynamics, Systems, Mission Planning, and Flight Operations teams, beginning about 18 months prior to launch. Plans existed to allow this under-fly for any possible launch date of Landsat 9. This included 16 ascent plans and 16 data acquisition plans, one for every day of the Landsat orbital repeat period, with a minimum of 5 days of useful coverage overlap between the sensors. After the Landsat 9 launch, the plan executed, and led to the acquisition of over 2000 partial to full overlapping scene pairs. Although containing less than the expected number of scenes, this dataset was larger than previous Landsat mission under-fly efforts and more than sufficient for performing cross-calibration of the Landsat 8 and Landsat 9 sensors. The details of the planning process and execution of this under-fly are presented.

Remote Sensing↗

Rapid movement of frozen debris-lobes: implications for permafrost degradation and slope instability in the south-central Brooks Range, Alaska

We present the results of a reconnaissance investigation of unusual debris mass-movement features on permafrost slopes that pose a potential infrastructure hazard in the south-central Brooks Range, Alaska. For the purpose of this paper, we describe these features as frozen debris-lobes. We focus on the characterisation of frozen debris-lobes as indicators of various movement processes using ground-based surveys, remote sensing, field and laboratory measurements, and time-lapse observations of frozen debris-lobe systems along the Dalton Highway. Currently, some frozen debris-lobes exceed 100 m in width, 20 m in height and 1000 m in length. Our results indicate that frozen debris-lobes have responded to climate change by becoming increasingly active during the last decades, resulting in rapid downslope movement. Movement indicators observed in the field include toppling trees, slumps and scarps, detachment slides, striation marks on frozen sediment slabs, recently buried trees and other vegetation, mudflows, and large cracks in the lobe surface. The type and diversity of observed indicators suggest that the lobes likely consist of a frozen debris core, are subject to creep, and seasonally unfrozen surface sediment is transported in warm seasons by creep, slumping, viscous flow, blockfall and leaching of fines, and in cold seasons by creep and sliding of frozen sediment slabs. Ground-based measurements on one frozen debris-lobe over three years (2008–2010) revealed average movement rates of approximately 1 cm day −1 , which is substantially larger than rates measured in historic aerial photography from the 1950s to 1980s. We discuss how climate change may further influence frozen debris-lobe dynamics, potentially accelerating their movement. We highlight the potential direct hazard that one of the studied frozen debris-lobes may pose in the coming years and decades to the nearby Trans Alaska Pipeline System and the Dalton Highway, the main artery for transportation between Interior Alaska and the North Slope.

Alaska↗

U–Pb, Rb–Sr, and U-series isotope geochemistry of rocks and fracture minerals from the Chalk River Laboratories site, Grenville Province, Ontario, Canada

As part of the Geologic Waste Management Facility feasibility study, Atomic Energy of Canada Ltd. (AECL) is evaluating the suitability of the Chalk River Laboratories (CRL) site in Ontario, situated in crystalline rock of the southwestern Grenville Province, for the possible development of an underground repository for low- and intermediate-level nuclear waste. This paper presents petrographic and trace element analyses, U–Pb zircon dating results, and Rb–Sr, U–Pb and U-series isotopic analyses of gneissic drill core samples from the deep CRG-series characterization boreholes at the CRL site. The main rock types intersected in the boreholes include hornblende–biotite (±pyroxene) gneisses of granitic to granodioritic composition, leucocratic granitic gneisses with sparse mafic minerals, and garnet-bearing gneisses with variable amounts of biotite and/or hornblende. The trace element data for whole-rock samples plot in the fields of within-plate, syn-collision, and volcanic arc-type granites in discrimination diagrams used for the tectonic interpretation of granitic rocks. Zircons separated from biotite gneiss and metagranite samples yielded SHRIMP-RG U–Pb ages of 1472 ± 14 (2 σ ) and 1045 ± 6 Ma, respectively, in very good agreement with widespread Early Mesoproterozoic plutonic ages and Ottawan orogeny ages in the Central Gneiss Belt. The Rb–Sr, U–Pb, and Pb–Pb whole-rock errorchron apparent ages of most of the CRL gneiss samples are consistent with zircon U–Pb age and do not indicate substantial large-scale preferential element mobility during superimposed metamorphic and water/rock interaction processes. This may confirm the integrity of the rock mass, which is a positive attribute for a potential nuclear waste repository. Most 234 U/ 238 U activity ratios (AR) in whole rock samples are within errors of the secular equilibrium value of one, indicating that the rocks have not experienced any appreciable U loss or gain within the past 1 Ma. However, 234 U/ 238 U AR in fracture mineral samples collected down to borehole lengths of about 740 m deviate from the secular equilibrium value and 234 U/ 238 U model ages calculated for fracture mineral samples showing excess 234 U range from 593 to 1415 ka, thus providing evidence of fracture flow in the associated bedrock during the past 1.5 Ma. Rare earth element patterns are variable in fracture-filling calcites and Fe oxides/hydroxides but are similar to those observed in associated whole-rock samples. The observed Ce anomalies are very small ( Ce N / Ce N ∗ ≈ 1 "> CeN/CeN∗≈1 ), do not vary with depth, and, therefore, do not contain evidence that the studied fracture minerals precipitated from oxidizing waters at the conceptual depth of a repository.

Ontario↗

K, Rb, Sr and 87Sr 86Sr in rocks from the Mid-Indian Oceanic Ridge

The 87 Sr/ 86 Sr ratios of Mid-Indian Oceanic Ridge (MIOR) basalts are nearly identical (0.7032 to 0.7035), with the exception of one more highly radiogenic sample (0.7043). These values are consistently higher than the strontium isotopic ratios of the ridge basalts from Atlantic and Pacific Oceans, suggesting that the source of the MIOR basalts was depleted in alkalies more recently and/or to a lesser degree than the basalts from other oceans.

Earth and Planetary Science Letters↗

Ancient blue oaks reveal human impact on San Francisco Bay salinity

San Francisco Bay is one of the most important estuaries on the west coast of the Americas. Its water quality is controlled primarily by streamflow from the Sacramento and San Joaquin rivers. In fact, freshwater inflow from the Sacramento-San Joaquin Delta explains 86% of the salinity variability at the mouth of the San Francisco Bay estuary [ Peterson et al ., 1989]. The massive diversion of streamflow by the California State Water Project and the Central Valley Project, part of the largest manmade water control system on Earth [ Reisner , 1988], has raised salinity in the estuary on daily, seasonal, and annual timescales [ Nichols et al., 1986; Peterson et al., 1989].

California↗

U-Pb, Re-Os, and Ar/Ar geochronology of rare earth element (REE)-rich breccia pipes and associated host rocks from the Mesoproterozoic Pea Ridge Fe-REE-Au deposit, St. Francois Mountains, Missouri

Rare earth element (REE)-rich breccia pipes (600,000 t @ 12% rare earth oxides) are preserved along the margins of the 136-million metric ton (Mt) Pea Ridge magnetite-apatite deposit, within Mesoproterozoic (~1.47 Ga) volcanic-plutonic rocks of the St. Francois Mountains terrane in southeastern Missouri, United States. The breccia pipes cut the rhyolite-hosted magnetite deposit and contain clasts of nearly all local bedrock and mineralized lithologies. Grains of monazite and xenotime were extracted from breccia pipe samples for SHRIMP U-Pb geochronology; both minerals were also dated in one polished thin section. Monazite forms two morphologies: (1) matrix granular grains composed of numerous small (<50 μ m) crystallites intergrown with rare xenotime, thorite, apatite, and magnetite; and (2) coarse euhedral, glassy, bright-yellow grains similar to typical igneous or metamorphic monazite. Trace element abundances (including REE patterns) were determined on selected grains of monazite (both morphologies) and xenotime. Zircon grains from two samples of host rhyolite and two late felsic dikes collected underground at Pea Ridge were also dated. Additional geochronology done on breccia pipe minerals includes Re-Os on fine-grained molybdenite and 40 Ar/ 39 Ar on muscovite, biotite, and K-feldspar. Ages (±2 σ errors) obtained by SHRIMP U-Pb analysis are as follows: (1) zircon from the two host rhyolite samples have ages of 1473.6 ± 8.0 and 1472.7 ± 5.6 Ma; most zircon in late felsic dikes is interpreted as xenocrystic (age range ca. 1522–1455 Ma); a population of rare spongy zircon is likely of igneous origin and yields an age of 1441 ± 9 Ma; (2) pale-yellow granular monazite—1464.9 ± 3.3 Ma (no dated xenotime); (3) reddish matrix granular monazite—1462.0 ± 3.5 Ma and associated xenotime—1453 ± 11 Ma; (4) coarse glassy-yellow monazite—1464.8 ± 2.1, 1461.7 ± 3.7 Ma, with rims at 1447.2 ± 4.7 Ma; and (5) matrix monazite (in situ)—1464.1 ± 3.6 and 1454.6 ± 9.6 Ma, and matrix xenotime (in situ)—1468.0 ± 8.0 Ma. Two slightly older ages of cores are about 1478 Ma. The young age of rims on the coarse glassy monazite coincides with an Re-Os age of 1440.6 ± 9.2 Ma determined in this study for molybdenite intergrown with quartz and allanite, and with the age of monazite inclusions in apatite from the magnetite ore ( Neymark et al., 2016 ). A 40 Ar/ 39 Ar age of 1473 ± 1 Ma was obtained for muscovite from a breccia pipe sample. Geochronology and trace element geochemical data suggest that the granular matrix monazite and xenotime (in polygonal texture), and cores of coarse glassy monazite precipitated from hydrothermal fluids during breccia pipes formation at about 1465 Ma. The second episode of mineral growth at ca. 1443 Ma may be related to faulting and fluid flow that rebrecciated the pipes. The ca. 10-m.y. gap between the ages of host volcanic rocks and breccia pipe monazite and xenotime suggests that breccia pipe mineral formation cannot be related to the felsic magmatism represented by the rhyolitic volcanic rocks, and hence is linked to a different magmatic-hydrothermal system.

Missouri↗

Progress report on lithium-related geologic investigations in Bolivia

The September 1, 1981, Samoa Islands Region earthquake occurred at the extreme northern end of the Tonga arc in a region where the Pacific plate may be disjointed along a hinge fault. In the last 50 years, magnitude 7 or greater earthquakes have occurred in this region on the average of once every six years, but four 7+ events have now occurred within the last six years. The mainshock was preceded about two hours earlier by a foreshock that was used as a calibration event for the Joint Epicenter Determination relocation of the mainshock and nearby seismicity occurring within a period seven months prior to and one week after the mainshock. The foreshock, better-located events of the prior seismicity, and most aftershocks are concentrated in a group near the mainshock epicenter, but several more distant aftershocks suggest that the aftershock zone may have been as large as 125 km in length and trended about S35?E. Identification of depth phases from a full suite of broadband records gives source depths of 25-3km for the mainshock and 29.5?3 km for the foreshock using a JB earth model. Source parameters were determined for the mainshock utilizing WWSSN analog and GDSN digital data. The preferred fault plane solution based on P-wave first motion data is a south by southwesterly steeply dipping normal fault, remarkably similar to the mechanism reported by Johnson and Molnar (1972) for the nearby earthquake of April 20, 196B. A waveform inversion technique described by Sipkin (1982), when applied to long-period P waveforms, gives an 'average' point source solution for a purely deviatoric moment rate tensor at a preferred source depth of 22 km. Very similar results were obtained from long-period GDSN body-wave and mantle-wave data using a centroid-moment tensor inversion technique described in Dziewonski, and others (1981). Both techniques provide solutions very close to a double couple source with a south by southwesterly shallow-dipping normal fault mechanism. To obtain the scalar mantle wave moment, GDSN vertical and transverse records 20,000 see in length were processed as described by Buland and Taggart (1981). Averaging all the data from Rayleigh and Love waves yields an estimate of 3.8 x 10^27 dyne-cm (as compared to about 1.9 x 10^27 from body-wave moment tensor inversions) or a moment magnitude (Mr) of 7.6. For the portion of the waveform analysed (50-5B sec), the body-wave inversion performed by Sipkin gives a source time function of duration approximately 28 sec with two peaks in activity. Simultaneous analysis of short-period records, and broadband ground displacements and velocities, a method described by Harvey and Choy (1982) and Choy and Boatwright (1981) revealed a complex rupture consisting of two subevents, of about the same moment, separated in time by about 25 sec, and with durations of about 25 sec each. The two peaks in activity resolved by the body-wave moment tensor inversion correspond to the first of these subevents.

Open-File Report↗

2022 Crustal Deformation Modeling Workshop Report

The 2022 Crustal Deformation Modeling Workshop was held June 20–24 at the Colorado School of Mines in Golden, Colorado. The workshop included two days of tutorials on the use of the open-source software PyLith for crustal deformation modeling followed by three days of science talks and discussions. The workshop focused on three primary themes: (1) Earthquake cycle modeling; (2) Inversions for fault slip; and (3) Faulting, fluids, and surface loading. The talks highlighted how computational models provide insight into intriguing observations of Earth and planetary behavior. These included (1) earthquake synchronization of rupture patches due to their close proximity to each other, (2) the influence of fault geometry and damage zones on hypocenter depth and rupture propagation, (3) a lack of steady-state faulting behavior due to long time scales for grain size evolution in the mid-crust, (4) crustal deformation due to tidal, hydrological, and atmospheric loads, and (5) plumes of gas and icy particles due to tidal driven faulting on Enceladus (one of Saturn’s moons). The talks also described new computational modeling capabilities for incorporating complex geologic structure into Bayesian inversions for fault slip and efficient implementation of earthquake cycle models using a symmetric interior discontinuous Galerkin method. The complete agenda is available on the Computational Infrastructure for Geodynamic (CIG) website.

Report↗

A synthesis of the basal thermal state of the Greenland Ice Sheet

The basal thermal state of an ice sheet (frozen or thawed) is an important control upon its evolution, dynamics and response to external forcings. However, this state can only be observed directly within sparse boreholes or inferred conclusively from the presence of subglacial lakes. Here we synthesize spatially extensive inferences of the basal thermal state of the Greenland Ice Sheet to better constrain this state. Existing inferences include outputs from the eight thermomechanical ice-flow models included in the SeaRISE effort. New remote-sensing inferences of the basal thermal state are derived from Holocene radiostratigraphy, modern surface velocity and MODIS imagery. Both thermomechanical modeling and remote inferences generally agree that the Northeast Greenland Ice Stream and large portions of the southwestern ice-drainage systems are thawed at the bed, whereas the bed beneath the central ice divides, particularly their west-facing slopes, is frozen. Elsewhere, there is poor agreement regarding the basal thermal state. Both models and remote inferences rarely represent the borehole-observed basal thermal state accurately near NorthGRIP and DYE-3. This synthesis identifies a large portion of the Greenland Ice Sheet (about one third by area) where additional observations would most improve knowledge of its overall basal thermal state.

Journal of Geophysical Research F: Earth Surface↗

Advances in the simulation and automated measurement of well-sorted granular material: 2. Direct measures of particle properties

1. In this, the second of a pair of papers on the structure of well-sorted natural granular material (sediment), new methods are described for automated measurements from images of sediment, of: 1) particle-size standard deviation (arithmetic sorting) with and without apparent void fraction; and 2) mean particle size in material with void fraction. A variety of simulations of granular material are used for testing purposes, in addition to images of natural sediment. Simulations are also used to establish that the effects on automated particle sizing of grains visible through the interstices of the grains at the very surface of a granular material continue to a depth of approximately 4 grain diameters and that this is independent of mean particle size. Ensemble root-mean squared error between observed and estimated arithmetic sorting coefficients for 262 images of natural silts, sands and gravels (drawn from 8 populations) is 31%, which reduces to 27% if adjusted for bias (slope correction between observed and estimated values). These methods allow non-intrusive and fully automated measurements of surfaces of unconsolidated granular material. With no tunable parameters or empirically derived coefficients, they should be broadly universal in appropriate applications. However, empirical corrections may need to be applied for the most accurate results. Finally, analytical formulas are derived for the one-step pore-particle transition probability matrix, estimated from the image's autocorrelogram, from which void fraction of a section of granular material can be estimated directly. This model gives excellent predictions of bulk void fraction yet imperfect predictions of pore-particle transitions.

Journal of Geophysical Research F: Earth Surface↗

Overview and USGS activities

The discipline of cartography is undergoing a number of profound changes that center on the emerging influence of digital manipulation and analysis of data. Perhaps the most fundamental distinction between the digital representation of cartographic data and the conventional printed graphic is the need to explicitly and unambiguously code the attributes and spatial relationships among the various data elements. It is also necessary to follow acceptable practices for automated data processing. These requirements have led to the development by the USGS National Mapping Division of several documents that establish in-house digital cartographic standards. In an effort to fulfill lead agency requirements for promulgation of Federal standards in the earth sciences, the documents have been assembled with explanatory text into a USGS Circular. This Circular describes some of the pertinent issues relative to digital cartographic data standards, documents the digital cartographic data standards currently in use within the USGS, and details the efforts of the USGS related to the definition of national digital cartographic data standards. It consists of several chapters; the first is a general overview, and each succeeding chapter is made up from documents that establish in-house standards for one of the various types of digital cartographic data currently produced.

Circular↗

Manganese: it turns iron into steel (and does so much more)

Manganese is a common ferrous metal with atomic weight of 25 and the chemical symbol Mn. It constitutes roughly 0.1 percent of the Earth’s crust, making it the 12th most abundant element. Its early uses were limited largely to pigments and oxidants in chemical processes and experiments, but the significance of manganese to human societies exploded with the development of modern steelmaking technology in the 1860s. U.S consumption of manganese is about 500,000 metric tons each year, predominantly by the steel industry. Because manganese is essential and irreplaceable in steelmaking and its global mining industry is dominated by just a few nations, it is considered one of the most critical mineral commodities for the United States.

Fact Sheet↗

Louisiana and Landsat

Louisiana holds a unique historical, economic, and cultural position in the national consciousness. Its off-shore oil operations help fuel the U.S. economy. The Port of South Louisiana is the busiest in the United States by cargo volume; the nearby Port of New Orleans is the sixth busiest. The former French and Spanish colony served as a key connection to the Caribbean long before U.S. independence, and Louisiana’s multinational effects soon melded into a Creole culture that had an outsized effect on America. That heritage remains a powerful draw for the tourism industry in Louisiana. Gulf Coast breezes carry the aromas of tropical flowers, sweet beignets, and savory crawfish through the 13 colorful blocks of New Orleans’ French Quarter. Interwoven with the sounds of jazz, rock, country music, and zydeco, the city’s charms delight more than 18 million visitors each year. The proximity of the Gulf Coast and the city’s elevation, however—just 6.5 feet above sea level—also offer an ominous warning of the ever-present threat of climate change and natural disaster. Hurricane Katrina battered New Orleans in 2005, an incident tied to more than 1,800 deaths that marks one of the most notorious U.S. weather-related tragedies in the 21st century. Environmental changes have amplified threats from tropical storms. Through more frequent and powerful storms, sea level rise threatens low-lying areas such as Lake Charles and creates unpredictable weather patterns that threaten the cities and agricultural operations to the north. Landsat data offer rich information that can aid in early warning, disaster response, and the monitoring of recovery from natural disasters. Its historic, unparalleled 50-year archive of repeat Earth observations also serves to guide resiliency plans and feeds modeling that can help States like Louisiana prepare for coming coastal and inland change. Here are just a few examples of how Landsat has been used to study and understand Louisiana.

Louisiana↗

hypoDD-A Program to Compute Double-Difference Hypocenter Locations

HypoDD is a Fortran computer program package for relocating earthquakes with the double-difference algorithm of Waldhauser and Ellsworth (2000). This document provides a brief introduction into how to run and use the programs ph2dt and hypoDD to compute double-difference (DD) hypocenter locations. It gives a short overview of the DD technique, discusses the data preprocessing using ph2dt, and leads through the earthquake relocation process using hypoDD. The appendices include the reference manuals for the two programs and a short description of auxiliary programs and example data. Some minor subroutines are presently in the c language, and future releases will be in c. Earthquake location algorithms are usually based on some form of Geiger’s method, the linearization of the travel time equation in a first order Taylor series that relates the difference between the observed and predicted travel time to unknown adjustments in the hypocentral coordinates through the partial derivatives of travel time with respect to the unknowns. Earthquakes can be located individually with this algorithm, or jointly when other unknowns link together the solutions to indivdual earthquakes, such as station corrections in the joint hypocenter determination (JHD) method, or the earth model in seismic tomography. The DD technique (described in detail in Waldhauser and Ellsworth, 2000) takes advantage of the fact that if the hypocentral separation between two earthquakes is small compared to the event-station distance and the scale length of velocity heterogeneity, then the ray paths between the source region and a common station are similar along almost the entire ray path (Fréchet, 1985; Got et al., 1994). In this case, the difference in travel times for two events observed at one station can be attributed to the spatial offset between the events with high accuracy. DD equations are built by differencing Geiger’s equation for earthquake location. In this way, the residual between observed and calculated travel-time difference (or double-difference) between two events at a common station are a related to adjustments in the relative position of the hypocenters and origin times through the partial derivatives of the travel times for each event with respect to the unknown. HypoDD calculates travel times in a layered velocity model (where velocity depends only on depth) for the current hypocenters at the station where the phase was recorded. The double-difference residuals for pairs of earthquakes at each station are minimized by weighted least squares using the method of singular value decomposition (SVD) or the conjugate gradients method (LSQR, Paige and Saunders, 1982). Solutions are found by iteratively adjusting the vector difference between nearby hypocentral pairs, with the locations and partial derivatives being updated after each iteration. Details about the algorithm can be found in Waldhauser and Ellsworth (2000). When the earthquake location problem is linearized using the double-difference equations, the common mode errors cancel, principally those related to the receiver-side structure. Thus we avoid the need for station corrections or high-accuracy of predicted travel times for the portion of the raypath that lies outside the focal volume. This approach is especially useful in regions with a dense distribution of seismicity, i.e. where distances between neighboring events are only a few hundred meters. The improvement of double-difference locations over ordinary JHD locations is shown in Figure 1 for about 10,000 earthquakes that occurred during the 1997 seismic crisis in the Long Valley caldera, California. While the JHD locations (left panel) show a diffuse picture of the seismicity, double-difference locations (right panel) bring structural details such as the location of active fault planes into sharp focus.

Open-File Report↗

Mineral Resource of the Month: Bromine

Bromine, along with mercury, is one of only two elements that are liquid at room temperature. Bromine is a highly volatile and corrosive reddish-brown liquid that evaporates easily and converts to a metal at extreme pressures — above about 540,000 times atmospheric pressure. Bromine occurs in seawater, evaporitic (salt) lakes and underground brines associated with petroleum deposits.

Earth↗