USGS Science⌕ Search

SEARCH · USGS Science

Search USGS Science

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

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79Linked to original sources
USGS2026

Migration water temperature and heat stress assessments in western Alaska Chinook salmon overlapping the 2019 heatwave

Chinook salmon population declines span their geographic range with climate hypothesized as a major driver. Concerns of warming freshwater temperatures in their northern range gained urgency during 2019 when a heatwave coincided with premature mortality. This study examined heat stress during the 2019 heatwave compared to subsequent years and described water temperatures in western Alaska to understand the degree to which freshwater temperatures may be a stressor. Heat stress was prevalent among Chinook salmon captured in the 2019 heatwave (Kuskokwim tributaries: 90% in Kwethluk and 63% Takotna river), and variable in subsequent years (∼8% to 60% across Kuskokwim tributaries and Norton Sound rivers). A review of water temperature data indicated that potentially stressful temperatures (≥18 °C) were most common and prolonged in the Yukon River, moderately common and prolonged in the Kuskokwim River, and relatively rare in the Norton Sound region. Water temperatures in 2019 broke several records for overall maximum and frequency of temperatures ≥ 18 °C. Migration water temperatures and heat stress in northern Pacific salmon habitats vary more widely than previously recognized (up to 25 °C).

Alaska↗
USGS2026

Water-budget simulations for selected watersheds in Cameron County, Texas, 2022–23

The U.S. Geological Survey, in cooperation with the City of Brownsville, Texas, configured and calibrated a set of hydrologic models for a 217-square-mile study area in Cameron County in south Texas during 2022–23. The models were used for estimating runoff and quantities of water diverted from the Rio Grande/Rio Bravo del Norte (hereinafter referred to as the “Rio Grande”) to maintain water-surface elevations in the canals and resacas (former distributary channels cut off from the main channel of the Rio Grande). Resacas provide habitat to aquatic species and help reduce the effects of flooding. Because of the large size of the study area and diversity of hydrologic conditions, the study area was divided into 11 watersheds, and separate hydrologic models were developed for 9 of the watersheds. Six of the nine modeled watersheds are drained mostly by canals (canal watersheds), and three of the modeled watersheds drain to resacas (resaca watersheds). The Hydrological Simulation Program—FORTRAN was selected for modeling the study area watersheds because it is flexible in simulating a wide variety of watershed conditions. The models were calibrated with streamflow data collected during 2022–23. The calibrated models were used to simulate water budgets (streamflow, evapotranspiration, water-storage volumes, and water diversions and withdrawals) during 2022–23. Model simulations showed that the resaca watersheds required more diversions from the Rio Grande and released less runoff than did the canal watersheds. Management practices maintaining resaca water levels constrained their runoff.

Texas↗
USGS2026

Flood-inundation maps for Río Grande De Loíza in and near Caguas, Puerto Rico, 2026

Digital flood-inundation maps for a 2.7-mile reach of Río Grande De Loíza in Caguas, Puerto Rico, were created by the U.S. Geological Survey. Water-surface profiles were computed for the stream reach by using a one-dimensional, steady-state, step-backwater model. The model was calibrated to the current (2025) stage-streamflow relation (rating curve) for the U.S. Geological Survey streamgage 50055000, Río Grande De Loíza, Puerto Rico. The resulting hydraulic model was then used to compute 16 water-surface profiles for water levels (flood stages) ranging from 19.00 to 34.00 feet at the streamgage; these flood stages range from “moderate flood stage” to above “major flood stage” as defined by the National Weather Service. The 34.00-foot stage exceeds the historical maximum peak stage of 33.20 feet, recorded at the streamgage in 1945. The simulated water-surface profiles were used in combination with a digital elevation model derived from light detection and ranging (lidar) data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood-mitigation strategies and to define flood-hazard areas to help protect life and property, to coordinate flood-response activities such as evacuations and road closures, and to aid post-flood recovery efforts.

Puerto Rico↗
USGS2026

Captivity may broadly impede bison ecological function: Space, season, and management

