USGS Science⌕ Search

Geology topics

Browse USGS reports

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

At least 865 records · Page 48Linked to original sources
USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗
USGS2026

Tropicalization of the temperate zone: Spatiotemporal variability of winter warming and declining freeze days across the United States

We investigate changes in cool-season and winter daily minimum (T min ) and maximum (T max ) temperatures, and the occurrence of freeze days, from 1952 to 2024 across the conterminous United States (CONUS). Emphasis is placed on the tropical-temperate transition zone (TTTz) in the southeastern CONUS. During winter, ~70% of the land area exhibited T min warming rates exceeding those of T max . The countywide coldest T min became milder across 57% of the CONUS, while the coldest T max showed little change and even cooled east of the Rocky Mountains in the central CONUS. Across the TTTz, 75% of freeze days occur within a ~25–100-day window, often fewer than 75 days in the southernmost areas. Approximately 80% of counties exhibited significant contractions in freeze-day concentration, with the largest and most spatially consistent changes occurring in the Southeast, primarily driven by later start dates. Roughly 85% of the CONUS experienced a significant decline in freeze days, with the largest relative declines in regions where average winter T min is above freezing, while parts of the Pacific Northwest showed no significant change. An analysis of freeze day isopleths (30, 45, 60 and 75 days) across 20-year periods showed that the mean latitude of freeze days has migrated poleward substantially. Between 101° W and 79° W in the TTTz, the 30 freeze-day isopleth for the late period (2005–2024) was, on average, 122 km (~1.1° latitude) farther north than in the early period (1952–1971). Generally, the largest latitudinal shifts and percentage losses in freeze days occurred across low-elevation, low-relief regions at lower latitudes (e.g., the Mississippi River Valley), with abrupt shifts occurring near topographic gradients. Regions with sharp elevational gradients (e.g., Balcones Escarpment, Ouachita Mountains and Tennessee Valley) exhibited smaller temporal changes, likely reflecting the barrier-like influence of higher terrain on the poleward retreat of freeze days.

conterminous United States↗
USGS2026

Accounting for emigration reveals high survival and bimodal size at departure from a loggerhead sea turtle (Caretta caretta) foraging area

The life history of hard-shelled sea turtles includes several ontogenetic shifts in habitat use and these complex permanent emigration patterns can impact estimates of stage-specific population rates, including survival. We developed several multistate mark recapture models to estimate survival of adult and juvenile loggerhead turtles from a coastal bay in the northern Gulf of America (also commonly referred to as the Gulf of Mexico) while, in some cases, accounting for permanent emigration and transient individuals. Our mark-recapture dataset consisted of 228 individual turtles with 37 total recaptures from 2011 to 2024. Of the models we fit, those that incorporated emigration produced higher estimates for annual survival than models that did not, and higher estimates than what is commonly seen in the literature for loggerheads. All models suggested a major permanent emigration pulse at the typical size of sexual maturity (70 cm straight carapace length) and another major pulse at > 90 cm. This bimodal pattern of departure may reflect differences in size at sexual maturity among loggerheads, possible genetic variability within the assemblage, or both. To assess the models’ ability to effectively recover true parameter values, we developed a simulation study of 50 randomly generated independent data sets under our specified models of similar sample size to our study dataset. Simulation results suggested that models that accounted for permanent emigration and transient individuals produced relatively unbiased estimates of survival, while models that did not often underestimated survival rates. Mark-recapture studies that may exhibit emigration and suffer from low recapture rates would benefit from auxiliary data collection such as acoustic telemetry detections to better estimate true rates of emigration and survival. Obtaining unbiased estimates of true survival by accounting for processes like emigration can support effective conservation of endangered long-lived species like loggerheads.

Florida↗
USGS2026

Patterns of floodplain forest mortality and recruitment along the Upper Mississippi and Illinois Rivers: Associations with forest fragmentation and flood inundation

