Oxygen Consumption and Sulfate Reduction in Vegetated Coastal Habitats: Effects of Physical Disturbance

anoxic conditions, especially in mangrove and seagrass sediments, as well as sediment acidiﬁcation with depth, likely decreased microbial remineralisation rates of sedimentary carbon. However, physical disturbance of sediments and thereby exposure of deeper sediments to O 2 seemed to stimulate aerobic metabolism in the exposed surface layers, likely reducing carbon stocks in VCHs. material increases microbial activity at the sediment surface, ensuring low O 2 penetration. (3) High carbon accumulation in deep, anoxic sediment (cm-m). (4) Low resuspension of sediment due to reduced ﬂow across sediment. (5) Sediment acidiﬁcation with depth potentially leading to inhibition of bacterial growth. Hence, VCHs, such as mangrove forests, salt marshes and seagrass meadows, are characterized by three important environmental properties: (i) high organic loading, (ii) high carbon burial rate, and (iii) efﬁcient protection of the buried carbon, as opposed to unvegetated sediment.

Vegetated coastal habitats (VCHs), such as mangrove forests, salt marshes and seagrass meadows, have the ability to capture and store carbon in the sediment for millennia, and thus have high potential for mitigating global carbon emissions. Carbon sequestration and storage is inherently linked to the geochemical conditions created by a variety of microbial metabolisms, where physical disturbance of sediments may expose previously anoxic sediment layers to oxygen (O 2 ), which could turn them into carbon sources instead of carbon sinks. Here, we used O 2 , hydrogen sulfide (H 2 S) and pH microsensors to determine how biogeochemical conditions, and thus aerobic and anaerobic metabolic pathways, vary across mangrove, salt marsh and seagrass sediments (case study from the Sydney area, Australia). We measured the biogeochemical conditions in the top 2.5 cm of surface (0-10 cm depth) and experimentally exposed deep sediments (>50 cm depth) to simulate undisturbed and physically exposed sediments, respectively, and how these conditions may affect carbon cycling processes. Mangrove surface sediment exhibited the highest rates of O 2 consumption and sulfate (SO 4 2− ) reduction based on detailed microsensor measurements, with a diffusive O 2 uptake rate of 102 mmol O 2 m −2 d −1 and estimated sulfate reduction rate of 57 mmol S tot 2− m −2 d −1 . Surface sediments (0-10 cm) across all the VCHs generally had higher O 2 consumption and estimated sulfate reduction rates than deeper layers (>50 cm depth). O 2 penetration was <4 mm for most sediments and only down to ∼1 mm depth in mangrove surface sediments, which correlated with a significantly higher percent organic carbon content (%C org ) within sediments originating from mangrove forests as compared to those from seagrass and salt marsh ecosystems. Additionally, pH dropped from 8.2 at the sediment/water interface to <7-7.5 within the first 20 mm of sediment within all ecosystems. Prevailing

