Sulfur Cycling in an Iron Oxide-Dominated, Dynamic Marine Depositional System: The Argentine Continental Margin

The interplay between sediment deposition patterns, organic matter type and the quantity and quality of reactive mineral phases determines the accumulation, speciation and isotope composition of pore water and solid phase sulfur constituents in marine sediments. Here, we present the sulfur geochemistry of siliciclastic sediments from two sites along the Argentine continental slope—a system characterized by dynamic deposition and reworking, which result in non-steady state conditions. The two investigated sites have different depositional histories but have in common that reactive iron phases are abundant and that organic matter is refractory—conditions that result in low organoclastic sulfate reduction rates. Deposition of reworked, isotopically light pyrite and sulfurized organic matter appear to be important contributors to the sulfur inventory, with only minor addition of pyrite from organoclastic sulfate reduction above the sulfate-methane transition (SMT). Pore-water sulfide is limited to a narrow zone at the SMT. The core of that zone is dominated by pyrite accumulation. Iron monosulfide and elemental sulfur accumulate above and below this zone. Iron monosulfide precipitation is driven by the reaction of low amounts of hydrogen sulfide with ferrous iron and is in competition with the oxidation of sulfide by iron (oxyhydr)oxides to form elemental sulfur. The intervals marked by precipitation of intermediate sulfur phases at the margin of the zone with free sulfide are bordered by two distinct peaks in total organic sulfur. Organic matter sulfurization appears to precede pyrite formation in the iron-dominated margins of the sulfide zone, potentially linked to the presence of polysulfides formed by reaction between dissolved sulfide and elemental sulfur. Thus, SMTs can be hotspots for organic matter sulfurization in sulfide-limited, reactive iron-rich marine sedimentary systems. Furthermore, existence of elemental sulfur and iron monosulfide phases meters below the SMT demonstrates that in sulfide-limited systems metastable sulfur constituents are not readily converted to pyrite but can be buried to deeper sediment depths. Our data show that in non-steady state systems, redox zones do not occur in sequence but can reappear or proceed in inverse sequence throughout the sediment column, causing similar mineral alteration processes to occur at the same time at different sediment depths.

The interplay between sediment deposition patterns, organic matter type and the quantity and quality of reactive mineral phases determines the accumulation, speciation, and isotope composition of pore water and solid phase sulfur constituents in marine sediments. Here, we present the sulfur geochemistry of siliciclastic sediments from two sites along the Argentine continental slope-a system characterized by dynamic deposition and reworking, which result in non-steady state conditions. The two investigated sites have different depositional histories but have in common that reactive iron phases are abundant and that organic matter is refractory-conditions that result in low organoclastic sulfate reduction rates (SRR). Deposition of reworked, isotopically light pyrite and sulfurized organic matter appear to be important contributors to the sulfur inventory, with only minor addition of pyrite from organoclastic sulfate reduction above the sulfate-methane transition (SMT). Pore-water sulfide is limited to a narrow zone at the SMT. The core of that zone is dominated by pyrite accumulation. Iron monosulfide and elemental sulfur accumulate above and below this zone. Iron monosulfide precipitation is driven by the reaction of low amounts of hydrogen sulfide with ferrous iron and is in competition with the oxidation of sulfide by iron (oxyhydr)oxides to form elemental sulfur. The intervals marked by precipitation of intermediate sulfur phases at the margin of the zone with free sulfide are bordered by two distinct peaks in total organic sulfur (TOS).
Organic matter sulfurization appears to precede pyrite formation in the iron-dominated margins of the sulfide zone, potentially linked to the presence of polysulfides formed by reaction between dissolved sulfide and elemental sulfur. Thus, SMTs can be hotspots for organic matter sulfurization in sulfide-limited, reactive iron-rich marine sedimentary systems. Furthermore, existence of elemental sulfur and iron monosulfide phases meters below the SMT demonstrates that in sulfide-limited systems metastable sulfur constituents are not readily converted to pyrite but can be buried to deeper INTRODUCTION Dynamic depositional systems strongly impact sedimentary geochemical processes; sediments formed under such geochemical regimes do not always show the typical sequence of terminal electron acceptor processes predicted for steady state environments (e.g., Claypool and Kaplan, 1974;Froelich et al., 1979;Berner, 1981). Depending on the availability and reactivity of the solid-phase electron acceptors, the redox sequence can even be reversed in some cases. For example, iron and manganese reduction can reappear or persist at depth below the zone dominated by sulfate reduction (e.g., Postma and Jakobsen, 1996;Kasten et al., 1998;Hensen et al., 2003;Riedinger et al., 2014;Treude et al., 2014). In particular, non-steady state depositional conditions can have a strong impact on the inventory of sulfide minerals in anoxic marine environments (e.g., Kasten et al., 1998Kasten et al., , 2003Aller et al., 2010;Borowski et al., 2013;Aller, 2014;Peketi et al., 2015). Iron sulfide-containing sediments that are eroded and/or transported down-slope in oxygenated seawater can be subject to reoxidation, leading to the conversion of iron monosulfides and fine-grained pyrite to (amorphous) ferric hydroxides (Luther et al., 1982;Morse, 1991). Shielded from sulfidic conditions in the upper sediment column due to rapid burial-and in the presence of mostly reworked, unreactive organic matter-those oxidized reactive iron phases are then preserved in deeper subsurface sediments März et al., 2008;Riedinger et al., 2014). The continental margin off Uruguay and Argentina is characterized by such highly dynamic depositional conditions (e.g., Riedinger et al., 2005;Henkel et al., 2011Henkel et al., , 2012Krastel et al., , 2013, and these locations are likely to be representative of environments that are common throughout the world along continental margins. Sediments along the continental margin off Uruguay and Argentina, at the western rim of the Argentine Basin, are not only characterized by redistribution/reworking, including recycling of organic matter, but also by an input of iron (as iron (oxyhydr)oxides). Those inputs of reactive iron greatly outpace sedimentary sulfide production, which results in an iron-dominated system (Haese et al., 2000;Hensen et al., 2003;Riedinger et al., 2005). This situation arises because residual, less reactive organic matter, which has already been degraded at the shallower sites of initial deposition, is reworked and transported and further remineralized along with sulfide minerals, which are oxidized to yield high amounts of Fe(III). Hensen et al. (2003) concluded that the non-steady state processes and associated iron oxidation in this dynamic system favors the retention of any reduced sulfur generated following redeposition, which was corroborated by the quantification of sulfur burial based on a transport and reaction model  and subsequent direct measurements of sulfur constituents . Although these studies demonstrated that sedimentary sulfur cycling in an iron oxide-dominated, dynamic marine system strongly affects sulfur sequestration and subseafloor microbial processes, some questions remain unresolved; in particular, it is unclear why a wide variety of different reduced sulfur compounds is formed and if there is a "classical" zonation-analogous to the standard redox towerfor diagenetic sulfur transformations. To study the fate of various sulfur species, their contents and isotopic signatures, as well as the iron inventory, in more detail in these dynamic depositional systems, we collected samples via gravity cores from the lower slope in the western Argentine Basin during the RV Meteor Expedition M78/3 in May-July 2009.

Sulfur Cycling in Marine Sediments
The main driver of the marine sedimentary sulfur cycle is the microbial reduction of sulfate (e.g., Goldhaber and Kaplan, 1974;Froelich et al., 1979;Jørgensen, 1982;Bowles et al., 2014). The two major catabolic microbial sulfate reduction pathways are organoclastic sulfate reduction and sulfate reduction coupled to anaerobic oxidation of methane (AOM); both processes release hydrogen sulfide to the pore water (e.g., Goldhaber and Kaplan, 1974;Jørgensen, 1982;Hoehler et al., 1994;Boetius et al., 2000). Microbial sulfate reduction discriminates against the heavier sulfur isotopes, leading to a relative depletion in 34 S in the produced hydrogen sulfide and a corresponding relative 34 S enrichment in the remaining sulfate pool (Jones and Starkey, 1957;Harrison and Thode, 1958;Thode et al., 1961;Kaplan and Rittenberg, 1964;Rees, 1970;Kaplan, 1974, 1980;Froelich et al., 1979;Bottrell and Raiswell, 2000;Brunner and Bernasconi, 2005;Sim et al., 2011;Wing and Halevy, 2014). Typically, with increasing sediment depth, as the sulfate pool becomes smaller, the remaining sulfate becomes isotopically heavier (e.g., Goldhaber and Kaplan, 1974;Torres et al., 1996). This trend is reflected in the production of hydrogen sulfide that also becomes isotopically heavier. As a consequence, with increasing sediment depth, any sulfur phase resulting from sulfide oxidation or directly precipitated iron sulfide will also become enriched in 34 S (e.g., Goldhaber and Kaplan, 1980). Furthermore, at the sulfate-methane transition (SMT) where sulfate is almost entirely consumed, iron sulfides are precipitated with the heaviest isotope composition compared to the upper sediments Borowski et al., 2013;Lin et al., 2016).

