Sediment and water samples were collected from three stations located in the Gdańsk Basin (southeastern Baltic Sea) at water depths of 50–103 m: ZGG, MET2, and MET1-MP (Fig. S1). Station MET1-MP represents an active deepwater pockmark influenced by both methane seepage and freshwater infiltration. Station MET2 is characterized by shallow methane accumulation without distinct seabed morphology and freshwater infiltration (Jas´niewicz et al., 2019; Brodecka-Goluch et al., 2022). Station ZGG, designated as the reference site, shows no evidence of methane occurrence or freshwater discharge. Sampling was conducted from RV Oceanograf in July 2023 and February 2024. Hydroacoustic surveys were performed prior to coring using a Reson Teledyne SeaBat 7125 multibeam echosounder, an EdgeTech SB-216S subbottom profiler, and a Simrad EK80 split-beam echosounder to determine precise sampling locations.

Sediment cores were collected with a Rumohr-Lot gravity corer equipped with Plexiglas liners (7.5 cm diameter, 150 cm length). At each station, seven cores (90–110 cm long) were collected during the July 2023 sampling campaign and additional seven cores in February 2024. These cores were used for analyses of: (1) total organic carbon (TOC), total nitrogen (TN), total sulfur (TS), stable isotopes of organic matter (δ13C, δ15N), loss on ignition (LOI), and water content (W); (2) CH₄, CO₂ and their stable isotopes (δ13C-CH₄, δ13C-CO₂, δ2H-CH₄), and oxidation–reduction potential (ORP); (3) pH and total alkalinity (TA); (4) dissolved inorganic carbon (DIC) and δ13C-DIC; (5) sulfate (SO42-), chloride (Cl⁻), calcium (Ca²⁺), magnesium (Mg²⁺), hydrogen sulfide (∑H₂S, the sum of H₂S, HS⁻ and S²⁻), ammonia (∑NH4+, the sum of NH4+ and NH₃), phosphate (PO43-), manganese (Mn²⁺) and iron (Fe²⁺) in pore water; (6) mineralogy and micromorphology; and (7) experimental rates of organoclastic sulfate reduction (ROSR) and anaerobic oxidation of methane (RAOM).

Pore waters were extracted from intact, sealed cores using Rhizon^® samplers (0.15 µm pore size). For each core liner, 4 mm holes matching the Rhizon^® sampler fitting were pre-drilled at 5 cm intervals and sealed with a tape. After cores retrieval aboard, the tape was pierced and samplers inserted; to minimize contact with air, cores remained sealed throughout porewater collection. The uppermost sampler (1–4 cm above the sediment surface, depending on station) collected bottom water overlying the core. Porewater was drawn into attached 20 cm³ plastic syringes. To minimize oxidation, all subsequent subsampling and treatment were performed immediately after porewater sampling. Immediately after extraction, TA, pH, ∑H₂S, ∑NH4+, and PO43- were measured onboard. Samples for major ions, Mn²⁺, and Fe²⁺ were acidified with HNO₃. DIC samples were poisoned with HgCl₂ and stored in sealed, nearly full tubes (< 1 % headspace) at 4 °C. Sampling for carbon isotopic composition of DIC (δ13C-DIC) was carried out in an N₂ glove box; porewater samples were collected into 2 cm³ septum vials that were filled by injecting porewater with one needle while displacing N₂ through a second needle; vials were completely filled with porewater and stored at 4 °C prior to analysis.

For CH₄ analysis, sediment was subsampled through 12 mm holes in the liners and transferred immediately to glass vials containing 2.5 % NaOH, which were sealed with butyl rubber stoppers and aluminum crimp caps (Jørgensen et al., 2001). NaOH was added to halt microbial activity, preventing biological production or consumption of CH₄. A headspace was equilibrated by shaking, and vials were stored caps-down at 4 °C until gas chromatographic analysis. Samples for carbon isotopic composition of methane and carbon dioxide (δ13C-CH₄, δ13C-CO₂) were collected in the same manner as the CH₄ samples but the procedure was carried out inside an N₂ glove box. Instead of adding NaOH, vials were pre-filled with Milli Q water.

Cores designated for solid-phase analyses were sectioned at 5 cm intervals using a PVC ring and plastic spatula; samples were taken from the central part of each section. Each sample was placed in a separate polyethylene bag. Samples for microbial analyses were collected from five depth intervals (0–20, 20–40, 40–60, 60–80, and 80–100 cm) using sterile tools. Samples for mineralogical analyses were taken from 10–15, 40–45, and 80–85 cm. All sediment samples were stored at -21 °C until analysis.

OM content is highest in surface sediments (0–20 cm) at all stations (Fig. 1). TOC correlates positively with TN (Fig. 1a, b), indicating that nitrogen is predominantly associated with organic matter. TOC ranges from 2.8 wt % to 10.7 wt %, TN from 0.2 wt % to 1.4 wt %, and TS from 0.8 wt % to 3.6 wt % (Fig. 1a–c). The molar TOC : TN ratio varies between 7.4 and 13.7 (Fig. 1d).

The highest OM contents occur at station MET1-MP. TOC ranges from 5.5 wt % to 10.7 wt %, and TN from 0.6 wt % to 1.4 wt %. This station also shows the strongest downcore variability. The lowest OM contents are observed at MET2, with similar values in both seasons. TOC ranges from 2.5 wt % to 5.0 wt % and TN from 0.2 wt % to 0.5 wt %. Station ZGG exhibits intermediate values. TOC ranges from 3.0 wt % to 6.3 wt % and TN from 0.3 wt % to 0.8 wt %.

TS contents exceed 1 wt % at most depths (Fig. 1c). The highest average TS (2.2 wt %) occurs at ZGG. Maximum values (≥ 3.2 wt %) are recorded at 70–75 cm in both seasons. The lowest average TS (1.2 wt %) occurs at MET2. At MET1-MP, mean TS is 1.7 wt %. The TOC : TN ratio at MET1-MP ranges from 7.4 to 9.1. At ZGG, values range from 8.7 to 11.9. The highest ratios occur at MET2 (9.3–13.7).

Sediment oxidation–reduction potential (ORP) ranges from -502 to -5 mV. The highest values occur in surface layers (Fig. S2). Loss on ignition (LOI) ranges from 5.5 wt % to 20.2 wt % (Fig. S2). Values are highest at MET1-MP (13.5 wt %–20.2 wt %), intermediate at ZGG (6.5 wt %–14.4 wt %), and lowest at MET2 (5.5 wt %–10.5 wt %). Water content (W) follows a similar pattern (Fig. S2). The highest values occur at MET1-MP (67.4 wt %–97.7 wt %). At ZGG, W ranges from 66.0 wt % to 90.6 wt %. At MET2, W ranges from 50.4 wt % to 79.5 wt %.

The isotopic composition of organic carbon (δ13C) ranges from -27.3 ‰ to -24.4 ‰, and that of nitrogen (δ15N) from 1.2 ‰ to 6.7 ‰ (Fig. 1e, f). The lowest δ13C values occur at station ZGG in both seasons. The highest values are observed at MET2 in summer. The largest δ13C variability (∼ 1.5 ‰) is recorded at ZGG (-27.3 ‰ to -25.5 ‰) and MET2 (-26.0 ‰ to -24.4 ‰). At MET1-MP, δ13C varies within a narrower range (< 1 ‰), from -26.2 ‰ to -25.2 ‰. The greatest δ15N variability occurs at MET2 (1.2 ‰–6.7 ‰). At this station, deeper sediments are depleted in ¹⁵N relative to surface layers. The depletion reaches ∼ 3 ‰ in the deepest intervals. At MET1-MP and ZGG, δ15N varies over a narrower range (2.0 ‰–4.5 ‰).

