USGS Science⌕ Search

SEARCH · USGS Science

Results for “Critical Reviews in Environmental Science and Technology”

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

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS2026

Indirect mineral import reliance and provenance

Mineral commodity supply chain analyses rely on international trade data reported by individual countries as quantities of a mineral commodity form imported from (or exported to) a partner. However, export quantities frequently exceed a country’s domestic production, or occur when no production data are reported, suggesting that the trade partner is merely an intermediary in a transshipment. These discrepancies can result in misleading conclusions regarding supply chain vulnerabilities and dependencies. We present a two-stage methodology to reconcile gaps between reported material sources and actual producers. First, we construct trade networks for specific mineral forms, treating production as a type of import to distinguish producing nations from entrepôts. By tracing flows through these networks, we attribute a target country’s imports to original producers via both direct (in a single trade link) and indirect (transferring through intermediaries) pathways. Second, these production-attributed flows are incorporated into multi-stage supply chains to determine the upstream provenance of feedstock for domestic refining and processing. This approach provides a more representative picture of trade reliance. For example, while the United States (U.S.) Geological Survey reports no imports of unwrought antimony metal from Russia in 2022 (U.S. Geological Survey (2025). Mineral Commodity Summaries 2025. 10.3133/mcs2025), our analysis reveals that over 16% of U.S. imports can be traced back to Russian mining through intermediate processing in countries such as China, India, and Vietnam. Additionally, our analysis of the aluminum supply chain shows that while the U.S. is reported as 52% net import reliant on aluminum materials in 2022, it is 100% reliant on foreign bauxite, 7% of which arrived indirectly. This unreported reliance, which is predominantly tied to bauxite mined in Brazil (43%) and Jamaica (28%), highlights our methods ability to capture the supply chain’s dependence on foreign feedstock that may be missing in single-stage trade data.

Mineral Economics↗
USGS2026

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River near Kansas City, Missouri, August 8–9, 2023

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 8 bridge crossings of the Missouri River near Kansas City, Missouri, on August 8–9, 2023. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches that extended about 1,550 to 1,640 feet longitudinally and generally extended laterally across the active channel from bank to bank during low floodflow to nonflood conditions. These surveys provided the channel geometry and hydraulic conditions of the river at the time of the surveys and provided characteristics of scour holes, which may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. The data collected from the surveys may also be useful to the Missouri Department of Transportation as a record of low floodflow conditions in regards to the stability and integrity of the bridges with respect to bridge scour. Bathymetric data were collected around every in-channel pier. Scour holes were at most piers where bathymetry could be obtained, except for those piers on banks or surrounded by riprap. All the bridge sites in this study were surveyed and documented in previous studies. The average difference between the bathymetric surfaces ranged from 0.07 to 4.16 feet higher in 2023 than 2019, which indicates overall deposition between the survey dates, as might be expected based purely on streamflow at the time of the survey. However, the average difference between the bathymetric surfaces ranged from 1.44 feet higher to 1.88 feet lower in 2023 than 2015, which indicates a dynamic equilibrium of scour and deposition overall between those surveys, despite the lower flow conditions in 2023. Similarly, the average difference between the bathymetric surfaces ranged from 3.18 feet higher to 5.19 feet lower in 2023 than 2011, which indicates a relative equilibrium between scour and deposition overall, albeit the trend was toward scour as might be expected because of the substantial flood event in 2011. Riprap blankets and alignment to flow had a substantial effect on the size of the scour hole for a given pier. Piers that were partially or fully surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers with no rock or riprap and effectively mitigated the scour holes historically observed at these piers. Several of the structures had piers that were skewed to primary approach flow. At most of the structures, the scour hole was deeper and longer on the side of the pier with impinging flow than the leeward side, with some amount of deposition on the leeward side, as typically observed at piers skewed to approach flow.

Kansas, Missouri↗
USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2025

Extremophile hotspots linked to containerized industrial waste dumping in a deep-sea basin

Decaying barrels on the seafloor linked to DDT contamination have raised concerns about the public health implications of decades old industrial waste dumped off the coast of Los Angeles. To explore their contents, we collected sediment cores perpendicular to five deep-sea barrels. The concentration of DDT and its breakdown products were highly elevated relative to control sites yet did not vary with distance from the barrels, suggesting that they were not associated with the contamination. Sediment cores collected through white halos surrounding three barrels were enriched in calcite and had elevated pH. The associated microbial communities were low diversity and dominated by alkalophilic bacteria with metagenome-assembled genomes adapted to high pH. A solid concretion sampled between a white halo and barrel was composed of brucite, a magnesium hydroxide mineral that forms at high pH. Based on these findings, we postulate that leakage of containerized alkaline waste triggered the formation of mineral concretions that are slowly dissolving and raising the pH of the surrounding sediment pore water. This selects for taxa adapted to extreme alkalinity and drives the precipitation of “anthropogenic” carbonates forming white halos, which serve as a visual identifier of barrels that contained alkaline waste. Remarkably, containerized alkaline waste discarded >50 years ago represents a persistent pollutant creating localized mineral formations and microbial communities that resemble those observed at some hydrothermal systems. These formations were observed at one-third of the visually identified barrels in the San Pedro Basin and have unforeseen, long-term consequences for benthic communities in the region.