American bison ( Bos bison ) were spared from numerical extinction but considered to be ecologically extinct, partially due to spatial restrictions in their former range. A key component of the ecological function of bison is their ability to adjust their spatial behaviors at multiple scales to match the seasonal distributions and availability of resources. To understand how management influences the amount and seasonality of space bison use at various scales, we estimated home range using autocorrelated kernel density estimators (95% utilization distribution) and foraging patch size by the maximum variance in first-passage time at four temporal scales (annual, growing, nongrowing, and focal seasons) for two free-ranging and three captive bison herds. The seasonal differences in space-use within management regimes and study areas were treated as nonparametric related measures while differences amongst management regimes and study areas were treated as nonparametric independent measures. Seasonal variations in space-use were detected in both management regimes, but only free-range bison depicted this pattern at the larger home range scale (p < 0.001). Free-range bison had significantly larger home ranges across seasons compared to captive bison (p < 0.001), but captive bison had larger foraging patch sizes during the growing season than free-range bison. Study areas with captive bison did not differ in their amount of space-use at any scale or season (p > 0.05). The strongest differences were apparent among study herds in annual home range size (p < 0.001) and growing season foraging patch size (p < 0.001). Our results demonstrate that captivity likely prevents bison from using their available landscape seasonally at broad scales, which could have ecological consequences for this historically nomadic herbivore. Although sociopolitical realities prevent bison restoration at a continental scale, increasing the amount of space available to bison where possible has the potential to replicate previous ecological processes. Further research can evaluate the influence of local environmental conditions on seasonality of bison space-use and can assess effects of scale on the ecological impacts of bison.

Colorado, Montana, Texas, Utah↗
USGS2026

Geodetic investigations of the Europa Clipper mission

The Europa Clipper mission will investigate the geophysical properties of Europa, one of Jupiter’s Galilean moons, to assess its habitability. Geodetic measurements will play a critical role in determining Europa’s internal structure, rotational state, and tidal deformation. The Geodesy Focus Group (GFG) coordinates cross-instrument efforts to measure Europa’s global shape, rotational parameters, gravity field, and degree-2 tidal Love numbers ( k 2 and h 2 ). The data provided by the Gravity/Radio Science (G/RS) investigation, Europa Imaging System (EIS), Radar for Europa Assessment and Sounding (REASON), and Europa Ultraviolet Spectrometer (UVS) will be used to infer geodetic constraints on the interior structure and construct a precise cartographic reference system. This combined dataset will provide new constraints on Europa’s tidal response, ice shell thickness, and the properties of its subsurface ocean. The resulting geodetic information will contribute to a deeper understanding of Europa’s internal dynamics and the potential habitability of its ocean.

Space Science Reviews↗
USGS2026

Large streamflow differences between forested and urbanized watersheds in the energy-limited eastern United States: The role of evapotranspiration and impervious surfaces

Urban forests and other green infrastructures have been viewed as part of the “Nature-based Solutions” (NbS) to mitigate emerging urban environmental change. This study focuses on the role of evapotranspiration (ET) in regulating water balances of small watersheds in the eastern United States. We compared streamflow and ET patterns at daily, monthly and annual scales and linked these hydrological variables to the physical properties of 11 paired watersheds dominated by forests (FW) or urban (UW) land covers. The annual precipitation ranged from 1028 mm to 1683 mm and potential ET (PET) from 815 mm to 1450 mm. The mean annual flow/precipitation (Q/P) ratios were 0.26 ± 0.13 and 0.41 ± 0.1 for FW and UW, respectively. Overall, UW had lower annual ET (772 mm in UW vs. 947 mm in FW), but higher mean annual and (∼58% higher), monthly water yield (17%–186% higher), and peakflow rates (up to 100 times higher) than FW. The streamflow differences between FW and UW were most pronounced during the growing season and early winter (June-November). The mean Q/P ratios for 30 large hurricane events (2016–2021) were 0.12 ± 0.11 and 0.38 ± 0.23 for FW and UW, respectively. The flow rates in the dormant season (around December-May) in UW were similar or lower than FW. We developed conceptual models to explain the seasonal and storm event streamflow differences using background climate (PET), ET, and land surface characteristics. Urban NbS designs should factor in strategies that maximize ET while minimizing impervious surfaces enhancing watershed “sponge” and “pump” functions.

eastern United States↗
USGS2026

Rapid hydrothermal triggering of induced seismicity at the Coso Geothermal Field

The long-term producing Coso Geothermal Field (CGF) in California operates over 100 wells tapping into a reservoir characterized by an extensive fracture network, complex fluid pathways, and regular seismic activity. Understanding the interaction between seismicity and injection can shed important light on the hydrothermal characteristics of the field. Here, we analyze 15 years of local seismic and daily operational data from the CGF, identifying a strong correlation between short-term increase in seismicity rate and seasonal volumetric and temperature variations in the reinjected fluid. Furthermore, the seismic footprint during peak injection of colder fluids reveals a near-instantaneous response up to 2 km away from the injection well, too rapid for pore pressure diffusion alone. This short-term and distant response is observed to have directional preference, indicating structural or permeability anisotropy within the reservoir. Additionally, the seismic response correlates with the initial volumetric increase of colder fluids, but also with temperature decrease during stable injection periods, suggesting thermal effects alone can play an important role in triggering distant seismicity.

