USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal of Coal Geology”

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.

1,660 records · Page 33Linked to original sources
USGS1976

Diamond drilling at the Ma'milah Mine, Kingdom of Saudi Arabia

The Ma'millah gold mine is about 90 km southeast of At Ta'if at lat 21°03'N., long 41°18'E., in southwestern Saudi Arabia. The deposit, which was worked extensively by ancient miners, consists of several veins and massive lenses of quartz along a regional fault. The extensive ancient mine workings, large size of some quartz outcrops, and significant gold content of several samples indicated that subsurface exploration was warranted. Accordingly, two holes were drilled into the deposit during the interval December 1972 to February 1973. The drill findings were disappointing; the grade of vein material was found to be too low to be of economic interest and the quartz bodies were found to pinch out both along strike and downdip. The drilling program therefore was terminated. The deposit is not worth further exploration

Ma'milah Mine↗
USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗
USGS1995

Water-quality assessment of the Ozark Plateaus study unit, Arkansas, Kansas, Missouri, and Oklahoma: Analysis of information on nutrients, suspended sediment, and suspended solids, 1970-92

Water-quality data collected during water years 1970-90 (October 1 to September 30) for 83 surface-water sites and during 1970-92 for 395 ground-water sites in the 48,000 square mile Ozark Plateaus study unit of the National Water Quality Assessment Program were analyzed using selected descriptive and statistical methods. The water- quality data include nutrient (nitrogen and phosphorus), suspended sediment, and suspended- solids data, and ancillary information such as fertilizer use, animal waste, sewage treatment plant, and land use. Statistically significant differences exist in surface-water quality that can be attributed to physiography, land use, and other effects. The sites that were considered to be substantially affected by sewage treatment plants had the largest concentrations of nutrients. Nutrient concentrations generally were larger at sites associated with agricultural basins than at sites associated with forested basins. Statistically significant differences existed in the quality of ground water that can be attributed to hydrogeologic and land-use effects. Nutrient concentrations generally were largest where the water source is indicated to be shallow in origin and where parts of the hydrogeologic units are in agricultural land-use areas. Water quality has changed at several surface-water sites since 1970. Nutrient concentrations appear to have increased at some sites and decreased at other sites. Causes of these apparent trends are not known, but many of the sites with apparent trends are in agricultural areas. Surface-water loads of nutrients and suspended sediment were affected by several factors including streamflow, climate, drainage area, reservoir operation, and inputs from point and nonpoint sources. Annual loads were largest in large basins, with large inputs of nutrients or sediment during periods of high streamflows at locations where reservoir operation effects are not substantial.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗
USGS1976

Guidelines for collection and field analysis of ground-water samples for selected unstable constituents

The unstable nature of many chemical and physical constituents in ground water requires special collection procedures and field analysis immediately after collection. This report describes the techniques and equipment commonly used in the collection and field analysis of samples for pH, temperature, carbonate, bicarbonate, specific conductance, Eh, and dissolved oxygen.

Open-File Report↗
USGS2024

Paleoseismology of the Sawtooth fault and implications for fault behavior in the epicentral region of the 2020 Mw 6.5 Stanley, Idaho, earthquake

The 2020 moment magnitude ( M w ⁠⁠ ) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans‐Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1δ ⁠ ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault ( M w ⁠ 6.8–7.4), (2) that shorter ruptures ( M w ≤ ⁠⁠6.9 ) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast‐trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.

Idaho↗
USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗
USGS2026

High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys. Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Nevada↗
USGS1983

Technique for estimating depth of floods in Tennessee

Estimates of flood depths are needed for design of roadways across flood plains and for other types of construction along streams. Equations for estimating flood depths in Tennessee were derived using data for 150 gaging stations. The equations are based on drainage basin size and can be used to estimate depths of the 10-year and 100-year floods for four hydrologic areas. A method also was developed for estimating depth of floods having recurrence intervals between 10 and 100 years. Standard errors range from 22 to 30 percent for the 10-year depth equations and from 23 to 30 percent for the 100-year depth equations.

Tennessee↗
USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗
USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗
USGS2025

Apatite geo-thermochronology and geochemistry constrain Oligocene-Miocene growth and geodynamics of the northeastern Tibetan Plateau

Understanding the geodynamics of plateau evolution requires examining the spatial and temporal aspects of mountain building in the northeastern Tibetan Plateau, which are still under debate. Here we integrate apatite geo-thermochronological and geochemical data from the Oligocene-Miocene succession of the Xunhua Basin to elucidate the evolution of the regional topography. The results suggest sediment provenance changes, at ca. 28, 20 and 12 Ma, indicating topographic growth of the West Qinling, Laji Shan, and Jishi Shan source areas, respectively. Our findings reveal Oligocene-Miocene stepwise deformation and middle Miocene stress reorganization within the northeastern Tibetan Plateau. We attribute this to Oligocene lithospheric foundering beneath the south-central Tibetan Plateau and the middle Miocene initiation of sinistral strike-slip faults due to the underthrusting of India and North China. This study highlights the roles of lithospheric removal and continental underthrusting in plateau growth, and the value of multi-proxy apatite analyses in provenance studies.