INTRODUCTION
Vegetated coastal habitats (VCHs) such as seagrass meadows, salt marshes and mangrove forests are increasingly recognized as high-value ecosystems (Costanza et al., 1997;Harborne et al., 2006;Barbier et al., 2011) due to their extraordinary high carbon capture (sequestration) and storage capacity in sediments (Duarte et al., 2005;McLeod et al., 2011;Fourqurean et al., 2012;Ricart et al., 2015;Atwood et al., 2017), frequently termed "blue carbon". VCHs also provide other important ecosystem services such as protecting coastal areas against erosion (Ward et al., 1984;Fonseca and Cahalan, 1992) and increasing marine biodiversity . Nevertheless, VCHs are being lost and existing in degraded conditions due to conversion to alternative land-uses (Alongi, 2002), large-scale physical damage [e.g., dredging, dragging and urban development (Fonseca et al., 1984;Peterson et al., 1988;Neckles et al., 2005)] and degradation [e.g., moorings ] and increased land run-off and sedimentation (Brodersen et al., 2017a), putting essential ecosystem services at risk. Sediments of VCHs are generally considered anoxic due to high organic carbon input and microbial activity Hansen et al., 2000;Nielsen et al., 2001;Holmer and Laursen, 2002) and slow diffusion of O 2 from seawater into the porewater. The high microbial activity is driven by large pools of organic material originating from the plants themselves (Miyajima et al., 1998) or as a result of increased sedimentation capturing organic material from the water column and adjacent land catchments (Burdige, 2005). Sedimentary carbon occurs in the form of both particulate organic carbon (e.g., plant material) and dissolved organic carbon (DOC) (e.g., root/rhizome exudates). Microbes in the aerobic surface sediment breakdown a large fraction of this carbon and ultimately convert it into dissolved inorganic carbon, i.e., CO 2 , HCO 3 − and CO 3 2− (Fenchel et al., 1998), termed microbial carbon remineralization. However, in the absence of O 2 introduction at depth by plant roots (termed radial O 2 loss) or burrowing animals, transport of gasses and solutes across the surface of cohesive sediments occurs through the diffusive boundary layer (DBL), a thin water layer in which molecular diffusion is the only means of transport (Jørgensen and Revsbech, 1985). Microbial respiration strongly affects the supply of O 2 to the sediment, particularly in finegrained sediment and low-energy environments where advective O 2 supply is less important (Burdige, 2011). Consequently, high rates of microbial activity in the sediment surface deplete O 2 within a few millimeters of the sediment surface (Jørgensen and Revsbech, 1985), where after the sediment become anoxic and highly reduced (Canfield et al., 1993). Respiratory products of anaerobic metabolism in the anoxic part of the sediment include sulfide (S 2− ), pyrite (FeS 2 ) or diagenetically produced organic sulfur, dinitrogen (N 2 ), ammonium (NH 4 + ), methane (CH 4 ), and other products (Burdige, 2011).
The characteristic physico-chemical micro-gradients of the uppermost sediment layers determines the metabolic capacity of the microbes inhabiting vegetated coastal sediments, which in turn can control the fate of carbon within the sediments (other factors likely controlling carbon sequestration include climate, bathymetry and flow regime; see Belshe et al., 2017). The presence or absence of O 2 directly determines which metabolic pathways and enzymes can be used for carbon degradation (Heider and Fuchs, 1997). In oxic environments, large recalcitrant molecules, such as lignin polymers, are degraded through hydrolysis whereby organic molecules are broken down to smaller, more readily-degradable molecules facilitated by a group of enzymes named hydrolytic mono-or di-oxygenases (Heider and Fuchs, 1997). In anoxic sediments, degradation of complex recalcitrant molecules, such as lignocellulose, usually occurs at much lower rates through various reductive processes (Evans, 1977;Benner et al., 1984;Heider and Fuchs, 1997;Boll et al., 2002). Oxygen exposure therefore strongly affects how much recalcitrant carbon gets degraded (Hartnett et al., 1998). In the thin oxic surface layer of the sediment, aerobic heterotrophs gain energy by oxidizing carbon with O 2 as the electron acceptor. Below this oxic layer, diverse microbial populations use a variety of alternative electron acceptors for catabolic reactions, such as NO 3 − , Mn 4+ , Fe 3+ , and SO 4 2− (Fenchel et al., 1998). Sulfate reduction is generally the dominant microbial process in anoxic, marine sediments (Jørgensen, 1982), where H 2 S is produced using SO 4 2− as an electron acceptor, and can thus be used as a proxy for the overall anaerobic microbial activity and carbon turnover rate in the sediment (Jørgensen, 1982;Canfield et al., 1993). This is also true for VCHs, such as seagrass sediments, where sulfate-reducing bacteria (SRB) are often the most abundant microbial functional group (Welsh, 2000;Jensen et al., 2007;Cúcio et al., 2016). The activity of SRB at and around the oxic/anoxic interfaces, such as in the rhizosphere of plants, results in steep gradients of phytotoxic H 2 S below the oxidized zone of the sediment (e.g., Brodersen et al., 2014;Schrameyer et al., 2018).
Many wetland plants, including seagrasses, salt marsh and mangrove plants, release DOC (Moriarty et al., 1986;Pollard and Moriarty, 1991) and O 2 into the rhizosphere from roots and the rhizome (Pedersen et al., 1998;Koren et al., 2015;Brodersen et al., 2015bBrodersen et al., , 2016Koop-Jakobsen et al., 2018). Exudate DOC can be easily degraded by microbes both within the oxic and anoxic sediments and thus lead to local stimulation of the microbial community (Isaksen and Finster, 1996;Hansen et al., 2000;Nielsen et al., 2001;Alongi, 2005;Lovell, 2005), which results in high microbial respiration rates in the surface sediment (and plant rhizosphere) and thus high sediment O 2 demand. In this way, VCHs can modulate their surrounding below-ground environment to solubilize nutrients (Alongi, 2005;Lovell, 2005;Long et al., 2008;Pagès et al., 2011;Brodersen et al., 2017b) and oxidize SRB-formed phytotoxic sulfide (Brodersen et al., 2015b(Brodersen et al., , 2017b. However, a cross-ecosystem comparison of the key biogeochemical processes and chemical conditions within sediments of VCHs in the context of carbon remineralization and sequestration, and how the sediment processes may be affected by physical disturbance, are rare despite being needed to understand organic carbon sequestration and cycling in VCHs (Belshe et al., 2017). Physical disturbance of VCHs can lead to a reduction in sediment organic carbon stocks (Bourque et al., 2015), likely via O 2 exposure of sub-surface sediment layers and the subsequent stimulation of microbial activity (Macreadie et al., 2015;Trevathan-Tackett et al., 2018a,b).
In this study, we measured O 2 consumption and sulfate reduction rates (i.e., estimates of the aerobic and anaerobic microbial activity) within sediments of seagrass, salt marsh and mangrove ecosystems from the Sydney area, New South Wales, Australia in order to characterize the chemical conditions that could ultimately influence the efficiency of carbon sequestration. Sediment samples were divided into "surface" cores and "deep" cores to compare the activity of the microbial community within the two different sediment layers after experimental exposure to O 2 , approximately simulating undisturbed and physically exposed sediments, respectively. We hypothesized that surface sediments will have higher microbial activity, and thereby higher rates of O 2 consumption and sulfate reduction, than deeper sediment layers due to differences in the availability of readily degradable organic carbon (Costa et al., 2007). In contrast, "deep" (previously anoxic) sediment layers will be susceptible to O 2 intrusion and thus dispersal and settlement of facultative (or obligate) aerobes; potentially increasing remineralization of deep VCH carbon stocks. Detailed information on key biogeochemical processes within surface and deep and newly exposed sediment layers is important for predicting the potential carbon sequestration capacity for the ecosystems in question and can provide means for better guidance of coastal management and protection Lovelock et al., 2017;Trevathan-Tackett et al., 2017).

Study Site and Sediment Core Sampling
Sediments from VCHs were collected from the sub-tropical east coast of New South Wales near Sydney, Australia in the summer in February 2014. Mangrove and salt marsh sediment cores were collected at Towra Point in Botany Bay, NSW, Australia (34.0064 • S, 151.1653 • E). Avicennia marina was the dominant species at the mangrove site, whereas the salt marsh site was a mixture of Juncus kraussii and Sarcocornia quinqueflora.
Towra Point consists of a sandy peninsular largely isolated from terrestrial influence, with a salinity of 32 ± 3 (mean ± SD) in nearby open water (Geoscience Australia, 2015). Seagrass sediment cores were collected at Fagans Bay, NSW, Australia (33.4345 • S,151.3223 • E), where the sole seagrass species was Zostera muelleri. Fagans Bay sits at the mouth of Narara Creek and typically has salinities ranging from 25 to 30, however, it is subject to salinity fluctuations during flooding or drought periods (Collier and Mackenzie, 2008). All sediment cores where collected intertidally at low tide (tidal range from ∼0.2 to 1.2 m) near the rhizomes and roots and no belowground plant material was removed. Sediments were collected in triplicate with openbarrel push cores made of PVC pipes (200 cm length, 5 cm internal diameter), which were pushed into the sediments until refusal was reached (Figure 1). Expanding rubber plugs were then inserted into the top-end of the PVC pipes to achieve internal suction while retrieving the sediment cores. Sediment cores were then returned to the shore and sliced intact into smaller 10 cm sections with depth ("surface" = 0-10 cm depth core section; "deep" = a 10 cm interval from a depth >50 cm) using a reciprocating saw. All core sections were capped with PVC caps within 60 s of sectioning to reduce O 2 exposure to sediments and returned to the laboratory for incubation.

Experimental Design
Laboratory incubations were established in mesocosms (32 L, rectangular tanks) at the University of Technology Sydney (UTS), Sydney, Australia following Macreadie et al., (unpublished Figure 1). Surface and deep sediment cores were incubated in the dark under ∼26 • C (resembling water temperatures around Sydney during summer-time), a salinity of ∼35 and 100% air saturation in the water column in the large seawater reservoirs for 90 days (flushed and refilled with fresh seawater twice daily and continuously monitored; Macreadie et al., unpublished). The deep cores were incubated this way to simulate physical exposure as a result of disturbance (e.g., a dredging operation and/or dragging event) and to expose the new sediment surface layers to oxic conditions for a longer time period to allow for potential dispersal and settlement of facultative aerobes. While exposing the sediments to different tidal regimes would be more representative of natural inundation conditions, simulated tidal conditions would preclude us from making cross-habitat comparisons on the effect of O 2 exposure; therefore, all sediment treatments were standardized to full inundation. Water column O 2 status was maintained by air bubbling using an air stone connected to an air pump, and continued water flow was ensured using water pumps.

Sediment Organic Carbon Content
Measurement of organic carbon content (C org ) within sediments was done following procedures described by Campbell et al. (2014). In brief, after the 90 days mesocosm incubation bulk sediment samples (10 cm thick) were dried at 60 • C to a constant weight. Cores were subsampled at 1 cm intervals to 2 cm depth and dried. Dry bulk density (g dry weight/volume) was calculated before the sediments were ground using a ball mill and split into two equal parts. A Costech 410 Elemental Analyzer was used to measure the total sediment carbon (TC sediment ) of the first set of duplicates. The second set of duplicates was ashed in a muffled furnace at 500 • C for 6 h before analysis on the elemental analyser to determine the inorganic carbon content of the ash (IC ash ). The inorganic carbon content of the original sediment (IC sediment ) was calculated using the IC ash value back-calculated from the original weight of the un-ashed sample. Finally, we calculated C org (% dry weight; mean ± SEM; n = 6) as the following:

Microsensor Profiling
Microsensor profiling characterized the chemical environment within the top 0-3 cm of three replicate "surface" and "deep" sediment cores (Figure 1). Prior to profiling, each sediment core was submerged in a small (2 L) aquarium with seawater, which itself was positioned in a water bath preheated to ∼26 • C via a submersible aquarium heater (Sonpar Aquarium Equipment, China). Water aeration in the aquarium was ensured and maintained via a continuous flow of air from a submerged air-stone connected to an air pump. Water movement was achieved with a small water pump attached to the aquarium wall and adjusted to a height that allowed for sufficient water flow across the sediment surface (i.e., continuous movement of organic particles, but no movement of larger sand particles). Each sediment core was incubated in the aquarium for at least 20 min at constant flow after moving from the incubation tank (Figure 1e) to the microprofiling station (Figure 1f) in order to settle before profiling commenced to ensure steady state profiles in the sediment (as confirmed by repeated microsensor profiles).
Oxygen (O 2 ), hydrogen sulfide (H 2 S) and pH microprofiles were measured in individual sediment cores with Clark-type O 2 microsensors (Revsbech, 1989) and H 2 S microsensors (Jeroschewski et al., 1996;Kühl et al., 1998) (tip diameter = 100 µm, 90% response time <8 and 6 s, respectively, stirring sensitivity <1%; Unisense A/S, Aarhus, Denmark) and a pH microelectrode (e.g., Kühl and Revsbech, 2001) (tip diameter = 500 µm, 90% response time <10 s; Unisense A/S, Denmark; linear range pH 4-9; signal to pH ratio >50 mV/pH unit). The pH microelectrode was used in combination with an external standard 2 mm reference electrode (Ionode LLC, Australia) immersed in the aquarium seawater. All three sensor types were connected to a 4-channel microsensor multimeter (Unisense A/S, Denmark), interfaced with a laptop computer running dedicated data acquisition and microsensor positioning software (SensorTrace Pro v.3.1.1, Unisense A/S, Denmark). Microsensors were mounted on a motorized micromanipulator (Unisense A/S, Denmark), one sensor at a time, and the micromanipulator itself was mounted on a heavy metal stand to avoid vibrations. Each microsensor was linearly calibrated from signal readings (2-to 4-point calibrations) at experimental temperature and salinity before microprofiling commenced (full calibration procedures are provided on the manufactures website 1 ). Briefly, the O 2 microsensors were calibrated in 0% (obtained by flushing seawater with dinitrogen gas and/or by adding an oxygen scavenger, such as sodium dithionite) and 100% air-saturated seawater. The H 2 S microsensors were calibrated in anoxic, acidic (pH 4) H 2 S solutions (i.e., 0, 50, 100, 250 µM H 2 S). The pH microelectrodes were calibrated in buffers with a known pH of 4, 7, and 10. Prior to calibration, the Clark-type O 2 microsensors were pre-contaminated with sulfide to avoid drifting signals during measurements (Brodersen et al., 2015a).
For each depth profile, the tip of the respective microsensor was positioned at the sediment surface (defined as 0 µm depth on graphs) by means of a manual micromanipulator, while observing the sensor tip and the sediment surface through a stereo-microscope (AmScope, Irvine, CA, United States) mounted on an articulating arm (SM-6TZ; Amscope, Irvine, CA, United States). Measurements were done starting at a set distance from the sediment surface and then measured in vertical depth increments of 100 µm for O 2 profiles, and 500 or 1000 µm for H 2 S and pH profiles, respectively, down to a maximum vertical depth of ∼25 mm depending on the sediment O 2 consumption and the depth of the sulfide front, controlled through the position and data acquisition software (SensorTrace Pro; Unisense A/S). Three replicate profiles were done with each microsensor at various random locations in each sediment core while avoiding the very edges to limit edge effects and areas of complex flow. Finally, the three depth profiles within each sediment core were averaged to produce one replicate depth profile for each replicate sediment core per microsensor (i.e., n = 2-3, sediment core replicates; n = 6-9, microsensor profile replicates; note, one sediment core replicate was excluded from the seagrass-inhabited "surface" sediment, because the sediment core had fallen over during incubation).

Flux Calculations
Net O 2 and total sulfide (S tot 2− ) fluxes were calculated via Fick's first law of diffusion based on the slope of the linear change in 1 www.unisense.com concentration with depth: where J is the flux (nmol cm −2 s −1 ) of the compound, (dimensionless) is the porosity of the medium (here seawater or sediment), D 0 (cm 2 s −1 ) is the molecular diffusion coefficient for the respective molecule in seawater at the experimental temperature and salinity, and ∂C (nmol cm −3 ) is the linear change in concentration over the distance ∂z (cm). Specifically, net O 2 fluxes [termed diffusive O 2 uptake (DOU)] were calculated from the slope of the linear change in O 2 concentration in the diffusive boundary layer just above the sediment surface (Jørgensen and Revsbech, 1985), with = 1 (water) and D 0,O2 = 2.29·10 −5 cm 2 s −1 at a temperature of 26 • C and a salinity of 35 (values obtained from Unisense gas table; Ramsing and Gundersen, 1994 "Seawater and gasses-Tabulated physical parameters of interest to people working with microsensors in marine systems" 1 ). The net fluxes of total sulfide were calculated from the bottom linear segment of the sulfide concentration profile in the sediment, with = 0.65 (mangrove), 0.60 (salt marsh), 0.51 (seagrass) in surface sediment cores and 0.59 (mangrove), 0.58 (salt marsh), 0.72 (seagrass) in deep sediment cores, and the apparent diffusivity of H 2 S molecules calculated as D 0,H2S = 0.7573 × D 0,O2 = 1.73 × 10 −5 cm 2 s −1 (conversion factor to calculate the diffusivity of H 2 S was taken from the Unisense gas table). As the H 2 S microsensor is only sensitive to the gaseous H 2 S component of the total sulfide pool (defined as S tot 2− = [H 2 S] + [HS − ] + [S 2− ]), the speciation of which is directly related to the porewater pH (Jeroschewski et al., 1996); where dissolved H 2 S gas predominates at porewater pH levels below 7 and HS − predominates at higher pH levels. The total sulfide concentration was calculated for each H 2 S measurement as described in Jeroschewski et al. (1996) based on the pH measured in the same vertical depth and region within the sediment. We only used the upward flux of total sulfide to calculate the estimated sulfate reduction rates inside the sediments (Kühl and Jørgensen, 1992); as the sediment cores were sealed at the bottom the downward flux of sulfide can be assumed insignificant. The average volume-specific O 2 consumption rate (R) was calculated by dividing the DOU with the measured O 2 penetration depth, i.e., R = DOU/O 2 penetration depth (assuming constant reactivity and diffusivity over this depth interval; Frederiksen and Glud, 2006). Calculations were done based on the averaged profiles (n = 3) from each sediment core, and the presented fluxes are the average (± 1 SE, n = 3) of the three replicate sediment cores. The data for the seagrass microprofiles at day 90 have been shared with another paper where they are used in a different context and only included in the supplementary material (Supplementary Figure S6 in Macreadie et al., unpublished).

Data Analysis
Two-Way Analysis of Variance (ANOVA) tests in conjunction with Tukey's post hoc tests were performed to assess the differences in O 2 and total sulfide flux and O 2 penetration depth between ecosystem type and sediment depth. Data were log transformed where necessary to meet the assumption of equal variances. In cases where transformation did not improve variance, one-way ANOVAs or independent t-tests were performed separately for ecosystem and depth factors.

Percent Organic Carbon in Surface and Deep Sediments
The percent organic carbon content (%C org ) of the mangrove, salt marsh and seagrass surface sediments were 4.50 ± 0.96, 1.61 ± 0.52, and 1.25 ± 0.26%, which was significantly higher than in the deep sediments (except for in seagrasses) where the %C org amounted to 1.91 ± 0.76, 0.34 ± 0.05, and 1.56 ± 0.09%, respectively (p < 0.05). Mangrove sediments had significantly higher %C org than seagrass and salt marsh sediments (p = 0.001). No significant difference was found between seagrass and salt marsh sediments.

Depth Microprofiles of O 2 , Total Sulfide and pH
O 2 was rapidly consumed in the mangrove and seagrass surface sediments (O 2 penetration depth of 0.8 ± 0.2 mm and 1.6 ± 0.4 mm, respectively), which was contrasted by the salt marsh sediments where O 2 penetration was ∼5 times deeper (vertical depth of 5.3 ± 1.1 mm) (Figure 2A and Table 1; depth x ecosystem interaction for O 2 penetration depth, F 2,11 = 6.482, p = 0.014; schematic diagram shown in Figure 2E). In the deep sediment, the O 2 penetration depth was very similar in all three ecosystems (∼3 mm) ( Figure 2B and Table 1); and ∼2-3 times deeper than in the surface sediment for the mangrove and seagrass sediment (Figures 2A,B).
The surface sediment sulfide production was highest in the mangrove samples, indicated by a significantly shallower sulfide front than the other ecosystem sediments (0.7 ± 0.2 mm; habitat main effect, F 2,16 = 10.979, p = 0.001) and high maximum total sulfide concentration (1000 ± 307 µmol S tot 2− L −1 ) ( Figure 2C and Table 1). Salt marsh surface sediments had the lowest sulfide production (sulfide front of 7.3 ± 1.6 mm, maximum total sulfide concentration of 23 ± 10 µmol S tot 2− L −1 ), whereas sulfide production in surface seagrass sediments was intermediate (variable sulfide front of 3 and 19 mm, maximum total sulfide concentration of 837 µmol S tot 2− L −1 ) ( Figure 2C and Table 1). In contrast to the surface sediments, the deep sediments showed a shift in the sulfide front and large reduction in maximum total sulfide concentrations for mangrove and seagrass sediments (3.2 ± 1.2 mm and 310 ± 109 µmol S tot 2− L −1 , and 6.3 ± 1.7 mm and 81 µmol S tot 2− L −1 , respectively; Figure 2D and Table 1). The sulfide production for deep salt marsh sediments was similar to surface sediments, with the sulfide front ∼3 mm deeper than in the surface sediments and a similar maximum total sulfide concentration (33 ± 29 µmol S tot 2− L −1 ; Figures 2C,D).
Depth microprofiles of pH showed a pronounced decrease in pH with depth in all ecosystems (Figure 3). In mangrove surface sediments, pH dropped from 8.2 at the surface to 6.5 at 4.5 mm depth ( Figure 3A). Similar sediment acidification was found in the deeper mangrove sediments with pH reaching 6.9 at 8.5 mm depth ( Figure 3B). In salt marsh sediments, the sediment acidification with depth was very similar between surface and deep sediment with pH reaching 6.7 at 17 mm depth and 6.6 at 20 mm depth, respectively (Figure 3). In seagrass surface sediments, pH reached 7.5 at 8 mm depth where after a small increase in pH with depth was observed and pH 7.1 at 25 mm depth in deep sediment (Figure 3).

Areal and Volumetric O 2 Consumption Rates
O 2 consumption across the sediment-water interface was significantly higher in the surface sediment as compared to the deep sediment in seagrass and mangrove sediments, as shown by 1.8-5.5-fold higher diffusive O 2 uptake (DOU) for surface sediments (Figure 4A). In contrast, no significant change was found between depths for salt marsh sediments ( Figure 4A). Comparing across ecosystem types, the DOU of the mangrove surface sediment (102 ± 17 mmol O 2 m −2 d −1 ) was significantly higher than both salt marsh (22 ± 7 mmol O 2 m −2 d −1 ) and seagrass (34 mmol O 2 m −2 d −1 ) sediments ( Figure 4A). In the deeper sediment, the DOU was very similar across all VCHs, approximately 20 mmol O 2 m −2 d −1 (Figure 4A).
The depth-independent volume-specific O 2 consumption rate (R) also showed a significantly higher activity at the sediment surface and at the oxic/anoxic interface below in the mangrove sediments (122 ± 21 mol O 2 m −3 d −1 ) as compared to the salt marsh (4 ± 1 mol O 2 m −3 d −1 ) and seagrass (22 mol O 2 m −3 d −1 ) sediments ( Figure 4B). As for the DOU, the volumespecific O 2 consumption rates of the deep sediment were similar across sediments from the different VCHs and amounted to approximately 7 mol O 2 m −3 d −1 (Figure 4B).

Estimated Sulfate Reduction Rates
The estimated sulfate reduction rates (SRR) were significantly different among ecosystems ( Figure 4C). Specifically, the SRR activity for mangrove surface sediment was 11-400 times higher and deep sediment 29-94 times higher than seagrass and salt marsh sediments, respectively ( Figure 4C). In salt marsh sediments SRR was negligible, for both depths, with maximum rates of ∼0.2 ± 0.2 mmol S tot 2− m −2 d −1 (Figure 4C). In contrast, SRR was higher in the surface sediments compared to deep sediments in both mangrove and seagrass sediments (∼3-9-fold difference; Figure 4C). The averaged estimated sulfate reduction rate of the mangrove surface sediment was 57 ± 9 mmol S tot 2− m −2 d −1 as compared to 5 mmol S tot 2− m −2 d −1 in seagrass surface sediment (Figure 4C). In the deep sediment, the estimated sulfate reduction rate amounted to 17 ± 8 mmol S tot 2− m −2 d −1 and 1 ± 0 mmol S tot 2− m −2 d −1 in mangrove and seagrass sediment, respectively ( Figure 4C). >50 cm depth). Microsensor profiling was performed within the top ∼0-2 cm depth of each "surface" and "deep" sediment core, depending on the sulfide concentration. The total sulfide concentration was calculated from the measured H 2 S and pH microprofiles, as shown previously by Kühl and Jørgensen (1992).  Surface, sediment surface (0-10 cm depth); deep, deeper sediment layers (10 cm core section from >50 cm depth). Values represent mean ± standard error of the mean (SEM, n = 2-3).