California↗
USGS2025

Comparison of creek and bay influences on salt marsh sediment budget and deposition patterns

The resilience of salt marshes with low organic production depends on their effective capture and retention of mineral sediment from adjacent waters. Little prior work has directly compared mechanisms of sediment import from wave-influenced marsh boundaries against those of tidal creeks. We used simultaneous deployment of net-deposition tiles and oceanographic sensors to identify the timing and magnitude of sediment import/export to, and redistribution within, a marsh in south San Francisco Bay. As the marsh has both an eroding bay-exposed scarp and a prominent tidal creek, we investigated the mechanisms and magnitudes of sediment import from the marsh-bay versus the marsh-creek interface. The strong daily sea breezes of the summer season produced most of the wave-driven erosion of the marsh scarp and controlled suspended sediment concentrations; the winter season had weaker winds punctuated by a few storms. A large seasonal difference in suspended sediment concentrations influenced both flood and ebb sediment fluxes to the marsh and led to much higher rates of import in the summer. Both bay-side and creek-side processes were important to total marsh sediment budget. Bay-side sediment contributions were more variable in time due to the bay-influenced environment, and creek-side contributions were overall larger, reflecting the large proportion of the marsh fed by creek water. Sediment was redistributed throughout the system, with erosion near the bay-edge, accretion near the creek-edge and slow import to the marsh interior. The marsh was net importing sediment in the summer and exporting in the winter from different rates of these processes; on an annual scale, the marsh was net importing despite rapid lateral marsh loss. These findings emphasize that a positive sediment budget does not imply a stable marsh and that both creek- and edge-side dynamics are important for marsh sedimentation and geomorphic trajectories. Further, we expand understandings of non-storm and seasonal controls on marsh sedimentation.

California↗
USGS2026

Evaluating the U.S. Geological Survey’s earthquake shaking hazard forecasts and their implications on seismic risk

We analyze the last six update cycles of the long-term probabilistic earthquake shaking hazard forecast from the U.S. Geological Survey (USGS) and discuss the changes in hazard estimates from the 1996 to the latest 2023 update of the National Seismic Hazard Model (NSHM) for the conterminous U.S. We summarize how our understanding of earthquake hazards has evolved over the last two decades and quantify the implications of changing hazards on people, buildings, and lifeline infrastructure. The net effect of changes in hazard estimates, with reference to the mean hazard from the 2023 NSHM update that has 2% probability of exceedance in 50 years, suggests an overall increase in total geographic areas facing very strong shaking levels (modified Mercalli Intensity (MMI) of VII or more) when compared with any of the previous cycles. This increase puts ~159.2 million people (an increase from 111.8 million when using the 2018 NSHM), ~47.3 million residential buildings (an increase from 32.5 million with the 2018 NSHM), and ~$25.6 trillion of economic exposure of all buildings (an increase from $17.4 trillion when using the 2018 NSHM) at risk from earthquakes. We demonstrate that small changes in hazard estimates do not necessarily imply small changes in risk estimates. Therefore, changes in risk estimates can be used to highlight key USGS NSHM updates for improving risk mitigation efforts.

Earthquake Spectra Journal↗
USGS2025

A 10-year continuous daily simulation of chloride flux from a suburban watershed in Fairfax County, Virginia, USA

