USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Waterway, Port, Coastal, and Ocean Engineering”

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,722 records · Page 70Linked to original sources
USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2026

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗
USGS2025

A crosswalk of the 2015 World Terrestrial Ecosystems to the International Union for the Conservation of Nature Global Ecosystem Typology Framework

To support ecosystem mapping and accounting applications, we aligned the 2015 U.S. Geological Survey/Esri/The Nature Conservancy-World Terrestrial Ecosystems (WTEs) with the International Union for Conservation of Nature Global Ecosystem Typology (GET) framework. This process, known as “crosswalking,” enabled the development of a global map of GET level 3 Ecosystem Functional Groups (EFGs) at a 250-meter spatial resolution. Crosswalking involved manually assigning 1,781 biogeographically stratified WTEs to their most probable EFG based on similarities in climate, terrain, vegetation, and geographic distribution. We compared attributes of the WTE dataset with summary characteristics of the EFGs. The resulting crosswalked global map of International Union for Conservation of Nature GET ecosystems is intended to be useful for standardizing ecosystem classification and reporting under frameworks such as the Kunming-Montreal Global Biodiversity Framework and the United Nations System of Environmental-Economic Accounting. We discuss key challenges in reconciling non-identical classifications, such as many-to-one relationships and variation in data quality.

Open-File Report↗
USGS2025

Foundational uncertainties in terminal Ediacaran chronostratigraphy revealed by high-precision zircon U-Pb geochronology of the Nama Group, Namibia

The Nama Group of southern Namibia and northwestern South Africa hosts the best-dated mixed carbonate-siliciclastic foreland basin succession of the terminal Ediacaran [ca. 551 million years (Ma) ago to <538 Ma] and is key for resolving the chronology of early metazoan evolution. Numerous silicified volcanic tuff interbeds are present, but differing interpretations regarding the fidelity of their ages lead to different regional stratigraphic correlations, especially for the Urusis Formation of the Schwarzrand Subgroup. An expanded record of the Urusis Formation is found in the Swartpunt area of southern Namibia, which has yielded an important metazoan biota. But the succession in this area is preserved as a series of thrusts at the leading edge of the Gariep orogenic belt and zircon U-Pb data show systematic age repetition. We use regional stratigraphic and structural mapping, integrated with carbonate carbon isotope (δ 13 C carb ) chemostratigraphy and high-precision radioisotope U-Pb zircon geochronology from outcrop and recently acquired drill core to develop a temporally calibrated basin-wide depositional model. This integrated dataset either reflects complex zircon reworking, inheritance, or potential analytical biases (Scenario 1) or the presence of a Gariep-related cryptic décollement within the Spitskop Member that has resulted in stratigraphic repetition (Scenario 2). We investigate the evidence for and against both scenarios and consider their implications for stratigraphic and δ 13 C carb correlations between the Swartpunt area and coeval exposures along the Orange River border with South Africa. Given that these issues are in an area that hosts numerous silicified ash beds and extensive exposure, an inability to confidently discount either scenario highlights a level of compounding uncertainty in zircon U-Pb geochronology that must be considered when attempting to build global chronostratigraphic frameworks. Scenario 1 implies that some of the weighted mean ages and Bayesian eruption ages from the Swartpunt area may be >1 Myr older than the depositional age of their respective ash beds when assuming existing stratigraphic correlations. If this scenario is preferred, then a cautious approach would be to consider all weighted mean zircon U-Pb ages from ash beds to reflect maximum depositional ages. Both scenarios support deposition of the Huns Member >540 Ma in the Swartpunt area if the oldest weighted mean age reported here represents a near-depositional age, which has significant implications for the temporal calibration of important terminal Ediacaran ichnofossil assemblages and future cyclostratigraphic studies. Stratigraphic correlations common to both scenarios allow us to temporally calibrate a basin evolution model for the Nama Group. Temporal trends in initial hafnium isotope (εHf) compositions of zircon grains from ash beds throughout the succession may support progressive crustal thickening associated with underplating of the Damara orogenic belt along the northern periphery of the Kalahari craton from ca. 547 Ma to ca. 538 Ma. The compilation of new and published zircon U-Pb ages may also imply that the locus of carbonate platform development migrated from north to south (present co-ordinates), tracking the migration of foredeep subsidence.