DISCUSSION
This study aimed to identify how biogeochemical conditions varied across seagrass, mangrove and salt marsh sediment in the Sydney area (NSW, Australia), in order to assess how the different conditions may affect organic carbon (C org ) cycling. The controlled laboratory experimental design allowed us to standardize O 2 exposure via full inundation and continued water flow. However, the power to control these conditions limits several naturally occurring variables that could also influence the sediment biogeochemistry, including tidal inundation, flow regime, daily or seasonal temperature, water-column O 2 and light shifts and the effects of plant exudation, bioturbation and infauna. Taking into consideration the experimental design and the subsequent constraints of the data interpretation, our results provide valuable insight into how micro-scale biogeochemical processes and chemical conditions can affect carbon sequestration and storage. First, we show that mangrove and seagrass sediments had higher microbial activity in the surface sediments than the salt marsh sediments. Additionally,  Diffusive O 2 uptake rates (DOU; in mmol m −2 d −1 ) in surface ("surface"; 0-10 cm) and deeper ("deep"; >50 cm depth) sediment layers of sediments originating from mangrove, salt marsh, and seagrass ecosystems (Sydney area, NSW, Australia). Rates are calculated from the measured O 2 microprofiles. (B) Volume-specific O 2 consumption rates (R; in mmol m −3 d −1 ) in surface ("surface"; 0-10 cm depth) and deeper ("deep"; >50 cm depth) sediment layers of sediments originating from mangrove, salt marsh, and seagrass ecosystems (Sydney area, NSW, Australia). Rates are calculated from the measured O 2 microprofiles divided by the O 2 penetration depths. (C) Estimated sulfate reduction rates (in mmol m −2 d −1 ) in surface ("surface"; 0-10 cm depth) and deeper ("deep"; >50 cm depth) sediment layers of sediments originating from mangrove, salt marsh, and seagrass ecosystems (Sydney area, NSW, Australia). The estimated sulfate reduction rates are calculated based on the flux of total sulfide (S tot 2− ) calculated from the measured H 2 S and pH microprofiles, as shown previously by Kühl and Jørgensen (1992). Rates are mean values ± SEM; n = 6-9, technical replicates (biological replication of 2-3). Letters indicate statistical difference between ecosystem type (main effect).
the microbes in surface sediments seemed more active compared to deeper sediments in both seagrass and mangrove ecosystems. In salt marsh sediments, however, metabolic activity seemed similar among depths. As a result of the rapid microbial consumption of O 2 within the first few millimeters of the undisturbed surface sediments and the reduced microbial activity in the deeper sediments, anoxic acidic sediment conditions could be ideal for accumulating and storing carbon in these ecosystems. However, exposure of deep sediments to O 2 due to physical disturbance can stimulate aerobic metabolism, which could potentially reduce the carbon stocks of these exposed sediment layers.