Results from summer and winter show that spatial differences among stations exceed seasonal variability (Figs. 2, 3). Depth profiles of major ions (Cl⁻, SO42-, Ca²⁺, Mg²⁺) are shown in Fig. 2. The concentration of Cl⁻ ranges from 70.30 to 180.40 mmol dm⁻³, SO42- from 0.03 to 9.23 mmol dm⁻³, Ca²⁺ from 2.52 to 4.84 mmol dm⁻³, while Mg²⁺ from 8.34 to 18.37 mmol dm⁻³.

Major-ion concentrations in bottom water and surface-most sediments reflect salinity. Salinity is highest at ZGG (12.2 in summer; 12.1 in winter) and lowest at MET2 (8.0 in summer; 7.5 in winter). At MET2, Cl⁻, Ca²⁺, and Mg²⁺ remain nearly constant with depth. In contrast, these ions decrease markedly with depth at MET1-MP. At ZGG, Ca²⁺ increases significantly downcore. Sulfate remains relatively high (≥ 3.5 mmol dm⁻³) throughout the ZGG profile. At MET1-MP and MET2, sulfate is restricted to the upper 30–40 cm.

∑H₂S reaches maximum concentrations at ZGG (up to 2026 µmol dm⁻³). At MET1-MP and MET2, ∑H₂S ranges from < LOQ to 1783 µmol dm⁻³ and is confined to discrete sediment intervals (Fig. 2e), corresponding to zones of sulfate-dependent methane oxidation. Bottom waters at MET2 are oxic (3.1 mg dm⁻³ O₂ in summer; 5.9 mg dm⁻³ in winter). At ZGG and MET1-MP, O₂ concentrations are < 0.2 in summer and < 1.7 mg dm⁻³ in winter.

Pore-water Fe²⁺ concentrations are generally < 30 µmol dm⁻³ (Fig. 2h). Elevated values occur only at MET1-MP, reaching 139.3 µmol dm⁻³ below 60 cm in summer. Mn²⁺ concentrations range from < LOQ to 51.6 µmol dm⁻³, with highest values at ZGG (Fig. 2g).

∑NH4+ ranges from 22.2 to 3794.9 µmol dm⁻³. PO43- ranges from 1.4 to 950.2 µmol dm⁻³. Both species attain maximum concentrations at MET1-MP (Fig. 3d, e). At ZGG and MET2, concentrations are lower, ranging from 22.2 to 2027.0 µmol dm⁻³ for ∑NH4+ and from 1.4 to 361.1 µmol dm⁻³ for PO43-.

The pH at all sampling sites was from 7.00 to 8.31, with an average of 7.57 ± 0.21 at station MET1-MP, 7.66 ± 0.26 at station MET2, and 7.57 ± 0.20 at station ZGG (Fig. S2).

The TA in the pore water varied in the range of 1911.3–32 747.2 µmol kg⁻¹. Overall, the highest TA in pore waters was found in the sediment of pockmark MET1-MP, while the lowest values were observed at ZGG, similarly to DIC, which is the main component of TA in seawater (Fig. 3a).

DIC concentrations range from 1851 to 33 323 µmol kg⁻¹ (Fig. 3b). The highest values occur at MET1-MP. The lowest concentrations and the smallest downcore variability are observed at the non-methane station ZGG.

The δ13C-DIC values range from -16.2 ‰ to 17.3 ‰ (Fig. 3c). At ZGG, δ13C-DIC remains negative throughout the profile. At MET2, negative values are restricted to the upper ∼ 25 cm, whereas deeper layers show positive values. At MET1-MP, δ13C-DIC exhibits limited variability. Negative values occur only in the surface layer (0–5 cm) during summer.

Methane concentrations at MET1-MP range from 0.17 to 3.24 mmol dm⁻³, and at MET2 from < LOQ to 8.49 mmol dm⁻³ (Fig. 2f). At MET1-MP, maximum CH₄ concentrations occur near the sediment surface (∼ 5 cm b.s.f.). Concentrations decrease with depth and stabilize at ∼ 2 mmol dm⁻³. At MET2, CH₄ increases progressively with depth and reaches a maximum at ≤ 80 cm b.s.f. Methane is not detected at ZGG.

Depth profiles of δ13C-CH₄, δ2H-CH₄, and δ13C-CO₂ at MET1-MP and MET2 are shown in Fig. S3. At MET1-MP, δ13C-CH₄ ranges from -62.6 ‰ to -55.8 ‰. At MET2, values are lower, from -76.6 ‰ to -69.6 ‰. The δ2H-CH₄ values range from -271 ‰ to -257 ‰ at MET1-MP and from -250 ‰ to -239 ‰ at MET2. At MET2, δ13C-CO₂ remains negative throughout the profile (-8.1 ‰ to -2.6 ‰). At MET1-MP, δ13C-CO₂ ranges from -1.0 ‰ to 5.7 ‰ and is predominantly positive, except in the surface layer and at ∼ 90 cm b.s.f.

The measured rates of organoclastic sulfate reduction (ROSR) and anaerobic oxidation of methane (RAOM) presented in Table 1 represent the maximum metabolic capacity of the microbial community (potential rates); incubations were conducted without substrate limitation and therefore reflect relative microbial activity between samples rather than absolute in situ rates.

ROSR ranges from 15 to 164 µmol dm⁻³ d⁻¹ (Table 1). At ZGG, values range from 24 to 46 µmol dm⁻³ d⁻¹. At MET2, rates vary between 21 and 48 µmol dm⁻³ d⁻¹. The highest variability occurs at MET1-MP (16–163 µmol dm⁻³ d⁻¹). At MET1-MP and MET2, maximum rates are observed in the 0–20 cm interval. At ZGG, the highest ROSR occurs at 20–40 cm depth.

RAOM ranges from 6 to 1077 µmol dm⁻³ d⁻¹. At MET1-MP, rates are elevated throughout the profile (59–1077 µmol dm⁻³ d⁻¹). Maximum values occur at 40–60 and 80–100 cm depth (Table 1). At MET2, RAOM is lower (6–33 µmol dm⁻³ d⁻¹) and is concentrated in surface sediments. At ZGG, rates are relatively uniform with depth (8–10 µmol dm⁻³ d⁻¹). Additional raw data concerning the determination of OSR and AOM are provided in the Supplement (Tables S2 and S3).

Depth-integrated ROSR and RAOM (0–100 cm) and the corresponding DIC production are summarized in Table 1. DIC production from OSR is comparable among stations. In contrast, MET1-MP exhibits exceptionally high RAOM and associated DIC production, exceeding values at the other stations by two orders of magnitude. The lowest process rates and DIC production occur at ZGG.

The studied sediments are fine-grained (Brodecka et al., 2013) and classified as silt (MET1-MP), sandy silt (MET2), and clayey silt (ZGG). They are dominated by a matrix composed of clay minerals, biogenic silica, and organic matter (Fig. S5). This matrix constitutes ∼ 20 vol %–50 vol % of the sediment and generally increases with depth. At MET1-MP it ranges from 30 % to 54 %, at MET2 from 20 % to 39 %, and at ZGG from 36 % to 53 %.

