Acoustic Evidence of Shallow Gas Occurrences in the Offshore Sin ú Fold Belt, Colombian Caribbean Sea

: High-resolution seismic analysis and bathymetry data, used in the Offshore Sin ú Fold Belt (OSFB), have revealed seabed and sub-surface anomalies, which were probably caused by the presence of shallow gas within the sedimentary records. Shallow gas is widely detected by the frequent presence of anomalous acoustic reﬂections including acoustic blanking, enhanced reﬂections, acoustic plumes, pockmarks, and dome structures. More than 30 anomalies that occur within a subsurface depth of ~65 m were acoustically detected within an area of 1000 km 2 on the continental shelf and upper continental slope, in water depths ranging from − 20 to − 1300 m. Moreover, a map with the spatial distribution of the gas occurrences is shown. A close relationship between the locally elevated seabed (dome structures), pockmarks, and acoustic blanking was found. Most of the active pockmarks may be closely related to the submarine path of the Uramita Fault, indicating that the gas occurrences are controlled by active faulting. The shallow gas occurrence was conﬁrmed by the generation of authigenic carbonate and the occurrence of chemosymbiotic biological communities sampled in the area. Although there is an admixture of biogenic gas, it is believed that many of the features observed relate to thermogenic gas. The identiﬁcation of these anomalies represents a useful basis for an assessment of marine geohazards and can serve as a hydrocarbon exploration tool.