Earth-Science Reviews↗
USGS2000

Ground-water quality atlas of Oakland County, Michigan

The U.S. Geological Survey (USGS), in cooperation with Oakland County Health Division (OCHD), collected 140 water samples from 38 wells in Oakland County during 1998 to better understand ground-water quality. OCHD had observed temporal variations in concentrations of various constituents, so two additional sets of samples were collected to evaluate potential short-term variability related to sample collection procedures and long-term seasonal variability. Replicate samples from 28 wells were analyzed in the Michigan Department of Environmental Quality (MDEQ) Drinking Water Laboratory to compare MDEQ&rsquo;s analytical results to those obtained from the USGS National Water Quality Laboratory. Several additional databases describing population, land use, water supply, soils, geology, and flows of ground water and surface water are presented in the first part of the report to assist in interpreting the water-quality data. Maps created from these databases are provided in the first portion of the report as an extension of the study-area description. The U.S. Environmental Protection Agency (USEPA) has established Maximum Contaminant Levels (MCL) and Secondary Maximum Contaminant Levels (SMCL) for which samples were analyzed in this study. Water from the 38 wells sampled by the USGS did not exceed the SMCL or MCL for sulfate, fluoride, or nitrite. However, water from 26 wells exceeded the SMCL for iron, water from 12 wells exceeded the SMCL for manganese, and water from 12 wells exceeded the SMCL for dissolved solids. Water from two wells exceeded the MCL for nitrate, although nitrate concentrations in water from most wells was below the detection limit. Water from seven wells exceeded the SMCL for chloride, and water from all wells contained detectable concentrations of chloride. Water from five wells exceeded the MCL for arsenic, and most of the wells sampled contained detectable concentrations of arsenic. These five wells were identified from previous MDEQ analyses to have elevated arsenic concentrations, and were sampled to obtain additional chemistry information. Replicate samples were collected from 26 of the 38 wells for analysis at the MDEQ Drinking-Water Laboratory to compare the results with the USGS National Water Quality Laboratory. The results of the replicate analyses indicate close agreement between the laboratories, with mean differences for nitrate, chloride, and arsenic of 0.10 milligrams per liter (mg/L) as nitrogen, 6.8 mg/L, and 0.0008 mg/L, respectively between the USGS and MDEQ analyses. Potential health effects associated with ingesting nitrate, chloride, and arsenic are provided with the water-quality data, along with references for further information.

Michigan↗
USGS2024

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2026

Albuquerque Seismological Laboratory strategic vision

This circular presents a strategic outlook for the U.S. Geological Survey (USGS) Albuquerque Seismological Laboratory (ASL) for the next 10 years (2026–36). The ASL is a USGS field office in the Geological Hazards Science Center that operates portions of the Advanced National Seismic System and the Global Seismographic Network and focuses on fundamental research for instrumentation testing and data quality. The strategic outlook is categorized into two types of tasks: “Foundational Tasks” and “Aspirational Tasks.” Foundational Tasks are those that maintain the laboratory’s basic operations and services, including regional and global seismic monitoring, improving data quality, and providing instrument testing and support. A suite of Aspirational Tasks is also articulated; these can be considered priority targets of ASL that could improve ASL’s seismic monitoring capabilities and mission. Such tasks include improvements to remote stations, testing capabilities of nonseismic geophysical instruments, detection threshold monitoring, rapid aftershock deployments, and expanding seismic monitoring networks. This report was written with input from the USGS Geological Hazards Science Center, the USGS Earthquake Hazards Program (EHP), and colleagues with an interest in the work done by the ASL. Although the details of these tasks may change, this document can provide guidance on the overarching tasks at the ASL from 2026 to 2036 and an overview of the various components of the ASL and how they fit into the EHP and the Global Seismographic Network Program.

Circular↗
USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗
USGS2025

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences↗
USGS1976

Non-radiometric borehole geophysical detection of geochemical halos surrounding sedimentary uranium deposits