Sedimentary Setting
The sedimentary environment of the Argentine Basin, including the continental margin off Argentina and Uruguay (Figure 1), is controlled by dynamic depositional processes, such as gravity-controlled sediment transport and strong current circulation (Ewing and Leonardi, 1971;Klaus and Ledbetter, 1988;Hernández-Molina et al., 2009;Preu et al., 2013). In the upper waters, the southward flowing Brazil Current and the northward flowing Malvinas (Falkland) Current meet in front of the Rio de la Plata (Peterson and Stramma, 1991). The Brazil Malvinas Confluence (BMC) leads to an increase in primary production over a distinct area related to the mixing of these tropical and Antarctic water masses, which results in strong gradients in nutrient, salinity, and temperature (Antoine et al., 1996;Behrenfeld and Falkowski, 1997;Chiessi et al., 2007) and elevated organic carbon input into the sediment along the shelf and upper slope. Southward-flowing North Atlantic Deep Water (NADW, 2,000 and 4,000 m) and northward-flowing Antarctic Bottom Water (AABW >4,000 m) are parallel to the continental margin; below 4,000 m water depth, the AABW transports mostly fine-grained sediment.
Predominantly terrigenous material, delivered from the numerous fluvial tributaries along the coast of Argentina and Uruguay (Iriondo, 1984;Piccolo and Perillo, 1999), is transported downslope from the continental shelf via gravity-controlled mass flows, including turbidity currents and density flows (Biscaye and Dasch, 1971;Klaus and Ledbetter, 1988;Sachs and Ellwood, 1988;Romero and Hensen, 2002;Hensen et al., 2003;Henkel et al., 2011Henkel et al., , 2012Gruetzner et al., 2012;Voigt et al., 2013). These mass flows also transport reworked organic matter further downslope, resulting in burial of refractory organic carbon (Hedges and Keil, 1995) in the deeper parts of the basin. Thus, high mean sedimentation rates, including mass transport deposits (MTD), combined with low levels of reactive organic carbon, leads to rapid burial of (highly) reactive iron minerals present as both primary terrigenous phases and reworked and oxidized authigenic components Riedinger et al., 2005Riedinger et al., , 2014.

