Revista Cubana de Meteorología Vol. 30, No. sp., 2024, ISSN: 2664-0880
Código QR
CU-ID: https://cu-id.com/2377/v30nspe07
Original Article

Assimilation of satellite observations in support of hydrocarbon drift forecasting in the Cuban EEZ

Asimilación de observaciones satelitales en apoyo al pronóstico de deriva de hidrocarburos en la ZEE de Cuba

iDAlejandro Rodríguez Pupo1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.*✉:alejandro.rodriguez@insmet.cu

iDAlexander Lobaina Lao2Center of Exact Sciences and Engineering (CUCEI). University of Guadalajara. Mexico.

iDDailín Reyes Perdomo1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

iDPaola Nicolosi3e-GEOS SpA, Italia.

iDMelany Abreu del Sol1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

iDCosimo Garbellano3e-GEOS SpA, Italia.

iDFilippo Daffinà3e-GEOS SpA, Italia.

iDOsvaldo E. Pérez López1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

iDFrank E. Ortega Pereira1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

iDMelissa Abreu del Sol1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

iDReinaldo Casals Taylor1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

iDAmilcar E. Calzada Estrada4Laboratory of Physical, Coastal and Estuarine Oceanography (LOFEC), Department of Oceanography (DOCEAN), Federal University of Pernambuco, Recife, PE, Brazil.


1Centre for Marine Meteorology. Institute of Meteorology of Cuba. Cuba.

2Center of Exact Sciences and Engineering (CUCEI). University of Guadalajara. Mexico.

3e-GEOS SpA, Italia.

4Laboratory of Physical, Coastal and Estuarine Oceanography (LOFEC), Department of Oceanography (DOCEAN), Federal University of Pernambuco, Recife, PE, Brazil.

 

*Corresponding author: Alejandro Rodríguez Pupo. E-mail: alejandro.rodriguez@insmet.cu

ABSTRACT

The initialization of the runs of the numerical model of lagrangian approach PETROMAR-3D, from SAR images, is shown. With this new feature, the model is prepared to give operational response to an oil pollution event in the seas around Cuba. The oil slick polygon outputs of the SEonSE Engine and SNAP-ESA satellite information processing tools were used, developing new functions and Python libraries in the PETROMAR-3D software. With this result, a semi-automatic oil slick surveillance and monitoring service was developed in the seas around Cuba, with the computer tools of satellite data processing and the PETROMAR-3D Lagrangian oil spill model. With this service, users receive the information through a first email from the moment the slick is detected, and then receive another one with the oil drift modeling map for 72 hours. For this purpose, a script was developed that runs automatically through a scheduled task in the operating system, and handles all the software applications of the system designed with new Python functions and libraries.

Key words: 
PETROMAR, SAR imagery, SEonSE, marine pollution
RESUMEN

Se muestra la inicialización de las corridas del modelo numérico de enfoque lagrangiano PETROMAR-3D, a partir de imágenes SAR. Con esta nueva característica, el modelo está preparado para dar respuesta operativa ante un evento de contaminación por petróleo en los mares alrededor de Cuba. Fueron empleadas las salidas de polígonos de manchas de hidrocarburos de las herramientas informáticas de procesamiento de información satelital SEonSE Engine y SNAP-ESA, desarrollando nuevas funciones y bibliotecas de Python en el software PETROMAR-3D. Con este resultado, se desarrolló un servicio semi-automático de vigilancia y monitoreo de manchas de petróleo en los mares alrededor de Cuba, con las herramientas informáticas del procesamiento de datos satelitales y el modelo lagrangiano de derrames de petróleo PETROMAR-3D. Con este servicio, los usuarios reciben la información mediante un primer email desde el momento en que se detecta la mancha, y luego recibe otro con el mapa de la modelación de la deriva de petróleo para las 72 horas. Para tal propósito, fue desarrollado un script que corre automáticamente mediante tarea programada en el sistema operativo, y maneja todas las aplicaciones informáticas del sistema diseñado con nuevas funciones y bibliotecas de Python.

Palabras clave: 
PETROMAR, imágenes SAR, SEonSE, contaminación marina

Received: 02/10/2023; Accepted: 12/12/2023

Alejandro Rodríguez Pupo. Centre for Marine Meteorology. Institute of Meteorology of Cuba.

Alexander Lobaina LaO. Center of Exact Sciences and Engineering (CUCEI). University of Guadalajara. Mexico. E-mail: alobainalao@gmail.com

Dailín Reyes Perdomo. Centre for Marine Meteorology. Institute of Meteorology of Cuba. E-mail: dailin.reyes@insmet.cu

Paola Nicolosi. e-GEOS SpA, Italia. E-mail: paola.nicolosi@e-geos.it

Melany Abreu del Sol. Centre for Marine Meteorology. Institute of Meteorology of Cuba. E-mail: melanyabreu1699@gmail.com

Cosimo Garbellano. e-GEOS SpA, Italia. E-mail: cosimo.Garbellano@e-geos.it

Filippo Daffinà. e-GEOS SpA, Italia. E-mail: filippo.daffina@e-geos.it

