On the modeling of the 2010 Gulf of Mexico Oil Spill

Abstract Two oil particle trajectory forecasting systems were developed and applied to the 2010 Deepwater Horizon Oil Spill in the Gulf of Mexico. Both systems use ocean current fields from high-resolution numerical ocean circulation model simulations, Lagrangian stochastic models to represent unresolved sub-grid scale variability to advect oil particles, and Monte Carlo-based schemes for representing uncertain biochemical and physical processes. The first system assumes two-dimensional particle motion at the ocean surface, the oil is in one state, and the particle removal is modeled as a Monte Carlo process parameterized by a one number removal rate. Oil particles are seeded using both initial conditions based on observations and particles released at the location of the Maconda well. The initial conditions (ICs) of oil particle location for the two-dimensional surface oil trajectory forecasts are based on a fusing of all available information including satellite-based analyses. The resulting oil map is digitized into a shape file within which a polygon filling software generates longitude and latitude with variable particle density depending on the amount of oil present in the observations for the IC. The more complex system assumes three (light, medium, heavy) states for the oil, each state has a different removal rate in the Monte Carlo process, three-dimensional particle motion, and a particle size-dependent oil mixing model. Simulations from the two-dimensional forecast system produced results that qualitatively agreed with the uncertain “truth” fields. These simulations validated the use of our Monte Carlo scheme for representing oil removal by evaporation and other weathering processes. Eulerian velocity fields for predicting particle motion from data-assimilative models produced better particle trajectory distributions than a free running model with no data assimilation. Monte Carlo simulations of the three-dimensional oil particle trajectory, whose ensembles were generated by perturbing the size of the oil particles and the fraction in a given size range that are released at depth, the two largest unknowns in this problem. 36 realizations of the model were run with only subsurface oil releases. An average of these results yields that after three months, about 25% of the oil remains in the water column and that most of the oil is below 800 m.

[1]  Chuanmin Hu,et al.  Tracking the Deepwater Horizon Oil Spill: A Modeling Perspective , 2011 .

[2]  Frank E. Hoge,et al.  Interaction of Hurricane Katrina With Optically Complex Water in the Gulf of Mexico: Interpretation Using Satellite-Derived Inherent Optical Properties and Chlorophyll Concentration , 2009, IEEE Geoscience and Remote Sensing Letters.

[3]  Arthur J. Mariano,et al.  A Particle Filter for Inverse Lagrangian Prediction Problems , 2010 .

[4]  Øistein Johansen,et al.  Statmap-a 3-dimensional model for oil spill risk assessment , 2004, Environ. Model. Softw..

[5]  P. D’haeseleer,et al.  Deep-Sea Oil Plume Enriches Indigenous Oil-Degrading Bacteria , 2010, Science.

[6]  T. Prasad,et al.  Upper‐ocean response to Hurricane Ivan in a 1/25° nested Gulf of Mexico HYCOM , 2007 .

[7]  James Cummings,et al.  A comparison of sequential assimilation schemes for ocean prediction with the HYbrid Coordinate Ocean Model (HYCOM): Twin experiments with static forecast error covariances , 2011 .

[8]  G. Halliwell,et al.  Evaluation of vertical coordinate and vertical mixing algorithms in the HYbrid-Coordinate Ocean Model (HYCOM) , 2004 .

[9]  Eric P. Chassignet,et al.  North Atlantic Simulations with the Hybrid Coordinate Ocean Model (HYCOM): Impact of the Vertical Coordinate Choice, Reference Pressure, and Thermobaricity , 2003 .

[10]  R. Bleck,et al.  Global Ocean Prediction with the Hybrid Coordinate Ocean Model , 2009 .

[11]  P. Niiler,et al.  Circulation over the continental shelf in the northern Gulf of Mexico , 2005 .

[12]  F. Muller‐Karger,et al.  Use of thermal infrared remote sensing data for fisheries, environmental monitoring, oil and gas exploration, and ship routing. , 2006 .

[13]  T. M. Chin,et al.  Lagrangian Analysis and Prediction of Coastal and Ocean Dynamics: Lagrangian data assimilation in ocean general circulation models , 2007 .

[14]  Poojitha D. Yapa,et al.  BEHAVIOR OF OIL AND GAS FROM DEEPWATER BLOWOUTS , 2004 .

[15]  A. Barth,et al.  Impact of GODAE products on nested HYCOM simulations of the West Florida Shelf , 2009 .

[16]  Annalisa Griffa,et al.  Applications of stochastic particle models to oceanographic problems , 1996 .

[17]  Raymond F. Kokaly,et al.  Estimated Minimum Discharge Rates of the Deepwater Horizon Spill-Interim Report to the Flow Rate Technical Group from the Mass Balance Team , 2010 .

[18]  W. Large,et al.  Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization , 1994 .

[19]  A. Budhiraja,et al.  Modified particle filter methods for assimilating Lagrangian data into a point-vortex model , 2008 .

[20]  D. Yoerger,et al.  Tracking Hydrocarbon Plume Transport and Biodegradation at Deepwater Horizon , 2010, Science.

[21]  J. Cummings,et al.  Operational multivariate ocean data assimilation , 2005 .

[22]  R. Weisberg,et al.  Evaluation of Global Ocean Data Assimilation Experiment products on South Florida nested simulations with the Hybrid Coordinate Ocean Model , 2009 .

[23]  Annalisa Griffa,et al.  Lagrangian Analysis and Predictability of Coastal and Ocean Dynamics 2000 , 2002 .

[24]  Omar M. Knio,et al.  Many Task Computing for modeling the fate of oil discharged from the Deep Water Horizon well blowout , 2010, 2010 3rd Workshop on Many-Task Computing on Grids and Supercomputers.

[25]  Remy Baraille,et al.  The HYCOM (HYbrid Coordinate Ocean Model) data assimilative system , 2007 .

[26]  P. Malanotte‐Rizzoli,et al.  Interhemispheric water exchange in the Atlantic Ocean , 2003 .

[27]  A. Solberg,et al.  Oil spill detection by satellite remote sensing , 2005 .

[28]  A. Mariano,et al.  Lagrangian Data in a High Resolution Numerical Simulation of the North Atlantic. I: Comparison with , 2001 .

[29]  Eric P. Chassignet,et al.  US GODAE: Global Ocean Prediction with the Hybrid Coordinate Ocean Model (HYCOM) , 2004 .