Sampling and Sample Processing
Gravity cores were collected east of the Rio de la Plata mouth from water depths of 3821 m and 3687 m (GeoB 13824-1 38 • 13.14 ′ S, 53 • 21.29 ′ W; GeoB 13863-1 39 • 18.70 ′ S, 53 • 57.16 ′ W; Figure 1) during R/V Meteor expedition M78/3 . The cores were taken at two sites from different depositional settings (Figure 2). GeoB 13824 was collected at the foot of the Mar del Plata Canyon , which is characterized by high accumulation of sediments that bypassed the shelf through the canyon, as well as potential erosion events. Core GeoB 13863 was recovered from a site on the lower slope south of the Rio de la Plata where recent sediment input is mainly influenced by currents and minor amounts of mass gravity flow. The site experienced profound variations in depositional conditions over glacial/interglacial timescalesmost likely related to changes in sea level .
After retrieval, the cores were cut immediately into 1-m segments on deck. For methane analysis, 3 cc syringe samples were taken from every core segment and transferred into 20 mL headspace vials pre-filled with 10 mL of a 5 M NaCl solution and stored at 4 • C. Additionally, samples for determinations of sulfate reduction rates (SRR) were taken from every core segment (for Hole GeoB 13824 only). Pore water samples were extracted on segment-halves in a cold room (∼4 • C) via the Rhizon method (Seeberg-Elverfeldt et al., 2005;Dickens et al., 2007). For sulfate and hydrogen sulfide concentrations and isotope analyses, 5 mL subsamples of pore water were added to a 2.5% zinc acetate (ZnAc) solution in order to fix all sulfide present as zinc sulfide (ZnS). The pH and Eh were measured using punch-in electrodes. Solid phase samples were taken at 20-30 cm intervals and placed and sealed under a nitrogen atmosphere in aluminum bags and stored frozen (−20 • C) until processed onshore. The geochemical data reported here can be accessed via the information system PANGAEA operated by the World Data Centers for Marine Environmental Sciences (https://doi.pangaea. de/10.1594/PANGAEA.856812).

Pore Water Analyses
Total dissolved sulfide concentrations ( H 2 S = H 2 S + HS − + S 2− ) were analyzed onboard the ship spectrophotometrically using the methylene blue method (Cline, 1969). GeoB 13824-1 samples for sulfate (SO 2− 4 ) concentrations were also analyzed onboard at 1:50 and 1:100 dilutions using a high performance liquid chromatography (HPLC) Sykam solvent delivery system coupled to a Waters 430 conductivity detector. Sulfate concentrations for Site GeoB 13863 were measured at the Max Planck Institute for Marine Microbiology (MPI-MM) in Bremen, Germany, via suppressed ion chromatography at a 1:100 dilution with double de-ionized water on a Metrohm 761 compact IC. Standard calibrations were performed using seawater provided by the International Association for the Physical Sciences of the Oceans (IAPSO) and in-house standards. The error of replicate analyses of sulfide and sulfate was <2 and 3%, respectively. Methane (CH 4 ) was measured with a Hewlett Packard 5890A gas chromatograph using a splitless injector, a stainless steel Porapak-Q column and a flame ionization detector at the MPI-MM. Chromatographic response on the GC instrument was calibrated against three different standards with variable concentrations of CH 4 . The measured concentrations were corrected for sediment porosity.
Potential polysulfide concentrations ( S 2− n -Calc., n = 2,...,8) were calculated using the thermodynamic constants of Kamyshny et al. (2007) with the program MinteQ R , assuming an activity of 1 for solid elemental sulfur and using the measured sulfide concentration. Values were calculated for seawater at FIGURE 2 | Multibeam bathymetric map of the study site area and seismic profiles across the investigated locations. Seismic profiles are displayed for Site GeoB 13824 from S to N, and for Site GeoB 13863 from NE to SW.
2 • C and in situ pH of 8. Compared to seawater at pH 8.2 at 25 • C, where 1.15 mol of S 0 can be solubilized per mol of sulfide (Kamyshny, 2009), only 0.21 mol of S 0 can be solubilized per mol of sulfide at 2 • C because of the lower dissolution of S 0 as polysulfides at lower temperature (Kamyshny et al., 2007). For a chosen sulfide concentration, the value for S 2− n -Calc. is lower than the value of solubilized S 0 because the latter calculates total concentration of zero-valent sulfur from polysulfide species with different sulfur chain length, while the former sums the concentration of polysulfides. For example, 0.21 mol of S 0 can be solubilized in presence of 0.6 mM sulfide at 2 • C, which corresponds to a total polysulfide concentration of 0.037 mM.

Sulfate Reduction Rates
Using the 35 S radiotracer method (Jørgensen, 1978), sulfate reduction rate (SRR) experiments were carried out onboard the ship in a refrigerated container (4 • C). Sediment slurries injected with 35 SO 2− 4 were incubated under in situ pressure (38 MPa; for further details regarding pressure incubations, see Vossmeyer et al., 2012). After the pressure was released, the incubation experiments were terminated by transferring the samples into 15 mL tubes containing a 20% ZnAc solution. These samples were kept frozen at −20 • C during transport and storage. TSRR were determined via the single-step cold chromium reduction (Kallmeyer et al., 2004), and measurements were carried out by scintillation counting at the Center for Geobiology, Aarhus University, Denmark. The detection limit for the SRR was between 0.2 and 1.3 pmol cm −3 d −1 .

Solid Phase Analyses
Multi-acid total digestions (hydrofluoric, HF; hydrochloric, HCl; and nitric, HNO 3 ) were performed on ∼50 mg of dry sediment sample using a microwave system (CEM Mars Xpress) at the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany (AWI). The accuracy of the measurements was verified using NIST SRM 2702 and inhouse (MAX) standards. Major elements were analyzed via inductively coupled plasma-atomic emission spectrometry (ICP-AES, Thermo Scientific IRIS Intrepid instrument). The reference standard contents were within the accepted analytical error for all elements measured. Contents of total carbon (TC) and total inorganic carbon (TIC) were determined by measuring dried and homogenized samples using an ELTRA CS 500 carbon sulfur analyzer equipped with acidification and furnace modules. The accuracy was ±3% and ±4%, respectively. The amount of total organic carbon (TOC) was calculated by subtracting the TIC fraction from TC.
Sequential iron extractions were carried out under anoxic conditions using frozen subsamples. Ascorbate, dithionite, and oxalate steps were applied on ca. 150-200 mg samples to determine, respectively, the fractions present as adsorbed ferrous iron and highly reactive/bioavailable ferric iron (Fe bio ), crystalline iron oxides such as goethite and hematite (Fe oxide ), and magnetite (Fe magn ;Ferdelman, 1988;Poulton and Canfield, 2005;Raiswell et al., 2010;Wehrmann et al., 2014;Henkel et al., 2016). The solutions were flushed with N 2 prior to extraction. All solutions were freshly prepared prior to extraction, and reagent blanks were taken. The iron concentration was analyzed by inductively coupled plasma-mass spectrometry (ICP-MS; Agilent 7500ce) after dilution in trace-metal grade 2% HNO 3 . Replicate sample extractions yielded reproducibility within 8%, and all iron concentrations of reagent blanks were below detection limit.
All solid sulfur phases were analyzed sequentially (for detailed method descriptions see Riedinger et al., 2014). Elemental sulfur (S 0 ) was extracted from approximately 2-3 g of wet (freshly thawed) sediment by shaking for ∼12 h in 10 mL pure methanol with a sample-to-extractant ratio of ∼1/10 (Zopfi et al., 2004). The headspace was flushed with N 2 to avoid oxidation of the reduced species. The concentration of S 0 was analyzed at the MPI-MM using a Sykam pump (S1100), a UV-Vis Detector (Sykam S3200), a Zorbax ODS-column (125 × 4 mm, 5 µm; Knauer, Germany) and 100% methanol (HPLC grade) at a flow rate of 1 mL per minute. Elemental sulfur was eluted after 3.5 min and detected at 265 nm, with a detection limit of about 1 µM and an analytical precision of ±0.5% SD. Based on replicate sample extractions and in-house standards, the precision and accuracy of the elemental sulfur measurements is better than 0.001 wt.%. Therefore, our samples with reported contents below 0.001 wt.% may represent elemental sulfur-free samples. This consideration is important in the discussion of potential presence or absence of polysulfide. For determination of acid volatile sulfide (AVS; mainly iron monosulfides-"FeS") and chromium reducible sulfur (CRS; which, due to the preceding removal of AVS and S 0 , corresponds mainly to pyrite), the samples were treated with the sequential, two-step acid/Cr(II) method (Fossing and Jørgensen, 1989). The sulfide produced in each step was trapped as ZnS in a 5% Zn-acetate solution and analyzed following dilution using the methylene blue method (Cline, 1969). Reproducibility was better than 7% based on an in-house standard. The fraction of total organic sulfur (TOS) was determined on the solid residue following the Cr(II) step (Werne et al., 2003). The samples were filtered and rinsed with double distilled water, dried and analyzed using a carbon-sulfur elemental analyzer (ELTRA CS 500). All solid-phase data are reported in dry weight units.
The degree of pyritization (DOP; Raiswell and Canfield, 1998) was calculated from the analyzed iron phases by dividing pyrite Fe (Fe CRS ; calculated from pyrite S) by the total highly reactive Fe present as iron oxides, iron monosulfide and pyrite Fe (Fe CRS /[Fe bio +Fe oxide +Fe mag +Fe AVS +Fe CRS ]). Highly reactive iron (Fe HR ) was defined according to Raiswell and Canfield (1998) as the sum of Fe bio , Fe oxide , Fe AVS , and Fe CRS .
Sediment porosity was determined according to standard IODP procedures (Blum, 1997) at the MARUM-the Center for Marine Environmental Sciences at the University of Bremen, Germany-using helium-displacement penta-pycnometers. The data were corrected for evaporated seawater, specifically the mass of precipitated salts, as described by .

Sulfur Isotopes
For isotope analyses of sulfate, filtered pore water aliquots were acidified, and sulfate was precipitated as barium sulfate (BaSO 4 ) by addition of barium chloride solution (BaCl 2 , 1 M). For determination of the sulfur isotope compositions of H 2 S, S 0 , AVS, and CRS, ZnS precipitates were converted to Ag 2 S by addition of AgNO 3 and subsequent washing with NH 4 OH to remove colloidal silver. Zinc sulfide was obtained from elemental sulfur by using the hot acid/Cr(II) distillation method. Sulfur isotope ratios of TOS were measured on bulk sediments (Werne et al., 2003). Isotope compositions of pore water sulfatesulfur were determined at the MPI-MM, all other isotope measurements were carried out at the University of California-Riverside (UCR). The measurements were performed by sample combustion with elemental analyzers that were connected via continuous helium flow to Thermo-Finnigan Delta R gas source isotope ratio mass spectrometers. All sulfur isotope measurements were calibrated with reference materials NBS 127 (δ 34 S = +21.1‰) and IAEA-SO-6 (δ 34 S = −34.1‰), and the standard error (1σ) of the measurements was <0.2‰ for δ 34 S.
The sulfur isotope composition is reported with respect to Vienna Canyon Diablo Troilite (V-CDT). In the case of the sulfur isotope measurements of total dissolved sulfide, the sample sizes were very small (∼1/10th typical analysis weight, 0.04 mg), which increases the uncertainty for the reported isotope composition. Based on repeated measurements of very small amounts of our laboratory standard, we estimate the standard error (1σ) for very small samples to be <2‰ for δ 34 S.

RESULTS
The sulfidic zone is defined as the interval where dissolved hydrogen sulfide accumulates in the pore water. The intervals above and below the sulfidic zone are defined as the postoxic and methanic zones, respectively (after Berner, 1981). Due to the typical loss of the top few centimeters of the core during recovery, the uppermost oxic zone was not sampled.

Pore Water
The pore water concentration profiles at both sites show a similar trend, with a linear decrease in sulfate concentrations from about 27.5 mM at the top of the cores to complete depletion at the SMT located at ∼5.5 mbsf at Site GeoB 13824 and ∼5 mbsf at Site GeoB 13863 (Figure 3). Below this transition, methane increases with depth at the two sites to concentrations of 6.1 mM and 9.9 mM, respectively. Pore water accumulations of free sulfide are restricted to a narrow interval close to the SMT-the sulfidic zone. At Site GeoB 13824, this sulfidic zone lies between 3.8 and 6.9 mbsf with maximum H 2 S of 668 µM; at Site GeoB 13863, the sulfidic zone is confined to 3.8-5.7 mbsf, and H 2 S reaches 597 µM. The dissolved iron (Fe 2+ ) concentrations at Site GeoB 13824 show a decrease from the top of the core (15.3 µM) down to about 3.5 mbsf (1.2 µM). Below this depth and throughout the sulfidic zone, dissolved iron was close to or below the detection limit (0.5 µM). Below the sulfidic zone, dissolved iron concentrations strongly increase with depth, with the exception of one small excursion at 8.5-9 mbsf (marked by a drop to 15 µM), with concentrations reaching 72.6 µM ( Figure 3A). The dissolved iron profile at Site GeoB 13863 shows a similar trend to that observed at Site 13824, with slightly elevated values in the upper sediment layers and a maximum concentration of 17.3 µM at 0.8 mbsf and Fe 2+ below detection in the sulfidic zone. Below this zone, dissolved iron concentrations increase again to a maximum of 50 µM before they decrease slightly with depth but remain above 25 µM ( Figure 3B; Riedinger et al., 2014).
At Site GeoB 13824, the redox potential (Eh) decreases linearly with depth, from 142 mV at the sediment surface to −288 mV at 4.45 mbsf. Below this depth, the redox potential increases slightly again with depth but stays below −120 mV. In contrast to the Eh profile at Site GeoB 13824, the redox potential at Site GeoB 13863 shows a strong correlation with dissolved sulfide concentration, including an Eh decrease to −244 mV in the sulfidic zone. Above and below this zone, Eh-values increase to 204 and 74 mV, respectively. The pH at this site does not show strong variation and stays in the range of 7.66-8.14, with slightly elevated values in the sulfidic zone. No pH measurements are available for Site GeoB 13824.

Sulfate Reduction Rates
Sulfate reduction rates (SRR) are only available for sediments from Site GeoB 13824. The rates are low (<12 pmol cm −3 d −1 ) throughout the core with values close to or below the detection limit. As a likely consequence of the low sampling resolution, the expected higher SRR at the sediment surface, as well as at the SMT as related to sulfate reduction coupled to anaerobic methane oxidation, were not captured.

Solid Phase
The investigated sediments are dominated by detrital material with low contents of calcium carbonate (CaCO 3 ) restricted to the upper 1.8 m at Site GeoB 13824 and the upper 0.5 m at Site GeoB 13863, along with high amounts of total iron ranging from 2.6 to 4.7 wt.% and aluminum from 6 to 9 wt.% (Figure 4). The Fe/Al ratios at both sites stay constant throughout the sediment column, with an average value of 0.5, similar to an average crustal value (Taylor and McLennan, 1985). At Site GeoB 13824, the CaCO 3 contents show a peak between 5 and 6 mbsf, corresponding to the depth of the SMT. The sediments at Site GeoB 13863 do not display such a peak likely because of our sampling protocol. Specifically, prior to sampling, discrete carbonate crystals were removed from about 5 mbsf, within the SMT, for separate analysis and are thus not included in the bulk element contents. Thus, for both sites, carbonate is present at the depth of the SMT, as depicted in the sediment log (Figure 4). TOC contents at Site GeoB 13824 show only slight variation throughout the core, with a range from 1.43 to 0.76 wt.%; sediments in the uppermost meter and at ∼2.5 mbsf show the highest contents (>1 wt.%). At Site GeoB 13863, the TOC contents range from 1.15 wt.% (in the uppermost 20 cm) to 0.44 wt.% and average 0.78 wt.% .
Total sulfur contents at Site GeoB 13824 are between 0.14 and 0.53 wt.% above and below the sulfidic zone. Within the sulfidic zone, values increase to up to 0.89 wt.%. Similarly, the total sulfur contents at Site GeoB 13863 show an increase within the sulfidic zone (up to 1.26 wt.%), while above and below this zone the values are low (0.09-0.42 wt.%). The iron oxide phases (Fe oxide ) at Site GeoB 13824, including labile iron (oxyhydr)oxides and crystalline phases (hematite and goethite), are in the range of 0.24 to 0.05 wt.% in the upper 6 m and strongly increase (up to 0.51 wt.%) in the deeper sediments-starting at the lower boundary of the sulfidic zone, with a drop at 8.15 mbsf to a value of 0.13 wt.%. The concentration profile for iron in magnetite (Fe magn ) resembles that for Fe oxide , with a less pronounced decrease in the sulfidic layer (0.13 wt.%) and a smaller increase below the sulfidic zone (up to 3.7 wt.%). The Fe oxide and Fe magn contents at Site GeoB 13863 scatter between 0.07 and 0.41 wt.%, with increased values above and below the sulfidic zone at 1.2-4 mbsf and 5.6-6.5 mbsf, respectively. While the ratio of highly reactive iron (Fe HR ) to total iron (Fe T ) is similar at both study sites, with values between 0.1 and 0.25, the DOP shows a strong difference: DOP at Site GeoB 13824 is elevated between 1 and 6.2 mbsf, with values reaching 0.73 (with one outlier of 0.18 at 1.5 mbsf). DOP at Site GeoB 13863 displays increased ratios up to only 0.63 in the sulfidic zone. Outside these intervals of increased DOP at both sites, the ratio stays low (<0.4) throughout the remaining sediment column.

Stable Sulfur Isotope Composition of Pore Water and Solid Phase Sulfur Compounds
The stable sulfur isotope composition of pore water sulfate (δ 34 S-SO 4 ) displays a typical enrichment in 34 S with depth-with a starting value of δ 34 S-SO 4 (+21.1‰) at the sediment surface at both sites equal to the global seawater values (+21.1‰ e.g., Rees, 1970;Böttcher et al., 2007). The δ 34 S-SO 4 values increase with depth to +64.7‰ at 5.2 mbsf and +54.0‰ at 4.5 mbsf at Sites GeoB 13824 and 13863, respectively (Figure 5). The few data points for the sulfur isotope composition of total free sulfide (δ 34 S-H 2 S) show no clear trend and scatter between +16.4 and +43.0‰ at both sites ( Figure 5). Some of the scatter in the δ 34 S-H 2 S data may be the result of very low concentrations (because of limited sample volume) and associated instrumental uncertainty at these levels. Consequently, interpretations of δ 34 S-H 2 S data should be viewed with caution.
Pyrite contents (CRS) at Site GeoB 13824 are elevated above and within the sulfidic zone, reaching values of up to 0.5 wt.%. In contrast, in the uppermost meter of the core and below the sulfidic zone, pyrite contents remain below 0.2 wt.% with one exception (0.27 wt.%) at 9.45 mbsf ( Figure 5A). The stable sulfur isotope composition of pyrite (δ 34 S-CRS) shows a trend similar to the δ 34 S-SO 4 profile, starting with −49.1‰ at the surface and showing the highest enrichment in 34 S (+29.0‰) at 6.5 mbsf-just below the SMT. Below the sulfidic zone, the δ 34 S-CRS data show a slight progressive depletion in 34 S with depth, from −9.8 to −39.3‰ in the lowermost sediment layer. In contrast to pyrite, the contents of iron monosulfide phases (AVS) are low (<0.016 wt.%) in the upper 7.5 meters at Site GeoB 13824 and increase below the sulfidic zone to a maximum of 0.042 wt.% at 8.45 mbsf (Figure 5A). The stable sulfur isotope composition of AVS (δ 34 S-AVS) at the same site varies between −33.0 and +32.4‰, with the highest values found below the sulfidic zone. The contents of S 0 show similar trends to AVS, with highest contents below the sulfidic zone reaching 0.164 wt.%. Although the stable sulfur isotope composition of elemental sulfur (δ 34 S-S 0 ) scatters throughout the sediment core, the heaviest values are observed below the sulfidic zone, with values of up to +46.5‰. TOS contents are low (<0.03 wt.%) in the upper 2.45 meters at this site and below the sulfidic zone ( Figure 5A). Distinct peaks in TOS are observed slightly above of the sulfidic zone, between 2.5 and 3.5 mbsf, at the top (3.5 to 4.5 mbsf) and within the sulfidic zone (5-6.5 mbsf). The stable sulfur isotope composition of TOS (δ 34 S-TOS) is more negative above the SMT (between −34.2 and −16.4‰) compared to values below of the SMT (−13.6 to +19.8‰).
The pyrite content at Site GeoB 13863 in the sulfidic zone ranges up to 0.72 wt.%, while below and above this zone, the pyrite contents stay below 0.34 wt.% (Figure 5B). Similar to the pyrite concentration profile, the δ 34 S-CRS data display a positive excursion in the sulfidic zone, reaching +50.9‰. Above and below the sulfidic zone, δ 34 S-CRS data stay below −11‰, with heavier values below and lighter values above. In contrast to pyrite contents, increased amounts of AVS are limited to the intervals above and below the sulfidic zone, with values of up to 0.012 wt.% ( Figure 5B). Outside of these intervals, AVS stays below 0.008 wt.%. Due to the low contents, no δ 34 S-AVS measurements could be carried out on samples from the sulfidic zone. Below the sulfidic zone, however, δ 34 S-AVS values are measurable and high (up to +39.5‰), while values are between −32.2 to −18.8‰ above the sulfidic zone. Profiles very similar to AVS are observed for elemental sulfur, with low contents (<0.015 wt.%) except for narrow intervals directly above and below the sulfidic zone where contents increase to up to 0.071 wt.%. δ 34 S-S 0 values are negative above the sulfidic zone (<−20‰), whereas elemental sulfur is enriched in 34 S below this zone with values reaching +44.8‰. Two TOS peaks are observed within the sulfidic zone, one at the upper and one at the lower boundary, with contents of up to 0.194 wt.% ( Figure 5B). Values above and below the sulfidic zone are constant at ∼0.02 wt.%. Above the SMT, δ 34 S-TOS is essentially constant around −10‰ with two exceptions, one at the surface (+3.0‰) and one at 3.85 mbsf (−27.3‰). Below the SMT, δ 34 S-TOS values are around 0‰ (from −3.7 to +2.8‰), except for the lower TOS peak (at the lower boundary of the sulfidic zone) where the 34 S is enriched (up to +45.3‰).

DISCUSSION
The investigated sediments are dominated by terrigenous inputs of silt and clay with high amounts of total iron, including abundant reactive ferric iron minerals (Figure 4). The observed concentration profiles of the various iron sulfide phases reflect the ongoing alteration of iron oxides via diverse reaction pathways several meters below the seafloor. The high amounts of reactive iron are the cause for the lack or low concentrations of dissolved hydrogen sulfide over most of the sediment column (for a detailed discussion on iron cycling in these sediments see Riedinger et al., 2014). Produced sulfide is immediately scavenged, either via oxidation coupled to the reduction of the iron oxides or by precipitation with dissolved ferrous iron produced by the reduction of ferric iron (Berner and Westrich, 1985;Hartgers et al., 1997;Riedinger et al., 2014). Sulfide can only build up in the pore water when the rates of sulfide formation exceed those of iron-phase alteration, leading to the establishment of a sulfidic zone (Figure 3). Because of the low reactivity of the organic matter , SRR are extremely low even in the presence of abundant sulfate in the upper meters ( Figure 3A), with average rates of 5.8 pmol cm −3 d −1 -about one order of magnitude lower than rates usually found in sediments from similar water depths (e.g., Fossing et al., 2000;Sawicka et al., 2012). Thus, sulfide production is mainly restricted to the SMT. The presence of free hydrogen sulfide only in the center of the SMT and the sulfide-limited conditions above and below result in a distinct sequence of metastable sulfur phases and pyrite across the sulfidic zone. Ultimately, this geochemical regime leads to the accumulation of metastable sulfur phases, such as authigenic monosulfides (Figures 5, 6) and elemental sulfur, as well as sulfurization of organic matter several meters below the sea floor.

Pyrite Formation at the Center of the Sulfidic Zone
Elemental sulfur and iron monosulfides form at the upper and lower boundary of the sulfidic zone. In contrast, elevated pyrite concentrations are found at the center of this zone. This relationship suggests that the pyrite is mainly formed in the presence of excess hydrogen sulfide. The relatively high (for an iron-dominated system) DOP ≥ 0.6 for the sediments within the sulfidic zone at both sites indicates almost complete alteration of the iron oxide phases, including magnetic minerals. Canfield and Berner (1987) postulated that magnetic minerals (such as magnetite) can be replaced by pyrite if sulfide concentrations remain high for relatively long periods of time (several hundreds of years), leading to an alteration (or loss) of the magnetic signal. The distribution of iron sulfide phases at Site GeoB 13863 agrees well with the results from numerical reactive-transport modeling for the depositional scenario specific to this location (see Riedinger et al., 2005). Specifically, the modeling predicts elevated pyrite contents that are limited to the sulfidic zone, coinciding with a strong decrease in magnetic susceptibility . Similar observations are also reported for other iron-dominated areas such as the Amazon Fan (Jørgensen and Kasten, 2006), Zambesi Fan (März et al., 2008) or cold-seep systems off southwestern Taiwan (Hsu et al., 2014).

Pyrite Outside of the Sulfidic Zone
Relatively high DOP values are also observed in the postoxic zone at Site GeoB 13824. This relationship could be attributed to sulfide accumulation from organoclastic sulfate reduction related to slightly higher TOC amounts and a longer duration of pore water steady-state conditions at this site compared to Site GeoB 13863. Alternatively, it is possible that at Site GeoB 13824, the SMT was previously located at a shallower location in the sediment. This, however, is not supported by the δ 34 S-CRS data, which remain low in the postoxic zone. A different possibility is that lower availability of reactive iron phases at Site GeoB 13824 relative to Site GeoB 13863 would also result in higher sulfide accumulation rates and thus more complete pyritization of the reactive phases.
In light of the extremely low measured SRR and the presence of high amounts of reactive iron phases, we are forced to question whether all of the pyrite in the surface sediments formed in situ. The highly 34 S-depleted pyrite isotope signature (∼ −30 to −50‰ in the upper 4 m; Figure 5) is very similar to other sulfur isotope signatures in shallow marine environments (e.g., Bottrell et al., 2009). We suggest that a major portion of the pyrite found in the surface sediments is derived from sediment reworking of shallow-water sediments, where organoclastic sulfate reduction and associated pyrite formation is more prevalent.

Formation of Iron Monosulfide and Elemental Sulfur
Intermediate sulfur phases are mainly observed below the sulfidic zone at Site GeoB 13824. At Site GeoB 13863, these intermediate sulfur phases also occur at the upper and the lower boundary of the sulfidic zone (Figure 5). Formation of intermediate sulfur phases at the upper and lower boundaries of the sulfidic zonethat is, simultaneous formation in sediments of different agesreflects the presence of high amounts of available iron (oxyhydr) oxide phases that react with sulfide to form zero-valent sulfur and ferrous iron (e.g., Aplin et al., 1993), while the precipitation of iron monosulfides is driven by the reaction of hydrogen sulfide with ferrous iron. The ferrous iron can result from iron (oxyhydr) oxide reduction coupled to in situ sulfide oxidation, but it can also diffuse to the sulfidic zone-originating from the upper sediment layers and associated organoclastic iron reduction and the lowermost sediments via iron-reduction coupled to AOM as discussed by Riedinger et al. (2014) for Site GeoB 13863. In absence of free sulfide, or at low pH (Kamyshny et al., 2004(Kamyshny et al., , 2007(Kamyshny et al., , 2008Kamyshny and Ferdelman, 2010), elemental sulfur is the stable zero-valent sulfur phase. At the investigated sites, the pH is high (7.66-8.14 in the sulfidic zone at Site GeoB 13863). This pH range makes it likely that elemental sulfur is replaced by polysulfides at the fringes of the sulfidic zone (Rickard and Luther, 2007;Kamyshny, 2009), which is corroborated by the near absence of elemental sulfur within the sulfidic zone (Figures 5, 6).
The co-occurrence of elevated elemental sulfur and AVS below the sulfidic zone indicates that in a sulfide-limited system, metastable minerals can persist and be buried into deeper sediment depths (Figure 5). This relationship indicates that intermediate sulfur phases can be more stable in situ than expected based on laboratory and/or modeling results. The availability of reactive iron oxide minerals in these deeper sediments can inhibit the transformation of iron monosulfide and elemental sulfur to pyrite, a reaction that likely only proceeds via an aqueous phase (e.g., polysulfides) (e.g., Rickard and Luther, 2007). The buried elemental sulfur might provide a deep-subsurface sulfur source for microbial communities, such as sulfur disproportionating organisms (e.g., Thamdrup et al., 1993), or for the reduction to sulfide in the methanic zone by archaea (Stetter and Gaag, 1983)-and may even lead to the accumulation of a deep sulfate pool (Riedinger et al., 2010;Treude et al., 2014).

Sulfurization of Organic Matter
The low δ 34 S-TOS values in the uppermost sediments at both sites ( Figure 5) are comparable to those in other continental margin settings (e.g., Bottrell et al., 2009). In the presence of excess hydrogen sulfide, sulfurization of selected labile organic compounds appears to precede the formation of pyrite in marine surface sediments (François, 1987;Ferdelman et al., 1991;Mossmann et al., 1991;Vairavamurthy et al., 1992Vairavamurthy et al., , 1995Amrani and Aizenshtat, 2004). The limited availability of such labile organic compounds in the investigated uppermost sediments concurrent with high amounts of reactive iron suggest, however, that a the marginal sulfurized organic compounds, depleted in 34 S, were mainly transported to this site, most likely from reworked shallow slope and/or shelf sediments (see also Riedinger et al., 2014).
Polysulfides are known catalysts for the sulfurization of organic compounds (e.g., François, 1987;Werne et al., 2000Werne et al., , 2008. Thus, enhanced sulfurization of organic matter can be expected to occur in geochemical zones where polysulfide formation is prevalent. Elevated TOS contents are not found within the entire sulfidic zone but rather occur in two distinct peaks. At site GeoB 13863, these two peaks align almost perfectly with the transition between the sulfide-free and sulfide-rich portions of the sediment, exactly at the position where one would expect to observe the accumulation of polysulfide-because of the peak elemental sulfur formation and subsequent reaction with sulfide-as shown by the calculated potential polysulfide concentrations ( Figure 6B). The δ 34 S-TOS data indicate different sources of sulfur for the upper and lower TOS peaks, likely reflecting the δ 34 S-H 2 S at those depths. There appears to be a third TOS peak in the sulfide-free postoxic zone at site GeoB 13824; however, the same overall picture emerges-an absence of elemental sulfur (and to a lesser degree absence of AVS) coincides with a peak in TOS ( Figure 5A). The TOS content of sediments located in the sulfidic zone between the two TOS enrichment peaks fall to almost background levels, likely indicating that sulfurization of particulate organic matter occurs at the sulfidic fringes of the sulfide zone but not in the center of the SMT where hydrogen sulfide displays the highest concentrations. Polysulfides have been shown to be stronger nucleophiles for the sulfurization of organic matter, and nucleophilic substitution has also been recognized to be the dominant sulfurization mechanism (Amrani, 2014, and references therein). If sulfurization of organic matter is indeed tied to polysulfides, our findings imply that the center of the SMT has lower polysulfide concentrations than the border region of the SMT. This could be explained by the absence of polysulfides due to a lack of oxidants.

IMPLICATIONS FOR THE GEOCHEMISTRY OF DYNAMIC, IRON-DOMINATED SEDIMENTARY SYSTEMS
The lack of dissolved hydrogen sulfide in the pore water of the surface sediments at our sites (Figure 3) indicates that the sulfide produced by organoclastic sulfate reduction reacts immediately with dissolved ferrous iron and thus little or no polysulfide is produced. In the absence of polysulfides, organic substances appear to be protected from sulfurization in a major portion of the sediment column. This observation is in good agreement with similar findings, for example, from the Cariaco Basin (Werne et al., 2003(Werne et al., , 2008. Since sulfide accumulation in the pore water is often restricted to a distinct zone, as produced by sulfate reduction coupled to anaerobic methane oxidation, the fringes of the SMT become hotspots for polysulfide formation and consequent organic matter sulfurization (Quijada et al., 2016). With the SMT being a zone of enhanced organic matter sulfurization, there is potential for enhanced preservation of organic molecules in this zone, which implies that biomarkers from organisms that thrive in such zones have an elevated preservation potential (Hebting et al., 2006).
The formation of S 0 at the SMT margins due to diffusion of sulfide into ferric-rich sediments, which then reacts with sulfide to form polysulfide, coincides with the occurrence of AVS that is the product of sulfide reacting with counter-diffusing ferrous iron. This raises the interesting question about the competition for polysulfide for AVS to pyrite transformation vs. organic sulfurization and how that plays out in the core of the SMT vs. the margins.

Competition between Sulfurization of Organic Matter and Pyrite Formation
There are two major pyrite formation mechanisms: the polysulfide and the hydrogen sulfide pathways (for a review, see Butler et al., 2004). The polysulfide pathway (Berner, 1970;Rickard, 1975;Sørensen and Jørgensen, 1984;Schoonen and Barnes, 1991;Canfield and Thamdrup, 1994;Rickard and Luther, 2007) consists of the transformation of solid iron monosulfide to an aqueous species, which can react with polysulfide to form pyrite, and a shorter-chained polysulfide species, according to The hydrogen sulfide pathway follows as similar pattern, with an initial transformation of solid iron monosulfide to an aqueous iron monosulfide, which then reacts with hydrogen sulfide, whereby the hydrogen ions are reduced to dihydrogen gas (e.g., Rickard and Luther, 1997;Butler and Rickard, 2000;Rickard and Luther, 2007), according to FeS (s) + H 2 S (aq) → FeS (aq) + H 2 S (aq) = FeS 2(s) + H 2(g) (2) The latter reaction only proceeds with H 2 S, and not bisulfide (HS − ). This means that for the reaction to occur, neutral or slightly acidic conditions are preferred; however, the reaction can also take place under slightly basic conditions (Rickard and Luther, 1997).
At our sites, sulfurization of reworked organic matter at the edges of the SMT appears to out-compete pyrite formation, as is evidenced by the peaks in TOS at the fringes of the sulfidic zonespecifically, organic matter is sulfurized, while AVS is preserved. In the center of the sulfidic zone pyrite predominates and no further generation of sulfurized organic matter occurs. Here, we consider two speculative explanations for this observation. The scenario presumes that polysulfide is also present in the sulfidic zone, but instead of being used for sulfurization of organic matter it is consumed in the formation of pyrite (Equation 1). In this scenario, one would have to speculate that under high sulfide concentrations, the competition between sulfurization of organic matter and pyrite formation is tilted in advantage of the latter.
Alternatively, one could consider the predominance of pyrite within the center of the sulfidic zone an indicator for absence of polysulfide, this would explain why no further generation of sulfurized organic matter occurs. One would then conclude that pyrite formation proceeds via the hydrogen sulfide pathway. It is important to note that in our sediments the pH is up to 8.14 in the center of the sulfidic zone-conditions that are not favorable for the hydrogen sulfide reaction (Equation 2) to occur, because most of the sulfide prevails as bisulfide. To solve this contradiction, one could speculate that the precipitation of authigenic calcite or dolomite from bicarbonate driven by the AOM coupled to sulfate reduction at the SMT (e.g., Kelts and McKenzie, 1982;Malone et al., 2002;Moore et al., 2004;Meister et al., 2007;Nöthen and Kasten, 2011;Wehrmann et al., 2011;Riedinger et al., 2014) liberates hydrogen ions, which pushes sulfide speciation slightly toward H 2 S. This sulfide then becomes available for pyrite formation (Equation 2). Thus, the precipitation of authigenic carbonate minerals, as observed at the SMT at both sites (Figure 4), could facilitate the formation of pyrite at the SMT via the hydrogen sulfide reaction (Equation 2). The combination of carbonate precipitation with pyrite formation balances hydrogen ion production and consumption, with no net change in pH, which remains high at a value around 8.

Interpretation of the Geochemical Evolution of Dynamic Sediments Based on Inventory and Isotope Composition of Sulfur Species
Sulfur and iron transformations in the sediments of the Argentine Basin occur in a dynamic sedimentary system. Rapid overall sedimentation is sometimes interrupted by phases of sediment winnowing caused by ocean currents. Episodic mass wasting processes can lead to instantaneous sediment accumulation (e.g., Hensen et al., 2003;Riedinger et al., 2005;Henkel et al., 2011Henkel et al., , 2012. These changes cause non-steady state conditions in the subsurface sediment/pore-water system (e.g., Kasten et al., 2003). A key question is if the analysis of the inventory of iron and sulfur species and the stable sulfur isotope composition of the sulfur constituents can fingerprint similar processes in other settings. To illustrate the potential of this approach, we present a three-stage geochemical scenario for the two investigated sites (Figure 7). The first stage is dominated by high sedimentation, as has been discussed previously by Hensen et al. (2003) and Riedinger et al. (2005) in the context of numerical modeling approaches. The high sedimentation rates could be due to one or multiple depositional events. During this time, reworked sediment from upslope settings accumulate, carrying TOS and pyrite with sulfur isotope signatures that correspond to these shallower settings ( Figure 7A). These phases are intermingled with refractory organic matter and ample amounts of iron (oxyhydr)oxides. A large portion of the rapidly buried organic matter reaches the methanic zone, resulting in an increase in the methane flux from below. During this phase, the sulfate concentration profile displays a "concave up" shape ( Figure 7A; Hensen et al., 2003;Kasten et al., 2003).
As the SMT migrates upward, the production of sulfide by sulfate reduction coupled to AOM sweeps across sediment FIGURE 7 | Schematic model displaying iron oxide alteration and sulfur phase distribution with its isotopic signature related to changes in depositional conditions (modified after Riedinger, 2005). (A) During high sedimentation (one or multiple depositional events) the SMT moves rapidly upward, a non-steady state condition that manifests itself as concave-up sulfate profile . Under these conditions, sulfide (H 2 S) produced by sulfate (SO 2− 4 ) reduction coupled to anaerobic oxidation of methane (CH 4 ) at the SMT is consumed faster than it can build up in the sediment by reacting with highly reactive iron oxide phases ( Fe oxides ; area shaded in light gray) and dissolved ferrous iron (Fe 2+ ). This results in the formation of intermediate sulfur phases (elemental sulfur, S 0 , and monosulfide phases, FeS) with distinct stable sulfur isotope signatures (δ 34 S). Only minor amounts of pyrite (FeS 2 ) are formed, the observed pyrite and TOS were mainly derived by sedimentation. (B) Once sedimentation rates decrease, H 2 S builds up at the SMT, which leads to the formation of pyrite in the sulfidic zone, as well as to the sulfurization of organic compounds (TOS) at the upper and lower rim of the sulfidic zone. Adjacent to the sulfidic zone, iron monosulfides and elemental sulfur are formed. All of these phases show stable sulfur isotope signatures that reflect the isotope composition of sulfide at the respective position in the sediment. These signatures are distinctively different from the isotope signatures of the same sulfur constituents that formed elsewhere in the sediment column, or were derived by sedimentation. (C) At low sedimentation rates, hydrogen sulfide can accumulate at the SMT and spread outwards. The result is a broadening of the sulfidic zone and consequently broader CRS and TOS peaks. that contains abundant reactive iron phases. The sulfide reacts immediately with these iron phases ( Fe oxides ; area shaded in light gray in Figure 7)-producing zero-valent sulfur phases (i.e., elemental sulfur) and ferrous iron, which in turn reacts with sulfide to form iron monosulfide. These newly formed sulfur phases carry the sulfur isotope signature of sulfide formed at the upward migrating SMT. Due to the refractory nature of the organic matter, organoclastic sulfate reduction is of minor importance, and only minor amounts of pyrite are formed, which may marginally alter (slight enrichment in 34 S) the sulfur isotope signature of the reworked pyrite delivered by sedimentation.
The second stage in our geochemical scenario is reached when sedimentation rates (strongly) decrease and the sulfate flux is primarily controlled by diffusion again, which causes a slowing down of the upward migration of the SMT. Under these conditions, the production of sulfide at the SMT is no longer out-competed by sulfide consumption through reaction with reactive iron phases. As a consequence, sulfide starts to build up in the pore water, and a sulfidic zone is established, which broadens over time. At the center of the sulfidic zone, pyrite is formed via the reaction of iron monosulfide with H 2 S, probably in conjunction with calcite precipitation. At the fringes of the sulfidic zone, iron monosulfides, and elemental sulfur are formed through the interaction between dissolved Fe 2+ and reactive iron oxide species and sulfide diffusing away from the center of the sulfidic zone. At these reaction fronts, found at the upper and lower boundaries of the sulfidic zone, additional hydrogen sulfide reacts with the produced elemental sulfur-resulting in the formation of dissolved polysulfides, which in turn sulfurize the organic matter. The sulfurization results in two TOS enrichment fronts, one near the upper and one at the lower boundary of the sulfidic zone, with less and more 34 S-enriched signatures, respectively (Figures 7B,C). This process also explains why TOS and elemental sulfur phases do not overlap.
The third stage in the scenario is reached when the sedimentation rates become very low (returning to mainly hemipelagic sedimentation). The system then adjusts to a quasisteady state in terms of the pore water, and the upward moving SMT is now "fixed" at a certain/relative depth (moving upward slowly at the same rate as sedimentation). Hydrogen sulfide accumulating at the SMT then spreads further out, resulting in a broadening of the sulfidic zone and broader CRS and TOS peaks ( Figure 7C). This is the stage that is captured in the investigated sediments, with steady state pore water profiles as displayed by sulfate, while the solid phases are integrated and mostly capture the preceding non-steady state conditions. Numerical modeling results for sediments from nearby locations suggest a timespan of several hundred to a few thousand years for the pore water sulfate concentration profile to regain the observed steady state conditions .
This three-stage scenario explains most of the observed patterns at Sites GeoB 13824 and 13863. Some refinements are required, however, to address phenomena specific to individual sites. For example, there is a third TOS peak in the postoxic zone of Site GeoB 13824, indicating that substantial polysulfide formation also took place at a shallower depth in the past. This relationship could be explained by an excursion of the SMT into shallower sediments, followed by a deepening of the SMT due to higher downward flux of sulfate into the SMT or a lower upward flux of methane. Such fluctuations can take place, for example, when sediment is removed by winnowing or if less sulfate is consumed by organocastic sulfate reduction during low sedimentation due to an increase in the refractory nature of organic matter. If the system experienced a constant low sedimentation rate over a long period of time, some of the observed isotope excursions could be overprinted. For example, the isotope signatures associated with the upper fringe of the SMT signal would be overprinted by the upward movement of the redox zones with a lower SMT fringe signature, which may add CRS and TOS with distinctly different isotope signatures to the bulk signal. In a different scenario, where high sedimentation rates follow a short period of low sedimentation, the CRS and TOS associated with the SMT would be buried, and the observed isotope excursions would be altered slowly over time if at all. The different isotope signatures, and the potential stacking of such signatures during shifts in methane or sulfate fluxes, opens the possibility of using high-resolution sulfur isotope signatures from organic matter to reconstruct past changes in the location of the SMT (Wehrmann et al., 2013).

SUMMARY AND CONCLUSIONS
Our study provides a refined view of the zonation of biogeochemical processes that occur within and in the vicinity of the sulfidic zone around the SMT in sedimentary systems dominated by reactive iron. Pyrite is formed in the center of the sulfidic zone. At the upper and lower fringe of the sulfidic zone, the interval where pyrite formation prevails is bounded by a zone in which sulfurization of organic matter dominates. The polysulfides are supplied by the reaction of free sulfide with zerovalent sulfur, which in turn is supplied by the oxidation of sulfide with iron (oxyhydr) oxides. Considering the refractory nature of the organic matter, it is interesting that sulfurization outcompetes pyrite formation. Above and below the sulfidic zone, elemental sulfur, and iron monosulfide phases build up, a process that is fueled by sulfide oxidation, which yields S 0 and ferrous iron and induces the precipitation of iron monosulfides through reaction with ferrous iron. Because of the distinct locations of these processes with respect to the center of the SMT, the produced sulfur phases record sulfur isotope signatures for sulfide that are representative for the isotope trends observed in the SMTthat is, a strong enrichment in 34 S from the top to the bottom of the SMT (e.g., Rudnicki et al., 2001;Brunner et al., 2016;Turchyn et al., 2016). These distinct isotope signatures allow the reconstruction of the complex history of biogeochemical sulfur cycling in dynamic sediments. The power of this approach has previously been demonstrated for pyrite enrichment-fronts that indicate the location of past (paleo or fossil) SMTs (Borowski et al., 2013;Lin et al., 2016) but can now be expanded and refined by the inclusion of the inventory and stable sulfur isotope composition of other sulfur constituents, such as the sulfur isotope composition of the organic phase (δ 34 S-TOS). We believe these analyses can be greatly augmented by employing highresolution SIMS techniques to identify the former presence of upper and lower fringes of SMTs, particularly in cases where such signatures were stacked on top of each other during fluctuations in the depth of the SMT. Such zones may have two isotopically distinct populations of pyrite and TOS, corresponding to the upper and lower SMT fringe. Techniques such as sulfur isotope analyses using SIMS have the potential to reveal such variations and clustering of isotope signatures that remain undetected in a bulk sample (Xiao et al., 2010;Bontognali et al., 2012;Farquhar et al., 2013;Fischer et al., 2014).
Our findings not only facilitate the reconstruction of the geochemical history of dynamic sediments, they also shed new light on several interlinked themes in the study of iron-sulfur cycling in marine sediments in iron dominated systems: (i) the occurrence of elemental sulfur and iron monosulfides at the upper and lower boundary of the sulfidic zone indicates that in a sulfide-limited system these metastable minerals are not (immediately) converted into pyrite and thus can be buried into deeper sediment depths where they can fuel deep biosphere processes (e.g., sulfate reducers via S 0 disproportionation); (ii) the sulfurization of organic compounds predominantly occurs when polysulfides become available in the margins of the sulfidic zone, leading to a preferential fixation and long-term preservation of biomarkers that occur in this zone. Our data provide evidence that organic matter sulfurization can occur at any sediment depth as long as hydrogen sulfide and oxidants are co-mingled to form dissolved polysulfides; (iii) pyrite is formed in the core of the sulfidic zone, likely via the reaction of iron monosulfide with H 2 S, a process that produces H 2 gas and that may require or induce precipitation of authigenic carbonate minerals and (iv) considerable amounts of reactive iron persist in sediment over long time due to rapid relocations of the SMT (thus short sulfide exposure time), which shifts hot spots of sulfide generation to specific intervals in the sediment column, where the reactive iron is more rapidly depleted due to the presence of free sulfide compared to the intervals below and above of this zone.
A further conclusion of our study is that the sulfur isotope signature of iron sulfide and organic sulfur phases in marine sedimentary systems is influenced by depositional as well as in situ geochemical processes-that is, there is a link between deposition in the shallow subsurface sediments and longterm signals being buried and preserved in the sedimentary record. Thus, dynamic depositional systems are characterized by non-traditional redox zonation, where inverse redox zonation sequences or overlap of specific zones can occur, resulting in multiple intervals of specific redox processes in different sediment strata. In other words, comparable geochemical reactions proceed in sediments of different sediment depth and ages.

AUTHOR CONTRIBUTIONS
All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.