Detrital aluminosilicates (mainly quartz, feldspars, and micas) occur as dispersed grains within this matrix. Their proportion remains relatively constant at MET2 (∼ 30 %–35 %) but decreases with depth at MET1-MP (13 % to 8 %) and ZGG (24 % to 13 %).

PXRD analyses of bulk samples (Fig. S6) and clay fractions (Fig. S7) indicate that phyllosilicates are dominated by mica/illite, chlorite, and kaolinite. A ∼ 17 Å reflection in glycolated samples that disappears upon heating indicates mixed-layer illite–smectite. The position of the (060) reflection is consistent with dioctahedral phyllosilicates (Moore and Reynolds, 1997). Microscopy also reveals glauconite, an Fe-rich mica-group mineral. It occurs sporadically at MET1-MP and MET2 and is more common at ZGG.

Authigenic carbonates and sulfides are minor components (< several wt %; Table 2). Pyrite (FeS₂) is the dominant sulfide. It occurs as isolated euhedral crystals (cubic or octahedral, ≤ 1 µm) or as framboids with diameters of a few to ∼ 20 µm (Fig. 4). Irregular framboidal aggregates also occur. Pyrite framboids are commonly associated with diatom frustules or organic-rich domains. Microscopy and PXRD indicate higher average pyrite contents at ZGG than at MET1-MP and MET2 (Table 2).

Dolomite [CaMg(CO₃)₂] is the predominant carbonate phase. Its presence is indicated by the 2.88 Å reflection in PXRD patterns (Fig. S6). Due to low concentrations, unambiguous interpretation is not possible because diffraction peaks (101), (015) and (021) indicative of space group of dolomite (Gregg et al., 2015; Lukoczki et al., 2026) are not apparent. However, SEM–EDS analyses (Fig. S8) show near-equimolar Ca and Mg, consistent with dolomite rather than high-Mg calcite. Also Raman microspectroscopy of carbonate crystals revealed typical dolomite (and not Mg-calcite) features (Fig. S9): strong sharp bands at 1100 (symmetric stretching), 305 (librational mode), and 180 (translational mode), as well as weak bands at 730 and 1445 cm⁻¹, related to in-plane bending and asymmetric stretching, respectively (Sun et al., 2014; Zheng et al., 2021; Alves et al., 2023). Comparison of the Raman spectra with that of proper dolomite (Rędziny quarry, Southern Poland, Kowalski et al., 1976) indicates that the band broadening (Sun et al., 2014; Zheng et al., 2021) that should be associated with structural disorder (expected in protodolomite and high-Mg calcite) does not occur, as the band half-widths are of the same order as for crystalline dolomite (Fig. S9b). Moreover, euhedral to subhedral rhombohedral habits are typical for authigenic dolomite (Fig. 4), and not nanocrystalline aggregates characteristic of protodolomite (Liu et al., 2020; Zheng et al., 2021). Therefore, despite the lack of unambiguous PXRD evidence, both the chemical composition, morphology, and Raman features suggest that the dominant carbonate in the sediments is dolomite.

Dolomite is most abundant at MET1-MP, slightly less abundant at MET2, and nearly absent at ZGG (Table 2). At MET1-MP, dolomite forms euhedral crystals several to ∼ 20 µm in size (Fig. 4). At MET2, both fine (< 20 µm) euhedral or subhedral crystals and larger subhedral to anhedral grains occur. The euhedral/subhedral morphology of many of the crystals indicates that they are authigenic. At ZGG, small (5–10 µm) subhedral to anhedral grains are restricted to the subsurface layer (10–15 cm).

All sediments contain trace amounts of calcite (CaCO₃), occurring as fine (< 10 µm) subhedral grains (Fig. 4h). Isolated euhedral siderite (FeCO₃) crystals were identified only at MET2 (Fig. 4i).

All Mössbauer spectra (Fig. S10; Table S4) are similar and dominated by Fe(II) components. The principal doublet, comprising ∼ 50 %–60 % of the spectral area, exhibits an isomer shift (IS) of ∼ 0.3 and quadrupole splitting (QS) of ∼ 0.6 mm s⁻¹. These parameters are characteristic of Fe(II) in pyrite (Byrne and Kappler, 2019). Two additional doublets contribute ∼ 25 %–30 % of the spectral area. They display high IS and QS values typical of Fe(II) in silicates (Dyar et al., 2006), although those of the minority doublet are also similar to Fe(II) in siderite (Rothwell et al., 2025). A further doublet, accounting for ∼ 10 %–28 % of the spectral area, shows relatively low IS (∼ 0.4–0.6 mm s⁻¹) and moderate QS (∼ 0.7–0.9 mm s⁻¹). These parameters indicate high-spin Fe(III) in octahedral coordination. The precise mineral host of this component cannot be resolved. The hyperfine parameters are consistent with paramagnetic Fe(III) oxyhydroxides at room temperature (Byrne and Kappler, 2019). However, Fe(III) incorporated in phyllosilicates such as micas/illite, chlorites, or glauconite is also plausible (De Grave et al., 1987).

Carbonate burial rates were calculated using dolomite contents, as dolomite is the dominant carbonate phase. Pyrite burial rates were also determined due to pronounced inter-station variability. The highest dolomite burial occurs at MET1-MP (Table 2). Rates are 484 ± 13 based on PXRD and 888 ± 87 µmol m⁻² d⁻¹ based on SEM–EDS. Intermediate values are observed at MET2 (254 ± 58 and 517 ± 7 µmol m⁻² d⁻¹, respectively). The lowest rates occur at ZGG (60 ± 35 from PXRD and 136 ± 110 µmol m⁻² d⁻¹ from SEM–EDS). In contrast, pyrite burial is highest at ZGG. Rates reach 959 ± 16 (PXRD) and 1341 ± 87 µmol m⁻² d⁻¹ (SEM–EDS). At MET1-MP, values are 556 ± 18 and 640 ± 153 µmol m⁻² d⁻¹, respectively. At MET2, pyrite burial rates are 341 ± 5 (PXRD) and 785 ± 4 µmol m⁻² d⁻¹ (SEM–EDS).

Nutrient enrichment associated with intensified fertilizer application has promoted widespread coastal eutrophication. Elevated nitrogen and phosphorus inputs stimulate primary production and increase the flux of OM to the seafloor. Enhanced organic carbon delivery raises sedimentary oxygen demand and promotes hypoxia and anoxia in bottom waters (Rabalais et al., 2010). These effects are amplified under strong density stratification, which limits vertical mixing and oxygen replenishment from surface layers (Carstensen et al., 2014). Under anoxic conditions, sulfate reduction becomes the dominant mineralization pathway, leading to hydrogen sulfide accumulation in bottom water.

Station ZGG, located in the western Gdańsk Deep in the southern Baltic Sea (Fig. S1), represents such an environment. The seafloor is covered by fine-grained sediments with mean water and OM contents of 76 % and 11 %, respectively, characteristic of an accumulation bottom type (Jönsson et al., 2005). TOC and TN values fall within the range typical for deep Baltic sediments (Us´cinowicz, 2011). Maximum concentrations occur in the uppermost centimeters (Fig. 1a, b), reflecting enhanced eutrophication.