California↗
USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗
USGS2026

Cotton farming affects ileal virome in a sedentary wild passerine

Although a few studies have focused on avian gut virome variation in response to environmental stressors, none have assessed virome in relation to the production of chemically intensive crop-based agriculture that alters food resources and detrimentally affects various aspects of avian health and fitness. In this study, we used shotgun metatranscriptomics to assess whether exposure to cotton ( Gossypium spp.) production had a deleterious effect on the ileal virome of sedentary northern mockingbirds ( Mimus polyglottos ) sampled from two cotton-producing areas (16 birds in total) and one uncultivated area (7 birds) in Texas, USA. We recovered 43 viruses representing 13 virus families, which included two viruses that appear to be potential vertebrate pathogens. Individual sample richness varied from 25 to 33 viruses. Both virome richness (Adj. r 2 = 0.247, F (2, 20) = 4.615, P = 0.022) and composition (r 2 = 0.370, F (2, 20) = 5.883, P = 0.001) differed among three sampling regions. Cotton production was associated with the increase of virome richness (Adj. r 2 = 0.283, df = 22, P = 0.005). Pesticide occurrence data collected using silicone bands at the three sites suggest that virome compositional changes are not only associated with total pesticide exposure but are also particularly sensitive to the pesticide combinations detected at each location.

Texas↗
USGS2026

Rupture into slow-slip fault regime during the 2018 Mw 6.9 Island of Hawaiʻi earthquake is followed by modest postseismic slip

On 4 May 2018, a M w 6.9 earthquake occurred on the south flank of Kīlauea, in the midst of an historic event that included a voluminous eruption from Kīlauea’s lower East Rift zone and caldera collapse at its summit. The earthquake was a consequence of both short‐ and long‐term stress buildup due to magmatic activity associated with the eruption and steady flank motion, respectively, and it revealed features of Kīlauea’s décollement fault that can inform understanding of future earthquake activity. We used geodetic data to determine the distributions of slip during the coseismic and postseismic periods and compared these with areas of known fault slip during past earthquakes and slow‐slip events (SSEs). The 2018 earthquake ruptured into an area of the décollement fault that was active during quasi‐regular SSEs that occurred in the two decades prior to 2018 but that have not been observed since. The coseismic slip model indicates that the amount of motion on the décollement fault was several times greater than what typically occurred during SSEs, suggesting that it may take decades for the fault to rebuild stress to the point at which SSEs will occur again. Postseismic afterslip also occurred in an area of the fault known to experience slow slip; however, unlike at other creeping faults, postseismic afterslip was rapid, being largely over within 2–3 days. The rapid nature and small magnitude of the postseismic afterslip may be due to the lack of a viscoelastic relaxation component, which is possibly a result of the shallow dip of the décollement fault not transferring stress efficiently into the lower crust.

Hawaii↗
USGS2026

Seasonal migrations and other movements

In the past 25 years new information has been obtained on the migrations and movements of mountain sheep (bighorn [Ovis canadensis], thinhorn [Ovis dalli]). This chapter provides a comprehensive overview of mountain sheep migration and other movements across their broad distribution in western North America. Across the range of mountain sheep, migrations and other seasonal movements define a complex movement portfolio that relates to and supports all aspects of mountain sheep ecology. Bighorn and thinhorn species have examples of migratory behaviors that span the continuum between annual residents and long-distance migrants. Migratory behaviors can be characterized as low- or high-elevation residents; elevational migrations and other variants; and geographic migrants. Native populations that have been extant on the landscape without notable human intervention have greater migratory propensity and more diverse migratory portfolios due to the maintenance of migration through cultural learning and social transmission. Restored and augmented populations, where the population-level knowledge of migration has been lost or greatly reduced, are largely nonmigratory, although translocations show some ability to restore short-distance elevational migrations. Seasonal spring and fall migrations are less common in desert bighorn sheep (O. canadensis spp.) or bighorn sheep living in canyon or prairie breaks landscapes. Mineral lick visitation is important and common across the range of mountain sheep. Managing for migratory diversity can help to sustain migratory behavior in the face of climate change and other anthropogenic pressures, which can limit landscape connectivity between seasonal ranges or alter the spatiotemporal dynamics of regional phenology with cascading effects to other biotic and abiotic interactions such as the need to balance forage and predation risk or the spatial refugia with increased temperatures.

