USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

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,666 records · Page 63Linked to original sources
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↗
USGS2022

U-Pb scheelite ages of tungsten and antimony mineralization in the Stibnite-Yellow Pine district, central Idaho

The Stibnite-Yellow Pine district contains the largest antimony resource in the United States, as well as significant gold, and is a historic producer of tungsten. Application of in situ laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) direct dating of scheelite from two Au-Sb-W ore deposits, Yellow Pine and Hangar Flats, yielded an older group of U-Pb ages in the range of 60.0 ± 2.8 to 57.0 ± 1.1 Ma and a younger U-Pb age for scheelite intergrown with stibnite of 47.4 ± 1.1 Ma. These in situ analyses were calibrated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) U-Pb lower intercept ages of two coarsely crystalline scheelite samples that yielded ages of 57.52 ± 0.22 and 56.62 ± 0.16 Ma. Scheelite of the latter age is of sufficient quality to serve as a primary reference material for LA-ICP-MS scheelite U-Pb geochronology. The group of older U-Pb scheelite ages agrees with 40 Ar/ 39 Ar ages of 56.9 ± 1.2 to 56.38 ± 0.54 Ma on adularia from Yellow Pine and Hangar Flats, whereas the younger U-Pb scheelite age is similar to an 40 Ar/ 39 Ar age of 46.00 ± 0.40 Ma on adularia from an epithermal gold-silver deposit in the adjacent Thunder Mountain caldera. Our results indicate that the main stage of tungsten mineralization occurred at ca. 57 Ma, whereas the main stage of antimony mineralization occurred at ca. 47 Ma—thereby providing first-time age constraints for antimony and tungsten mineralization in the Stibnite-Yellow Pine district.

Idaho↗
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↗
USGS2024

Deep structure of Siletzia in the Puget Lowland: Imaging an obducted plateau and accretionary thrust belt with potential fields

Detailed understanding of crustal components and tectonic history of forearcs is important due to their geological complexity and high seismic hazard. The principal component of the Cascadia forearc is Siletzia, a composite basaltic terrane of oceanic origin. Much is known about the lithology and age of the province. However, glacial sediments blanketing the Puget Lowland obscure its lateral extent and internal structure, hindering our ability to fully understand its tectonic history and its influence on modern deformation. In this study, we apply map-view interpretation and two-dimensional modeling of aeromagnetic and gravity data to the magnetically stratified Siletzia terrane revealing its internal structure and characterizing its eastern boundary. These analyses suggest the contact between Siletzia (Crescent Formation) and the Eocene accretionary prism trends northward under Lake Washington. North of Seattle, this boundary dips east where it crosses the Kingston arch, whereas south of Seattle the contact dips west where it crosses the Seattle uplift (SU). This westward dip is opposite the dip of the Eocene subduction interface, implying obduction of Siletzia upper crust at this southern location. Elongate pairs of high and low magnetic anomalies over the SU suggest imbrication of steeply-dipping, deeply rooted slices of Crescent Formation within Siletzia. We hypothesize these features result from duplication of Crescent Formation in an accretionary fold-thrust belt during the Eocene. The active Seattle fault divides this Eocene fold-thrust belt into two zones with different structural trends and opposite frontal ramp dips, suggesting the Seattle fault may have originated as a tear fault during accretion.

Washington↗
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↗
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↗
USGS1999

Characterization and simulation of the quantity and quality of water in the Highland Lakes, Texas, 1983-92

