Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Interpretation Research”

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 145 records · Page 8Linked to original sources
USGS2003

Location of long-period events below Kilauea Volcano using seismic amplitudes and accurate relative relocation

We present methods for improving the location of long-period (LP) events, deep and shallow, recorded below Kilauea Volcano by the permanent seismic network. LP events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport. However, it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals. At Kilauea, similar LP waveform signatures suggest the existence of LP multiplets. The waveform similarity suggests spatially close sources, while catalog solutions using arrival time estimates are widely scattered beneath Kilauea's summit caldera. In order to improve estimates of absolute LP location, we use the distribution of seismic amplitudes corrected for station site effects. The decay of the amplitude as a function of hypocentral distance is used for inferring LP location. In a second stage, we use the similarity of the events to calculate their relative positions. The analysis of the entire LP seismicity recorded between January 1997 and December 1999 suggests that a very large part of the LP event population, both deep and shallow, is generated by a small number of compact sources. Deep events are systematically composed of a weak high-frequency onset followed by a low-frequency wave train. Aligning the low-frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated. This observation favors an interpretation in terms of triggering and resonance of a magmatic conduit. Instead of defining fault planes, the precise relocation of similar LP events, based on the alignment of the high-energy low-frequency wave trains, defines limited size volumes. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Crustal structure of accreted terranes in southern Alaska, Chugach Mountains and Copper River Basin, from seismic refraction results

Seismic refraction data were collected along a 320-km-long "transect' line in southern Alaska, crossing the Prince William, Chugach, Peninsular, and Wrangellia terranes, and along several shorter lines within individual terranes. Velocity structure in the upper crust (less than 9-km depth) differs among the four terranes. In contrast, layers in the middle crust (9- to 25-km depth) in some case extend across projected terrane boundaries. The top of a gently north dipping sequence of low- and high-velocity layers (5.7-7.8 km/s), more than 10 km thick, extends from near the surface in the southern Chugach terrane to more than 20-km depth beneath the southern Peninsular terrane. This sequence, truncated by the suture between the Prince William and Chugach terranes, is interpreted to be an underplated "terrane' made up of fragments of the Kula plate and its sedimentary overburden that were accreted during subduction in the late Mesozoic and/or early Tertiary, during or between times of accretion of the Prince William and Chugach terranes.

Alaska↗
USGS2012

Geology of the Mount St. Helens area: Record of discontinuous volcanic and plutonic activity in the Cascade Arc of southern Washington

The Quaternary edifice of Mount St. Helens volcano was built upon a deeply eroded terrane of gently folded and altered volcanic and plutonic rocks that represent the core of the Tertiary Cascade magmatic arc. These rocks constitute an east dipping homoclinal sequence, several kilometers thick, of subaerially erupted mafic to silicic flows and volcaniclastic strata; K‐Ar ages from this section range from about 28 to 23 Ma (late Oligocene and earliest Miocene), which corresponds to an apparent lull in Cascade volcanism to the north of Mount Rainier. Volcanism was essentially continuous during this period of time, and neither a well‐defined base nor top of the section is exposed within the mapped area. Basalt and basaltic andesite dominate the lower part of the mid‐Tertiary section, whereas andesitic and dacitic rocks comprise most of the upper part. This section was intruded by numerous mafic to silicic dikes, sills, and irregular plutonic bodies, most no more than a few million years younger than their host rocks, and subjected to pervasive burial metamorphism and widespread hydrothermal alteration. Large areas of hornfelsed rock surrounding even relatively small intrusions indicate that the proportion of plutonic rock becomes significantly greater at shallow depth beneath the existing erosion surface. A large granitic pluton intruded the mid‐Tertiary section north of Spirit Lake at about 21 Ma. The Earl porphyry copper deposit occurs within the pluton but appears too young (17 Ma) to be genetically related to it. In contrast to the rather continuous and voluminous Oligocene to early Miocene activity, volcanism since then in the Mount St. Helens area has been localized and Volumetrically minor. Products of three younger eruptive periods have been recognized: a sequence of 15 m.y. old pyroxene andesite flows resting unconformably on mid‐Tertiary strata south of Mount St. Helens, widespread shallow dikes and sills of pyroxene andesite between 10 and 8 m.y. old, and compositionally diverse rocks erupted during the past 3 m.y. The Quaternary lavas are more potassic than the Tertiary lavas and typically contain phenocrysts of hornblende and biotite, which are absent from the older rocks. A number of Tertiary structures define a broad NNE trending zone that may reflect a deep‐seated lithospheric flaw that has controlled the locus of Cascade magmatism in southern Washington for the past 25 m.y. Mount St. Helens lies within this zone at the intersection of the NNW striking St. Helens seismic zone (SHZ) and an ENE trending alignment of Pleistocene silicic plug‐domes. No surface breakage has been detected along the SHZ, which is apparently very young. The linear zone of silicic vents is probably controlled by a fault that has been interpreted from seismic records to occur directly beneath the volcano. This zone parallels the directions of regional maximum horizontal compressive stress and North America/Juan de Fuca plate convergence. Mount St. Helens is an example of a low‐volume tectonically controlled magmatic system in an early stage of development.