Osvaldo E. Pérez López. Centre for Marine Meteorology. Institute of Meteorology of Cuba. E-mail: osvaldo.perez@insmet.cu

Frank E. Ortega Pereira. Centre for Marine Meteorology. Institute of Meteorology of Cuba. E-mail: frank.ortega@insmet.cu

Melissa Abreu del Sol. Centre for Marine Meteorology. Institute of Meteorology of Cuba. E-mail: melissa.abreu98@icloud.com

Reinaldo Casals Taylor. Centre for Marine Meteorology. Institute of Meteorology of Cuba. E-mail: reinaldo.casals@insmet.cu

Amilcar E. Calzada Estrada. Laboratory of Physical, Coastal and Estuarine Oceanography (LOFEC), Department of Oceanography (DOCEAN), Federal University of Pernambuco, Recife, PE, Brazil. E-mail: amilcar.calzada@ufpe.br.

Conflict of interest: We declare that we have no conflict of interest.

Author roles: Alejandro Rodríguez Pupo: Investigation, Software, Writing - Original Draft Preparation. Alexander Lobaina LaO: Investigation, Software. Dailín Reyes Perdomo: Data Curation, Formal Analysis, Investigation. Paola Nicolosi: Methodology, Investigation, Administration, supervision, Review & Editing. Melany Abreu del Sol: Data Curation, Formal Analysis, Investigation. Cosimo Garbellano: Data Curation, Software, Review & Editing. Filippo Daffinà: Conceptualization, Formal Analysis, Methodology, Software, Review & Editing. Osvaldo E. Pérez López: Data Curation, Investigation. Frank E. Ortega Pereira: Data Curation, Investigation. Melissa Abreu del Sol: Data Curation, Investigation. Reinaldo Casals Taylor: Data Curation, Formal Analysis, Investigation. Amilcar E. Calzada Estrada: Conceptualization, Formal Analysis, supervisión.

CONTENT

INTRODUCTION

 