The Highland Lakes, located in central Texas, are a series of seven reservoirs on the Colorado River (Lake Buchanan, Inks Lake, Lake Lyndon B. Johnson, Lake Marble Falls, Lake Travis, Lake Austin, and Town Lake). The reservoirs provide hydroelectric power for the area. In addition, Lake Austin and Town Lake also provide the public water supply for the Austin metropolitan area. Saline water released from Natural Dam Salt Lake during 1987–89 caused increased concern among water managers that high-salinity water entering the Highland Lakes could result in waterquality problems, necessitating additional treatment of the water. The maximum dissolved solids concentrations for the reservoirs after the saline inflow were about two to three times the average concentrations before the inflow. The maximum concentrations of chloride and sulfate after the inflow were about three to five times the average concentrations before the inflow. The concentrations of dissolved solids, chloride, and sulfate in Lake Buchanan, Inks Lake, Lake Lyndon B. Johnson, and Lake Marble Falls were less than the concentrations of the applicable water-quality standards by the end of 1990. Concentrations of these constituents in Lake Travis, Lake Austin, and Town Lake did not decrease to previous levels, which were less than the concentrations of the applicable waterquality standards, until the end of 1991. Constituent concentrations for Lake Buchanan and Inks Lake; for Lake Lyndon B. Johnson and Lake Marble Falls; and for Lake Travis, Lake Austin, and Town Lake were similar because of the relative storage capacities and location of tributary inflows. From the initial increase in constituent concentrations in Lake Buchanan (summer 1987) in response to the saline inflow, the high-salinity water passed through the entire Highland Lakes in about 3.5 years. A mathematical mass-balance model was used to simulate the input and movement of highsalinity water through the Highland Lakes and to estimate monthly mean concentrations of dissolved solids, chloride, and sulfate for wet, average, and dry hydrologic conditions. The simulated median monthly concentrations during the 10-year simulation period for each reservoir generally are larger for the average condition than for the wet condition and generally are larger for the dry condition than for the average condition. The simulated concentrations of dissolved solids, chloride, and sulfate decreased to levels less than the concentrations of the applicable water-quality standards in about 2 to 5 years after the saline water inflow of 1987–89 was simulated for the three hydrologic conditions. Results from the simulations indicate that saline inflows to the Highland Lakes similar to those of the releases from Natural Dam Salt Lake during 1987–89 are unlikely to cause large increases in future concentrations of dissolved solids, chloride, and sulfate in the Highland Lakes. The results also indicate that high-salinity water will continue to be diluted as it is transported downstream through the Highland Lakes, even during extended dry periods.

Texas↗
USGS1999

Estimation of magnitude and frequency of floods for streams in Puerto Rico: New empirical models

Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.

Puerto Rico↗
USGS2024

Managing basin-scale carbon sequestration: A tragedy of the commons approach

The Tragedy of the Commons is a well studied problem in the literature of ecology, economics, and environmental policy which illustrates the deleterious consequences of managing common pool resources when individual and social incentives are misaligned. In this work, we apply a simple model of carbon sequestration in a deep saline aquifer by two neighboring geologic CO 2 storage (GCS) operators to begin investigating if a Tragedy of the Commons framework applies to GCS. Specifically, we consider the pressure space as a “commons” because the injection by each firm at its own well increases the downhole injection pressure at both wells. We assume that a firm will decrease its injection rate if the downhole pressure at its well exceeds a predefined maximum (i.e., exceeds the “pressure limit”). With this assumption in place, we find that the same injection flowrates are optimal for both wells, regardless of whether they are owned by the same firm or competing firms. This suggests that GCS may not be best represented by a pure Tragedy of the Commons framework under our initial assumptions. However, there could be economic incentives or contractual obligations that may result in either or both GCS operators being unwilling to reduce their injection rates. Thus, we conclude the conference paper with a discussion of future extensions of our approach that may demonstrate closer alignment with the Tragedy of the Commons, including explicit definitions of pore-space rights, firm uncertainty regarding the parameters of the Theis equation, and the potential role of unitization.

Conference Paper↗
USGS2025

Hydrogeology of unconsolidated and bedrock aquifers along the Salmon River, including Malone, Franklin County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, investigated the hydrogeology of the unconsolidated and bedrock aquifers along the Salmon River corridor in northern Franklin County, New York. The study area covers roughly 147 square miles and includes the village of Malone and parts of the Towns of Malone, Bellmont, Burke, Constable, Westville, Bangor, Duane, and Franklin, New York. Groundwater is the primary source for water supply within the study area. Eighty-three percent of all residents use public water supplied from two production wells that draw water from a thick, highly productive sand-and-gravel aquifer likely receiving induced infiltration from the Salmon River. Twenty-four percent of the 187 verified domestic wells in the study area outside of the production well service area boundary are screened in typically discontinuous deposits of stratified sand and sand and gravel. Seventy-six percent of the wells are completed in bedrock aquifers including the Potsdam Sandstone, metamorphic rocks, Theresa Formation, and unknown bedrock. Characterizing and understanding potential groundwater resources is critical for protecting the quality of the groundwater. The information in this report may be used to guide delineation of groundwater contributing areas, assess potential threats to aquifers from both point and nonpoint sources, respond to contamination from spills or leaks from underground storage facilities or other sources, and to support assessments for future development of municipal water supplies.

New York↗