Context Different rates of floodplain forest recruitment and mortality can reveal important changes in ecosystem processes that drive forest dynamics, resulting in net changes in forest cover, thereby influencing a wide range of river habitat and morphological characteristics. Objectives We evaluated characteristics of forest change areas in the Upper Mississippi River System. Methods An overlay technique was used to map patches of forest loss, gain, and persistence between 2010 and 2020 in relation to a series of explanatory variables. Results We quantified a net decline in forest cover ranging from 3.2 to 16.8% in the uppermost five study reaches, and a net increase in forest cover ranging from 0.5 to 4.6% in the southernmost three reaches. Patches of forest loss and persistence were similarly tall (> 15 m), dense (> 90% cover), silver maple ( Acer saccharinum ) dominated forests, whereas forest gain patches were short (< 15 m), less dense (< 66% cover) and more likely to be dominated by willow ( Salix) species. Both forest loss and gain patches were smaller than forest persistence patches and were typically found in areas with low neighborhood forest density (< 50% forested 10 ha neighborhood). Areas that experienced more than three flood events per growing season, more than 100 consecutive days of inundation during a single flood event, and more than 60 mean total days of inundation per growing season from 2011 to 2020 showed a net loss of forest cover in all study reaches. In contrast, net increases in forest cover were restricted to areas that experienced less than a single flood event per growing season, less than 40 consecutive days of inundation during a single flood event and less than 30 mean total days of inundation per growing season from 2011 to 2020. Conclusions Forest mortality along these river reaches is associated with forest fragmentation and an increasingly wetter hydrological regime.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗
USGS2026

Refinement of a framework for Moving Aircraft River Velocimetry (MARV) and application to particle tracking along Alaskan rivers

Information on river velocities enhances understanding flood hazards, evaluating habitat conditions, and predicting the transport of floating materials. In this follow-up study, we used data from two new sites, one with a more complex morphology and the other with a lower suspended sediment concentration, to provide further evidence that Moving Aircraft River Velocimetry (MARV) can yield accurate velocity estimates ( R 2 up to 0.87 when compared to field measurements) for long segments of large, turbid rivers. The MARV workflow is packaged in freely available software and is robust to implementation details; neither buffering to mitigate edge effects nor a new approach to aggregating velocity vectors improved performance. MARV was not sensitive to parameters used to establish overlapping image sequences, but combining a long window with a short jump between consecutive windows was the optimal configuration. Although accuracy varied from one cross section to the next, agreement between remotely sensed velocities and those measured in the field was independent of position within a frame range. As an initial step toward application of the approach to help address practical problems, we showed how MARV can drive particle tracking models. Our first-order simulations suggest that channel morphology and flow velocity are the primary controls on travel time and particle fate, with diffusive processes playing a lesser role. Although MARV can be used to characterize an instantaneous flow field, a more comprehensive framework that accounts for other physical processes would be required to model specific types of events like oil spills.

Alaska↗
USGS2026

Meso-scale pressure reactor demonstrates biostimulation of coal-dependent methanogenesis

Replicating deep subsurface conditions remains a major challenge for advancing subterranean biotechnologies. Microbially enhanced coal bed methane production offers a promising approach to increase biogenic gas recovery from coal seams, yet experimental progress has been limited by the difficulty of simulating in situ conditions. A key question is whether nutrient amendments stimulate coal-dependent methanogenesis or simply convert injected nutrients to methane. Here, we report the use of a large-scale (325 L) reactor to investigate 13 C-labeled microalgae amendments on microbially-enhanced coal bed methane production at in situ pressure. Labeled methane was detected only during the initial stimulation phase, while most of the methane produced over five months originated from the coal. These field-relevant results demonstrate that microalgae can accelerate early methanogenesis and enhance coal-derived methane production over time. This study provides a critical proof-of-concept with broad implications for advancing subsurface biotechnologies toward field-scale deployment.

Communications Earth & Environment↗
USGS2026

Variability and consistency in wildfire susceptibility: Insights from a national compilation

Background Wildfire risk in the United States is rising and remains a land management priority. The quantitative wildfire risk assessment (QWRA) framework integrates fuels, topography, weather and values at risk to estimate the potential change in value from wildfire. Within this, response functions (RFs) represent how values respond to fire intensity. These are often based on expert judgment, but variation across assessments is unclear. Aims This study uses data from the US Geological Survey (USGS) Wildfire Hazard and Risk Assessment Clearinghouse to characterize consistency and variation across categories and contexts. Methods We applied descriptive statistics to summarize RFs, using tables, box-and-whisker plots and heat maps stratified by highly valued resource or asset (HVRA) category and spatial scale. Key results RFs and value definitions vary, especially for ecosystem-related resources. Some functions, such as for buildings in the wildland–urban interface (WUI), translate well across contexts, while others require more input. Conclusions Some functions are broadly transferable, while others need customization. This analysis provides references and starting points for improvement to RFs in QWRAs. Interpretations Expanding the clearinghouse and dataset and building more transparency in expert elicitation can build trust among communities, agencies and end-users, and can support efficient use of limited resources to mitigate wildfire risk.

International Journal of Wildland Fire↗
USGS2026

Oak masting breaks down in a highly fragmented, pollen limited landscape