O 2 , Sulfide and pH Distributions
O 2 was rapidly consumed within the first millimeters of the surface sediment of mangrove and seagrass ecosystems, and generally did not penetrate below 7 mm depth in any ecosystem. Sedimentary O 2 is consumed by microbial respiration and through chemical oxidation of reduced solutes that diffuse up from anoxic sediment layers below (Rasmussen and Jørgensen, 1992). The prevailing anoxic sediment conditions likely represents ideal chemical conditions for carbon burial as O 2 availability increases the degradation rate of recalcitrant carbon due to, e.g., the activation of aggressive mono-and dioxygenases enzymes in aerobic bacteria (Heider and Fuchs, 1997). This likely contributes to the high carbon sequestration rate observed in VCHs, where below-ground tissues are growing in largely anoxic sediment and therefore seldom exposed to oxic conditions during degradation (e.g., Frederiksen and Glud, 2006;Kristensen et al., 2008;Colmer et al., 2013;Brodersen et al., 2014Brodersen et al., , 2015aBrodersen et al., ,b, 2017bKoop-Jakobsen et al., 2017). These findings also correlate with the significantly higher percent organic carbon content (%C org ) found in the mangrove sediments compared to the seagrass and salt marsh sediments in this study. Hence, burial of carbon in the deeper sediment layers of VCHs is most likely a result of the aerobic respiration processes in the surface sediment occurring at a slower rate than the high inputs of organic matter to the system (carbon in > carbon out; see also, e.g., Kristensen et al., 2008 and references therein). The surface sediment of mangrove and seagrass ecosystems were characterized by high total sulfide concentrations that increased with depth below the oxic/anoxic interface and were similar to levels found in densely vegetated seagrass meadows (Borum et al., 2005). However, in some scenarios, the concentrations found in this study were close to reported phytotoxic levels that cause reduced aboveground production or death for VCHs (≥500 µmol L −1 for A. marina mangroves and Salicornia salt marsh species and >600 µmol L −1 for Zostera seagrass species; Lamers et al., 2013 and references herein). In naturally vegetated areas, the O 2 from the plant roots would help prevent sulfide intrusion into the plant, but this defense against phytotoxicity is regulated by the daylight hours, hypersalinity and the water-column O 2 availability (Pedersen et al., 2004;Brodersen et al., 2015bBrodersen et al., , 2017bKoren et al., 2015;Johnson et al., 2018). There was a pronounced reduction in the maximum total sulfide concentration in the deep sediment of mangrove and seagrass ecosystems by a factor of 4 to 24, respectively, which can be explained by less easily degradable carbon in the deeper sediment layers (Burdige, 2007;Costa et al., 2007). In the salt marsh sediments, the maximum total sulfide concentration only reached 33 µmol S tot 2− L −1 , which is far below the reported phytotoxic levels for salt marsh plants (≥500 µmol L −1 ; Lamers et al., 2013 and references herein) and below what generally is reported for salt marsh ecosystems (Howarth and Teal, 1979;Howarth and Giblin, 1983;Koch et al., 1990). It is therefore unlikely that the measured sediment sulfide levels would have led to reduced photosynthesis (i.e., lower DOC), growth and/or enhanced mortality rates [i.e., higher particulate organic carbon (POC)] potentially affecting the total organic matter input in any of the investigated coastal ecosystems.
We observed pronounced sediment acidification with depth in all investigated VCHs and sediment depths; with the lowest pH values (pH of <7) recorded in the mangrove and salt marsh sediments. This acidification of the sediment is likely due to O 2 -driven chemical oxidation of sediment-produced H 2 S, as the chemical reaction between O 2 and H 2 S generates protons (Nielsen et al., 2006) in combination with CO 2 release from microbial respiration (Soetaert et al., 2007). The low pH in the surface layers (down to pH of 6.5) just below the oxic/anoxic interface of the investigated mangrove and salt marsh sediments may inhibit the growth of some microbial taxa or groups, such as the bacterial phylum Bacteroidetes (Lauber et al., 2009;Thomas et al., 2011), with a decrease of up to ∼50% in relative abundance as compared to at a pH of 8.2 [based on Figure 5 in Lauber et al. (2009)]. Bacteroidetes is known to be important for the degradation of complex organic matter in marine environments, i.e., high molecular weight compounds, such as proteins and polysaccharides (Rosselló-Mora et al., 1999;Thomas et al., 2011). On the other hand, low pH conditions have been shown to benefit sulfate reducers, such as Desulfovibrio spp., in seagrass rhizospheres (Nielsen et al., 1999). In laboratory assays some sulfate reducers have the capacity to degrade recalcitrant organic matter (e.g., lignin; Goldfarb et al., 2011), however, isolates of Desulfovibrio from coastal sediments utilized primarily organic acids and some sugars (Nielsen et al., 1999). We therefore suggest that local sediment acidification can be an important driver of microbial community structure and function in vegetated coastal sediments, and thereby may also influence organic matter cycling within the sediments. However, such pH-induced effects on the microbial community composition in the sediments were beyond the scope of this paper and need to be further investigated.