Introduction
The study of submarine, coastal, and insular areas in the world has been widely developed in recent decades [1][2][3][4][5].This is particularly thanks to the development of technologies and methods in the acoustic, geodesic, geophysical, and remote sensing fields, as well as the improvement in geospatial data management [6][7][8][9][10][11][12].The incorporation of new methods for the exploration of deep environments has allowed researchers to detect and map shallow gas accumulations in marine deposits [13][14][15].Following Fleischer et al. [15], shallow gas is defined here as free gas trapped in sub-bottom sediments within a depth of 1000 m; this is the same for Holocene muddy sediments [16], which are considered the most important pathway for the emission of methane from deeper marine sediments to the atmosphere [17], originating in specific structures on the seabed and inducing acoustic features in the sedimentary record [13,15,[18][19][20].Therefore, the detection and mapping of shallow gas is usually performed through the analysis of high-resolution acoustic tools and techniques [13,[21][22][23].
Shallow gas in near-surface sediments causes numerous anomalies in seismic data due to acoustic impedance between gas-bearing and gas-free sediments [13,21,[24][25][26][27].These anomalies include: (1) Acoustic turbidity caused by the scattering of acoustic energy due to the presence of gases in the sediment pore space, revealing a poor amplitude and low continuity of reflectors [13,[28][29][30].Acoustically turbid zones represent gas-charged sediments [31].According to Fanin [32], this effect is already generated with only 1% of gas present in sediments.(2) Acoustic blanking zones related to the absence of sediment layering due the complete absorption of seismic signals in gas-charged sediments [13,29,33].The most evident acoustic anomaly is blanking, appearing as diffuse and chaotic seismic facies, masking all other reflections [34].The difference between acoustic turbidity and blanking is explained by the amount of free fluid accumulated in the sediment [35].
(3) Enhanced reflections caused by reflectors that are discontinuous, with high reflectivity or amplitude due to gas content within the shallow sediments [36].(4) Bright spots resulting from the increase in the acoustic impedance contrast between sediments in which fluid is present and their surrounding areas.This anomaly indicates the top of a fluid-charged sediment with a high-amplitude reverse polarity [37].The other acoustic features that may indicate the occurrence of shallow gas in seismic data include gas seeps and pockmarks.The acoustic blankings are explained as vertical disturbances of the sediment layer due the acoustic dispersion effect of the upward migration of gases or fluids in pores [36,38].The vertical disturbances may seem to suggest the presence of small gas seeps [28].The pockmarks are defined as crater-shaped seabed depressions that are generally related to gas and/or fluid discharge in shallow sediments [39].Also, large-scale gas leaks can produce numerous bubbles throughout the water column that are detected, using highfrequency echograms, as "gas plumes" or "hydroacoustic flares" [40].In addition, fluid/gas movement indicators can cause seafloor irregularities such as diapirism, mud volcanism, subsidence and/or landslides, and other indicators [41,42].
The Colombian Caribbean region has been a target for energetic resource explorations since the early 1950s.Since the discovery of the Chuchupa and Ballenas fields in the 1970s [43,44], the region has been considered a gas province originating mainly from microbial activity [45].Also, the latest discoveries in areas located below 2000 m with the drilling of the Orca (2014), Kronos (2015), Purple Angel (2017), and Gorgon (2022) exploratory wells confirm the presence of an emerging basin dominated by methane [46].Despite these discoveries, the knowledge of shallow gas-bearing sediments is still limited and the shelf remains poorly understood.In this study, we analyzed multibeam bathymetry and seismic sub-bottom profilers (SBPs) in the Offshore Sinú Fold Belt (OSFB), in the Colombian Caribbean Sea, to describe and map the acoustic anomalies related to the occurrence of gas.The study was largely focused on the shallow gassy sediments.The knowledge obtained will be helpful analysis as a proxy for seep distribution and seep activity; it can serve as an indicator for deeper hydrocarbon reserves and thus as an exploration tool.

Geological Configuration
The study area covers an area of 6700 km 2 and extends along the shelf and continental slope with a maximum depth of 2870 m (Figure 1).It is located on a wide platform, generally with muddy bottoms, with turbid, warm, and low-salinity environments [47,48].Sediment supply is influenced by the deltaic systems of four rivers: Atrato, Mulatos, Sinú, and Magdalena (Figure 1B).The Atrato River presents a moderate annual sediment flow of 11.26 × 10 6 ton yr −1 [49], most of which is deposited within the gulf.The Mulatos River discharges around 0.21 × 10 6 ton yr −1 , and the Sinú River discharges 6.1 × 10 6 ton yr −1 .Finally, the Magdalena River discharges 143.9 × 10 6 ton yr −1 [49][50][51].[52]), onshore geology, and the main tectonic elements (based on [53]).It covers an area of 6700 km 2 that extends along the shelf and part of the continental slope, with a maximum depth of 2870 m. (B).Regional maps showing the location of the main rivers that contribute sediments: 1. Atrato River; 2. Mulatos River; 3. Sinú River; and 4. Magdalena River.The location of the NPDB and SCDB is based on [52].The bathymetric data and drainage are based on [50,53] respectively.Abbreviations: NPDB: North Panama Deformed Belt; SCDB: South Caribbean Deformed Belt.(C).Seismic section across the Offshore Sinú Fold Belt, showing thrust, normal faults, and diapirs as the main structures in the continental shelf and continental slope.Adapted from [54][55][56].
The study area in the OSFB is the frontal part of the southwestern portion of an accretionary wedge called the South Caribbean Deformed Belt (SCDB), formed from the subduction of the Caribbean Plate under the South American Plate [55].The Sinú Fold Belt is separated from the North Panama Deformed Belt by the NNW-SSE-trending Uramita Fault (Figure 1A).The belt is configured over a Mesozoic basement [56,57] and involves a sedimentary cover with a thickness of approximately 8 km from the Cenozoic age [58].The study area is mostly established in the continental shelf and only in the upper continental slope.According to Figure 1C [54], the vertical thickness of the wedge is greater at the deformation front, decreasing towards the shelf break.Beneath the deformation front, the depth of detachment is intermediate, and beneath the shelf break it is the lowest.The continental shelf that forms the inner part of the belt is characterized by large normal faults, thrust and related folds [52,[59][60][61], mud diapirs, and mud volcanoes [52,[62][63][64][65][66].The source of shale diapirism, considered the main detachment, is located within the Paleogene shales [52,[63][64][65]67].The external part related to the continental slope is dominated by the NW-trending imbricate thrust, piggy-back basins, thrusts and related folds, tear faults subperpendicular to the thrusts, and shallow normal faults [52,56,59,68].Figure 1C shows a structural-style interpretation across the OSFB in the continental shelf and continental slope, exhibiting a wedge-shaped geometry gently inclined to the NW, and gently dipping due to a basal detachment to the SE [69].

Methods
For this work, the seabed and sub-bottom anomalies in the OSFB were analyzed using high-resolution marine geophysical data.Multibeam bathymetry and SBP data were acquired and used to characterize bottom geomorphology, map the distribution of submarine landslides, examine shallow sub-bottom structures, and discuss the nature of gas occurrences in this region.The acquisition of the information used was carried out in 2018  [52]), onshore geology, and the main tectonic elements (based on [53]).It covers an area of 6700 km 2 that extends along the shelf and part of the continental slope, with a maximum depth of 2870 m. (B).Regional maps showing the location of the main rivers that contribute sediments: 1. Atrato River; 2. Mulatos River; 3. Sinú River; and 4. Magdalena River.The location of the NPDB and SCDB is based on [52].The bathymetric data and drainage are based on [50,53] respectively.Abbreviations: NPDB: North Panama Deformed Belt; SCDB: South Caribbean Deformed Belt.(C).Seismic section across the Offshore Sinú Fold Belt, showing thrust, normal faults, and diapirs as the main structures in the continental shelf and continental slope.Adapted from [54][55][56].
The study area in the OSFB is the frontal part of the southwestern portion of an accretionary wedge called the South Caribbean Deformed Belt (SCDB), formed from the subduction of the Caribbean Plate under the South American Plate [55].The Sinú Fold Belt is separated from the North Panama Deformed Belt by the NNW-SSE-trending Uramita Fault (Figure 1A).The belt is configured over a Mesozoic basement [56,57] and involves a sedimentary cover with a thickness of approximately 8 km from the Cenozoic age [58].The study area is mostly established in the continental shelf and only in the upper continental slope.According to Figure 1C [54], the vertical thickness of the wedge is greater at the deformation front, decreasing towards the shelf break.Beneath the deformation front, the depth of detachment is intermediate, and beneath the shelf break it is the lowest.The continental shelf that forms the inner part of the belt is characterized by large normal faults, thrust and related folds [52,[59][60][61], mud diapirs, and mud volcanoes [52,[62][63][64][65][66].The source of shale diapirism, considered the main detachment, is located within the Paleogene shales [52,[63][64][65]67].The external part related to the continental slope is dominated by the NW-trending imbricate thrust, piggy-back basins, thrusts and related folds, tear faults subperpendicular to the thrusts, and shallow normal faults [52,56,59,68].Figure 1C shows a structural-style interpretation across the OSFB in the continental shelf and continental slope, exhibiting a wedge-shaped geometry gently inclined to the NW, and gently dipping due to a basal detachment to the SE [69].

Methods
For this work, the seabed and sub-bottom anomalies in the OSFB were analyzed using high-resolution marine geophysical data.Multibeam bathymetry and SBP data were acquired and used to characterize bottom geomorphology, map the distribution of submarine landslides, examine shallow sub-bottom structures, and discuss the nature of gas occurrences in this region.The acquisition of the information used was carried out in 2018 aboard the ARC Roncador Research Vessel operated by the Caribbean Oceanographic and Hydrographic Research Center (CIOH).
The Kongsberg EM-302 Multibeam Echosounder (MBES) was used to obtain the bathymetric data.This MBES was designed to carry out soundings up to 7000 m deep, with openings up to 8 km and a frequency of 30 kHz.For the survey, an opening angle of 70 • and a 30% overlap were applied, which ensured 100% seabed coverage.The software used for data acquisition was the Seafloor Information System (SIS), which allows for viewing the topography of the seabed and the depth profile in real time.In addition, SIS constantly monitors the input data to ensure their quality.For positioning, the SEAPATH 330 differential positioning system was deployed, which guarantees centimeter accuracy in the bathymetric survey and uses an MRU-5 motion sensor.The equipment implemented to measure the speed of sound in the water column was a VALEPORT SVS.The CARIS HIPS and SIPS v3.35 software were used to post-process data, with the objective of controlling the quality of the information, eliminating erroneous depth soundings and applying corrections.Finally, a digital elevation model of the area was created covering a surface area of 6700 km 2 with a spatial resolution of 30 m.
High-resolution seismic data were acquired simultaneously via bathymetry, using a SBP300 (Kongsberg, Norway) sub-bottom profiler emitting FM (frequency-modulated) signals.This tool operates with a frequency range of 5 kHz, reaching a maximum vertical resolution of 15 to 25 cm.Eleven high-resolution seismic-acoustic profiles were obtained with an approximate total length of 100 km (Figure 1A).The profiles were processed and analyzed using the Sonar Wiz V5 software.Proper bottom alignment was performed, and chirp removal and gain correction calculations were applied to improve signal intensity and optimize acoustic facies' interpretations.Data generated via the SBP 300 was recorded in the SEG-Y format, which allows for post-processing using some standard seismic processing packages.The processed SBP data were used to identify the shallow deposits.The interpretations were carried out taking into account the tectonic environment, making it possible to visualize bottom and sub-bottom variations associated with the structural style of the area.The sediment cores reported in this work to discuss vent-related d13C values were acquired with a gravity piston corer on board the R/V ARC-Quindío in December 2009 [70].

Seafloor Morphology
The region covers the transition zone between the shelf and the upper continental slope.Unlike the other zones of the Sinú Fold Belt, only the upper continental slope occurs in this region.The bathymetry of the area is the most regular of the belt, suggesting almost no recent tectonic activity.The water depth varies from 18 to 2870 m and the shelf break occurs at a depth of 125 m.The continental shelf exhibits a relatively flat relief with depths between 18 and 120 m.Taking the coastline as a reference, the platform has an amplitude of 37 km towards the northwest.Its geometry shows a thickening towards the Gulf of Urabá that reaches approximately 100 km.Locally, it is interrupted by reef structures, dome structures, and some depressions in the terrain morphology.The former is interpreted in this study as related to occasional emanations of fluids or pockmarks (Figure 2A).
Starting at approximately 125 m, the continental slope is found up to a maximum depth of 2750 m.The morphology of the slope in the area shows profiles of inclined and irregular topography (Figure 2C).Numerous systems of wide channels and turbidity channels (canyons) are formed on the upper continental slope, with steep slopes indicating incision.These systems are slightly sinuous and reach lengths between 5 and 20 km, and do not show a direct connection with an important river on the continent.Cross-sectional topographic profiles along the channels are mainly U-shaped, with walls that reach 7 • (Figure 2C).Some sections of the slope show numerous topside canyons with no evidence of significant mass loss, while other sections are characterized by large-scale mass loss with no canyons or gullies.7° (Figure 2C).Some sections of the slope show numerous topside canyons with no evidence of significant mass loss, while other sections are characterized by large-scale mass loss with no canyons or gullies.The continental slope is highly deformed by mass wasting and shows features indicating faulting and a complex interplay of sediment transport processes, directed downslope along the continental slope and deep into the basin.Figure 2B shows numerous seafloor features generated primarily by these turbidity currents and landslides.It comprises a variety of phenomena that include landslides, rockfalls, and debris flows.We also identified a series of mass-transport complexes (MTCs), mass failures, and slide scars (Figure 2B).MTCs are particularly attached to the continental slope associated with the continental shelf break, with areas of up to 100 km 2 where the deformation associated with mud diapirism overlaps; see Figure 2B.A total of five probable mass failures were identified based on the occurrence of headwall scarps (red dotted lines in Figure 2B), which, together with MTCs, evacuated sediments from the continental shelf and slope.At the base of the slope and on the deep basin seafloor below each of the escarpments, there are large areas of debris fans of positive relief.This morphology reveals the emplacement of rough rubble lobes with blocky, mounded, and irregular seabed morphologies.
Based on the bathymetric model, in the Sinú Fold Belt area, two different shallow gas-related features were identified: (1) dome structures and (2) active and inactive pockmarks.

Dome Structures
Based on bathymetric data, twenty-three submarine dome structures (DS) were identified in the study area.The domes were found at depths between 25 and 1300 m, and some of these volcanoes can be seen in detail (Figure 3).The structures are distributed on The continental slope is highly deformed by mass wasting and shows features indicating faulting and a complex interplay of sediment transport processes, directed downslope along the continental slope and deep into the basin.Figure 2B shows numerous seafloor features generated primarily by these turbidity currents and landslides.It comprises a variety of phenomena that include landslides, rockfalls, and debris flows.We also identified a series of mass-transport complexes (MTCs), mass failures, and slide scars (Figure 2B).MTCs are particularly attached to the continental slope associated with the continental shelf break, with areas of up to 100 km 2 where the deformation associated with mud diapirism overlaps; see Figure 2B.A total of five probable mass failures were identified based on the occurrence of headwall scarps (red dotted lines in Figure 2B), which, together with MTCs, evacuated sediments from the continental shelf and slope.At the base of the slope and on the deep basin seafloor below each of the escarpments, there are large areas of debris fans of positive relief.This morphology reveals the emplacement of rough rubble lobes with blocky, mounded, and irregular seabed morphologies.
Based on the bathymetric model, in the Sinú Fold Belt area, two different shallow gasrelated features were identified: (1) dome structures and (2) active and inactive pockmarks.

Dome Structures
Based on bathymetric data, twenty-three submarine dome structures (DS) were identified in the study area.The domes were found at depths between 25 and 1300 m, and some of these volcanoes can be seen in detail (Figure 3).The structures are distributed on the continental shelf and slope and appear as a dome-conical building that rises above the seabed.The heights oscillate between 15 and 45 m and the diameters at the base oscillate between 250 m and 2 km.Some of the structures are calderas due to the extension of their crater (e.g., DS3, DS4, and DS5).Most of them are isolated and some are grouped (e.g., DS3 and DS7).Some show some type of plume or gas leak above the water column, revealing areas of central vents that probably constitute the source of fluid for the formation of the dome structures, e.g., at SBP1 (Figure 4) and SBP5 (Figure 5). the continental shelf and slope and appear as a dome-conical building that rises above the seabed.The heights oscillate between 15 and 45 m and the diameters at the base oscillate between 250 m and 2 km.Some of the structures are calderas due to the extension of their crater (e.g., DS3, DS4, and DS5).Most of them are isolated and some are grouped (e.g., DS3 and DS7).Some show some type of plume or gas leak above the water column, revealing areas of central vents that probably constitute the source of fluid for the formation of the dome structures, e.g., at SBP1 (Figure 4) and SBP5 (Figure 5).

Pockmarks
The pockmarks are mostly located on the continental shelf.Due to their expression, it is very difficult to differentiate them in the slope area.They appear as local crater-shaped depressions of the seabed with depths between 3 and 15 m and with diameters that reach 50 to 170 m (Figure 2A).These pockmarks indicate abrupt and sudden discharges of gas/fluids from deeper zones towards the sea floor and the water column.An acoustic blanking signal below the pockmark is usually associated with active gas seeps, as SBP7 (Figure 6) and SBP9 (Figure 7) show.The feeder channel is displayed as an acoustic turbid zone in SBPs.This phenomenon of seismic suppression could be due to the existence of gas.

Pockmarks
The pockmarks are mostly located on the continental shelf.Due to their express it is very difficult to differentiate them in the slope area.They appear as local crater-sha depressions of the seabed with depths between 3 and 15 m and with diameters that r 50 to 170 m (Figure 2A).These pockmarks indicate abrupt and sudden discharge gas/fluids from deeper zones towards the sea floor and the water column.An acou blanking signal below the pockmark is usually associated with active gas seeps, as S (Figure 6) and SBP9 (Figure 7) show.The feeder channel is displayed as an acoustic tu zone in SBPs.This phenomenon of seismic suppression could be due to the existenc gas.

Sub-Seafloor Morphology
Seismic anomalies related to gas occurrences in superficial marine sediments of the OSFB have been clearly identified in the SBPs (Figures 4-8).Acoustically transparent zones on the sub-bottom profile sections were detected between 25 and 65 m below the seabed and extend across the entire shelf to the upper continental slope.The occurrences have been classified according to their specific seismic signatures, dimensions, and geometry.
Acoustic blanking (AB) zones appear as a reflection of great amplitude with extensions between 350 and 700 m.Their upper surface has a flat geometry that completely masks the entire underlying sequence due to the total absorption of seismic energy (Figures [4][5][6][7][8].This is the reason why AB shows up as an acoustically transparent zone in many sections.This phenomenon is one of the most common seismic anomalies produced by the occurrence of shallow gas in the area.AB also shows a diffuse disturbance due to the dispersion of various amplitudes, displaying less influence on the sedimentary record.Its

Sub-Seafloor Morphology
Seismic anomalies related to gas occurrences in superficial marine sediments of the OSFB have been clearly identified in the SBPs (Figures 4-8).Acoustically transparent zones on the sub-bottom profile sections were detected between 25 and 65 m below the seabed and extend across the entire shelf to the upper continental slope.The occurrences have been classified according to their specific seismic signatures, dimensions, and geometry.
effect is related to disturbed or slightly discontinuous reflectors whose superficial degassing has led to the blurring of the layer.In some cases, the continuity of the reflectors can be observed, although the generating source is not differentiated.According to Fannin [32], this effect is already generated with only 1% of gas present in sediments.These areas are characterized by chaotic sequences with variable extensions, whose superficial outgassing leads to the aforementioned blurring of the layer (Figures 4-8).
Acoustic blanking is also identified as a vertical disturbance of the sediment layers due to the acoustic dispersion effect of the upward migration of gases or fluids in pores with small amplitude They can be observed in SBP7 (Figure 6) and SBP9 (Figure 7), which reveal amplitudes of approximately 3 m and heights of 20 m (these "heights" are clearly identified as the elevation of the acoustic plume in the water column).They can be found isolated, as in SBP4 (Figure 5), aligned ,as in SBP8 and SBP9 (Figure 7), or in groups, as in SBP7 (Figure 6).They are also seen immediately below the pockmarks, indicating a feeder channel, as in SBP7 (Figure 6) and SBP9 (Figure 7).The origin of these acoustics is located between −50 and −60 m below the seabed.
Enhancement reflectors (ER) show the discontinuous horizons of high acoustic reflection with lengths reaching 1 to 2 km.They are found on both the shelf and the continental slope up to depths of 1500 m.The enhanced reflectors are shown in SBP1 (Figure 4) and SBP6 (Figure 5).Some sub-bottom records reveal expulsions or plumes in the water column, e.g., at SBP1 (Figure 4) and SBP5 (Figure 5).They are located on dome structures, in areas of pockmarks and acoustic blanking, with heights that reach 15 m.Sometimes they are individual forms, and in others, they appear in high density, forming fields of plumes (Figure 7).Some type of acoustic blanking was observed under these forms of gas escape at SBP7 (Figure 6) and SBP9 (Figure 7).
The pockmarks and dome structures can also be recognized better on the SBPs.The pockmarks appear as depressions in the seabed and can be recognized on the SBP records as highly reflective areas.Pockmarks are shown in SBP4, SBP6 (Figure 5), SBP7 (Figure 6), and SBP8 (Figure 7).
Mud volcanoes in the SBP are characterized by acoustic blanking, interpreted as sediments disturbed by fluid migration, by either recent or past activity.They are shown in SBP1, SBP2 (Figure 4), SBP3, SBP5 (Figure 5), SBP10, and SBP11 (Figure 8).

Discussion
The high-resolution seismic analysis together with bathymetry data allowed us to identify and map different acoustic anomalies due to the occurrences of shallow gas in the Offshore Sinú Fold Belt (OSFB).Shallow gassy sediments can particularly be recognized by acoustic anomalies such as attenuation and scattering of the seismic signal [18,71,72].Acoustic blanking (AB) zones appear as a reflection of great amplitude with extensions between 350 and 700 m.Their upper surface has a flat geometry that completely masks the entire underlying sequence due to the total absorption of seismic energy (Figures 4-8).This is the reason why AB shows up as an acoustically transparent zone in many sections.This phenomenon is one of the most common seismic anomalies produced by the occurrence of shallow gas in the area.AB also shows a diffuse disturbance due to the dispersion of various amplitudes, displaying less influence on the sedimentary record.Its effect is related to disturbed or slightly discontinuous reflectors whose superficial degassing has led to the blurring of the layer.In some cases, the continuity of the reflectors can be observed, although the generating source is not differentiated.According to Fannin [32], this effect is already generated with only 1% of gas present in sediments.These areas are characterized by chaotic sequences with variable extensions, whose superficial outgassing leads to the aforementioned blurring of the layer (Figures 4-8).
Acoustic blanking is also identified as a vertical disturbance of the sediment layers due to the acoustic dispersion effect of the upward migration of gases or fluids in pores with small amplitude They can be observed in SBP7 (Figure 6) and SBP9 (Figure 7), which reveal amplitudes of approximately 3 m and heights of 20 m (these "heights" are clearly identified as the elevation of the acoustic plume in the water column).They can be found isolated, as in SBP4 (Figure 5), aligned, as in SBP8 and SBP9 (Figure 7), or in groups, as in SBP7 (Figure 6).They are also seen immediately below the pockmarks, indicating a feeder channel, as in SBP7 (Figure 6) and SBP9 (Figure 7).The origin of these acoustics is located between −50 and −60 m below the seabed.
Enhancement reflectors (ER) show the discontinuous horizons of high acoustic reflection with lengths reaching 1 to 2 km.They are found on both the shelf and the continental slope up to depths of 1500 m.The enhanced reflectors are shown in SBP1 (Figure 4) and SBP6 (Figure 5).Some sub-bottom records reveal expulsions or plumes in the water column, e.g., at SBP1 (Figure 4) and SBP5 (Figure 5).They are located on dome structures, in areas of pockmarks and acoustic blanking, with heights that reach 15 m.Sometimes they are individual forms, and in others, they appear in high density, forming fields of plumes (Figure 7).Some type of acoustic blanking was observed under these forms of gas escape at SBP7 (Figure 6) and SBP9 (Figure 7).
The pockmarks and dome structures can also be recognized better on the SBPs.The pockmarks appear as depressions in the seabed and can be recognized on the SBP records as highly reflective areas.Pockmarks are shown in SBP4, SBP6 (Figure 5), SBP7 (Figure 6), and SBP8 (Figure 7).
Mud volcanoes in the SBP are characterized by acoustic blanking, interpreted as sediments disturbed by fluid migration, by either recent or past activity.They are shown in SBP1, SBP2 (Figure 4), SBP3, SBP5 (Figure 5), SBP10, and SBP11 (Figure 8).

Discussion
The high-resolution seismic analysis together with bathymetry data allowed us to identify and map different acoustic anomalies due to the occurrences of shallow gas in the Offshore Sinú Fold Belt (OSFB).Shallow gassy sediments can particularly be recognized by acoustic anomalies such as attenuation and scattering of the seismic signal [18,71,72].In high-resolution seismic profiles, gas presence in the sediments produces the absorption of most of the high-frequency acoustic energy, masking the deeper horizons.Transparent zones (or anomalous zones without reflections) are found in the seismic profiles producing acoustic blanking, turbidity zones, acoustic voids, acoustic blanking, acoustic turbid layers, and bounded gas upthrust structures [13,25,73,74].These occurrences are common in the OSFB sedimentary record (Figures 4-8).
The predominant acoustic occurrence in the OSFB is the acoustic blanking in the sediments.An area of acoustic blanking in the platform and slope area was analyzed.The profiles also reveal high-frequency, high-amplitude areas such as enhanced reflectors.These higher intensity reflectors indicate variations in acoustic impedance due to high gas flow [75][76][77].They are commonly found in the vicinity of pockmarks, as in SBP6 (Figure 5).This type of event reveals a local gas occurrence in the upper part of the superficial sediments (in the vicinity of the seabed) at depths of 25 to 65 m below the seafloor.The presence of acoustic blanking with a columnar shape is commonly observed in high-resolution profiles, showing a columnar disturbance in the OSFB surface sediments.These vertical anomalies, generally located less than 50 m below the seabed, imply that gas has propagated upward from the deep subsurface to surface sediments [13,76,[78][79][80], indicating important migration pathways [81].
The most common surface expressions associated with gas seeps on the continental margins are pockmarks and dome structures [82][83][84].The pockmarks' origin is related to the expulsion of shallow gas pockets [15,36].The mechanism of the formation of pockmarks is strongly controlled by local sedimentation processes, such as sediment type, thickness, permeability, pore pressure, and fluid content [20].Their development requires the presence of fine granular sediments to support their structure and generate sediment removal, local subsidence, and escape failures [20].The results allow us to establish that pockmarks are very common characteristics in the OSFB.Pockmarks were observed on the continental slope in water depths between 20 and 2000 m.These pockmarks generate sinkholes that reach a depth of 15 m and are commonly found in the vicinity of dome structures.These depressions are produced by gas migrations, which displace part of the sediment and generate conical depressions [85,86].The presence of these structures in the area indicates that the seafloor morphology may be formed by silty sediment packages and fine-grained clays (as suggested by [47,83,84]), which is mainly related to the sediment supply of the Atrato, Sinú, and Magdalena rivers.Bottom profiles reveal inactive pockmarks as well as pockmarks associated with active gas seeps in the OSFB.Most pockmarks reveal the presence of zones of acoustic turbidity and enhanced reflectors below their base, suggesting a continuous supply of gases [87][88][89].Some acoustic blanking also occurs under pockmarks, indicating active vertical gas leaks [20].Although some pockmarks do not show clear evidence of recent activity, their formation suggests that they are remnants of past active degassing [36].
Dome structures have been reported in many places around the world [90][91][92][93][94][95][96][97][98][99][100].They have been widely studied given their close relationship with the occurrence of hydrocarbons, and they also constitute the most important migration pathways of methane and CO 2 from marine sediments to the atmosphere [17,91,95].Around 23 dome structures have been identified in water depths between 20 and 2000 m in the OSFB.They are characterized by internal dispersion and acoustic blanking in seismic records.These behaviors are interpreted as sequences of shallow sediments that have been disturbed by the upward migration of fluids that arrived near the seabed.
Various acoustic plumes over dome structures or pockmarks were observed in the seismic records.Their presence is interpreted as gas bubbles escaping into the water column, indicating active degassing, given the abundance of gas in subsoil sediments at high pressure conditions [18,101].If the free gas in the sediments cannot be held below the seabed, it will eventually escape and form gas plumes [102].The abundance of observed gas plumes in the bottom profiles indicates that gas is abundant in the surface sediments.
The origin of gas in shallow marine sediments is associated with thermogenic and/or biogenic processes [15,22,38,73,103].In both cases, the gases result from organic matter via thermogenic processes that depend on pressure and temperature, or as a result of bacterial activity.Since the discovery of the Chuchupa and Ballenas fields in the 1970s [43,44], the offshore Colombian Caribbean has been considered a gas province originating mainly from microbial activity [45].The latest discoveries in areas located below 2000 m with the drilling of the Orca (2014), Kronos (2015), Purple Angel (2017) and Gorgon (2022) exploratory wells confirm the presence of an emerging basin dominated by methane [46].This suggests occurrences not only of biogenic origin but also the thermogenic generation of hydrocarbons [46,104].
Rua et al. (2016) [70] analyzed sediment cores in the southern Offshore Sinú Fold Belt, collected via a gravity piston at a site that showed acoustic anomalies in high-resolution subsoil profiles (see Figure 1 for location).The carbon isotopic signature measurements on six samples revealed low δ13C values (as low as 30%).These results suggest that the origin of the superficial gas occurrences is a deep thermogenic source [73,103,105].Gracia et al. (2012) [106] reported, during the INVEMAR-ANH I project carried out in 2008, the discovery of new chemosymbiotic species in the south of the Sinú Fold Belt.Several fragments and shells corresponding to families of chemosymbiotic bivalves were collected, including two morphotypes of the family Lucinidae (Graecina colombiensis) and three unidentified species of Lucinoma.The occurrence of these biological communities dominated by chemosymbiotic metazoans is highly endemic to reducing conditions, such as those observed in environments where there are seeps of methane or other hydrocarbons [107,108].Other indicators of seepage activity, such as pieces of authigenic carbonates, were also sampled.This is a typical feature of sites where methane is expelled from the seabed.The data presented above [70,106], including our results, biologically and geologically support the presence of gas seeps on the seabed in the area of the Offshore Sinú Fold Belt.
Natural gas is considered one of the main sources of energy in the world [109,110].
In Colombia, it is the second-most used energy source.The expansion of this area could significantly contribute to supporting the country's energy security and open the possibility of future gas exports, offering a fundamental energy source for the energy transition that the Colombian government has been developing.

Conclusions
One of the primary sources of energy across the globe is natural gas, which is the reason that numerous explorations with different methodologies are being carried out worldwide.
The presence of several acoustic anomalies in the Offshore Sinú Fold Belt (OSFB) was detected via the combined use of high-resolution seismic analysis and bathymetry data.This method made possible the identification and mapping of several shallow gassy sediments.
Gas-related features observed on the very high-resolution acoustic surveys include blanking, enhanced reflections, acoustic plumes, pockmarks, and dome structures.
Seismic records show a variety of sonic plumes over dome structures or pockmarks.Moreover, minor gas escape was observed throughout the water column as gas plumes.
It was also detected that, in some places, shallow gas continues under the sea floor at about -60 m depth, often marked by step-like offsets, indicating different levels of gas penetration, possibly due to interbedded deposits of fine sand and mud.
Finally, the results indicate that gas is abundant in the surface sediments (as evidenced by the frequency of gas plumes detected in the bottom profiles) and, thus, the use of more expensive exploration methods is not only justified but recommended.

Figure 1 .
Figure 1.(A).Location of the study area showing the digital elevation model (based on[52]), onshore geology, and the main tectonic elements (based on[53]).It covers an area of 6700 km 2 that extends along the shelf and part of the continental slope, with a maximum depth of 2870 m. (B).Regional maps showing the location of the main rivers that contribute sediments: 1. Atrato River; 2. Mulatos River; 3. Sinú River; and 4. Magdalena River.The location of the NPDB and SCDB is based on[52].The bathymetric data and drainage are based on[50,53] respectively.Abbreviations: NPDB: North Panama Deformed Belt; SCDB: South Caribbean Deformed Belt.(C).Seismic section across the Offshore Sinú Fold Belt, showing thrust, normal faults, and diapirs as the main structures in the continental shelf and continental slope.Adapted from[54][55][56].

Figure 1 .
Figure 1.(A).Location of the study area showing the digital elevation model (based on[52]), onshore geology, and the main tectonic elements (based on[53]).It covers an area of 6700 km 2 that extends along the shelf and part of the continental slope, with a maximum depth of 2870 m. (B).Regional maps showing the location of the main rivers that contribute sediments: 1. Atrato River; 2. Mulatos River; 3. Sinú River; and 4. Magdalena River.The location of the NPDB and SCDB is based on[52].The bathymetric data and drainage are based on[50,53] respectively.Abbreviations: NPDB: North Panama Deformed Belt; SCDB: South Caribbean Deformed Belt.(C).Seismic section across the Offshore Sinú Fold Belt, showing thrust, normal faults, and diapirs as the main structures in the continental shelf and continental slope.Adapted from[54][55][56].

Figure 2 .
Figure 2. (A).Multibeam bathymetric relief of the study area with spatial distribution of the gas occurrence features; (B).3D view of the multibeam bathymetry.Continental slope has been modified by mass wasting events; (C).Profiles across the continental slope; (D).Cross-sectional topographic profile along a channel dike complex.

MTC 3 Figure 2 .
Figure 2. (A).Multibeam bathymetric relief of the study area with spatial distribution of the gas occurrence features; (B).3D view of the multibeam bathymetry.Continental slope has been modified by mass wasting events; (C).Profiles across the continental slope; (D).Cross-sectional topographic profile along a channel dike complex.

Figure 3 .
Figure 3. Morphology of submarine volcanoes identified in the area of the shelf and the continental slope.Most of them appear isolated and some appear in groups (See Figure 2 for locations).Figure 3. Morphology of submarine volcanoes identified in the area of the shelf and the continental slope.Most of them appear isolated and some appear in groups (See Figure 2 for locations).

Figure 3 . 16 Figure 4 .
Figure 3. Morphology of submarine volcanoes identified in the area of the shelf and the continental slope.Most of them appear isolated and some appear in groups (See Figure 2 for locations).Figure 3. Morphology of submarine volcanoes identified in the area of the shelf and the continental slope.Most of them appear isolated and some appear in groups (See Figure 2 for locations).J. Mar.Sci.Eng.2023, 11, x FOR PEER REVIEW 7 of 16

Figure 6 .
Figure 6.Sub-bottom profiler 7 presenting acoustic blanking (AB).Several plumes are observed the pockmarks (P) and active pockmarks (AP) that reveal active escapes of gas from the soil tow the water.

Figure 6 .
Figure 6.Sub-bottom profiler 7 presenting acoustic blanking (AB).Several plumes are observed over the pockmarks (P) and active pockmarks (AP) that reveal active escapes of gas from the soil towards the water.

Figure 6 .
Figure6.Sub-bottom profiler 7 presenting acoustic blanking (AB).Several plumes are observed over the pockmarks (P) and active pockmarks (AP) that reveal active escapes of gas from the soil towards the water.

Figure 7 .
Figure 7. Sub-bottom profilers 8 and 9 showing numerous pockmarks (P) and acoustic blanking (AB) in the area of the shelf break.

Figure 7 .
Figure 7. Sub-bottom profilers 8 and 9 showing numerous pockmarks (P) and acoustic blanking (AB) in the area of the shelf break.