Masting, the production of variable and synchronous seed crops, is common among wind pollinated trees, and periodic pollen limitation is thought to be an important component of the process of masting. However, little is known either about the dynamics of mast seed production or pollen limitation in isolated trees. We tracked seed production of valley oak ( Quercus lobata ) over 11 years in a highly fragmented landscape in California’s Central Valley in which many trees lacked conspecific neighbors within 500 m. We found that oaks with closer neighbors produced more acorns on average than oaks with fewer neighbors. Consistent with theories linking mast seeding and pollination, the isolation effect was most apparent in the highest seed production years. While individual trees showed substantial interannual variation in seed production, we observed little synchrony in seed production among individuals compared to intact oak populations throughout California. Our results indicate that highly fragmented populations suffer signs of pollen limitation and do not produce synchronous seed crops. Fragmentation, alongside climate change, is a dimension of global change that interrupts the typical seed production dynamics of trees.

California↗
USGS2026

PDCD-DAT: A global database of pyroclastic density current deposit field data

Pyroclastic density currents represent one of the deadliest hazards posed by active volcanoes. Analysis of their deposits provides valuable insights into their internal dynamics and informs numerical simulations of pyroclastic density currents which underpin many volcanic hazard assessments. We present PDCD-DAT, a global database of pyroclastic density current deposit characteristics compiled from peer-reviewed literature. The database includes both quantitative datasets (e.g., grain size, density, bedform dimensions, thickness) and qualitative descriptors (e.g., sedimentary structures, lithofacies). PDCD-DAT includes data from 85 source publications, covering 97 eruptions or eruptive phases, and 214 individual depositional units from 55 globally distributed volcanoes. Eruptions recorded in the database range from VEI 1–8. We highlight examples of potential applications of the database, which include (i) comparison of single deposit case studies to global datasets, (ii) informing input parameters and conditions for numerical and analogue models of pyroclastic density currents, (iii) validation of numerical and analogue models against a wide variety of natural case study deposit architectures, and (iv) estimating hazard impact metrics of pyroclastic density currents from past eruptions. We show that the database represents a useful tool for improving our ability to model pyroclastic density currents, predict their associated hazards, and understand the relationships between the internal dynamics of pyroclastic density currents and the properties of their deposits. PDCD-DAT is integrated with the FlowDat Mass Flow Database, which provides a sustainable platform for the database. We aim for PDCD-DAT to be expanded in the future through addition of pyroclastic density current deposit datasets from new field studies conducted by the volcanology research community.

Journal of Applied Volcanology↗
USGS2026

Not so fatal attraction: Captive female Burmese Python lures do not improve wild python detection

The ongoing invasion of Python bivittatus (Burmese Python; henceforth, Python) across the Greater Everglades Ecosystem (GEE) has led to near total collapse of the affected mammal community over the past few decades. Management efforts to eliminate Pythons and control their spread have been hampered by the Python's low detectability, which may be improved by using a lure. In controlled settings, male Pythons show an attraction to the scent from reproductively active females. To test the effectiveness of using reproductively active female Pythons as a lure for attracting wild male Pythons in the field, we conducted a paired experiment with wild-caught female Pythons in pens and empty control pens. We monitored Python visitation at all sites using camera traps, which resulted in >3,000,000 photographs that we filtered to 4 independent detections of Pythons using AI software. Python detection was low at sites with female Pythons (3 observations) and control sites (1 observation) over 90 days at 12 sites. Stress associated with captivity may have halted reproductive females from producing pheromones, eliminating the chemosensory cue that lures males. Identifying and implementing husbandry techniques to reduce stress in wild-caught female Pythons could improve the effectiveness of this technique. Little is currently known about the chemical ecology of Pythons, and pheromonal communication in particular, and further research in this area could aid in the identification and production of effective, low-cost lures to increase detection and removal of this invasive species.

Florida↗
USGS2026

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗
USGS2026

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S. Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

New Hampshire↗
USGS2026

Hydrogeologic framework and conceptual groundwater-flow model of the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, 1998–2022

