USGS Science⌕ Search

SEARCH · USGS Science

Results for “Limnology and Oceanography”

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,656 records · Page 92Linked to original sources
USGS2024

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education↗
USGS2024

Shallow lake, strong shake: Record of seismically triggered lacustrine sedimentation from the 1959 M7.3 Hebgen Lake earthquake within Henrys Lake, Idaho

We investigate a shallow lake basin for evidence of a large historic intraplate earthquake in western North America. Henrys Lake, Idaho is an atypical candidate for a lacustrine paleoseismic study given its shallow depth (~7 m) and low relief (≤2° slopes ). Here, we test the earthquake-recording capacity of this basin type by showing sedimentological evidence of the 1959 M7.3 Hebgen Lake earthquake within sediment cores, using anthropogenically produced 137 Cs activity to constrain timing. In addition to expanding the morphologic range of basins targeted for lacustrine paleoseismic studies, this work has implications for sediment response in dam-enhanced basins. Lack of sedimentological evidence for other earthquakes coupled with radiocarbon chronology reveals that the 1959 event is the only clearly recorded earthquake within Henrys Lake since the mid-Holocene. Henrys Lake offers a proxy for paleo-earthquake signatures within similar lacustrine environments and underscores the importance of further paleoseismic studies in the region.

Idaho↗
USGS2026

Effect of mineral deposit data on predictions from the three-part approach to quantitative mineral resource assessment—A study of 16 previous U.S. Geological Survey assessments

The three-part approach to quantitative mineral resource assessment requires information about the properties of undiscovered mineral deposits in an assessment area. These properties are unknown, so the properties of discovered mineral deposits of the same mineral deposit type are used instead. In the three-part approach, these discovered mineral deposits come from around the world, and their properties constitute the pooled data for that mineral deposit type. Alternatively, these discovered mineral deposits could come from the assessment area, and their properties constitute the tract data for that mineral deposit type. Tract data may be more representative of the undiscovered mineral deposits in the assessment area than the pooled data. The goal of this study was to determine whether resource predictions using pooled data are equivalent to resource predictions using tract data. To this end, 16 previous U.S Geological Survey assessments were studied. For each assessment, resources were predicted for one undiscovered mineral deposit in the assessment area. One set of predictions used pooled data, and another used tract data. The two sets of predictions were compared with an equivalence test, using the six assessment statistics that are commonly reported for mineral resource assessments. Practical equivalence is the condition that two corresponding assessment statistics are within a factor of 1.5 of one another. For each of 2 assessments, all 6 assessment statistics were practically equivalent. For both assessments, the assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 1.30 times smaller to 1.03 times larger. For each of 14 assessments, 1 or more of the 6 assessment statistics were not practically equivalent. The assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 26.6 times smaller to 5.53 times larger. The use of pooled data has been a standard procedure in the three-part approach since at least 1986. The 16 assessments in this study are not a representative sample of those prior assessments that used pooled data. So, it is inappropriate to use the study results to infer whether pooled data affected the resource predictions for those prior assessments.

Scientific Investigations Report↗
USGS2024

Rock sample photogrammetry

This step-by-step protocol describes the photogrammetry process used by the U.S. Geological Survey Spokane Imaging Lab (SPIMG) lab to create 3D models of geologic samples. Steps related to photographing small objects are applicable to photogrammetry in general, however, SPIMG-specific steps involving lab hardware and software may not be.

Protocols.io↗
USGS2025

Haunted Summerville: Ghostly lights or earthquake lights?

Among the colorful local lore in the Charleston, South Carolina, area, are a number of ghost stories, shared not only over campfires but also in published books. Among the most well-known of the stories is the tale of the Summerville Light. Local lore holds that a strange light sometimes seen in a remote area is a lantern carried by the ghost of a local woman who once waited hours for her husband, who turned out to have been decapitated earlier that day in a train accident (DePoppe, 2023). Extant sources suggest the ghost stories began to circulate in the 1950s to 1960s. So pervasive was the lore that (Old) Sheep Island Road became known among local residents as Light Road, with a local stretch of road known today as Old Light Road. Reviewing the location where the lights appear as well as the nature of accounts, I suggest that many if not all of the anecdotal observations can be most readily attributed to natural phenomena, including earthquake lights from earthquakes that were too small to be felt. Accounts of lights near Summerville cluster in proximity to the generally accepted epicenter of the 1886 Charleston, South Carolina, earthquake, where foreshocks to the 1886 mainshock were apparently concentrated, and within a few kilometers of three M3.5 – 4.4 earthquakes in 1959 and 1960.

South Carolina↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗
USGS2025

Borehole geophysical time-series logging to monitor passive ISCO treatment of residual chlorinated-ethenes in a confining bed, NAS Pensacola, Florida