Northeastern Tibetan Plateau↗
USGS2025

Geochemistry and thermometry of magnetite veins and replacements in iron ore deposits from the Iron Springs district, SW Utah, USA: Relation to magmatic and hydrothermal processes

Ore-genesis, texture, chemical compositions, and thermometry of magnetite were used to better evaluate the origins of iron ore deposits across the Iron Spring district, Utah. Trace element compositions, determined by LA-ICP-MS, were used to calculate temperatures based on X Mg numbers and characterize the trace element signatures of magmatic and hydrothermal magnetite. Magnetite occurs as coarser crystalline grains within veins, breccias, and cavities along with finer crystalline replacements. Magnetite is most commonly intergrown with apatite, phlogopite, and carbonate in veins, breccias, and cavities where apatite occurs separately from the other gangue minerals. The results using empirical trace element discrimination diagrams generally suggest that the magnetite is overall of hydrothermal origin regardless of textural variations or field relationships. The X Mg numbers of the magnetite in the magnetite-apatite ores, however, are suggestive of high temperature magmatic conditions; whereas those of vein, breccia, and replacement occurrences containing other gangue minerals are suggestive of lower temperature hydrothermal processes.

Utah↗
USGS2025

Geochemical evidence for the origin of late Quaternary loess, Seward Peninsula, Alaska

Loess is the most widespread surficial deposit in the state of Alaska. Although loess of last glacial age is common in mid-continental North America, records of last glacial loess in Alaska have been elusive. Here we report a record of last glacial loess on the Seward Peninsula, along with studies of a transect of loess localities also thought to date to the last glacial period. Highly diverse hypotheses on the origin of loess on the Seward Peninsula include local volcanic rocks and mountain ranges to the north and south, as well as offshore deposits. We present particle size and mineralogical and geochemical data that indicate that the most likely sources for loess in this part of eastern Beringia are glaciogenic silt particles of the Kobuk River and Noatak River, both of which drain the Brooks Range to the northeast. Northeasterly paleowinds were therefore responsible for loess on the Seward Peninsula. Such winds probably occurred during late summer or early autumn, when river discharge was at a minimum but temperatures were still above freezing. Northeasterly winds during the last glacial period are in good agreement with a recent paleoclimate model.

Alaska↗
USGS2024

The impact of source time function complexity on stress drop estimates

Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time‐ and frequency‐domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016 ) source time function (STF) database and analyze how rupture complexity drives stress‐drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress‐drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.

Bulletin of the Seismological Society of America↗
USGS2026

Unraveling protracted modification of Archean and Paleoproterozoic crust in central Laurentia, Penokean orogen, with garnet and accessory mineral geochronology and microstructural analysis

Proterozoic metamorphism and deformation of the southern margin of the Superior craton in the Lake Superior region is attributed to the Penokean orogeny (1890−1830 Ma). This model includes a period of crustal inversion in which Archean basement blocks were exhumed through overlying Paleoproterozoic strata, producing the corridor of gneiss domes that parallels the trend of the Penokean orogen across the northern Midcontinent, USA. However, recent geologic mapping and 40 Ar/ 39 Ar geochronology challenge this interpretation, suggesting instead that the gneiss dome structures reflect younger episodes of tectonic activity along the southern margin of Laurentia. In absence of integrated pressure-temperature-time-deformation constraints for these rocks, interpretations are largely limited to their final cooling history, making it difficult to both identify the tectonic forces that shaped the architecture of the Penokean orogenic belt and assess the extent to which later Proterozoic tectonism modified the southern Superior craton. We address this problem with an approach joining thermodynamic modeling, garnet and accessory mineral geochronology, and microstructural analysis for several metamorphic rocks across the gneiss dome corridor. The U-Pb ages of titanite reveal that the Proterozoic geometries of exhumed basement gneiss domes are governed by preexisting Archean structures. Garnet Lu-Hf geochronology constrains the timing of prograde-to-peak metamorphism in the Penokean orogenic belt. Granulite facies metamorphism is related to the final stages of the Penokean orogeny at 1837 Ma and localized in a belt of high-grade rocks near a major Penokean suture. Garnet Lu-Hf ages of samples adjacent to gneiss domes reflect regional metamorphism following the accretionary phase of the Penokean orogeny, between 1825 Ma and 1782 Ma, which we suggest reflects continued crustal thickening related to convergence farther south during this time interval. Combination of garnet microstructures and Sm-Nd ages reflects later exhumation of gneiss domes and buried metasedimentary rocks by ca. 1750 Ma, consistent with previously published 40 Ar/ 39 Ar cooling ages across the region. Reset Lu-Hf and Sm-Nd garnet ages and U-Pb ages of syn-kinematic titanite reflect reactivation of primary Penokean structures during this period of basement uplift. These data document significant modification of the Penokean orogen and the Archean crust of the southern Superior province between 1800 Ma and 1700 Ma. Tectonic activity during this interval coincides with collisional events recognized in western Laurentia, suggesting that the period immediately following the Penokean orogeny may be a broadly important time for crustal growth and modification in proto-North America.

Michigan, Wisconsin↗