The molar TOC : TN ratio ranges from 8.5 to 11.8, exceeding typical marine values of 4–8 and indicating a contribution of terrestrial OM, for which TOC : TN commonly exceeds 12 (Meyers, 1997; Carneiro et al., 2021; Yadav et al., 2025). Carbon isotopic compositions further support this interpretation. δ13C values at ZGG range from -27.3 ‰ to -25.5 ‰ (Fig. 1e). These values are more negative than those of marine phytoplankton (-18 ‰ to -24 ‰; Yadav et al., 2025) and typical mid-latitude marine OM (-22 ‰ to -20 ‰; Meyers, 1994), and fall within the range characteristic of terrestrial C₃ plants (-35 ‰ to -25 ‰; average ∼-27 ‰). The combined TOC : TN and δ13C data indicate a significant input of vascular plant-derived organic carbon, which is carbon-rich and isotopically depleted relative to marine primary production (Meyers, 1994, 1997), in all the studied sediments. Riverine discharge is a major source of terrestrial carbon, contributing up to 10 %–30 % of the total carbon pool in Baltic sediments (Miltner and Emeis, 2001). Nitrogen isotopic compositions (δ15N=2.0 ‰–3.6 ‰; Fig. 1f) are at or slightly below the typical range for marine OM (3 ‰–12 ‰; Maksymowska et al., 2000; Wada and Hattori, 1991). In contrast, terrestrial C₃ plants exhibit a broader range (-10 ‰ to +10 ‰). The relatively low δ15N values at ZGG therefore further support a substantial allochthonous contribution. Increased biomass development may also lead to a modification of the isotopic composition of both nitrogen and carbon in sediments, the effect of which may also be visible in the case described. For example, a decrease in the δ15N value of OM can be caused by the addition of isotopically light nitrogen to the N pool, bound by diazotrophic cyanobacteria, as in the case described by Struck et al. (2000) for the Gotland Basin. In turn, increased primary production can lead to an increase in ¹³C content in sediment organic matter. For example, in the Oder Haff an increase of δ13C in OM of 1.7 ‰ was found, while in the Gotland Basin it was 2.6 ‰ (Struck et al., 2000).

Sediments at ZGG are also characterized by elevated TS contents (1.59 wt %–3.58 wt %; Fig. 1c). The TOC : TS ratio provides insight into depositional redox and salinity conditions (Berner et al., 1979; Berner and Raiswell, 1984). A ratio of 2.8 ± 0.4 is typical of normal marine sediments underlying oxygenated waters (Berner, 1982). Lower values indicate enhanced sulfide formation under anoxic conditions, whereas TOC : TS > 2.8 is characteristic of freshwater environments. The mean TOC : TS ratio at ZGG (2.5 ± 0.6) is slightly lower below the marine oxic benchmark. This suggests suboxic to anoxic depositional conditions, controlled by substantial water depth, persistent stratification, and bottom-water oxygen depletion.

In the uppermost centimeters, pore-water composition closely resembles that of bottom water (Fig. 2). With increasing depth, diagenetic processes progressively modify the chemical composition. Sulfate concentrations decrease rapidly within the upper ∼ 30 cm and stabilize at ∼ 3.5 mmol dm⁻³ in deeper layers (Fig. 2d). In parallel, H₂S concentrations increase (Fig. 2e), reflecting active microbial sulfate reduction. Concentrations of NH4+, PO43-, DIC, and TA also increase with depth (Fig. 3), indicating that a common mineralization process supplies these components to pore waters. Sulfate can be consumed during both OSR (Table S1: R5) and AOM (Table S1: R14). Although methanogenesis and AOM may occur deeper in the sediment at ZGG, sulfate remains relatively abundant throughout the studied profile, and CH₄ concentrations are below LOQ (Fig. 2f). These observations indicate that sampling captured the sulfate reduction zone above the sulfate–methane transition (SMT). This interpretation is consistent with the inhibitory effect of sulfate on methanogenesis, as sulfate-reducing microorganisms outcompete methanogens for H₂ and acetate (Reeburgh, 2007). The δ13C-DIC profile further supports this conclusion. δ13C-DIC values decrease with depth (Fig. 3c), whereas methanogenesis typically increases δ13C-DIC due to preferential removal of ¹²C into CH₄ (Whiticar, 1999). Moreover, DIC and δ13C-DIC exhibit a strong linear correlation (R2= 0.98), and no δ13C-DIC minimum characteristic of the SMT is observed (Malinverno and Pohlman, 2011). The slope of the regression corresponds to δ13C values typical of marine OM (-22 ‰ to -20 ‰; Meyers, 1994), indicating that DIC is derived predominantly from organic carbon oxidation (Chen et al., 2024).

These data demonstrate that OSR is the principal process consuming SO42- and controlling pore-water DIC and TA at ZGG. According to stoichiometry, OSR produces two moles of DIC per mole of sulfate reduced (Table S1: R5), whereas sulfate-driven AOM yields a 1 : 1 ratio (Table S1: R14). The relationship between the increase in DIC (ΔDIC) and the decrease in sulfate (ΔSO42-) therefore provides a diagnostic criterion for distinguishing these processes (Masuzawa et al., 1992; Rassmann et al., 2016). At ZGG, the ΔDIC : ΔSO42- ratio is 1.85 : 1 (Fig. S4), close to the theoretical 2 : 1 value expected for OSR. Minor deviations from the theoretical ratio may reflect partial DIC removal by authigenic carbonate precipitation (Akam et al., 2020). Chemoautotrophic CO₂ fixation may also contribute to DIC consumption (Treude et al., 2007). In addition, a fraction of sediment-produced DIC diffuses into the overlying water. Assuming oxidation of H₂S at the sediment–water interface, the estimated net diffusive benthic flux of carbonate alkalinity ([HCO3-] + 2[CO32-]) and DIC at ZGG is 983 ± 40 and 942 µmol m⁻² d⁻¹, respectively (Łukawska-Matuszewska and Kiełczewska, 2016; Łukawska-Matuszewska and Dwornik, 2025).

The ROSR values in ZGG sediments (24–46 µmol dm⁻³ d⁻¹) are consistent with rates reported from other Baltic regions using radioisotope methods. These include 0.2–42.3 in the Curonian and Vistula Lagoons (Pimenov et al., 2013), 0–43 in the Gdańsk Basin (Pimenov et al., 2010), and < 1.0–21.4 µmol dm⁻³ d⁻¹ in Kattegat (Iversen and Jørgensen, 1985). The measured values also fall within the lower range of rates reported from Eckernförde Bay (10–465 µmol dm⁻³ d⁻¹; Treude et al., 2005), Himmerfjärden (0.2–338 µmol dm⁻³ d⁻¹; Sawicka and Brüchert, 2017), and the Gotland Deep (< 250 µmol dm⁻³ d⁻¹; Piker et al., 1998). Incubation-derived potential ROSR values are several times higher than rates previously calculated for the same area from sulfate pore-water gradients (0.2–4.6 µmol dm⁻³ d⁻¹; Łukawska-Matuszewska and Dwornik, 2025). This discrepancy is expected, as incubation experiments determine gross sulfate reduction, whereas pore-water profiles reflect net sulfate consumption (Iversen and Jørgensen, 1985). Experiments with labeled substrates in the present study indicate the potential for AOM at ZGG (Table 1), with rates up to ∼ 10 µmol dm⁻³ d⁻¹. This suggests that methane could be oxidized if supplied from deeper layers. However, geochemical evidence indicates that the sampled interval corresponds to the sulfate reduction zone above the SMT. Consequently, AOM does not significantly contribute to DIC and TA production in the studied sediments. It should be also mentioned that experimentally-derived AOM rates represent the potential activity of the microbial community under substrate-replete conditions rather than the actual in situ process rates. Therefore, the similar experimentally-derived rates at MET-2 and ZGG indicate a comparable microbial potential for AOM at both sites. In contrast, the porewater profiles reflect in situ conditions, where methane availability differs substantially between the sites. At MET-2, methane is present and supports active AOM, which is reflected in the porewater profiles. At ZGG, methane concentrations are negligible, so AOM appears less significant in situ despite the presence of a microbial community capable of performing this process if methane became available.