O 2 Consumption and Estimated Sulfate Reduction
The higher O 2 consumption rates and estimated sulfate reduction rates in the mangrove and seagrass surface sediments as compared to salt marsh surface sediments, indicates higher overall microbial activity (maybe as a result of higher labile carbon availability), and thus results in a shallower oxidized sediment zone. The diffusive O 2 uptake and volume-specific O 2 consumption rates of the seagrass surface sediments (for direct comparison, R = 256 µmol m −3 s −1 ) were similar to O 2 consumption rates previously measured in seagrass meadows from the same region (Brodersen et al., 2017a). Estimated sulfate reduction rates and O 2 consumption were also higher in surface sediment than deeper sediment layers for both seagrass and mangrove ecosystems, indicating higher organic matter content in the surface sediments, which aligns with previous results (Saintilan et al., 2013). No significant sulfate reduction was found in salt marsh sediment, which suggests low overall microbial activity. However, salt marsh ecosystems have been generally reported to have high sulfate reduction rates (e.g., Howarth and Teal, 1979). We hypothesize that the low microbial activity, and the deviation to previous findings, could be due to a few causes. First, a higher fraction of recalcitrant or stable carbon could be deposited in the salt marsh sediments that we sampled compared to the other habitats, which could have reduced the ability of the microbial communities to utilize the sediment carbon (Unger et al., 2016;Leorri et al., 2018). Second, there could have been high precipitation rates of FeS below the oxic sediment zone (Thamdrup et al., 1994) that prevented H 2 S from accumulating in the sediment porewater. However, low H 2 S accumulation rates do not necessarily indicate lower sulfate reduction rates in the sediment. The precipitation of FeS is driven by H 2 S produced by sulfate reducers and the availability of Fe in the sediment, which is beyond what can be detected by the microsensor measurements in this study. Additionally, the low anaerobic metabolic rates in the salt marsh sediment could be due to their position in the upper intertidal zone, whereby the sediment surface is frequently exposed to air during low tide (i.e., molecular diffusion of O 2 is 10,000-fold faster in air as compared to water and 1 L of air holds 30-fold more O 2 than 1 L of water; Armstrong, 1979). The resulting large fluctuations in porewater salinity, temperature and O 2 penetration depth in the sediment are, in turn, likely to affect microbial metabolism.
Lastly, the sulfate reducers may be inhabiting sediment layers below the maximum microsensor profiling depth (i.e., down to ∼2 cm depth) used in this study and thus may not have been detected.