Increasing levels of chloride in surface water are associated with detrimental effects on water quality, aquatic ecosystems, infrastructure, and human health. Numerous mass-balance studies have inferred watershed transport processes by interpreting chloride inputs and outputs, but few represent internal dynamics explicitly. We constructed a coupled water/chloride mass balance model to gain insights into storage, residence time, and transport processes in a 10-km 2 urban watershed. The model, which operates over a 10-year period at a daily time scale, represents storage in a dynamic soil-moisture reservoir, quick-flow runoff from storm events, and slow-flow runoff that sustains streamflow in dry weather. The calibrated model accurately represented (a)the observed transition from a streamflow enrichment regime in cold months to a dilution regime in warmer months, (b) the observed tendency for late-summer concentrations to be higher after winters with heavy snowfall, and (c) a period-of-record downward trend in chloride concentration likely associated with a downward trend in annual snowfall. Estimated chloride inputs averaged 195 metric tons per year, while the average output was 270 metric tons per year. In contrast, estimated storage was only 107 metric tons. The estimated mean residence time in groundwater was 1.27 years. This short residence time indicates that efforts to reduce inputs will manifest as decreased concentrations in streamflow on a management-relevant time scale of several years. The coupled mass balance model yielded insights into internal watershed dynamics that would not be possible from simple input/output analysis; such models can be useful tools for gaining insight into small watershed hydrology and pollutant transport.

Virginia↗
USGS2026

Measuring stress In high pressure deformation experiments with high speed fiber-optics

High‐pressure, high‐temperature rock deformation experiments are essential for understanding deep Earth processes, but accurately measuring stress is challenging due to the inaccuracy introduced by seal friction within the apparatus and large inertia, which affects earthquake process measurements made far from the experimental fault. To overcome these limitations, we developed and implemented a simple, compact fiber‐optic sensor based on an External Cavity Fabry–Perot Interferometer for in situ load measurement on a piston inside the pressure vessel of a Griggs‐type high‐pressure apparatus. The sensor can be used at sample temperatures up to 800°C for both slow creep and fast rupture testing with bandwidth capability from DC to 6 MHz. Two important experimental results are described: (1) the first direct measurements of seal friction in this apparatus type; these measurements record seal friction approximately one‐tenth of the confining pressure under both low (150 MPa) and high (1 GPa) pressure conditions; and (2) high‐temporal resolution capture of dynamic stress drops during stick‐slip faulting events at high‐confining pressure, revealing near‐fault wave propagation details, high‐frequency oscillations, and implying high‐transient slip rates (6–12 m/s) previously inaccessible with external sensors.

The Seismic Record↗
USGS2024

Geologic map of the northwest flank of Mauna Loa volcano, Island of Hawai‘i, Hawaii

Mauna Loa, the largest active volcano on Earth, has erupted 34 times since written descriptions became available in A.D. 1832. The most recent eruption of Mauna Loa occurred on November 27, 2022, after a 38 year hiatus; it lasted for 12 days. Some eruptions began with only brief seismic unrest, whereas others followed several months to a year of increased seismicity. Once underway, Mauna Loa’s eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows that erupted from the Southwest Rift Zone in 1950 advanced at an average rate of 9.3 kilometers per hour (5.8 miles per hour); all three lobes reached the ocean within ~24 hours. Near the eruptive vents, the flows likely traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa has great volcanic-hazard potential for the Island of Hawai‘i. Volcanic hazards on Mauna Loa can be anticipated, and risk substantially mitigated, by documenting its past activity to refine our knowledge of the hazards, and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. The map of the north and west flanks of Mauna Loa shows the distribution and relation of volcanic and surficial sedimentary deposits. It incorporates previously reported work published as generalized small-scale maps and a more detailed map. Within the mapped area, lava has flowed from three different source regions: the Northeast Rift Zone (22 percent), the summit (64 percent), and radial vents (14 percent). All three have different points of origin which, in turn, affect the flow characteristics and periodicity of activity. The map area includes the uppermost part of the NERZ and extends from the highest elevation––13,040 feet at the south end of the Kokoolau quadrangle, just below the summit caldera––to the sea northwest and west of the summit. Lava that erupts from the north and west flanks typically flows to the west, northwest, or north, depending on the vent location. Both morphologic lava flow types—‘a‘ā and pāhoehoe—are present. Pāhoehoe units tend to spread out or widen in low-slope regions, such as in the saddle regions between Mauna Loa and Mauna Kea or between Mauna Loa and Hualālai. In comparison, ʻaʻā flows generally produce narrower flow lobes that have higher relief. This map is the fifth in a series of five maps that will cover Mauna Loa volcano. NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii↗
USGS2023

Connecting flood-related fluvial erosion and deposition with vulnerable downstream road-stream crossings