In-situ chemical oxidation (ISCO) is a common method to remediate chlorinated ethene contaminants in groundwater. Monitoring the effectiveness of ISCO can be hindered because of insufficient observations to assess oxidant delivery. Advantageously, potassium permanganate, one type of oxidant, provides the opportunity to use its strong electrical signal as a surrogate to track oxidant delivery using time-series borehole geophysical methods, like electromagnetic (EM) induction logging. Here we report a passive ISCO (P-ISCO) experiment, using potassium permanganate cylinders emplaced in boreholes, at a chlorinated ethene contamination site, Naval Air Station Pensacola, Florida. The contaminants are found primarily at the base of a shallow sandy aquifer in contact with an underlying silty-clay confining bed. We used results of the time-series borehole logging collected between 2017 and 2022 in 4 monitoring wells to track oxidant delivery. The EM-induction logs from the monitoring wells showed an increase in EM response primarily along the contact, likely from pooling of the oxidant, during P-ISCO treatment in 2021. Interestingly, concurrent natural gamma-ray (NGR) logging showed a decrease in NGR response at 3 of the 4 wells possibly from the formation of manganese precipitates coating sediments. The coupling of time-series logging and well-chemistry data allowed for an improved assessment of passive ISCO treatment effectiveness.

Florida↗
USGS2025

Remote single-station seismic monitoring of the July–October 2022 earthquake swarm at Tau volcano, American Samoa

From July to October 2022, a non-eruptive volcanic earthquake swarm occurred within ~15 km of Taʻū Island, located in eastern American Samoa. Felt reports from local residents were the only available information about the swarm when it started, as American Samoa lacked a seismic monitoring network. We developed a consistent single-station catalog for the entire swarm, using seismic data from the nearest station IU.AFI, ~250 km away. We applied the EQTransformer deep-learning model (Mousavi et al., 2020), automatically picking Pn and Sn arrivals on IU.AFI continuous data. We retained only events with Sn-Pn times of 22.5–25 seconds, consistent with the expected locations based on felt reports, then detected smaller swarm events with subsequent template-matching. This single-station catalog characterized the swarm’s onset and escalation to peak activity before a multi-agency field response team installed a local seismic network in mid-August 2022. This permanent seismic network captured the swarm’s decline. EQTransformer identified short S-P times on the first two locally deployed seismometers, both Raspberry Shake sensors, to constrain the swarm’s distance from Taʻū Island. Modern seismological processing methods, combined with basic observations such as felt reports, can quickly contribute useful information during an earthquake response in a poorly monitored region.

American Somoa↗
USGS1998

Satellite imagery of the north-central Gulf of Mexico, June 1993-May 1997

This CD contains imagery collected by the Advanced Very High Resolution Radiometer (AVHRR) on the NOAA polar-orbiting weather satellites. The AVHRR provides almost daily coverage of a site at a resolution of approximately 1 km. The data types included are sea surface temperature (SST), water reflectance (REF), and a false-color infrared overview (FCI).

Alabama, Florida, Louisiana, Mississippi↗
USGS2025

Reflections on a trio of North American earthquakes in 1925

In 1925, three moderately large damaging earthquakes occurred in North America over four months: the 28 February (local time; LT) M 6.2 Charlevoix, 27 June (LT) M 6.6 Montana, and 29 June M 6.5 Santa Barbara earthquakes. The centennial anniversaries of these events motivated this retrospective consideration focused on the ground motions generated by the three events, including a reconsideration of early intensity assignments for the Montana earthquake. At the time, these three earthquakes appeared to support the arguments of some geologists who downplayed the severity of seismic hazard in southern California relative to other parts of the country. Some of the arguments advanced at that time, for example that Los Angeles “has the least to fear from ‘Acts of God’ of any city under the American flag,” ( Hill, 1928 ) sound naïve if not laughable now, but a comparison of well‐constrained shaking distributions for the three earthquakes reveals the dramatic difference in wave propagation efficiency in western versus eastern North America (ENAM), which leads to moderate ENAM events being felt to much larger distances. At M 6.2, the 1925 Charlevoix earthquake was a notably large event in ENAM. This earthquake was the largest event in eastern Canada since 1870 and caused damage in the epicentral region in addition to towns as far away as 200 km, with felt shaking extending over 1000 km. In contrast, felt shaking from the Santa Barbara earthquake barely extended beyond ∼200 km. Compiling published intensity distributions for larger ENAM earthquakes, we show that perceptible earthquake shaking is not uncommon in ENAM over century time scales, but experience with weakly felt shaking may incline people to downplay potential earthquake risk.

California, Montana, Quebec↗
USGS2026

U.S. Geological Survey Monitoring Milestones—Oe-151 at Woodgate, NY (433112075091501)

On July 9, 1926, monitoring well Oe-151 at Woodgate, New York (USGS ID 433112075091501) recorded its first groundwater data. Since then, the well has provided water data nearly continuously and has now reached a 100-year milestone for data collection. The well is part of the U.S. Geological Survey (USGS) Climate Response Network (CRN), which is a national network of wells selected to monitor natural groundwater conditions. Well Oe-151 is the first well in the network to reach a 100-year Centennial milestone.

New York↗
USGS2025

U.S. Geological Survey National Groundwater Climate Response Network

What is the U.S. Geological Survey National Groundwater Climate Response Network? The U.S. Geological Survey’s (USGS) National Water Monitoring Network is a network of networks that includes the Climate Response Network (CRN). The CRN is a network of wells selected to monitor the effects of climate variability, including droughts, on groundwater systems within the United States, Puerto Rico, and the U.S. Virgin Islands. The primary purpose of the CRN is to monitor these effects on groundwater levels in unconfined or near-surface confined aquifers that are minimally affected by pumping or other human-induced stresses.