Book chapter↗
USGS2026

New methods provide a 300–year perspective on modern area burned in two wilderness areas of the southwest United States

Climate change, expanding human ignitions, and increased fuels from fire exclusion are driving increases in area burned and fire severity in dry conifer forests of the western United States. Increasing area burned is occurring against the backdrop of a large fire deficit caused by over a century of fire exclusion. A key land management question is whether historically frequent fire regimes can be restored. Accurate estimates of historical annual area burned (prior to circa 1900) are necessary to evaluate modern area burned (after circa 1900), but are difficult to derive, and have rarely been calibrated or validated against modern fires, leaving their accuracy uncertain. We developed new methods to use tree-ring fire scars to reconstruct historical annual area burned and compare it to modern annual area burned. We focused on two southwestern US wilderness areas—Saguaro National Park (SAGU) and the Gila Wilderness (GILA)—that have a long history of using prescribed and managed fires. The abundant modern low- and moderate-severity fires allowed us to (1) calibrate and validate the fire-scar models against mapped fires to derive the first uncertainty estimates of reconstructed annual area burned and (2) test whether active fire management can help restore annual area burned to historical levels. A multi-model ensemble consisting of 10 individual member models accurately estimated area burned of mapped modern fires with no consistent biases. Each member model had distinct strengths and assumptions that made them suitable for specific applications (e.g., the synchrony model is easily applied, and Thiessen polygons provide spatially explicit area burned estimates). The accurate reconstruction of modern area burned from relatively sparse fire-scar data at GILA suggests that dense grids may not be necessary for accurate reconstructions. Our findings reveal that despite the near absence of fire in the early 20th century, both annual and 20-year sums of area burned in recent decades are back within historical levels at GILA, and trending toward historical levels at SAGU. These results demonstrate that fire management can help restore the historically prevalent, ecologically important process of widespread, frequent, low-to-moderate-severity fire in dry conifer forests.

Arizona, New Mexico↗
USGS2026

Assimilation efficiency of rainbow trout fed natural diets

Assimilation efficiency is a critical assumption of stable isotope mixing models and bioenergetics models, yet few studies examine how assimilation efficiency influences modeling inferences. We conducted laboratory experiments to determine rainbow trout ( Oncorhynchus mykiss ) assimilation efficiencies. Assimilation efficiency averaged 55.8% (SE ± 0.90) and 64.5% (SE ± 1.98) at the 10% and 25% ration levels, respectively, and did not differ significantly. Caloric energy egested in feces was not significantly different between ration levels. Caloric energy excreted in ammonia and urea was significantly different between the ration levels, with a higher amount at the 25% ration level. Absorption efficiency was significantly greater at a higher ration level. Percent energy egested in feces was significantly greater at lower ration levels. Percent energy excreted as ammonia and urea was high compared to previous studies and did not differ significantly by ration level. Our estimates of assimilation efficiency of rainbow trout (56% and 65%) were lower than some previously reported estimates, and our estimates of energy losses (feces and ammonia) were higher than some previous estimates. As knowledge of species-specific assimilation efficiencies increases, our ability to draw strong inferences and improve the accuracy of model predictions will improve.

Environmental Biology of Fishes↗
USGS2026

Invasive grass influences on the fire cycle and treatment effectiveness to control their abundance in the Intermountain West, USA

Many non-native invasive grass species increase wildfire activity and regenerate more quickly than native species. This invasive grass–fire cycle has severe negative consequences for ecosystems, creating a need to understand how different invasive grass species alter fuel characteristics and fire behavior, as well as effective treatments to control their abundance. To address these needs and increase fire and natural resource management preparedness, we performed a review and meta-analysis of recent (1985 to 2023) scientific literature. We focused on the Intermountain West, USA, where six dominant invasive grass species have already transformed ecosystems, including winter annuals—cheatgrass ( Bromus tectorum L.), medusahead [ Taeniatherum caput-medusae (L.) Nevski], red brome ( Bromus rubens L.), and Mediterranean grass [ Schismus arabicus Nees and Schismus barbatus (Loefl. ex L.) Thell]; and summer perennials—buffelgrass [ Pennisetum ciliare (L.) Link] and Lehmann’s lovegrass ( Eragrostis lehmanniana Nees). Within the 204 selected articles, B. tectorum was the most well-studied species, treatment effectiveness was the most common study type, and more studies addressed fuel accumulation than fire characteristics. While initial reductions in B. tectorum following wildfire were followed by large increases, P. ciliare initially increased and then steadily declined, and other invasive grass species had no significant post-fire changes over time. Chemical treatments were more effective than other treatments for B. tectorum , P. ciliare , and Schismus spp., although T. caput-medusae had a greater reduction with chemical treatments compared with the other species. In many cases, treatment effectiveness was enhanced when treatment types were combined or repeat treatments were conducted. Both B. tectorum and T. caput-medusae increased to pretreatment conditions within 3 and 5 yr, respectively, although there were no detectable trends for other species. Our results provide comprehensive comparisons of the effect of invasive grass species on fuel and fire characteristics and much needed insight on effective strategies for reducing invasive grass impacts to ecosystems.