It is estimated that 53% of the oil that reaches the marine environment is produced by man as the maritime transport industry is one of the main contributors to marine pollution. Spills from supertanker accidents represent only 8% of the total oil found in the oceans; However, both large spills (Perdomo et al., 2021Perdomo, D. R., Estrada, A. E. C., Hidalgo, D. C., Pupo, A. R., Domínguez, D. C., Secades, J. H., Infante, J. C., & LaO, A. L. (2021). Main terms and chronology of marine oil spill in the last 20 years. Revista Cubana de Meteorología, 27(4), e08. https://cu-id.com/2377/v27n4e08; Pupo et al., 2021Pupo, A. R., Estrada, A. E. C., Domínguez, D. C., Hidalgo, D. C., Perdomo, D. R., Veguería, P. J., LaO, A. L., Infante, J. C., & Secades, J. H. (2021). El modelo PETROMAR-3D bajo el paradigma del sistema de innovación. Revista Cubana de Meteorología, 27(3). http://rcm.insmet.cu/index.php/rcm/article/view/569), and small spills are the focus of attention (Polinov et al., 2021Polinov, S., Bookman, R., & Levin, N. (2021). Spatial and temporal assessment of oil spills in the Mediterranean Sea. Marine Pollution Bulletin, 167, 112338. https://doi.org/10.1016/j.marpolbul.2021.112338).

The transportation of hydrocarbons, via tankers, is a complex task that has given rise to oil spills as a result of collisions, groundings, fires and other accident events ranging from minor incidents to catastrophic spills. (Kark et al., 2015Kark, S., Brokovich, E., Mazor, T., & Levin, N. (2015). Emerging conservation challenges and prospects in an era of offshore hydrocarbon exploration and exploitation. Conservation Biology, 29(6), 1573–1585. https://doi.org/10.1111/cobi.12562; Perdomo et al., 2021Perdomo, D. R., Estrada, A. E. C., Hidalgo, D. C., Pupo, A. R., Domínguez, D. C., Secades, J. H., Infante, J. C., & LaO, A. L. (2021). Main terms and chronology of marine oil spill in the last 20 years. Revista Cubana de Meteorología, 27(4), e08. https://cu-id.com/2377/v27n4e08). Despite this, large oil spills have decreased considerably over recent decades (ITOPF, 2023ITOPF. (2023). ITOPF Handbook. ITOPF. https://www.itopf.org/fileadmin/uploads/itopf/data/Documents/Company_Lit/ITOPF_Handbook_2023). This has been achieved through the combined efforts of the oil-shipping industry and governments, as well as the actions of the International Maritime Organization (IMO) (Polinov et al., 2021Polinov, S., Bookman, R., & Levin, N. (2021). Spatial and temporal assessment of oil spills in the Mediterranean Sea. Marine Pollution Bulletin, 167, 112338. https://doi.org/10.1016/j.marpolbul.2021.112338).

However, another important source of oil pollution from ships is the deliberate discharge of hydrocarbon waste and ballast into the sea. These indisciplines constitute a constant threat to seabirds, benthic organisms and human health; and can cause a considerable number of victims over time, cited by (Polinov et al., 2021Polinov, S., Bookman, R., & Levin, N. (2021). Spatial and temporal assessment of oil spills in the Mediterranean Sea. Marine Pollution Bulletin, 167, 112338. https://doi.org/10.1016/j.marpolbul.2021.112338). These bad practices can be eliminated by strict implementation of existing regulations, control, monitoring and surveillance of maritime traffic.

PETROMAR-3D model (Calzada et al., 2021Calzada, A. E., Delgado, I., Ramos, C., Pérez, F., Reyes, D., Carracedo, D., Rodríguez, A., Chang, D., Cabrales, J., & Lobaina, A. (2021). Lagrangian Model PETROMAR-3D to Describe Complex Processes in Marine Oil Spills. Open Journal of Marine Science, 11(01), 17–40. https://doi.org/10.4236/ojms.2021.111002) has been developed at the Center for Marine Meteorology (CMM) of the Institute of Meteorology of Cuba (INSMET). This model is the result of the study of the characteristics of numerous oil spill models that operate in different regions of the world, being a technical solution according to the organizational conditions and resources available in the country. (Pupo et al., 2021Pupo, A. R., Estrada, A. E. C., Domínguez, D. C., Hidalgo, D. C., Perdomo, D. R., Veguería, P. J., LaO, A. L., Infante, J. C., & Secades, J. H. (2021). El modelo PETROMAR-3D bajo el paradigma del sistema de innovación. Revista Cubana de Meteorología, 27(3). http://rcm.insmet.cu/index.php/rcm/article/view/569). In addition to covering the simulations in the seas adjacent to Cuba, a work domain has recently been designed to run the model in the northeast of Brazil (Varona et al., 2024Varona, H. L., Noriega, C., Calzada, A. E., Medeiros, C., Lobaina, A., Rodriguez, A., Chang, D., Reyes, D., Araujo, J., Silva, M., Das Chagas Moura, M., & Araujo, M. (2024). Effects of meteo-oceanographic conditions on the weathering processes of oil spills in northeastern Brazil. Marine Pollution Bulletin, 198, 115828. https://doi.org/10.1016/j.marpolbul.2023.115828).

With PETROMAR-3D version 2.1 (Calzada et al., 2021Calzada, A. E., Delgado, I., Ramos, C., Pérez, F., Reyes, D., Carracedo, D., Rodríguez, A., Chang, D., Cabrales, J., & Lobaina, A. (2021). Lagrangian Model PETROMAR-3D to Describe Complex Processes in Marine Oil Spills. Open Journal of Marine Science, 11(01), 17–40. https://doi.org/10.4236/ojms.2021.111002), an oil slick fate forecasting service via email was developed operationally. Users had to send an email with a text file containing the model input file; but, if this file was not sent, the application had a default configuration.

Despite the improvement in operational performance in PETROMAR-3D version 2.1, the oil drift and weathering forecast lacked efficiency because it needed information on initial conditions only from the configuration file; and it did not assimilate a polygon with the morphology of the oil slick, but instead calculated the spread of the slick by applying the Fay equation (Fay, 1971)Fay, J. A. (1971). Physical processes in the spread of oil on a water surface. International Oil Spill Conference Proceedings, 1971(1), 463–467. https://doi.org/10.7901/2169-3358-1971-1-463 assuming a circle. Consequently, these specialized services to the entities in charge of maritime safety and the mitigation of pollution at sea did not constitute an efficient tool contributing to the Cuban Marine Pollution Early Warning System.

Starting in 2022, with the international collaboration project e-GEOS-INSMET: “Strengthening the Cuban marine system”, with alias “Marine Surveillance”; The acquisition of SAR images was achieved continuously and with full coverage of the territorial seas of the Cuban archipelago every 48 hours for 9 months. This process, completely unprecedented in Cuba, made it possible to download and digitally process SAR images from the COSMO-SkyMED satellite constellation using the SAR digital image processing computer tool SEonSE Engine (Oddone et al., 2021Oddone, A., Daffina, F., Quattrociocchi, D., Francioni, E., Gentile, V., Pieralice, F., Cesarano, L., Britti, F., Pietranera, L., Minati, F., Luzietti, L., Grandoni, D., Cardillo, P. F., Biscardi, M. A., & Bonetto, F. (2021). Satellite Sar Data Exploitation: Innovations and Potentialities for Earth Observation Advanced Monitoring and Servicing. In IGARSS 2021—2021 IEEE International Geoscience and Remote Sensing Symposium (pp. 1887–1890). IEEE. https://doi.org/10.1109/IGARSS47720.2021.9554026), provided by company e-GEOS SpA, in accordance with the project contract.

In this way, the conditions were created to develop the initialization of the PETROMAR - 3D lagrangian model from satellite information and introduce surveillance, monitoring and forecasting services for drift and weathering of oil slicks with the incorporated technology. So this article aims to show the development of an improved version of the PETROMAR-3D model, which is based on the information resulting from satellite image processing.

After the Introduction, section 2 describes the computer tools used for the digital processing of satellite images, the PETROMAR-3D model and the data used in the operational work of oil slick monitoring, and some considerations that were taken into account to carry out the modeling. Section 3 describes the results of the PETROMAR-3D assimilation of information from SAR image processing and, subsequently, the development of a semi-automatic surveillance and monitoring system for oil slicks at sea. Finally, conclusions and recommendations are presented in Section 4.

2. METHODS AND MATERIALS

 

Oil is depicted on SAR images by long stripes on the sea surface. The low water dissolution capacity of oil results in high shortwave damping and expansion in surface tension, decreasing wind erosion, appearing in SAR images as a dark region (Gee et al., 2016Gee, D., Sowter, A., Novellino, A., Marsh, S., & Gluyas, J. (2016). Monitoring land motion due to natural gas extraction: Validation of the Intermittent SBAS (ISBAS) DInSAR algorithm over gas fields of North Holland, the Netherlands. Marine and Petroleum Geology, 77, 1338–1354. https://doi.org/10.1016/j.marpetgeo.2016.08.014) that is detected in wind speed ranges between 3 and 12 m/s (Ramsey III et al., 2013Ramsey, E., Rangoonwala, A., & Bannister, T. (2013). Coastal Flood Inundation Monitoring with Satellite C‐band and L‐band Synthetic Aperture Radar Data. JAWRA Journal of the American Water Resources Association, 49(6), 1239–1260. https://doi.org/10.1111/jawr.12082).

A hydrocarbon spill into the sea depends on components such as its volume, its physical properties, and meteo-ocean conditions. The different physical-chemical-oceanographic processes in oil affect the sensitivity of SAR images over time (Chaturvedi et al., 2020Chaturvedi, S. K., Banerjee, S., & Lele, S. (2020). An assessment of oil spill detection using Sentinel 1 SAR-C images. Journal of Ocean Engineering and Science, 5(2), 116–135. https://doi.org/10.1016/j.joes.2019.09.004). The compilations of repeated SAR images, in the same event, are very useful for monitoring the evolution of an oil spill (Liu et al., 2013Liu, Y., MacFadyen, A., Ji, Z.-G., & Weisberg, R. H. (2013). Monitoring and modeling the Deepwater Horizon oil spill: A record breaking enterprise. John Wiley & Sons. https://doi.org/10.1002/2013EO200006).

2.1 SAR image processing tools for detection of oil slicks in the sea

 

With the e-GEOS-INSMET collaboration, two very important tools for SAR image processing were installed at the CMM: SEonSE and SNAP-ESA. SAR images of the COSMO-SkyMed constellation of the ScanSAR Huge type were acquired by Matera Remote Sensing Centre of e-GEOS in Italy and then downloaded by INSMET, comprising an area per image of 200 x 200 km with 100 m spatial resolution, covering 350,000.00 km2 of the Cuban archipelago in 48 hours, reaching 80% coverage.

SEonSE Engine

 

This is a client-server architecture application, implemented as a QGIS plugin, to detect and characterize ships and oil spills from SAR images acquired by satellite, helping the user to perform the processing. This software can work with COSMO-SkyMed, SAR RADARSAT2, and Sentinel-1 images without requiring pre-processing because the products are read in native format.

Being an application that focuses on maritime surveillance, it automatically loads data from AIS (Automatic Identification System) satellite tracking systems with the tracks and details of the ships present in the analysis area, in addition to wind measurements and waves. Possible polluters of the detected oil spills can then be derived with SAR image processing, among the closest vessels (Oddone et al., 2021Oddone, A., Daffina, F., Quattrociocchi, D., Francioni, E., Gentile, V., Pieralice, F., Cesarano, L., Britti, F., Pietranera, L., Minati, F., Luzietti, L., Grandoni, D., Cardillo, P. F., Biscardi, M. A., & Bonetto, F. (2021). Satellite Sar Data Exploitation: Innovations and Potentialities for Earth Observation Advanced Monitoring and Servicing. In IGARSS 2021—2021 IEEE International Geoscience and Remote Sensing Symposium (pp. 1887–1890). IEEE. https://doi.org/10.1109/IGARSS47720.2021.9554026).

SNAP-ESA

 

SNAP (Sentinel Application Platform) responds to a free program offered by the European Space Agency to process and analyze satellite images from the Sentinel fleet of satellites. The application is developed jointly by Brockmann Consult GmbH and has specific tools to work with the images based on the satellite model (Sentinel Booxes); be multiband images (Sentinel 2, Sentinel 3), multispectral images (Envisat, Landsat, MODIS), and SAR images.

The strength of the tool is the strategic combination of bands that allows working on various aspects of territorial issues. To do this, RGB combinations can be made with the satellite bands that are being used to compose images, both in real color and false color, for the treatment of territorial aspects based on agriculture, land use, urban planning, bodies of water, forest management among others.

As part of the SNAP tools are intended for the management of standardized Sentinel products through the folder and metadata structure of the Copernicus platform, the work band structure available for treatment and management can be used in a single step. of Sentinel satellite images, uploading an XML file for bulk incorporation of all image data in a single step (Foumelis et al., 2018Foumelis, M., Delgado Blasco, J. M., Desnos, Y.-L., Engdahl, M., Fernandez, D., Veci, L., Lu, J., & Wong, C. (2018). Esa Snap—Stamps Integrated Processing for Sentinel-1 Persistent Scatterer Interferometry. In IGARSS 2018—2018 IEEE International Geoscience and Remote Sensing Symposium (pp. 1364–1367). IEEE. https://doi.org/10.1109/IGARSS.2018.8519545).

2.2 The PETROMAR-3D oil spill model

 

The PETROMAR - 3D model (Calzada et al., 2021Calzada, A. E., Delgado, I., Ramos, C., Pérez, F., Reyes, D., Carracedo, D., Rodríguez, A., Chang, D., Cabrales, J., & Lobaina, A. (2021). Lagrangian Model PETROMAR-3D to Describe Complex Processes in Marine Oil Spills. Open Journal of Marine Science, 11(01), 17–40. https://doi.org/10.4236/ojms.2021.111002) is the modeling tool used in this work. It is a Lagrangian scheme to calculate oil spill trajectories, created at the Center for Marine Meteorology of the Institute of Meteorology. This model incorporates numerical forecast results from various atmospheric, hydrodynamic and wave models, as well as other input data for drift and weathering calculations. The data used in the modeling came from forecast outputs of atmospheric and oceanographic models detailed below:

  • GFS (Global Forecast System), a global model for numerical weather prediction owned by the National Oceanic and Atmospheric Administration (NOAA). The time domain of the model is presented in two parts: the first has a higher resolution and extends up to 180 hours (7 days), while the second part runs from 180 to 384 hours (16 days) at a lower resolution. . It is used as the spatial domain of the model in a horizontal direction that divides the surface of the earth into cells of 0.5 ° x 0.5 °, although there are outputs of lower resolution (1.0 ° x 1.0 ° and 2.5 ° x 2.5 °). Wind direction and speed variables are used (Ncdc, 2018Ncdc. (2018). Global Forecast System (GFS) | National Centers for Environmental Information (NCEI) Formerly Known as National Climatic Data Center (NCDC). https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs) with spatial and temporal resolution of 28 km and 3 hours respectively.

  • NAM-CARIBE (North American Mesoscale-Caribbean Forecast System). The domain considers the entire Greater Caribbean and presents a spatial resolution of 12 km and a time of 3 hours (Nam, 2017Nam. (2017). North America Mesoscale Forecasting System | National Center for Environmental Information (NCEI) formerly known as National Climatic Data Center (NCDC). https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/north-american-mesoscale-forecast-system-nam). The download site is: http://nomads.ncep.noaa.gov/pub/data.

  • cmems_mod_glo_phy_anfc_0.083deg_PT1H-m, belonging to the Copernicus product GLOBAL_ANALYSISFORECAST_PHY_001_024. The Operational Mercator global ocean analysis and forecast system at 1/12 degree is providing 10 days of 3D global ocean forecasts updated daily. This product includes hourly mean surface fields for sea level height, temperature and currents. The global ocean output files are displayed with a 1/12 degree horizontal resolution with regular longitude/latitude equirectangular projection. 50 vertical levels are ranging from 0 to 5500 meters (Lellouche et al., 2023Lellouche, J. M., Galloudec, O. L., Regnier, C., Gennip, S. V., Chune, S. L., Levier, B., Greiner, E., Drevillon, M., & Szczypta, C. (2023). Global Production Centre (Quality Information Document GLO_ANALYSISFORECAST_PHY_001_024; p. 75). Copernicus Marine Service. https://catalogue.marine.copernicus.eu/documents/QUID/CMEMS-GLO-QUID-001-024.pdf). The download is available at: https://data.marine.copernicus.eu/product/GLOBAL_ANALYSISFORECAST_PHY_001_024/description

  • NCOM (Navy Coastal Ocean Model) US Navy operational global ocean model (Martin et al., 2009Martin, P. J., Barron, C. N., Smedstad, L. F., Campbell, Ti. J., Wallcraft, A. J., Rhodes, R. C., Rowley, C., Townsend, T. L., & Carroll, S. N. (2009). Navy Coastal Ocean Model (NCOM) Version 4.0 (User’s Manual): Defense Technical Information Center. https://doi.org/10.21236/ADA494652; Ncom, 2017Ncom. (2017). Naval Oceanographic Office Regional Navy Coastal Ocean Model (NCOM) | National Centers for Environmental Information (NCEI) formerly known as National Climatic Data Center (NCDC). https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/navoceano-ncom-reg), which includes several domains including the Gulf of Mexico and the Caribbean Sea. The download website is: http://nomads.ncep.noaa.gov/pub/data

2.3 Considerations in modeling with PETROMAR-3D from SAR images

 

Operational monitoring and forecasting of oil slicks at sea is very complex. In most cases, when a marine oil spill event occurs, the properties of the hydrocarbons are unknown. The volume dumped, its density, the date and time of occurrence are examples of these properties; which are necessary for numerical modeling and forecasting of the fate and modification of the properties of the slick over time.

Also, it may be the case (and it happens most of the time) that the acquisition of the satellite image in the area of the oil spill occurs several hours after the event began, which makes the task even more complex for modelers. Taking into account these factors, some considerations have been taken into account to make the operational forecast of destination and changes in properties of the oil slicks detected through satellite images:

  • 1. Do not include the spreading calculation in the modeling in operational mode.

The PETROMAR-3D model performs the spreading calculation using the Fay equations (Fay, 1971)Fay, J. A. (1971). Physical processes in the spread of oil on a water surface. International Oil Spill Conference Proceedings, 1971(1), 463–467. https://doi.org/10.7901/2169-3358-1971-1-463. This physical process occurs in the first hours of the spill, and its duration depends on the volume spilled (Perdomo et al., 2021Perdomo, D. R., Estrada, A. E. C., Hidalgo, D. C., Pupo, A. R., Domínguez, D. C., Secades, J. H., Infante, J. C., & LaO, A. L. (2021). Main terms and chronology of marine oil spill in the last 20 years. Revista Cubana de Meteorología, 27(4), e08. https://cu-id.com/2377/v27n4e08). These data are unknown if the modeling is initialized from the processing of a satellite image.

  • 2. Carry out numerical modeling with an average density of the spilled oil.

Table 1 shows the classification of hydrocarbons according to density (Mohammadi et al., 2020Mohammadi, M., Khanmohammadi Khorrami, M., Vatani, A., Ghasemzadeh, H., Vatanparast, H., Bahramian, A., & Fallah, A. (2020). Rapid determination and classification of crude oils by ATR-FTIR spectroscopy and chemometric methods. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 232, 118157. https://doi.org/10.1016/j.saa.2020.118157).

Table 1.  Classification of hydrocarbons according to density
Hydrocarbon classification API Degree Density [g/cm3]
Ligth > 31.1 <0.87
Medium 22.3-31.1 0.92-0.87
Heavy 10-22.3 1.0-0.92

To carry out the numerical modeling, since the initial value of the density of the spilled oil is unknown, a density value has been considered according to the case. For light oil (in API degrees): 33, 25 for medium and 20 for heavy.

  • 3. Estimate the volume of the slick.

The coordinates of the oil slick polygon can be provided to the PETROMAR-3D model in two ways:

  • a) With the SNAP-ESA polygon creation tool. The polygons are created manually from both SAR and optical images.

In this case, a division of the polygon into pixels is assumed (square particles of 0.001 degrees, equivalent to 12 m2), and the area of the slick is calculated by the sum of the total area of the calculated pixels.

  • b) With the generation of a Shapefile file from SEonSE Engine. The polygons are created from the spill report in shape format and the data on the spill area is included in this file.

Once the area is obtained by any of the methods, the estimated volume of the slick is calculated by multiplying the area by the thickness (assuming that the thickness of the slick is homogeneous). In the same way as the density, 3 average values of slick thickness were established according to the magnitude of the spill. In this way, for hydrocarbon spills: heavy 10-2 m, medium 10-3 m and light 10-4 m (Perdomo et al., 2021Perdomo, D. R., Estrada, A. E. C., Hidalgo, D. C., Pupo, A. R., Domínguez, D. C., Secades, J. H., Infante, J. C., & LaO, A. L. (2021). Main terms and chronology of marine oil spill in the last 20 years. Revista Cubana de Meteorología, 27(4), e08. https://cu-id.com/2377/v27n4e08).

These considerations do not constitute a factor that undermines the modeling precision because the forecast is updated with the acquisition of a subsequent satellite image; and furthermore, in the hours following the start of the event, data on the properties of the hydrocarbons can be obtained.

3. RESULTS

 

3.1 Assimilation of satellite information

 

With the SNAP-ESA and SEonSE Engine software, a polygon output is generated using a file with the coordinate values of the oil slick. This file is assimilated by the PETROMAR-3D model through a Python script that is responsible for locating the polygon in the model's working domain.

The functions developed in the new PETROMAR-3D script are executed with pre-established parameters such as the date, the API degree value and the polygon area. These new features are very advantageous for the operational work of the model because previously the software config file was configured manually.

There has already been the case in which several oil slicks appear in the same SAR image. With the new features of PETROMAR-3D, several slicks can be modeled in the same scenario, each of which are introduced from their respective files (figure 1).

Figure 1.  SAR image from October 5, 2022, detected with the SEonSE Engine app.

The coordinates of the polygon enter as parameters into the model, which creates the Lagrangian particles that must be processed in the calculations of the functions responsible for modeling.

This is achieved by choosing the smallest rectangular domain from the maximum and minimum longitudes and latitudes that include the spot; and a square grid with a given dxdy is constructed by applying the polygon mask to discriminate the interior nodes. In each interior node an element is planted and its drift is modeled (figure 2).

Figure 2.  The image shows the process from which the polygon of the slick is assimilated, the grid points that are inside the slick are obtained, and the output of the plotting with PETROMAR-3D.

3.2 Development of a monitoring system for oil slicks in the sea

 

Based on the outputs of the SEonSE Engine, SNAP-ESA applications, and the initialization of PETROMAR-3D with the generated polygons, a Surveillance and Monitoring System for oil slicks around Cuba was developed. This task was carried out by means of a new Python script running automatically on a server.

This script uses Python libraries for email messaging such as: poplib, smtplib, email.mime.multipart, email.mime.text, email.mime.base, email.mime.image. It has a defined function that is responsible for sending notifications to users in case of oil slick detections through emails. It was scheduled as an automatic task to be executed every 1 hour, although it depends on the pre-processing of the SAR images by the SEonSE Engine operator on duty. This feature makes the oil slick surveillance and monitoring system semi-automatic.

Figure 3 shows a diagram of the semi-automatic surveillance process of the developed system. In it you can see the hardware and software tools that make it up.

The system is made up of a server, with the PETROMAR-3D model installed and the Python script developed in this research (run_petromar.py); and 2 NAS servers (network-connected data devices) storing the data from the outputs of the numerical models downloaded from the different regional and global metocean centers, such as: data on ocean currents, wind and waves; and the processed SAR images, polygons of the oil slicks and corresponding png images, the PETROMAR-3D scenario and trajectory outputs, and the “run_petromar.py” configuration files respectively.

In a folder created for the configuration files of “run_petromar.py”, there are files that contain the mailing list of the system users, the texts of the notification messages and the oil slick fate forecast for the next 72 hours, and a log file that saves the information of each service provided to clients.

So by running the script automatically every 1 hour, it creates a list with the date and time of the last 24 hours, checking the processed images folder to see if they contain any polygon files within the time interval of the previously created list. . The configuration log file that stores the information on the spills that have been simulated is reviewed, and the data is filtered by selecting the new spills detected to execute the modeling.

For each case of interest, the script sends a notification email to the previously configured addresses, attaching the png image associated with the spill; and then models the behavior of the oil slick for the next 72 hours by starting the PETROMAR-3D model run. At the end of the run, a second email is sent to the clients with the graphic outputs. The process ends with writing the run data to the log file.

This new service constitutes a new approach in marine surveillance in Cuba. With the previous approach, the user had to provide information on the occurrence of the spill, sent the data on the initial conditions and did not know, in most cases, the cause of the damage. From now on, the monitoring system is the one that discovers the oil slicks, and sends the notification and forecasts to the end users. This constitutes a new paradigm in marine surveillance in our geographical area.

Figure 3.  Diagram showing the design of the semi-automatic oil spill surveillance and monitoring system.

The challenge that this research team faces from now on is the sustainability of this service; given that there is the possibility of obtaining SAR Sentinel 1 images of the Copernicus site daily for free, but with spatial coverage limitations. This challenge can be solved by incorporating into the developed semi-automatic system the oil slick detection methodologies used in satellites with optical sensors. With this technology, free images are also available and several criteria are applied to detect oil slicks in the sea. Only in this case, there will be no automatic information regarding the vessels sailing through the scene of the incident.

CONCLUSION

 

At the end of this investigation, the following results were achieved:

  1. New functions were programmed into the PETROMAR-3D model, creating the capabilities for the initialization of the runs from the SAR images processed with the SNAP-ESA and SEonSE Engine tools.

  2. A semi-automatic system for surveillance and monitoring of hydrocarbon slicks in the seas of the EEZ of Cuba was developed.

RECOMMENDATIONS

 

The authors of this research recommend continuing to delve deeper into:

  1. The analysis of the spread of the oil spill from obtaining the polygon from the satellite image.

  2. Estimation of the thickness of oil slicks drifting on the marine surface, using SAR technology.

  3. Work on the sustainability of this service by incorporating optical satellites and their respective digital oil slick processing techniques.

  4. Continue working on the promotion of this service, as well as environmental education to preserve our marine environment.

REFERENCES

 

Calzada, A. E., Delgado, I., Ramos, C., Pérez, F., Reyes, D., Carracedo, D., Rodríguez, A., Chang, D., Cabrales, J., & Lobaina, A. (2021). Lagrangian Model PETROMAR-3D to Describe Complex Processes in Marine Oil Spills. Open Journal of Marine Science, 11(01), 17–40. https://doi.org/10.4236/ojms.2021.111002

Chaturvedi, S. K., Banerjee, S., & Lele, S. (2020). An assessment of oil spill detection using Sentinel 1 SAR-C images. Journal of Ocean Engineering and Science, 5(2), 116–135. https://doi.org/10.1016/j.joes.2019.09.004

Fay, J. A. (1971). Physical processes in the spread of oil on a water surface. International Oil Spill Conference Proceedings, 1971(1), 463–467. https://doi.org/10.7901/2169-3358-1971-1-463

Foumelis, M., Delgado Blasco, J. M., Desnos, Y.-L., Engdahl, M., Fernandez, D., Veci, L., Lu, J., & Wong, C. (2018). Esa Snap—Stamps Integrated Processing for Sentinel-1 Persistent Scatterer Interferometry. In IGARSS 2018—2018 IEEE International Geoscience and Remote Sensing Symposium (pp. 1364–1367). IEEE. https://doi.org/10.1109/IGARSS.2018.8519545

Gee, D., Sowter, A., Novellino, A., Marsh, S., & Gluyas, J. (2016). Monitoring land motion due to natural gas extraction: Validation of the Intermittent SBAS (ISBAS) DInSAR algorithm over gas fields of North Holland, the Netherlands. Marine and Petroleum Geology, 77, 1338–1354. https://doi.org/10.1016/j.marpetgeo.2016.08.014

ITOPF. (2023). ITOPF Handbook. ITOPF. https://www.itopf.org/fileadmin/uploads/itopf/data/Documents/Company_Lit/ITOPF_Handbook_2023

Kark, S., Brokovich, E., Mazor, T., & Levin, N. (2015). Emerging conservation challenges and prospects in an era of offshore hydrocarbon exploration and exploitation. Conservation Biology, 29(6), 1573–1585. https://doi.org/10.1111/cobi.12562

Lellouche, J. M., Galloudec, O. L., Regnier, C., Gennip, S. V., Chune, S. L., Levier, B., Greiner, E., Drevillon, M., & Szczypta, C. (2023). Global Production Centre (Quality Information Document GLO_ANALYSISFORECAST_PHY_001_024; p. 75). Copernicus Marine Service. https://catalogue.marine.copernicus.eu/documents/QUID/CMEMS-GLO-QUID-001-024.pdf

Liu, Y., MacFadyen, A., Ji, Z.-G., & Weisberg, R. H. (2013). Monitoring and modeling the Deepwater Horizon oil spill: A record breaking enterprise. John Wiley & Sons. https://doi.org/10.1002/2013EO200006

Martin, P. J., Barron, C. N., Smedstad, L. F., Campbell, Ti. J., Wallcraft, A. J., Rhodes, R. C., Rowley, C., Townsend, T. L., & Carroll, S. N. (2009). Navy Coastal Ocean Model (NCOM) Version 4.0 (User’s Manual): Defense Technical Information Center. https://doi.org/10.21236/ADA494652

Mohammadi, M., Khanmohammadi Khorrami, M., Vatani, A., Ghasemzadeh, H., Vatanparast, H., Bahramian, A., & Fallah, A. (2020). Rapid determination and classification of crude oils by ATR-FTIR spectroscopy and chemometric methods. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 232, 118157. https://doi.org/10.1016/j.saa.2020.118157

Nam. (2017). North America Mesoscale Forecasting System | National Center for Environmental Information (NCEI) formerly known as National Climatic Data Center (NCDC). https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/north-american-mesoscale-forecast-system-nam

Ncdc. (2018). Global Forecast System (GFS) | National Centers for Environmental Information (NCEI) Formerly Known as National Climatic Data Center (NCDC). https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs

Ncom. (2017). Naval Oceanographic Office Regional Navy Coastal Ocean Model (NCOM) | National Centers for Environmental Information (NCEI) formerly known as National Climatic Data Center (NCDC). https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/navoceano-ncom-reg

Oddone, A., Daffina, F., Quattrociocchi, D., Francioni, E., Gentile, V., Pieralice, F., Cesarano, L., Britti, F., Pietranera, L., Minati, F., Luzietti, L., Grandoni, D., Cardillo, P. F., Biscardi, M. A., & Bonetto, F. (2021). Satellite Sar Data Exploitation: Innovations and Potentialities for Earth Observation Advanced Monitoring and Servicing. In IGARSS 2021—2021 IEEE International Geoscience and Remote Sensing Symposium (pp. 1887–1890). IEEE. https://doi.org/10.1109/IGARSS47720.2021.9554026

Perdomo, D. R., Estrada, A. E. C., Hidalgo, D. C., Pupo, A. R., Domínguez, D. C., Secades, J. H., Infante, J. C., & LaO, A. L. (2021). Main terms and chronology of marine oil spill in the last 20 years. Revista Cubana de Meteorología, 27(4), e08. https://cu-id.com/2377/v27n4e08

Polinov, S., Bookman, R., & Levin, N. (2021). Spatial and temporal assessment of oil spills in the Mediterranean Sea. Marine Pollution Bulletin, 167, 112338. https://doi.org/10.1016/j.marpolbul.2021.112338

Pupo, A. R., Estrada, A. E. C., Domínguez, D. C., Hidalgo, D. C., Perdomo, D. R., Veguería, P. J., LaO, A. L., Infante, J. C., & Secades, J. H. (2021). El modelo PETROMAR-3D bajo el paradigma del sistema de innovación. Revista Cubana de Meteorología, 27(3). http://rcm.insmet.cu/index.php/rcm/article/view/569

Ramsey, E., Rangoonwala, A., & Bannister, T. (2013). Coastal Flood Inundation Monitoring with Satellite C‐band and L‐band Synthetic Aperture Radar Data. JAWRA Journal of the American Water Resources Association, 49(6), 1239–1260. https://doi.org/10.1111/jawr.12082

Varona, H. L., Noriega, C., Calzada, A. E., Medeiros, C., Lobaina, A., Rodriguez, A., Chang, D., Reyes, D., Araujo, J., Silva, M., Das Chagas Moura, M., & Araujo, M. (2024). Effects of meteo-oceanographic conditions on the weathering processes of oil spills in northeastern Brazil. Marine Pollution Bulletin, 198, 115828. https://doi.org/10.1016/j.marpolbul.2023.115828