Roll-type uranium deposits are formed by the concentration of uranium by ground water in geochemical cells. Non-uranium minerals having different solubilities may be deposited ahead of or behind the uranium minerals, forming halos that surround the ore. In addition, oxidizing and reducing environmental conditions may cause zones of mineral alteration to develop beyond the limits of the uranium deposit. Certain physical property anomalies that are commonly associated with halos can be detected by relatively fast and inexpensive borehole geophysical measurements made either in individual holes, or between two adjacent holes. Borehole measurements that have been found to be useful include electrical resistivity, induced polarization, and magnetic susceptibility. Electrical resistivity is increased by the presence of calcite and other cementing minerals that sometimes create permeability barriers in the neighborhood of uranium deposits. Induced polarization (IP) response is increased by sulfide and clay minerals that are commonly found in anomalous concentrations near roll-type deposits. Magnetic susceptibility is usually decreased by the oxidation of magnetite to hematite or limonite in the zone of chemical alteration that is left as a trail behind roll fronts. Borehole measurements of electrical resistivity, induced polarization and magnetic susceptibility were made in the vicinity of a uranium roll-type deposit in south Texas. Results indicate that mineral halos can be detected by borehole measurements made. in wide-spaced drill holes, and that the total amount of drilling needed to find a deposit can be reduced substantially by this exploration approach.

Open-File Report↗
USGS2026

Perspectives and insights gained from applying landscape ecological approaches to the Upper Mississippi riverscape

Advances in remote sensing, spatial data collection, hydraulic and geospatial models, and the maturation of landscape ecology during the early 2000s spurred the concept of rivers as landscapes or ‘riverscapes’. Since that time, the Upper Mississippi River Restoration (UMRR) Program - one of the world’s longest running large-river ecosystem restoration and monitoring programs - has applied landscape ecological principles to the Upper Mississippi River System (UMRS). I examine spatial patterns of aquatic nutrient concentrations, submersed aquatic vegetation, and fish and mussel communities, finding that within-river patterns are often patchy in nature, likely resulting from how nutrient processes and community traits interact with spatial variability in hydraulic exchange patterns. Local-scale patchy spatial patterns in the UMRS are embedded in system-scale longitudinal gradients and a hierarchical view of the river has provided a foundation for strategic river restoration planning. Key Words: Aquatic Habitat, hydraulic connectivity, patch, restoration, riverscape

BioScience↗
USGS1976

Distribution functions for statistics derived from bivariate normal and bivariate two-parameter log-normal populations

The distribution functions for statistics that may be used to assess the significance of differences between sample means, standard deviations, coefficients of skewness, and coefficients of variation are obtained by Monte Carlo experiments. The distributions are expressed as functions of sample size, cross correlation, and skewness. In general, the distributions are more sensitive to cross correlation than to skewness. As sample size increases, however, the distributions tend to become more sensitive to skewness.

Open-File Report↗
USGS2025

A global database of soil microbial phospholipid fatty acids and enzyme activities

Soil microbes drive ecosystem function and play a critical role in how ecosystems respond to global change. Research surrounding soil microbial communities has rapidly increased in recent decades, and substantial data relating to phospholipid fatty acids (PLFAs) and potential enzyme activity have been collected and analysed. However, studies have mostly been restricted to local and regional scales, and their accuracy and usefulness are limited by the extent of accessible data. Here we aim to improve data availability by collating a global database of soil PLFA and potential enzyme activity measurements from 12,258 georeferenced samples located across all continents, 5.1% of which have not previously been published. The database contains data relating to 113 PLFAs and 26 enzyme activities, and includes metadata such as sampling date, sample depth, and soil pH, total carbon, and total nitrogen. This database will help researchers in conducting both global- and local-scale studies to better understand soil microbial biomass and function.

Scientific Data↗
USGS2021

Factors influencing distribution of Coccidioides immitis in soil, Washington State, 2016

Coccidioides immitis and Coccidioides posadasii are causative agents of Valley fever, a serious fungal disease endemic to regions with hot, arid climate in the United States, Mexico, and Central and South America. The environmental niche of Coccidioide s spp. is not well defined, and it remains unknown whether these fungi are primarily associated with rodents or grow as saprotrophs in soil. To better understand the environmental reservoir of these pathogens, we used a systematic soil sampling approach, quantitative PCR (qPCR), culture, whole-genome sequencing, and soil chemical analysis to identify factors associated with the presence of C. immitis at a known colonization site in Washington State linked to a human case in 2010. We found that the same strain colonized an area of over 46,000 m 2 and persisted in soil for over 6 years. No association with rodent burrows was observed, as C. immitis DNA was as likely to be detected inside rodent holes as it was in the surrounding soil. In addition, the presence of C. immitis DNA in soil was correlated with elevated levels of boron, calcium, magnesium, sodium, and silicon in soil leachates. We also observed differences in the microbial communities between C. immitis -positive and -negative soils. Our artificial soil inoculation experiments demonstrated that C. immitis can use soil as a sole source of nutrients. Taken together, these results suggest that soil parameters need to be considered when modeling the distribution of this fungus in the environment.