Fact Sheet↗
USGS2024

Getting to the surface of the problem: A dynamic rupture benchmark for shallowly dipping faults near Earth’s surface

The Statewide California Earthquake Center (SCEC) Dynamic Rupture Group workshop “Getting to the Surface of the Problem: A Dynamic Rupture Benchmark for Shallowly-Dipping Faults Near Earth’s Surface” was convened virtually, in an online meeting on November 4, 2024. 52 people participated, including two who sent pre-recorded lightning talks earlier, but were unable to join us that day due to the international time-zone difference. One-half of our participants were students or postdocs, and the group included scientists from 11 countries (USA, Canada, China, Czech Republic, France, Germany, Iceland, Japan, New Zealand, Saudi Arabia, Switzerland).

Final Report↗
USGS2025

U.S. Geological Survey geomagnetic variometer data: Capitalizing on seismic infrastructure

The U.S. Geological Survey’s Geomagnetism Program is collaborating with the Earthquake Hazards Program and Global Seismographic Network Program to densify magnetic field observations. This collaboration focuses on the installation of magnetometers, or magnetic variometers, at existing seismic stations. Along with improving the density of space weather observations for hazard monitoring, these data can be used to correct colocated magnetic field induced noise in seismic data. Such corrections are especially useful during time periods of large magnetic storms where the magnetic field‐induced instrument noise can be of similar amplitude to earthquake ground‐motion records.

contiguous United States↗
USGS2024

Evaluating the impact of uncertainty in ground motion forecasts for post-earthquake impact modeling applications

The US Geological Survey’s (USGS) ShakeMap system provides a rapid characterization of strong ground shaking in areas directly affected by an earthquake. This study focuses on studying the aggregate effects of macroseismic shaking estimates from ShakeMap, expressed in terms of modified Mercalli intensity (MMI), when accounting for the uncertainty in forecasted ground motions. We use a Monte Carlo approach to generate numerous spatially correlated realizations of ground motions by utilizing a combination of circulant embedding and kriging techniques for efficiently handling the correlations. We then assessed the aggregate effects of shaking by looking at bin counts across these realizations. We demonstrate that the aggregate shaking regarding the mean macroseismic intensity estimates (from the ShakeMap output) is a biased representation of the aggregate shaking when shaking uncertainty is included. Incorporating shaking uncertainty can help to improve various downstream earthquake impact applications, such as the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) overall earthquake fatality distribution or estimates of shaking-induced ground failure impacts from consequential earthquakes.

Earthquake Spectra Journal↗
USGS1976

Digital telemetry system network installation for the National Earthquake Information Service

A Digital Telemetry System Network, developed by the Albuquerque Seismological Laboratory, was installed for the USGS National Earthquake Information Service (NEIS) during January 1975. NEIS is located in Golden, Colorado. The initial Digital Telemetry System Network consists of four independent digital telemetry systems that provide wide area coverage over most of the United States. The four systems were installed between: Albuquerque, New Mexico to NEIS; St. Louis, Missouri to NEIS; Salt Lake City, Utah to NEIS; and Menlo Park, California to NEIS. Each digital telemetry system will transmit fourteen channels of high accuracy seismic data in a digital format over a standard C-1 grade telephone line. Digitized data words will have a resolution of eleven digital bits plus one sign bit. Each digital telemetry receiver, located at NEIS, will have the capacity to provide fourteen channels of seismic data in a digital format suitable for either on-line digital computer analysis or magnetic tape recording. In addition the data can be presented in an analog format suitable for recording on helicorders or develocorders. This network is the first large scale digital telemetry system used by the U.S. Geological Survey for the transmission of real time seismic data.

Open-File Report↗
USGS2026

Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations

Earthquake ground‐motion intensities, such as pseudospectral accelerations (SAs), can vary significantly with horizontal orientation. However, conducting probabilistic seismic hazard analysis (PSHA) for each horizontal orientation is challenging because current ground‐motion models used in PSHA consider only a single horizontal intensity value, usually the median across all orientations, known as RotD50. To address this limitation, we employ physics‐based simulations for PSHA, which contain full waveforms from which ground‐motion intensities can be computed for all horizontal orientations to study directional seismic hazard. We apply our approach to the latest CyberShake study of the Greater Los Angeles metropolitan area, developed by the Statewide California Earthquake Center, finding that seismic hazard at a 2475‐yr return period, a common value used for earthquake‐resistant design, varies significantly with horizontal orientation. For instance, for SAs at 3 s, the maximum seismic hazard across all horizontal orientations is, on average, 15% higher than the median RotD50 hazard, with these differences becoming more pronounced at longer periods. These observed variations can generally be attributed to physical mechanisms that polarize seismic waves, such as the radiation pattern of the earthquake source and the influence of the subsurface structure. These results may have important implications for earthquake engineering applications, particularly for long‐period structures in areas with substantial horizontal variations in seismic hazard.

California↗