This study was conducted by the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, to update the hydrogeologic framework and conceptual flow model for the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, which together compose the Ogallala aquifer focus area. The study included the construction of a potentiometric surface, and available geologic and hydrologic data were used to evaluate saturated thickness of the aquifer. The water budget for the updated conceptual groundwater-​flow model was based on estimated inflows and outflows for the 1998–2022 study period. Saturated thickness of the Ogallala aquifer averaged 127 and 116 feet for the panhandle and northwest parts, respectively. Groundwater withdrawals from the Ogallala aquifer for 1998–2022 averaged 422,054 and 39,645 acre-feet per year (acre-ft/yr) for the panhandle and northwest parts, respectively. Recharge, the primary inflow, was estimated at 0.63 inch per year for the 1998–2022 study period, with the panhandle part of the Ogallala aquifer receiving 175,068 acre-ft/yr and the northwest part of the Ogallala aquifer receiving 49,376 acre-ft/yr. Additional inflows included irrigation return flows, estimated at 8,111 and 642 acre-ft/yr for the panhandle and northwest parts, respectively, of the Ogallala aquifer. Net lateral groundwater flows, considered to be aquifer outflows, were estimated to account for 31,908 acre-ft/yr for the Ogallala aquifer focus area. Streambed seepage, which was an outflow of 5,535 acre-ft/yr, was only present in the northwest part of the Ogallala aquifer. Vertical leakage and saturated-zone evapotranspiration were considered negligible outflows. These findings provide a revised conceptual groundwater-flow model water budget for the Ogallala aquifer focus area in Oklahoma.

Colorado, Kansas, Oklahoma, Texas↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗
USGS2026

Space between houses influences movement and habitat selection of ungulates: Width as a novel metric of development

Wildlife often lose access to habitat due to housing development. The magnitude of indirect habitat loss can be conditional on the configuration of individual houses, but commonly used metrics (i.e. density or distance) can overlook the configuration of development. We introduce a novel framework to index the configuration of development based on the width of space between houses and associated structures. We use resource selection functions to assess the degree that GPS-collared elk ( Cervus canadensis ), mule deer ( Odocoileus hemionus ), pronghorn and moose ( Alces alces ) on winter range and on migration routes use space between houses within northwest Wyoming, USA, near the towns of Cody and Jackson. Further, to help inform regulations aiming to promote wildlife movement across a gradient of land uses, we differentiated between individuals residing in primarily rural and exurban areas. Rural populations of elk, mule deer, pronghorn and moose avoided spaces narrower than 2 km and never used spaces narrower than 50 m between houses, whereas exurban populations of elk, mule deer and moose selected for spaces narrower than 2 km but avoided spaces narrower than 50 m. We identified cutoffs in rural and exurban areas where space may become too narrow for most animals to use. Through this metric, managers and policy makers can inform the necessary width to maintain wildlife movement through corridors. Our width metric can be applied to other systems and our workflow is publicly available ( https://wildlifemovetools.org/width-calculator ) so users can estimate the width of space between structures in their focal areas. Synthesis and applications . In rural areas, maintaining spaces between houses >500 m will likely facilitate ungulate movement. To fully conserve functional habitat, such as unimpeded habitat use along migration corridors, maintaining spaces approximately >2.5 km between houses will likely be necessary.

Wyoming↗
USGS2026

Lava flow emplacement dynamics during the Mauna Loa 2022 eruption

On November 27, 2022, Mauna Loa erupted for the first time in 38 years, threatening a major highway on the Island of Hawaiʻi. The opening phases of the eruption were focused at the volcano’s summit, before magma propagated into the upper Northeast Rift Zone to erupt at four fissure segments that localized at fissure 3 by December 2. This fissure produced a braided ʻaʻā flow that traversed the steep north flank at >1 km/hr, before slowing as slopes shallowed. Effusion rates peaked in the early hours of the eruption and then steadied at around 100 m3/s at fissure 3. Lava effusion became variable and dropped abruptly on December 7. The eruption ended on December 10 with a final flow length of 18 km. We tracked lava effusion rates, flow advance, channel dynamics, and 3-D flow field development with field and remote sensing observations throughout this eruption to inform lava flow forecasts and investigate ʻaʻā flow dynamics. Flows were mapped using optical and infrared imagery, uncrewed aircraft systems, and airborne single-pass InSAR. Video- and time-lapse cameras were deployed to record lava fountaining and channelized flow at frequent intervals, supplemented by field measurements. Analyses of these data built a detailed record of effusion rates and associated flow dynamics. We combine these high frequency observations of the proximal flow field with broader thermal and 3-D mapping of the evolution in the medial and distal flow field to examine ʻaʻā flow dynamics from source to toe. Within the steep, medial section of the flow field, channel blockages caused back-ups and overflows as the eruption progressed. Disrupted lava supply to the distal flow field and low slopes resulted in stalling, widening, and thickening at the toe. Overall, 142 million m3 of lava covered 36 km2 over 13 days. Our results inform the controls on development and evolution of lava flow fields, with implications for forecasting primary flow advance and secondary flow field growth.

Hawaii↗