Based on stoichiometry (Table S1: R5) and measured ROSR, OSR contributes 49–92 µmol dm⁻³ d⁻¹ of DIC to ZGG pore waters. The depth-integrated DIC production from OSR within the 0–100 cm interval equals 54 mmol m⁻² d⁻¹ (Table 1). Reduction of Mn(IV) during OM degradation may represent an additional DIC source, as indicated by the presence of Mn²⁺ in pore waters (Fig. 2g). However, Mn²⁺ concentrations are low relative to H₂S, the byproduct of sulfate reduction. Therefore, Mn(IV) reduction likely contributes only marginally to DIC production and is not included in the quantitative assessment. Besides Mn(IV) reduction, organic matter oxidation in marine sediments can proceed via denitrification. However, in the reducing sediments examined here nitrate was consistently below detection limit and nitrite occurred only at trace levels (data not shown), so denitrification was not quantified. While denitrification can produce DIC, we have no evidence for substantial availability of oxidized nitrogen species in these samples and therefore expect its contribution to be negligible.

The progressive increase in DIC with depth leads to pore-water supersaturation with respect to carbonate minerals, primarily calcite and dolomite, throughout most of the sediment column, except for the upper 5–10 cm (Fig. S11). Despite thermodynamic supersaturation, carbonate contents remain low (Table 2). Calcite occurs only in trace amounts, and dolomite was identified exclusively in the subsurface layer during summer. Moreover, anhedral morphology of many grains raises uncertainty as to whether all these carbonates are entirely authigenic. The dolomite burial rate is substantially lower than in methane-bearing sediments (Table 2) and does not exceed 271 µmol m⁻² d⁻¹. For comparison, burial rates reach 525 at MET2 and 994 µmol m⁻² d⁻¹ at MET1-MP. Comparable rates to those recorded at the ZGG – 60 ± 35 from PXRD and 136 ± 110 µmol m⁻² d⁻¹ from SEM–EDS – were reported by Zha et al. (2022) for sediments on the northern slope of the South China Sea. Using reactive transport modeling, they estimated the burial rate of Ca-Mg carbonates to be 65 in the upper 30 cm of sediment at one site and 79 µmol⁻² d⁻¹ in the upper 10 cm at another.

The DIC : TA ratio is a simple diagnostic of the relative contribution of non-carbonate bases to total alkalinity. In seawater, DIC defined as the sum of [CO2*] + [HCO3-] + [CO32-] with [CO2*] denoting the sum of dissolved [CO₂] and [H₂CO₃] (Dickson et al., 2007), typically represents over 95 % of TA (Emerson and Hedges, 2008). In this concept, lower DIC : TA values signal a larger noncarbonate alkalinity fraction (e.g., HS⁻, NH₃), whereas values near unity indicate TA is dominated by carbonate species. The mean DIC : TA ratio in ZGG pore waters is 0.75. This value indicates that approximately 25 % of total alkalinity derives from non-carbonate bases, primarily H₂S. Hydrogen sulfide constitutes an independent acid–base system capable of proton exchange (Hu et al., 2010) and is therefore included in TA calculations. Upward diffusion of H₂S toward the sediment–water interface followed by oxidation leads to partial consumption of alkalinity (Table S1: R10) generated during anaerobic mineralization (Krumins et al., 2013). In contrast, precipitation of sulfides within the sediment (Table S1: R15) prevents complete reoxidation and results in a net increase in pore-water TA.

Geochemical modeling indicates that pore waters at ZGG are supersaturated with respect to iron sulfides, including greigite, mackinawite, and pyrite, throughout the sediment column (Fig. S11). Increasing TS content with depth and relatively low TOC : TS ratios further support enhanced sulfide accumulation compared to the methane-affected stations. Mineralogical data confirm that pyrite is more abundant at ZGG than at MET1-MP and MET2. Its burial rate ranges from 940 to 1447 µmol m⁻² d⁻¹, several times higher than at the other sites, where OSR rates and pore-water H₂S concentrations are lower. These values are high and comparable to the upper range reported by Jørgensen et al. (1990) for sediments from the Baltic–North Sea transition zone (50–1240 µmol m⁻² d⁻¹).

The rate of sulfate reduction in the studied sediments likely exceeds the rate of reaction between H₂S and reactive iron phases. This imbalance promotes the accumulation of dissolved H₂S in pore waters (Boesen and Postma, 1988; Canfield et al., 1992). Accordingly, H₂S concentrations remain high throughout the sediment profile, whereas Fe²⁺ is largely depleted from pore waters (Fig. 2). Mössbauer spectroscopy indicates that iron is predominantly divalent (Fig. S10). Thus, the availability of reactive Fe(III) phases likely limits pyrite formation. Such limitation is typical of anoxic to euxinic marine environments (Raiswell and Berner, 1985; Boesen and Postma, 1988) and has previously been documented for sediments of the Gdańsk Deep (Łukawska-Matuszewska et al., 2019). Incomplete removal of H₂S as solid sulfides and its reoxidation (Table S1: R10) may contribute to pore-water acidification. The observed increase in Ca²⁺ concentration with depth (Fig. 2) may therefore reflect carbonate dissolution driven by lowered pH.