Intermountain West↗
USGS2026

Simulated ground motion dataset in the Azores Plateau, Portugal, on rock and soil sites

Building on a previously developed bedrock dataset, this study extends the Azores Plateau ground motion simulations to include soil-amplified records and introduces a comprehensive validation framework. Soil amplification is modeled using one-dimensional soil profiles. A stochastic source-based approach is employed to generate the dataset, incorporating randomization of input-model parameters to account for the aleatory uncertainty in seismic activity. The accuracy of the dataset is verified through a comprehensive validation framework, showing that the randomization effectively captures variance and inter-period correlation observed in records. This work provides a robust dataset for advancing seismic hazard and risk assessment in the Azores Plateau.

central and eastern Azores islands↗
USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗
USGS2026

Statistical facilitation in environmental science: Integrating results from complementary statistical analyses can improve ecological interpretations

Professionals working in biological conservation seek to understand, manage, and restore populations of native organisms using many techniques. A common approach for this discipline is using long-term data collections to inform decision making. However, several quantitative issues complicate statistical analysis of monitoring datasets and can reduce the utility of results for conservation decision making. Integrating results from multiple analyses applied to the same dataset (i.e., approaching the same biological problem using different techniques) is one way to address concerns related to field data that violate statistical assumptions. This process allows data analysts, researchers, and managers to assemble insights based on the weight of evidence. Here we tested whether three different statistical techniques [(1) multiple logistic regression on original data, (2) multiple logistic regression on standardized data (i.e., mean of 0 and standard deviation of 1), and (3) random forest analysis] identified a similar hierarchy for selecting natural and anthropogenic habitat regressors. Our examination of how environmental variables affected Plains Minnow ( Hybognathus placitus ), a state-threatened fish, is relevant to other taxa and locations. We gained useful information from redundancies (i.e., agreements across analyses). New directions also emerged by addressing ambiguities (i.e., disagreements among results across analyses). When multiple analyses were integrated into one ecological story, a clearer interpretation emerged. Viewing different statistical tests as facilitators that provide mutual advantages can advance the understanding and application of statistical analyses applied to non-experimental field datasets.

Kansas↗
USGS2026

Origins, evolutions, and future directions of Landsat science products for advancing global inland water and coastal ocean observations

In April 2020, the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center introduced a Level 2 provisional Aquatic Reflectance (AR) product for the Landsat 8 Operational Land Imager (OLI), marking the initial phase in developing a standardized global product for Landsat-derived surface water measurements. The goal of USGS EROS aquatic product research and development is to prepare for an operational processing architecture for Landsat Collection 3 in the late 2020s that will enable use of quality-controlled data for emerging Landsat aquatic science applications. To achieve this, we released a subset of the Landsat 8/9 provisional AR products (Crawford et al., 2025, https://doi.org/10.5066/P14MBBRM) and examined its general performance through the Science Algorithms to Operations (SATO) framework alongside quantitative assessment using community made inland water data records (GLObal Reflectance community dataset for Imaging and optical sensing of Aquatic environments, GLORIA) and radiometric coastal validation platforms (NASA’s Ocean Color component of the Aerosol Robotic Network, AERONET-OC). Variability within the validation datasets indicate that the performance of the Landsat 8/9 provisional AR retrieval is highly context-dependent; errors are minimal in optically simple waters (e.g., clear to moderately turbid coastal waters) but increase considerably in optically complex waters where factors such as elevated levels of turbidity, chlorophyll (Chl a ) concentrations, or colored dissolved organic matter (CDOM) dominate the water column. Additionally, this paper examines key algorithmic considerations for atmospheric correction, highlighting factors that influence accuracy, scalability, and computational efficiency necessary for collection processing in the operational Landsat Product Generation System (LPGS). This paper is intended to communicate with aquatic scientists, satellite oceanographers, and the broader Earth observation community on the origins, requirements, challenges, successes, and future objectives for operationalizing global AR data products for Landsat satellite missions.

Earth System Science Data↗