Fluvial erosion is increasingly responsible for infrastructure and building damages associated with floods as the intensity of extreme rainfalls hit rural and urban rivers in a variety of climate settings across the United States. Extreme floods in 2016 and 2018 caused widespread culvert blockages and road failures, including extensive damage along steep tributaries and ravines in the Marengo River, Wisconsin, watershed during 2016 and 2018. A study conducted by the U.S. Geological Survey (USGS), Wisconsin Wetlands Association (WWA), Ashland County, and the Northwest Wisconsin Regional Planning Commission (NWRPC) investigated the special concern of fluvial erosion hazards (FEHs) associated with gullying, streamside landslides, and the loss of wetland storage in headwaters. In 2019, a pilot study was begun to map and classify ephemeral and perennial streams and wetlands in terms of their sensitivity to FEHs. This study combined data from field-based rapid geomorphic assessments (RGAs) coupled with a stream network-wide geographic information system (GIS) approach for mapping stream segments, referred to as fluvial process zones (FPZ), sensitive to erosion, deposition, and channel change. The GIS approach used nationally available 10-meter (m) resolution topology and an extended stream network to map FPZs based on Strahler stream order, stream power, channel slope, presence of adjacent steep valley sides and headwater flats, and adjacent landform setting. Bankfull channel widths derived from RGA-based hydraulic geometry curves combined with drainage areas, an estimate of bankfull flow, and channel slope were used to calculate specific stream power for the FPZs. Lastly, the FPZs were characterized by their location within three major landform settings that affect erosion potential. The resulting vulnerability maps provided a screening framework to identify FPZs that are sensitive to incision, gullying and mass wasting along steep headwater ephemeral channels, as well as downstream perennial channels that have the potential for valley-side landslides, coarse sediment deposition, and channel change. Lastly, each FPZ was characterized in terms of hydrologic alteration associated with ditching. The vulnerability mapping products and rankings of sensitivity of FPZs will ultimately be used by Ashland County and their collaborators to prioritize natural flood management projects that mitigate FEHs, restore hydrology, and reconnect channels with adjacent wetlands and floodplains.

Wisconsin↗
USGS2026

A deep dive into subduction zones and the Japan Trench

Deep beneath the ocean, giant slabs of rock slide beneath one another in areas called subduction zones. These powerful movements shape the planet by building mountains, creating fiery volcanoes, and even causing huge earthquakes. Scientists study subduction zones to understand why these natural events happen and where and how they might occur so we can be better prepared. One important subduction zone is the Japan Trench, where a big earthquake happened in 2011. To learn more about subduction zones, scientists go on ocean drilling expeditions. They sail on research ships, drill into the seafloor, and collect rock samples. These samples help scientists understand what happened in the past and how subduction zones work. In turn, this information helps communities better understand and prepare for future earthquakes.

Frontiers for Young Minds↗
USGS2026

Key minerals in energy transmission

The electrical grid relies on minerals with unique properties to transform, transmit, and convert electricity. Transmission needs are growing—to support the expansion of data centers and link diverse energy resources into the grid. This expansion will require greater quantities of minerals. Although some minerals are produced in the United States, the United States relies on imports for others. This infographic shows the key minerals that support energy transmission, and how reliant the United States is on imports for these minerals.

General Information Product↗
USGS2026

Apatite (U-Th)/He thermochronology from Marsh Creek anticline reconciles Cenozoic and Holocene strain patterns and elucidates the Cenozoic canning displacement zone of Arctic Alaska

The northeastern Brooks Range is a vast fold-thrust belt that records the northernmost expression of Cenozoic deformation in Alaska, likely with complex dynamic linkages to distant tectonic elements. The principal focus of this manuscript is the Marsh Creek anticline (MCA), which is a large and recently active structure within the broader northeastern Brooks Range. The MCA is a complex, polyphase structure mostly buried beneath the Arctic coastal plain. We present apatite (U-Th)/He (AHe) dates from six outcrop samples, all near the structural crest of the anticline, as well as supporting detrital zircon U/Pb (ZUPb) age spectra from five outcrop samples. The AHe sample transect spans ~3 km of structural relief. The four shallowest AHe samples are not reset. However, the two deepest samples exhibit low intrasample variability among replicates and strong date versus grain radius correlations, suggesting post-depositional burial heating to within the AHe partial retention zone. Inverse thermal history modeling of these data resolves a protracted phase of exhumational cooling that began around 33–28 Ma. Detrital ZUPb age spectra from Paleogene strata support an interpretation of diverse inherited thermal histories recorded by the AHe data, and they are also consistent with extant reconstructions of Paleogene sediment dispersal. Placed into regional context, the results constrain a phase of increasing structural relief in the northeastern Brooks Range fold-thrust belt that began by the Eocene. The relief developed in a zone characterized by diminishing structural relief to the west called the Canning displacement zone. Although this zone was active since the Eocene, several structures including the studied part of the MCA were rapidly exhumed in the Oligocene. Lastly, we argue for kinematic compatibility between Cenozoic contraction and active strike-slip in the northeastern Brooks Range, with the temporal change in strain style likely reflecting the increase in structural relief and topography.