High OM loading from local primary production and terrigenous sources increases oxygen demand in coastal sediments and frequently exceeds oxygen supply (Hedges et al., 1997; Schlünz and Schneider, 2000). Oxygen penetration is therefore shallow. It is typically < 1 mm in muddy and silty sediments and only a few centimeters in sands (Glud, 2008), owing to limited transport and rapid consumption. As a consequence, anaerobic mineralization dominates in most coastal settings (Soetaert et al., 1996), with sulfate serving as the principal electron acceptor, as observed at station ZGG. When organic carbon input exceeds sulfate availability, mineralization proceeds via methanogenesis (Froelich et al., 1979). Diagenetic zonation in coastal sediments differs markedly from that in open-ocean environments. Redox zones are thinner and shifted closer to the sediment surface (D'Hondt et al., 2004; Egger et al., 2018). In low-salinity settings, sulfate reduction, the SMT, and methanogenesis commonly occur at shallow depths (Albert et al., 1998; Thang et al., 2013; Egger et al., 2015). Station MET2 is located in a nearshore zone (Fig. S1) with bottom-water salinity of ∼ 7–8. The molar TOC : TN ratio (9.3–13.7; Fig. 1d) indicates a substantial contribution of terrestrial OM, which is relatively nitrogen-poor compared to protein-rich marine phytoplankton. Carbon isotopic compositions support this interpretation. δ13C values range from -26.0 ‰ to -24.4 ‰, i.e., below -24 ‰, consistent with terrestrial input (Yadav et al., 2025). Nitrogen isotopic values (δ15N = 1.21 ‰–6.66 ‰) show the greatest ¹⁵N enrichment in the upper 0–20 cm, coinciding with elevated TOC and TN contents. Terrestrial OM is generally less degradable than marine OM due to the abundance of cellulose and lignin in vascular plant tissues (Dickens et al., 2006). The observed isotopic and compositional variations therefore reflect preferential degradation of labile phytoplankton-derived compounds and the relative accumulation of more refractory terrestrial material in the sediment.

Because MET2 is located in a shallower, transport-dominated setting, the sediment consists of sandy silt and contains lower TOC and TN than the silty and clayey sediments at the other stations (Fig. 1a, b). Owing to the relatively low salinity, major ion concentrations in bottom and pore waters are the lowest among all sites. These concentrations remain nearly constant with depth, except for SO42- (Fig. 2a–d). Sulfate decreases rapidly within the upper 0–40 cm due to microbial reduction and falls below LOQ at greater depths. In the upper ∼ 20 cm b.s.f., the ΔDIC : ΔSO42- ratio ranges from 1.0 to 2.1, consistent with contributions from both sulfate-driven AOM and OSR. The mean ratio of 1.51 : 1 indicates a mixed origin of DIC. Below this interval, DIC production becomes disproportionately high relative to sulfate consumption, and the ΔDIC : ΔSO42- ratio increases to ∼ 4.0 (Fig. S4). Such excess DIC may reflect processes that produce DIC but do not consume sulfate, including methanogenesis, or alternative AOM pathways coupled to Mn(IV) or Fe(III) reduction (Meister et al., 2019). Elevated Mn²⁺ concentrations between 30 and 80 cm and Fe²⁺ below 80 cm (Fig. 2g, h) may suggest involvement of these pathways. However, the absolute concentrations of Mn²⁺ and Fe²⁺ are low, indicating that metal-driven AOM likely contributes only marginally to DIC and TA production and is not considered further.

Methane concentrations at MET2 reach 6.5 in summer and 8.45 mmol dm⁻³ in winter (Fig. 2f). Stable isotope data for δ13C-CH₄, δ2H-CH₄, and δ13C-CO₂ (Fig. 5) indicate a microbiological origin, predominantly via hydrogenotrophic methanogenesis (Whiticar, 1999). Methane concentrations decrease upward from maxima in deeper layers to values below LOQ in the uppermost centimeters. Methane is absent from the overlying water column, demonstrating that the sediment acts as an efficient methane filter despite the shallow SMT at ∼ 5–13 cm b.s.f. In Baltic Sea sediments, SMT depths vary widely, ranging from 80–160 cm in the Danish Straits (Iversen and Jørgensen, 1985), ∼ 150 cm in Aarhus Bay (Thomsen et al., 2001), and 30–40 cm in the Gotland Deep (Piker et al., 1998). SMT depths comparable to those observed here have been reported from Storfjärden and Pojo Bay (Myllykangas et al., 2020; Jilbert et al., 2021), the western Gotland Basin (Ketzer et al., 2024), and Eckernförde Bay in the southwestern Baltic Sea (Treude et al., 2005; Maltby et al., 2018).

The occurrence of sulfate-dependent AOM within the 5–13 cm b.s.f. interval is indicated by the most negative δ13C-DIC values (Fig. 3c) (Malinverno and Pohlman, 2011). In the overlying sediment, OM is oxidized primarily via OSR. However, separation of OSR and AOM zones based solely on isotopic data is difficult, as both processes produce ¹³C-depleted DIC. OSR shifts δ13C-DIC from values near 0 ‰, typical of seawater, toward approximately -21 ‰, characteristic of marine OM. Subsequent sulfate-driven AOM may further deplete DIC in ¹³C through oxidation of isotopically light CH₄ (Meister and Reyes, 2019). At MET2, δ13C-DIC values range from -12.7 ‰ to -1.2 ‰. These values are markedly less negative than expected for only AOM or OSR end-members. If DIC generated by sulfate-AOM inherited the isotopic signature of methane entering the SMT, δ13C-DIC values would approach ∼-76 ‰. Conversely, DIC derived from OSR should approximate the δ13C of sedimentary OM (-25.6 ‰ at MET2). The relatively enriched δ13C-DIC values observed in both AOM and OSR intervals suggest mixing with ¹³C-enriched DIC diffusing upward from the methanogenic zone (Malinverno and Pohlman, 2011; Meister and Reyes, 2019). Overall, pore-water DIC at MET2 reflects a mixture of inorganic carbon produced by overlapping diagenetic processes. A comparable situation has been reported from Storfjärden, where OSR, sulfate-driven AOM, and methanogenesis co-occur within compressed redox zones (Jilbert et al., 2021).

The ROSR values at MET2 range from 21 to 48 µmol dm⁻³ d⁻¹. These rates are comparable to those measured at the other stations (Table 2) and fall within ranges reported in the literature (Sect. 4.1). The RAOM values (6–33 µmol dm⁻³ d⁻¹) are also consistent with data from other Baltic regions, including 0–50 in the Gdańsk Basin (Pimenov et al., 2010), 0.03–14 in Eckernförde Bay (Treude et al., 2005), and 0.002–6.2 µmol dm⁻³ d⁻¹ in Kattegat and Skagerrak (Iversen and Jørgensen, 1985). The highest RAOM occurs within the 0–20 cm layer. Rates in this interval are up to six times higher than those in the 20–40 cm layer. This distribution corresponds to geochemical indicators, including the steepest SO42- gradient, maximum H₂S concentrations, ΔDIC : ΔSO42- ratios lower than expected exclusively for OSR, and negative δ13C-DIC values in the upper sediment. At MET2, ROSR contributes 41–96 µmol DIC dm⁻³ d⁻¹ to pore waters. The depth-integrated DIC production from OSR within the 0–100 cm interval equals 59 mmol m⁻² d⁻¹, similar to values obtained for ZGG (Table 1). In addition, AOM contributes 6–33 µmol dm⁻³ d⁻¹ of DIC, corresponding to an integrated production rate of 9 mmol m⁻² d⁻¹. The pronounced chemical gradients at the sediment–water interface correspond to a diffusive benthic flux of carbonate alkalinity ([HCO3-] + 2[CO32-]) estimated at 1498 ± 80 µmol m⁻² d⁻¹ (Łukawska-Matuszewska and Dwornik, 2025). The mean DIC : TA ratio in MET2 pore waters is 0.97. This value indicates a relatively minor contribution of non-carbonate bases to total alkalinity compared to ZGG. The dominance of carbonate alkalinity suggests that alkalinity production at MET2 is less reversible. In contrast, at ZGG, substantial H₂S production and its subsequent oxidation release H⁺ and partially consume alkalinity (Table S1: R10).

Pore waters at MET2 exhibit higher saturation indices with respect to calcite and dolomite than those at methane-free ZGG (Fig. S11). Intense diagenesis leads to a strong increase in DIC, which raises porewater saturation with respect to carbonates and thereby promotes precipitation of authigenic phases. This is consistent with the greater abundance of dolomite throughout the sediment column, as documented by PXRD and microscopic observations (Table 2). The dolomite burial rate at MET2 (184–525 µmol m⁻² d⁻¹) substantially exceeds that at ZGG (1–271 µmol m⁻² d⁻¹). Taken together, the porewater chemistry and burial rates suggest that methane-bearing sediments in our study area sequester more inorganic carbon than methane-free, OSR-dominated sediments (ZGG) – likely because methane-related processes generate additional DIC.

In deeper layers at MET2, pore waters are also supersaturated with respect to siderite. In contrast, siderite is undersaturated throughout the ZGG profile (Fig. S11). This difference likely reflects Fe²⁺ availability. At ZGG, Fe²⁺ is efficiently removed from pore waters through reaction with abundant H₂S and subsequent sulfide precipitation. At MET2, H₂S is restricted to a narrow interval, with maximum concentrations at ∼ 13 cm b.s.f. (Fig. 2e). Below this depth, H₂S concentrations fall below LOQ. Under these conditions, a fraction of dissolved Fe²⁺ may precipitate as siderite rather than as iron sulfides.

In MET2 sediments, the mean TOC : TS ratio is 2.8. This value corresponds to typical marine sediments deposited beneath oxygenated bottom waters (Berner, 1982). Under such conditions, H₂S required for pyrite formation is produced within the anaerobic sediment during OM degradation. The rate of FeS₂ formation is controlled primarily by sedimentation rate and organic matter supply (Berner, 1982).

In contrast, at ZGG, H₂S frequently occurs in the overlying water column. This allows pyrite to form not only within the sediment but also in the water column (Raiswell and Berner, 1985) and, as a result, pyrite is more abundant at ZGG than at MET2 (Table 2). At MET2, the pyrite burial rate ranges from 334 to 789 µmol m⁻² d⁻¹, the lowest among the studied sites. This reflects oxygenated bottom waters and comparatively low H₂S concentrations in pore waters. These conditions likely explain the Mössbauer results, which show a slightly lower proportion of Fe(II) associated with pyrite in the upper sediment layers at MET2 compared to deeper layers and to the other stations (Table S4).

Sedimentation rate and organic carbon burial strongly regulate sulfate consumption and the depth of the SMT in marine sediments (Canfield, 1991; Egger et al., 2018). In coastal settings, submarine groundwater discharge (SGD) constitutes an additional controlling factor. SGD involves the exchange of terrestrial groundwater mixed with recirculated seawater that has infiltrated coastal aquifers. Its chemical composition often deviates from conservative mixing due to biogeochemical reactions within the aquifer (Moore, 2010). In the Gdańsk Basin, SGD occurs in both coastal and deeper offshore areas. It is expressed by anomalous salinity distributions, with lower salinity in bottom and pore waters relative to the overlying water column (Piekarek-Jankowska et al., 1994; Falkowska and Piekarek-Jankowska, 1999; Ehlert von Ahn et al., 2024; Matciak et al., 2024). In nearshore zones, SGD significantly affects water quality by supplying carbon, nutrients, and contaminants (Szymczycha et al., 2012, 2013, 2023). However, station MET1-MP, located in central Gulf of Gdańsk ∼ 80 m water depth, represents a deepwater pockmark where freshwater originates from the Upper Cretaceous aquifer beneath the seafloor (Kozerski et al., 1987). Fresh groundwater infiltration into the pockmark has been documented previously (Brodecka-Goluch et al., 2022) and is confirmed here by the Cl⁻ profile in pore waters (Fig. 2a). The inflow of freshwater dilutes pore waters and reduces SO42- availability for OM and methane oxidation. This alteration of sulfate supply modifies early diagenetic processes and promotes shallow methanogenesis (Pimenov et al., 2010; Ehlert von Ahn et al., 2024).

As a seabed depression formed by gas expulsion, this pockmark acts as a natural sediment trap. It accumulates both autochthonous and allochthonous OM. Consequently, silty sediments at MET1-MP exhibit the highest TOC and TN concentrations among the studied sites (Fig. 1a, b). The TOC : TN ratio (7.8–8.7) approaches that typical of marine OM. However, the isotopic composition of carbon and nitrogen in the organic fraction (Fig. 1e, f) indicates an additional contribution of terrestrial material.

The pockmark's physical and sedimentary characteristics concentrate OM, raising demand for electron acceptors and fueling intense diagenesis; as a result, porewater sulfate is rapidly consumed in the uppermost decimeters and the sulfate–methane transition shifts toward the sediment surface. Methanogenesis and AOM supply additional DIC alongside OSR, causing elevated porewater DIC through intense diagenesis. Pore waters at MET1-MP are also characterized by very high concentrations of NH4+ and PO43- (Fig. 3d, e). These elevated nutrient levels further indicate intense degradation of the substantial OM inventory accumulated within the pockmark sediments.

The substantial water depth and persistent stratification at MET1-MP further constrain oxygen supply to the sediment, similar to conditions at ZGG. Under these conditions, SO42- is rapidly depleted from pore waters. The SMT is extremely shallow. The δ13C-DIC profile (Fig. 3c) indicates that the SMT center is located immediately below the sediment surface, likely within 3–5 cm b.s.f. Within the upper ∼ 8 cm, ΔDIC : ΔSO42- ratios reflect contributions from OSR, AOM, or both processes. The mean ratio of 1.24 : 1 suggests that sulfate-driven AOM plays a dominant role in SO42- consumption in this interval (Miller et al., 2017). At greater depths, ΔDIC values exceed those expected from OSR and AOM, indicating an additional DIC source unrelated to SO42- consumption. With increasing depth, δ13C-DIC shifts toward positive values, reaching up to +16.4 ‰ (Fig. 3c). Similar values (up to +14 ‰) and comparable pore-water profiles were reported by Ehlert von Ahn et al. (2024) from a station located in similar area. This enrichment likely reflects the upward diffusion of ¹³C-enriched DIC generated during methanogenesis, as observed at MET2. However, groundwater infiltration may also contribute.

Groundwaters from the Upper Cretaceous aquifer in the Gdańsk region are predominantly of the HCO₃–Na and HCO₃–Ca types (Sadurski, 1985). They are alkaline (pH 7.5–8.3) and exhibit variable dissolved solute concentrations, with DIC values reaching ∼ 11 mmol dm⁻³ (Kozerski et al., 1987). These waters circulate through glauconitic sands forming the principal aquifer and migrate upward through carbonate and siliceous strata, including carbonates, gaizes, and marls (Sadurski, 1985). Infiltrating groundwater can dilute pore waters and reduce DIC concentrations in deeper sediment layers (Fig. 3b). Moreover, groundwater DIC is commonly enriched in ¹³C due to carbonate dissolution and isotopic exchange with CaCO₃, for which δ13C values typically exceed 0 ‰ (Walter et al., 2007; Han et al., 2014; Campeau et al., 2017). Consequently, both methanogenesis and SGD influence the δ13C-DIC signature in pore waters of the MET1-MP pockmark.

Stable isotope data indicate that CH₄ at MET1-MP is of microbiological origin (Fig. 5). Methane is generated not only by CO₂ reduction but also by acetate fermentation, a pathway commonly associated with freshwater environments (Whiticar, 1999; Batther et al., 2025). This dual signature further demonstrates the influence of freshwater infiltration on sedimentary microbial processes. In contrast to MET2, where CH₄ concentrations increase progressively with depth, the MET1-MP profile exhibits a subsurface or intermediate maximum (Fig. 2f). The extremely shallow SMT and the occurrence of methane close to the sediment surface may enhance the potential for gas escape into the water column. Lapham et al. (2024) showed that methane leakage across the SMT is widespread in organic-rich Baltic Sea sediments, with a highly variable amount of subseafloor gas escaping oxidation (0 %–100 %). Previous investigations documented gas emissions from multiple locations within the MET1-MP pockmark (Brodecka et al., 2013; Majewski and Klusek, 2014; Łukawska-Matuszewska et al., 2025). The pockmark is a large structure (∼ 1200 m long and ∼ 500 m wide), and gas discharge is spatially heterogeneous. Fluid migration pathways and rates vary considerably. During the present study, CH₄ concentrations in bottom water were below LOQ, whereas surface sediments contained ∼ 0.3 mmol dm⁻³ CH₄. Thus, diffusion into bottom water followed by aerobic oxidation cannot be excluded.

The ROSR attains a maximum of 164 µmol dm⁻³ d⁻¹ in the 0–20 cm layer. This peak corresponds to steep SO42- depletion and elevated H₂S concentrations, reflecting active sulfate reduction. Below 20 cm, ROSR decreases markedly and remains relatively uniform across deeper layers (Table 2). Excluding the surface interval, ROSR is lower than at the other stations. In contrast, RAOM is substantially higher, ranging from 59 to 1077 µmol dm⁻³ d⁻¹. These RAOM values exceed most rates reported for the Baltic Sea. Rates measured in the 0–20, 20–40, and 60–80 cm layers are comparable to those reported by Pimenov et al. (2008), who observed a maximum of 80.6 µmol dm⁻³ d⁻¹ in another gas-saturated pockmark in the Gdańsk Basin, among the highest previously documented in the region. By comparison, AOM rates in Eckernförde Bay range from 0.03 to 14 µmol dm⁻³ d⁻¹ (Treude et al., 2005), in the Curonian and Vistula Lagoons from 0.07 to 2.89 µmol dm⁻³ d⁻¹ (Pimenov et al., 2013), and in Kattegat and Skagerrak from 0.002 to 6.2 µmol dm⁻³ d⁻¹ (Iversen and Jørgensen, 1985). The highest RAOM values at MET1-MP approach those reported for gassy sediments at a cold seep in Monterey Bay (280 ± 40 µmol dm⁻³ d⁻¹) and fall within the lower range measured at a cold seep in the Gulf of Mexico, where rates reached 4800 ± 400 µmol dm⁻³ d⁻¹ (Bowles et al., 2019).

Direct comparison of microbial rates among studies remains challenging: even when identical radiotracer protocols are used, variability in microbial community structure, sediment properties and incubation conditions can produce substantially different rate estimates. At MET1-MP, RAOM varies by up to two orders of magnitude between sediment layers. Comparable intra-profile variability has been documented elsewhere in the Baltic Sea (e.g., Iversen and Jørgensen, 1985; Pimenov et al., 2008; Pimenov et al., 2013), where rates differed by as much as 300 % among layers. Nevertheless, the MET1-MP pockmark exhibits a markedly greater AOM capacity than the other sites investigated in this study and than most previously studied Baltic sediments. Depth-integrated calculations demonstrate that AOM is the dominant DIC source at MET1-MP. OSR contributes 58 mmol m⁻² d⁻¹ of DIC, whereas AOM supplies 331 mmol m⁻² d⁻¹, more than five times greater. Consistently, MET1-MP represents the strongest benthic DIC source among the studied sites. The estimated net diffusive flux of carbonate alkalinity equals 2017 ± 47 µmol m⁻² d⁻¹ (Łukawska-Matuszewska and Dwornik, 2025), exceeding fluxes from the other sediments by 50 %–100 %.

Intense DIC production within the MET1-MP pockmark results in markedly elevated DIC and TA concentrations relative to the other stations (Fig. 3a, b). In contrast to MET2 and ZGG, however, DIC does not increase progressively with depth. Below ∼ 30 cm, DIC concentrations decline. This pattern may partly reflect dilution by infiltrating freshwater. Because the exact chemical composition of the discharging waters is not fully constrained, their quantitative influence on DIC variability cannot be assessed. Notably, the decrease in DIC coincides with pronounced declines in Mg²⁺ and Ca²⁺ concentrations (Fig. 2b, c). This relationship suggests that authigenic carbonate precipitation represents an additional controlling process. Comparable pore-water trends have been reported from methane-affected settings, including shallow sediments of the western slope of the Mid-Okinawa Trough (Xu et al., 2018), pockmark sediments of the Northern Congo Fan (Nöthen and Kasten, 2011), and pockmarks on the inner shelf of the East China Sea (Wang et al., 2024).

Among the studied sites, MET1-MP affected by methanogenesis and freshwater infiltration exhibits the highest average dolomite content. Euhedral crystals are common, supporting an authigenic origin. In addition to dolomite and minor calcite, traces of siderite were identified. Authigenic carbonate formation at MET1-MP therefore constitutes an efficient inorganic carbon sink. This observation is consistent with Rzepa et al. (2026), who found that freshwater input in another deepwater pockmark in the Gdańsk Basin lowers porewater sulfate and promotes methanogenesis and Fe-mediated AOM, increasing dissolved inorganic carbon and thereby favoring authigenic carbonate formation. The dolomite burial rate ranges from 467 to 994 µmol m⁻² d⁻¹, substantially exceeding values at ZGG (1–271 µmol m⁻² d⁻¹) and MET2 (184–525 µmol m⁻² d⁻¹). Similarly, Chen et al. (2024) found that at the Makran continental margin, sediments with a high methane flux exhibited authigenic carbonate precipitation rates that were at least double those of the background area.

In the MET1-MP pockmark influenced by freshwater infiltration, the effect of CH₄ on redox-sensitive elements (Fe and Mn) is more pronounced than at MET2. Freshening and reduced sulfate availability enhance the relative importance of alternative methane oxidation pathways (Egger et al., 2015). The observed increase in Fe²⁺ and Mn²⁺ concentrations below the sulfate reduction zone likely reflects reductive dissolution of Fe(III) and Mn(IV) phases coupled to AOM. The measured RAOM values probably include contributions from both sulfate-dependent AOM and metal-dependent AOM, provided that microorganisms capable of these pathways were present. This distinction does not affect DIC production estimates, as oxidation of 1 mole of CH₄ yields 1 mol of HCO3- irrespective of the terminal electron acceptor (Table S1: R12–R14). MET1-MP exhibits the highest Fe²⁺ concentrations among the studied stations. This contrasts with ZGG, where Fe²⁺ is almost entirely removed from pore waters through reaction with abundant H₂S. The total sulfur content in MET1-MP sediments is lower than at ZGG, and the TOC : TS ratio is several times higher (Fig. 1). These relations indicate that H₂S availability limits pyrite formation in sediments influenced by freshwater infiltration. Despite this limitation, pyrite remains a common authigenic phase at MET1-MP. However, its mean concentration is significantly lower than at ZGG. The pyrite burial rate ranges from 453 to 828 µmol m⁻² d⁻¹, comparable to values at MET2. These rates are approximately half those observed at ZGG, where intense sulfate reduction promotes high H₂S production and enhanced iron sulfide precipitation.