USGS Science⌕ Search

USGS · 70178095

Integrated modeling approach for fate and transport of submerged oil and oil-particle aggregates in a freshwater riverine environment

Abstract

The Enbridge Line 6B pipeline release of diluted bitumen into the Kalamazoo River downstream of Marshall, Michigan, U.S.A., in July 2010 was one of the largest oil spills into freshwater in North American history. A portion of the oil interacted with river sediment and submerged requiring the development and implementation of new approaches for detection and recovery of oil mixed with river sediment. Hydrodynamic and sediment transport modeling became an integral part of containment and recovery operations for decision support about the potential fate and migration of submerged oil and oiled sediment. Three models were developed for the U.S. Environmental Protection Agency to cover a range of spatial scales of interest to onsite operations. Two-dimensional (2D) hydrodynamic and sediment transport models from the Environmental Fluid Dynamics Code and the sediment bed model SEDZLJ1 were used to simulate potential resuspension, migration, and deposition of submerged oil and oiled sediment along a 38-mile reach of the Kalamazoo River affected by the oil from Marshall to Kalamazoo. An algorithm was added to SEDZLJ to represent three additional particle size classes of oilparticle aggregates (OPAs) with a range of sizes, specific gravities, and settling velocities. Field and laboratory experiments and flume tests were done to support the numerical modeling of OPAs. A three-dimensional hydrodynamic model was developed to simulate hydrodynamics and OPA tracking through Morrow Lake, the most downstream impoundment. This model incorporated wind and dam operations into high and low flow, lake drawdown, and containment simulations. Finally, a 2D unstructured grid model, HydroSed2D, was used to simulate flows and sediment transport along 1- to 2-mile segments of the Kalamazoo River around islands and through side channels and backwater areas that are particularly prone to submerged oil deposition. Integrated models could be developed quickly due to the availability of information and services combined with spill response operations that included: bathymetry and topography data, fieldbased geomorphic mapping of submerged oil, and discharge measured at U.S. Geological Survey streamflow gauges. Modeling results were included in a multiple-lines-of-evidence approach that was used by the Federal On-Scene Coordinator and operations staff for decision-making related to assessment and recovery of submerged oil, as well as net environmental benefit analysis. Similar modeling approaches will likely be useful for future oil spills in riverine environments.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Faith A. Fitzpatrick, Rex Johnson, Zhenduo Zhu, David Waterman, Richard D. McCulloch, Earl Hayter, Marcelo H. Garcia, Michel C. Boufadel, Timothy Dekker, Jacob S. Hassan, David T. Soong, Christopher J. Hoard, Kenneth Lee. 2016. Integrated modeling approach for fate and transport of submerged oil and oil-particle aggregates in a freshwater riverine environment. https://pubs.usgs.gov/publication/70178095

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Authoritative knowledge infrastructure for GeoAI: Implications for scientific spatial data providers

Geospatial Artificial Intelligence (GeoAI) increasingly relies on automated analysis of authoritative geospatial data, exposing limitations in data infrastructures originally designed for simple dissemination rather than for knowledge-driven use. From the perspective of scientific spatial data providers, the challenge is enabling GeoAI responsibly while retaining semantic clarity, provenance, reproducibility, and governance. This short paper (i) examines how GeoAI exposed limitations within existing geospatial data infrastructures, and (ii) identifies provider-side requirements for an authoritative knowledge infrastructure. Drawing on GIScience, scientific workflow, and operational AI literature, the paper outlines implications for national mapping organizations that seek to support advanced analytics without jeopardizing scientific rigor or stewardship responsibilities.

Conference Paper↗

Springs of the Catskill Mountains: A conduit flow system in predominantly siliciclastic strata

The Catskill Mountains of New York state are composed of thick sequences of interbedded mudstone, sandstone, and siltstone that are intercalated with thin beds of paludal limestone and calcareous mudstone. Sandstones often contain rip-up clasts of calcareous material at their bases and contain minor amounts of carbonate cement throughout their mass. This minor amount of calcareous material weathers out to produce thin zones of solutional porosity and permeability integrated within perched aquifers that feed perennial hillslope springs. Solutional features are evident within thin carbonate lenses and along fractures in sandstones, occasionally resulting in explorable cave-sized conduits. This study focuses on the behavior of bedrock springs and seeps in the headwaters of the Neversink River basin above the Neversink Reservoir. Bedrock springs are numerous throughout the Catskill Mountains and serve to provide baseflow to the headwater streams that tends to be more persistent during dry times than other baseflow sources. We observed dissolution-influenced conduits within the bedrock, and suffosional conduits within unconsolidated materials. Continuous water temperature at six springs was monitored at hourly intervals from May to October 2022 to investigate the thermal response of groundwater to daily and seasonal changes in air temperature along with event-based meteorological events. Temperature monitoring was augmented at some springs with time-lapse photography to evaluate persistence of flow. The thermal responses of springs to precipitation events were highly variable across sites, with some springs exhibiting flashy responses indicating a strong connection to the surface, while other springs showed muted or no thermal response to events. The seasonal thermal response also varied across sites having lag times of 23-48 days between peak air temperature and peak spring water temperature. The water temperatures and geochemical compositions of these bedrock springs support a conceptual model of multiple residence times, showing both event and seasonal groundwater storage within the conduit flow system that feeds the springs of the Catskill Mountains.

New York↗

Rapid deglaciation of the Tuolumne Meadows region linked to glacier hypsometry

During the LGM deglaciation (19–11 ka), episodes of glacier retreat and stabilization in the western United States are often attributed to abrupt climate change events via changes in Atlantic Meridional Overturning Circulation (e.g., Becker et al. 2025; McGee et al. 2018; Phillips 2016). However, there is considerable variability in the magnitude and timing of glacier retreat across the western US that is difficult to reconcile with a uniform climate forcing (Marcott et al. 2019), indicating either heterogeneous climate change or the influence of local non-climatic variables, such as glacier hypsometry (Young et al. 2011). Prior work from basal lake radiocarbon ages (Gillespie and Clark 2011) and cosmogenic nuclide dating (Phillips 2016; Phillips et al. 2009; Rood et al. 2011) established that Sierra Nevada glaciers receded from LGM moraines to north-facing cirques between 18 and 15 ka, with most deglaciation occurring at 15.8 ± 0.5 ka. As discussed above, geomorphic evidence supports abrupt glacier stagnation and in situ melting in the Tuolumne Meadows region.

Calfiornia↗