Potential Influence of Physical Disturbance on VCH Carbon Sequestration
High carbon accumulation rates due to particle trapping coupled with low to no occurrence of aerobic remineralization beyond the top mm of the sediment surface or the few mm surrounding the root tips (Pedersen et al., 1998;Frederiksen and Glud, 2006;Brodersen et al., 2014), supports the sequestration and preservation of sediment organic carbon. However, VCHs can be physically disturbed (e.g., boating activity, conversion, bioturbation by invertebrates) that would introduce O 2 to previously anoxic layers (e.g., Macreadie et al., 2015). Exposing deep sediments to the surface provided insight into how disturbances and subsequence exposure to O 2 could affect carbon stocks . Oxygen uptake of the deeper sediments indicated the activation of aerobes or facultative anaerobes present in the previously anoxic layers and/or transport and dispersal of facultative (or obligate) aerobes via the water column to the now exposed deeper sediment layers. Both aerobic and anaerobic metabolism in mangrove and seagrass sediments was generally lower for the deep sediments compared to surface sediments, suggesting either lower carbon substrate availability (Burdige, 2007) or lower enzyme utilization or efficiency (Arnosti, 2011). Regardless, O 2 exposure did elicit aerobic microbial metabolism, a response that may represent a significant risk to carbon that has accumulated on decadal to millennial timescales. (1) High sedimentation rate due to presence of plants capturing particles from the water column, increases burial of surface applied organic matter. (2) High input of organic material increases microbial activity at the sediment surface, ensuring low O 2 penetration. (3) High carbon accumulation in deep, anoxic sediment (cm-m). (4) Low resuspension of sediment due to reduced flow across sediment. (5) Sediment acidification with depth potentially leading to inhibition of bacterial growth. Hence, VCHs, such as mangrove forests, salt marshes and seagrass meadows, are characterized by three important environmental properties: (i) high organic loading, (ii) high carbon burial rate, and (iii) efficient protection of the buried carbon, as opposed to unvegetated sediment.
Our study used >50 cm carbon stock (i.e., a sediment core from >50 cm depth) and likely represents extreme case of disturbance, e.g., dredging and shrimp aquaculture. However, with the frequency of natural (bioturbation) and anthropogenic (habitat clearing, boating scars) physical disturbances in VCHs (Teal et al., 2008;Pendleton et al., 2012), these results highlight the potential risk of increased vulnerability for carbon stock losses via microbial remineralization as a result of O 2 exposure.