Washington↗
USGS1997

Hydrologic and water-quality conditions in the Horse Creek basin, west-central Florida, October 1992-February 1995

A baseline study of the 241-square-mile Horse Creek basin was undertaken from October 1992 to February 1995 to assess the hydrologic and water-quality conditions of one of the last remaining undeveloped basins in west-central Florida. During the period of the study, much of the basin remained in a natural state, except for limited areas of cattle and citrus production and phosphate mining. Rainfall in 1993 and 1994 in the Horse Creek basin was 8 and 31 percent, respectively, above the 30-year long-term average. The lowest and highest maximum instantaneous peak discharge of the six daily discharge stations occurred at the Buzzard Roost Branch and the Horse Creek near Arcadia stations with 185 to 4,180 cubic feet per second, respectively. The Horse Creek near Arcadia station had the lowest number of no-flow days with zero days and the Brushy Creek station had the highest number with 113 days. During the study, the West Fork Horse Creek subbasin had the highest daily mean discharge per square mile with 30.6 cubic feet per second per square mile, and the largest runoff coefficient of 43.7 percent. The Buzzard Roost Branch subbasin had the lowest daily mean discharge per square mile with 5.05 cubic feet per second per square mile, and Brushy Creek and Brandy Branch shared the lowest runoff coefficient of 0.6 percent. Brandy Branch had the highest monthly mean runoff in both 1993 and 1994 with 11.48 and 19.28 inches, respectively. During the high-baseflow seepage run, seepage gains were 8.87 cubic feet per second along the 43-mile Horse Creek channel. However, during the low-baseflow seepage run, seepage losses were 0.88 cubic foot per second. Three methods were used to estimate average annual ground-water recharge in the Horse Creek basin: (1) well hydrograph, (2) chloride mass balance, and (3) streamflow hydrograph. Estimated average annual recharge using these three methods ranged from 3.6 to 8.7 inches. The high percentage of carbonate plus bicarbonate analyzed at the Carlton surficial aquifer well could indicate an upward ground-water flow from the underlying intermediate aquifer system. Based on constituent concentrations in water samples from the six daily discharge stations, concentrations generally are lower in the upper three subbasins, West Fork Horse Creek, Upper Horse Creek, and Brushy Creek than in the lower three subbasins. Typically, concentrations were highest for major ions at Buzzard Roost Branch and nutrients at Brushy Creek.

Florida↗
USGS2025

From hydrated silica to quartz: Potential hydrothermal precipitates found in Jezero crater, Mars

On Earth, silica-rich phases from opal to quartz are important indicators and tracers of geological processes. Hydrated silica, such as opal, is a particularly good matrix for the preservation of molecular and macroscopic biosignatures. Cherts, a type of silica-dominated rocks, provide a unique archive of ancient terrestrial life while quartz is the emblematic mineral of the Earth's continental crust. On Mars, hydrated silica has been detected in several locations based on remote sensing and rover-based studies. In the present article we report on the detection of cobbles made of hydrated silica (opal or chalcedony), as well as well-crystallized quartz. These detections were made with the SuperCam instrument onboard Perseverance (Mars 2020 mission), using a combination of LIBS, infrared and Raman spectroscopy. Quartz-dominated stones are detected unambiguously for the first time on the Martian surface, and based on grain size and crystallinity are proposed to be of hydrothermal origin. Although these rocks were all found as float, we propose that these detections are part of a common hydrothermal system, and represent different depths / temperatures of precipitation. This attests that hydrothermal processes were active in and around Jezero crater, possibly triggered by the Jezero crater-forming impact. These silica-rich rocks, in particular opaline silica, are very promising targets for sampling and return to Earth given their high biosignature preservation potential.

Earth and Planetary Science Letters↗