Alaska↗
USGS2024

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano↗
USGS2024

Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network

The Everglades Depth Estimation Network interpolates water-level gage data to produce daily water-level elevations for the Everglades in south Florida. These elevations were used to estimate flow vectors (gradients and directions) and volumetric flow rates using the Gradient-Related Discharge in an Everglades Network (GARDEN) application developed by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Flow rates in both the east-west and north-south directions were computed on a 400-meter square grid using modified parameters in the Manning’s equation. The frictional resistance parameter in the Manning’s equation was calibrated to measured flow rates at coastal creeks fed by Everglades Depth Estimation Network boundary flows. Levees and other features that act as barriers to flow were defined as “no-flow” grid cells where vectors were set to zero. The flow volume magnitudes were calibrated with 2020 daily values of coastal river flows, and verification was performed using 2021 data. Within a given day, the measured coastal river flows fluctuate more than the GARDEN boundary flows because of tidal and wind forcings. Because the GARDEN boundary flows were the upstream water source for the coastal rivers, calibration focused on matching average daily flow volumes rather than daily fluctuations. The Pearson’s correlation coefficient is 0.766 for the 2020 calibration period and 0.566 for the 2021 verification period. Applying GARDEN to periods with hydraulic-control-structure releases allows the propagation of structure flows to be seen in the daily flow-vector maps along with the multiday response of flows farther downgradient. Flow vectors may be overestimated near control structures because of difficulties in resolving the water gradient downstream from the structure. Flow vectors farther from the structure are more accurate than those near the structure.

Florida↗
USGS2024

Effects of drought and cloud-water interception on groundwater recharge and wildfire hazard for recent and future climate conditions, Kauaʻi, Oʻahu, Molokaʻi, Maui, and the Island of Hawaiʻi

The Water-budget Accounting for Tropical Regions Model (WATRMod) code was used for Kauaʻi, Oʻahu, Molokaʻi, Maui, and the Island of Hawaiʻi to estimate the spatial distribution of groundwater recharge, soil moisture, evapotranspiration, and climatic water deficit for a set of water-budget scenarios. The scenarios included historical and future drought conditions, and a land-cover condition where shrubland and forest within the cloud zone were converted to grassland. For the historical drought condition, island-wide mean annual recharge estimates range from a decrease of 30 percent (239 million gallons per day [Mgal/d]) for Kauaʻi to a decrease of 39 percent (2,706 Mgal/d) for the Island of Hawaiʻi, relative to the reference condition consisting of 1978–2007 rainfall and 2020 land cover. For the future drought condition, estimates of island-wide mean annual recharge range from a decrease of 40 percent (477 Mgal/d) on Maui to a decrease of 51 percent (116 Mgal/day) on Molokaʻi. Complete conversion of all shrubland and forest within the cloud zone to grassland for each drought condition produces estimated land-cover-related decreases in island-wide mean annual recharge (in addition to the drought-related decreases) of 11–12 Mgal/d on Oʻahu, 119–135 Mgal/day on Maui, and 689–849 Mgal/d on the Island of Hawaiʻi. The spatial distributions of increases in conditions indicative of moisture stress and potential wildfire hazard were quantified using the relative frequency of soil moisture less than a selected threshold value (monthly mean soil moisture less than 0.074, expressed as a fraction of available water capacity), evapotranspiration less than a selected threshold value (monthly evapotranspiration less than 0.96 inches), and climatic water deficit greater than a selected threshold value (monthly climatic water deficit greater than 0.77, expressed as fraction of potential evapotranspiration). For the historical drought condition, the greatest increases in the relative frequency for the moisture-stress indicators occur across parts of east and southwest Kauaʻi; central, east, and west Oʻahu; central Molokaʻi; central Maui and low- to mid-altitude parts of West Maui volcano; and the northwestern and southern parts of the Island of Hawaiʻi. For the future drought condition, the greatest increases in the relative frequency of the moisture-stress indicators occur across parts of west Kauaʻi; central and west Oʻahu and Molokaʻi; a band of mid-altitude area on the southern slope of West Maui volcano and across the southwestern slope of Haleakalā; and mid-altitude areas of the northwestern and southern parts of the Island of Hawaiʻi. Complete conversion of all shrubland and forest within the cloud zone to grassland for each drought condition results in land-cover-related increases in the relative frequency of moisture-stress indicators around Kaʻala in the Waiʻanae Range and the southeastern part of the Koʻolau Range on Oʻahu, the southern part of West Maui volcano and the southwestern slope of Haleakalā on Maui, and the upland regions of the western and southern parts of the Island of Hawaiʻi.

Hawaii↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