CONCLUSION
Our results showed that in the absence of biological O 2 introduction mangrove and seagrass surface sediments had higher rates of O 2 consumption and estimated sulfate reduction than salt marsh surface sediments in the Sydney area (NSW, Australia), with highest rates measured in the sediment surface of mangroves. We also found that surface of mangrove and seagrass sediments had higher O 2 consumption and estimated sulfate reduction rates than deeper sediments, which is indicative of higher microbial activity and higher carbon resources. In addition, depth microprofiles of pH and O 2 revealed pronounced sediment acidification with depth and low O 2 penetration in all VCH sediments. The highest percent organic carbon content was found in sediments originating from mangrove forests, which also exhibited lowest O 2 penetration depths and strongest sediment acidification with depth. The results from this case study suggests that the efficiency of the VCHs that we studied to store and bury carbon may in part be explained by (i) the prevailing anoxic conditions in the bulk sediment in concert with (ii) rapid sediment acidification with depth. Both of these conditions have the potential to favor slower anaerobic microbial remineralization processes, thus adding to other important factors (e.g., high organic matter accumulation and low sediment resuspension) that influence organic carbon accumulation in VCHs (Figure 5). Conversely in scenarios of exposure of deep sediments to the surface (e.g., as a result of a dredging and/or dragging events), the subsequence exposure to O 2 and likely activation of aerobic metabolism is hypothesized to be a mechanism for reducing carbon stocks in VCHs. Future research should thus be aiming at providing supporting evidence for many of the arguments made here, such as DIC measurements and the sources and quality of organic carbon, as well as determination of the iron content and/or iron speciation of the sediments (e.g., iron oxides, pyrite and FeS), as these geochemical variables can affect the calculations of the microbial activity determined from the detailed microsensor measurements (i.e., Belshe et al., 2017).