Washington↗
USGS1978

Morphology of chasma walls, Mars

The landforms developed on the walls of the Valles Marineris system of chasmas are of three major types, which are locally transitional. The most common type is composed of steep spurs and gullies. The dominant process in the formation or modification of this type appears to be the downslope movement of. material under the influence of gravity, resulting in the accumulation of extensive talus deposits. The type is morphologically similar to high, steep terrestrial scarps in desert or alpine environments. The second morphologic type consists of walls dissected by tributary canyons with characteristic V-shaped cross profiles and blunt canyon heads that locally contain lobate deposits. The tributary canyons may be relict features of the time, when the existence of running water was possible on the surface of Mars. The third morphologic type consists of landslide scars forming broad curved or straight recessed sections of chasma wall. This type is accompanied by landslide deposits that form hummocky floors at the base of the recessed sections. The landslides developed at the expense of other wall morphologies. Chains of rimless depressions and craters that parallel the main structural trends of the chasmas are best interpreted as collapse holes. The origin of the chasmas on Mars is conjectural and may have been structural (grabens), but, on the basis of morphologic studies of their walls, it is suggested that most of the present wall configuration is the result of erosional scarp retreat, where erosion follows preestablished structural planes of weakness.

Journal of Research of the U.S. Geological Survey↗
USGS1977

Chemical structure of humic acids - Part 2, the molecular aggregation of some humic acid fractions in N, N-dimethylformamide

Humic acid fractions form molecular aggregates in solution. In previous studies we have shown by small angle X-ray scattering that the size of these aggregates is a function of pH. In this study we have found that the size of the aggregates of two humic acid fractions in water and buffers and in dimethylformamide solutions can be changed by oxidation with molecular oxygen and air. These results cast new light on the bonding mechanisms that cause aggregation of the humic acid particles in solution. We have interpreted the changes in aggregation sizes as being brought about by changes in intermolecular and intramolecular hydrogen bonding of the humic particles. Solvation of the humic molecules by dimethylformamide interferes with some of the hydrogen bonding reactions between proton donor and acceptor groups on the same humic acid molecules or on different molecules.

Journal of Research of the U.S. Geological Survey↗
USGS1978

Potassium-argon ages of basement rocks from Saint George Island, Alaska

Potassium-argon ages from basement granitic rocks that intrude serpentinized peridotite on St. George Island, Alaska, range from 50 to 57 million years, with an initial argon isochron age. interpreted as a minimum figure, of 52 ± 2 m.y. The age of the granitic rocks and their association with serpentinized peridotite near the continental margin, which apparently marks the site of a Mesozoic subduction zone, suggest that the position of the subduction zone shifted southward to the Aleutian Trench in early Tertiary time.

Alaska↗
USGS1980

Summaries of data on and lists of references to metallic and selected nonmetallic mineral deposits in the McCarthy Quadrangle, Alaska

These summaries of data on metallic and selected nonmetallic mineral occurrences and lists of selected references to .them in Geological Survey, U.S. Bureau of Mines, and State of Alaska Division of Geological and Geophysical Surveys (and predecessor State and Territorial agencies) reports and maps and a very.few journal articles are designed to aid in library research on the mineral resources of the McCarthy quadrangle, Alaska. The references listed are selected in the sense that mainly statistical reports such as the annual Minerals Yearbook of the U.S. Bureau of Mines and many annual and biennial reports of the Alaska Division of Geological and Geophysical Surveys and its predecessor agencies are not included. Also not included are data on claims about which little more than their locations is known (for example, localities 126- 130 in MacKevett and Holloway, 1977 (OF 77-169A). These omissions should not be interpreted as a judgement that the claims are not on valid mineral occurrences, but only chat there are insufficient data to describe any mineral deposits that might be present. Geochemical anomalies determined by analyses of rock and stream-sediment samples in which no metallic mineral was indentified are also omitted.

Alaska↗
USGS2012

Rapid extension recorded by cooling‐age patterns and brittle deformation, Naxos, Greece

The metamorphic core complex exposed as the island of Naxos in the Aegean Sea records an unusually complete sequence of structures developed as a result of continental extension. The structures formed during Miocene rise and cooling from ductile, upper amphibolite facies and anatectic conditions to brittle near‐surface conditions beneath the Naxos detachment fault. Top‐to‐the‐north ductile fabrics in the footwall, which initially developed during amphibolite facies prograde metamorphism, were overprinted by a succession of north directed, normal sense lower‐temperature brittle structures as the footwall was tectonically unloaded and unroofed. Pseudotachylite and cataclasite formation, brittle faulting, alteration, and erosion of the footwall occurred during continued slip and tectonic denudation. Neogene conglomerate and megabreccia, in part derived from exhumation of this footwall, lie structurally above the peripheral Naxos fault. Published K‐Ar and 40 Ar/ 39 Ar ages for hornblende, white mica, and biotite in the footwall decrease northwestward; apparent ages 17–50 Ma in the southeast correspond to areas of low metamorphic grade where preextension argon was partially retained. Published ages 16–10 Ma in higher‐grade rocks of the domal core in the north are cooling ages that for each of the three minerals show a component of younging in the NNE direction of extension. Assuming this is the direction of unroofing, we interpret the rate of this younging as the fault slip rate as the footwall rocks moved >20 km SSW relative to their hanging wall along the base of the Naxos detachment fault. The calculated rates of slip average 5–8 mm/yr, comparable to maximum rates reported in the Basin and Range province.

Naxos↗
USGS2020

Paired air-water annual temperature patterns reveal hydrogeological controls on stream thermal regimes at watershed to continental scales

Despite decades of research into air and stream temperature dynamics, paired air-water annual temperature signals have been underutilized to characterize watershed processes. Annual stream temperature dynamics are useful in classifying fundamental thermal regimes and can enhance process-based interpretation of stream temperature controls, including deep and shallow groundwater discharge, when paired with air signals. In this study, we investigated multi-scale variability in annual paired air-water temperature patterns using sine-wave linear regressions of multi-year daily temperature data from streams of various sizes. A total of 311 sites from two spatially intensive regional datasets (Shenandoah National Park and Olympic Experimental State Forest) and a spatially extensive national dataset spanning the contiguous United States (U.S. Geological Survey gages) were evaluated. We calculated three annual air-water thermal metrics (mean ratio, phase lag, and amplitude ratio) to deduce the influence of groundwater and other watershed processes on stream thermal regimes at multiple spatial scales. Site-specific values of the three annual air-water thermal metrics ranged from 0.69 to 5.29 (mean ratio), −9 to 40 days (phase lag), and 0.29 to 1.12 (amplitude ratio). Regional patterns in the annual thermal metrics revealed persistent yet spatially variable influences of shallow groundwater discharge and high levels of thermal variability within watersheds, indicating the importance of local hydrogeological controls on stream temperature. Furthermore, annual thermal metric patterns from the regional datasets were generally concordant with the national dataset suggesting the utility of these annual thermal metrics for analysis at multiple scales. Analysis of the national dataset showed that previously defined thermal regimes based on water temperature alone could be further refined using air-water metrics and these metrics were related to physiographic watershed characteristics such as contributing area, elevation, and slope. This research demonstrates the importance of spatial scale and heterogeneity for inferring hydrological process in streams and provides guidance for the interpretation of annual air-water temperature metrics that can be efficiently applied to the growing database of multi-year temperature records. Results from this research can aid in the prediction of future thermal habitat suitability for coldwater-adapted species at ecologically and management-relevant spatial scales with readily available data.

Continental United States↗
USGS2010

Seismic characterization of the fall 2007 eruptive sequence at Bezymianny Volcano, Russia

We examine an eruptive sequence in late 2007 at Bezymianny Volcano to characterize the magmatic plumbing system and eruption-related seismicity. Earthquake locations reveal seismicity below and offset to the north of the volcano along a tectonic fault. Based on historical seismicity, the magma chamber is postulated to have a top at about 6 km depth. Minor dome explosions, large sub-plinian eruptions and dome collapses are analyzed using an automated event classification scheme. Low-frequency tremor, interpreted as gas escape, and low-frequency earthquakes are a dominant proportion of the energy released. We also examine multiplet earthquakes whose behavior during the study period changed significantly and systematically before the largest eruption, demonstrating the potential of tracking multiplets to assess changing conditions with the conduit.

Bezymianny Volcano↗
USGS1984

Distribution of Eurasian watermilfoil, Myriophyllum spicatum , in the St. Clair-Detroit River system in 1978

Submersed macrophytes were surveyed at 595 stations located throughout the St. Clair-Detroit River system between Lakes Huron and Erie, 23 August to 13 October 1978. Eurasian watermilfoil ( Myriophyllum spicatum ), first recorded in the system in 1974, became the fourth most common submersed macrophyte in the system by 1978. However, it has not been reported as a widespread nuisance in this system as it has in many other large water bodies in the United States. Observations made during the present study, and interpretation of an aerial photograph, suggest that M. spicatum was a minor nuisance to small boat navigation in portions of the system. Information presented in this study provides a baseline against which future changes in the occurrence of M. spicatum in the St. Clair-Detroit River system can be measured.

Journal of Great Lakes Research↗
USGS2006

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Thermal regime of the State 2-14 well, Salton Sea Scientific Drilling Project

Temperature logs were made repeatedly during breaks in drilling and both during and after flow tests in the Salton Sea Scientific Drilling Project well (State 2–14). The purpose of these logs was to assist in identifying zones of fluid loss or grain and to characterize reservoir temperatures. At the conclusion of the active phase of the project, a series of logs was begun in an attempt to establish the equilibrium temperature profile. Initially, we were able to log to depths below 3 km, but beginning in late May of 1986, it was impossible to log below about 1.8 km owing to casing failure. Our best estimates of formation temperature below 1.8 km are 305° ± 5°C at 1890 m and 355° ± 10°C at 3170 m. For the upper 1.8 km the latest temperature log (October 24, 1986), using a digital “slickline” (heat-shielded downhole recording) device, was within a few degrees Celsius of equilibrium, as confirmed by a more recent log (July 31, 1987) to a depth of ∼ 1 km. As in most other wells in the Salton Sea geothermal field, there is an impermeable, thermally conductive “cap” on the hydrothermal system; this cap extends to a depth of more than 900 m at the State 2–14 well. Thermal conductivities of 19 samples of drill cuttings from this interval were measured at room temperature. The conductivity values were corrected for in situ porosity as determined from geophysical logs and for the effects of elevated temperature. Thermal gradients decrease from about 250 mK m −1 (same as degrees Celsius per kilometer) in the upper few hundred meters to just below 200 mK m −1 near the base of the conductive cap. Using one interpretation, thermal conductivities increase with depth (mainly because of decreasing porosity), resulting in component heat flows that agree reasonably well with the mean of about 450 m W m −2 . This value agrees well with heat flow data from shallow wells within the Salton Sea geothermal field. A second interpretation, in which measured temperature coefficients of quartz- and carbonate-rich rocks are used to correct thermal conductivity, results in lower mean conductivities that are roughly constant with depth and, consequently, systematically decreasing heat flux averaging about 350 mW m −2 below 300 m. This interpretation is consistent with the inference (from fluid inclusion studies) that the rocks in this part of the field were once several tens of degrees Celsius hotter than they are now. The age of this possible disturbance is estimated at a few thousand years.

Journal of Geophysical Research Solid Earth↗
USGS2002

Measurement error associated with surveys of fish abundance in Lake Michigan

In fisheries, imprecise measurements in catch data from surveys add uncertainty to the results of fishery stock assessments. The USGS Great Lakes Science Center (GLSC) began to survey the fall fish community of Lake Michigan in 1962 with bottom trawls. The measurement error was evaluated at the level of individual tows for nine fish species collected in this survey by applying a measurementerror regression model to replicated trawl data. It was found that the estimates of measurement-error variance ranged from 0.37 (deepwater sculpin, Myoxocephalus thompsoni ) to 1.23 (alewife, Alosa pseudoharengus ) on a logarithmic scale corresponding to a coefficient of variation = 66 to 156%. The estimates appeared to increase with the range of temperature occupied by the fish species. This association may be a result of the variability in the fall thermal structure of the lake. The estimates may also be influenced by other factors, such as pelagic behavior and schooling. Measurement error might be reduced by surveying the fish community during other seasons and/or by using additional technologies, such as acoustics. Measurement-error estimates should be considered when interpreting results of assessments that use abundance information from USGS-GLSC surveys of Lake Michigan and could be used if the survey design was altered. This study is the first to report estimates of measurement-error variance associated with this survey.

Lake Michigan↗
USGS2008

Exponential decline of aftershocks of the M7.9 1868 great Kau earthquake, Hawaii, through the 20th century

The remarkable catalog of Hawaiian earthquakes going back to the 1820s is based on missionary diaries, newspaper accounts, and instrumental records and spans the great M 7.9 Kau earthquake of April 1868 and its aftershock sequence. The earthquake record since 1868 defines a smooth curve complete to M 5.2 of the declining rate into the 21st century, after five short volcanic swarms are removed. A single aftershock curve fits the earthquake record, even with numerous M 6 and 7 main shocks and eruptions. The timing of some moderate earthquakes may be controlled by magmatic stresses, but their overall long‐term rate reflects one of aftershocks of the Kau earthquake. The 1868 earthquake is, therefore, the largest and most controlling stress event in the 19th and 20th centuries. We fit both the modified Omori (power law) and stretched exponential (SE) functions to the earthquakes. We found that the modified Omori law is a good fit to the M ≥ 5.2 earthquake rate for the first 10 years or so and the more rapidly declining SE function fits better thereafter, as supported by three statistical tests. The switch to exponential decay suggests that a possible change in aftershock physics may occur from rate and state fault friction, with no change in the stress rate, to viscoelastic stress relaxation. The 61‐year exponential decay constant is at the upper end of the range of geodetic relaxation times seen after other global earthquakes. Modeling deformation in Hawaii is beyond the scope of this paper, but a simple interpretation of the decay suggests an effective viscosity of 10 19 to 10 20 Pa s pertains in the volcanic spreading of Hawaii's flanks. The rapid decline in earthquake rate poses questions for seismic hazard estimates in an area that is cited as one of the most hazardous in the United States.

Hawaii↗
USGS2004

Strain accumulation across the Coast Ranges at the latitude of San Francisco, 1994-2000

A 66-monument geodetic array spanning the Coast Ranges near San Francisco has been surveyed more than eight times by GIPS between late 1993 and early 2001. The measured horizontal velocities of the monuments are well represented by uniform, right-lateral, simple shear parallel to N29°W. (The local strike of the San Andreas Fault is ∼N34°W.) The observed areal dilatation rate of 6.9 ± 10.0 nstrain yr −1 (quoted uncertainty is one standard deviation and extension is reckoned positive) is not significantly different from zero, which implies that the observed strain accumulation could be released by strike-slip faulting alone. Our results are consistent with the slip rates assigned by the Working Group on California Earthquake Probabilities [2003] to the principal faults (San Gregorio, San Andreas, Hayward-Rodgers Creek, Calaveras-Concord-Green Valley, and Greenville Faults) cutting across the GPS array. The vector sum of those slip rates is is 39.8 ± 2.6 mm yr −1 N29.8°W ± 2.8°, whereas the motion across the GPS array (breadth 120 km) inferred from the uniform strain rate approximation is 38.7 ± 1.2 mm yr −1 N29.0°W ± 0.9° right-lateral shear and 0.4 ± 0.9 mm yr −1 N61°E ± 0.9° extension. We interpret the near coincidence of these rates and the absence of significant accumulation of areal dilatation to imply that right-lateral slip on the principal faults can release the accumulating strain; major strain release on reverse faults subparallel to the San Andreas Fault within the Coast Ranges is not required.

California↗
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗
USGS1977

A comparison of some analytical techniques for determining uranium, thorium, and potassium in granitic rocks

Geochemical exploration for uranium requires accurate and precise determinations of low-level concentrations. We have used seven different techniques and four different treatments of the fluorometric method to analyze for uranium in granitic rocks. In addition we have used four analytical techniques for thorium and three analytical techniques for potassium, two elements that are commonly present in anomalous amounts within uranium provinces. Our results show that commonly used techniques for thorium and potassium determinations are both adequately precise and accurate, but that many techniques used for uranium determinations lack the necessary precision or accuracy for complete geochemical prospecting. We suggest that a combination of delayed-neutron determinations for uranium and γ -ray spectrometric analyses for radium equivalent uranium, thorium, and potassium provides the best data base for geochemical exploration for uranium. If more detailed interpretations are desired, the combination of γ -ray spectrometry and α -spectrometry may be best. Carefully done fluorometric analyses should be adequate for water, ore, mineralized rock, and other applications where high precision and accuracy are not required.

Journal of Research of the U.S. Geological Survey↗