Biogeochemical cycling in an organic-rich coastal marine basin. VI - Temporal and spatial variations in sulfate reduction rates
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Acetate is important in anaerobic metabolism of non-vegetated sediments but its role in salt marsh soils was not investigated thoroughly. Acetate concentrations, oxidation (C-14) and SO4(2-) reduction (S-35) were measured in S. alterniflora soils in NH and MA. Pore water from cores contained greater than 0.1 mM acetate and in some instances greater than 1.0 mM. Non-destructive samples contained less than 0.01 mM. Acetate was associated with roots and concentrations were highest during vegetative growth and varied with changes in plant physiology. Acetate turnover was very low whether whole core or slurry incubations were used. Radiotracers injected directly into soils yielded rates of SO4(2-) reduction and acetate oxidation not significantly different from core incubation techniques. Regardless of incubation method, acetate oxidation did not account for a significant percentage of SO4(2-) reduction. These results differ markedly from data for non-vegetated coastal sediments where acetate levels are low, oxidation rate constants are high and acetate oxidation rates greatly exceed rates of SO4(2-) reduction. The discrepancy between rates of acetate oxidation and SO4(2-) reduction in marsh soils may be due either to the utilization of substrates other than acetate by SO4(2-) reducers or artifacts associated with measurements of organic utilization by rhizosphere bacteria.
We have combined several different methodologies to quantify rates of organic carbon mineralization by the various electron acceptors in sediments from the coast of Denmark and Norway. Rates of NH4+ and Sigma CO2 liberation sediment incubations were used with O2 penetration depths to conclude that O2 respiration accounted for only between 3.6-17.4% of the total organic carbon oxidation. Dentrification was limited to a narrow zone just below the depth of O2 penetration, and was not a major carbon oxidation pathway. The processes of Fe reduction, Mn reduction and sulfate reduction dominated organic carbon mineralization, but their relative significance varied depending on the sediment. Where high concentrations of Mn-oxide were found (3-4 wt% Mn), only Mn reduction occurred. With lower Mn oxide concentrations more typical of coastal sediments, Fe reduction and sulfate reduction were most important and of a similar magnitude. Overall, most of the measured O2 flux into the sediment was used to oxidized reduced inorganic species and not organic carbon. We suspect that the importance of O2 respiration in many coastal sediments has been overestimated, whereas metal oxide reduction (both Fe and Mn reduction) has probably been well underestimated.
Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.
To investigate the relationship between sigma(CO2) delta(C-13) values and rates of the dominant remineralization processes at the organic-rich field site of Cape Lookout Bight, NC, the isotopic composition of porewater sigma(CO2) was measured on a seasonal basis. The sigma(CO2) delta(C-13) values varies seasonally in response to changes in rates of sulfate reduction and methanogenesis, the dominant remineralization processes at this site. A tube incubation experiment was also performed to determine the isotopic signature of the sigma(CO2) produced by sulfate reduction and methanogenesis. The delta(C-13) of the sigma(CO2) produced in the sulfate reduction zone determined from the tube incubation was -14.3 plus or minus 1.9, a value enriched in C-13 relative to the labile organic fraction. The C-13-enrichment may be caused by low rates of methanogenesis occurring in the sulfate reduction zone. The delta(C-13) of the sigma(CO2) produced in the methanogenic zone was estimated to be +44 per mil, whereas the co-produced methane was -65 per mil. The fractionation factor for CO2 reduction was calculated to be 1.055, a value in agreement with previous estimates at this site. The measured concentration and delta(C-13) of the sigma(CO2) at Cape Lookout was closely reproduced by a diagenetic model using the measured rates of sulfate reduction and sigma(CO2) production, and the isotopic signature of the sigma(CO2) production in the two biogeochemical zones.
Acetate concentrations and rates of acetate oxidation and sulfate reduction were measured in S. alterniflora sediments in New Hampshire and Massachusetts. Pore water extracted from cores by squeezing or centrifugation contained in greater than 0.1 mM acetate and, in some instances, greater than 1.0 mM. Pore water sampled nondestructively contained much less acetate, often less than 0.01 mM. Acetate was associated with roots, and concentrations varied with changes in plant physiology. Acetate turnover was very low whether whole core or slurry incubations were used. Radiotracers injected directly into soils yielded rates of sulfate reduction and acetate oxidation not significantly different from core incubation techniques. Regardless of incubation method, acetate oxidation did not account for a substantial percentage of sulfate reduction. These results differ markedly from data for unvegetated coastal sediments where acetate levels are low, oxidation rate constants are high, and acetate oxication rates greatly exceed rates of sulfate reduction. The discrepancy between rates of acetate oxidation and sulfate reduction in these marsh soils may be due either to the utilization of substrates other than acetate by sulfate reducers or artifacts associated with measurements of organic utilization by rhizosphere bacteria. Care must be taken when interpreting data from salt marsh sediments since the release of material from roots during coring may affect the concentrations of certain compounds as well as influencing results obtained when sediment incubations are employed.
As marine sediments are buried, microbial communities transition from sulfate-reduction to methane-production after sulfate is depleted. When this biogenic methane diffuses into the overlying sulfate-rich sediments, it forms a sulfate-methane transition zone (SMTZ) because sulfate reducers deplete hydrogen concentrations and make hydrogenotrophic methanogenesis exergonic in the reverse direction, a process called the anaerobic oxidation of methane (AOM). Microbial participation in these processes is often inferred from geochemistry, genes, and gene expression changes with sediment depth, using sedimentation rates to convert depth to time. Less is known about how natural sediments transition through these geochemical states transition in real-time. We examined 16S rRNA gene amplicon libraries and metatranscriptomes in microcosms of anoxic sediment from the White Oak River estuary, NC, with three destructively sampled replicates with methane added (586-day incubations) and three re-sampled un-amended replicates (895-day incubations). Sulfate dropped to a low value (∼0.3 mM) on similar days for both experiments (312 and 320 days, respectively), followed by a peak in hydrogen, intermittent increases in methane-cycling archaea starting on days 375 and 362 (mostly Methanolinea spp. and Methanosaeta spp., and Methanococcoides sp. ANME-3), and a methane peak 1 month later. However, methane δ 13 C values only show net methanogenesis 6 months after methane-cycling archaea increase and 4 months after the methane peak, when sulfate is consistently below 0.1 mM and hydrogen increases to a stable 0.61 ± 0.13 nM (days 553–586, n = 9). Sulfate-reducing bacteria (mostly Desulfatiglans spp. and Desulfosarcina sp. SEEP-SRB1) increase in relative abundance only during this period of net methane production, suggesting syntrophy with methanogens in the absence of sulfate. The transition from sulfate reduction to methane production in marine sediments occurs through a prolonged period of methane-cycling by methanogens at low sulfate concentrations, and steady growth of sulfate reducers along with methanogens after sulfate is depleted.
Biological sulfate reduction (BSR) represents a promising strategy for bioremediation of sulfate-rich waste streams, yet the impact of metabolic interactions on performance is largely unexplored. Here, genome-resolved metagenomics was used to characterize 17 microbial communities in reactors treating synthetic sulfate-contaminated solutions. Reactors were supplemented with lactate or acetate and a small amount of fermentable substrate. Of the 163 genomes representing all the abundant bacteria, 130 encode 321 NiFe and FeFe hydrogenases and all genomes of the 22 sulfate-reducing microorganisms (SRM) encode genes for H 2 uptake. We observed lactate oxidation solely in the first packed bed reactor zone, with propionate and acetate oxidation in the middle and predominantly acetate oxidation in the effluent zone. The energetics of these reactions are very different, yet sulfate reduction kinetics were unaffected by the type of electron donor available. We hypothesize that the comparable rates, despite the typically slow growth of SRM on acetate, are a result of the consumption of H 2 generated by fermentation. This is supported by the sustained performance of a predominantly acetate-supplemented stirred tank reactor dominated by diverse fermentative bacteria encoding FeFe hydrogenase genes and SRM capable of acetate and hydrogen consumption and CO 2 assimilation. Thus, addition of fermentable substrates to stimulate syntrophic relationships may improve the performance of BSR reactors supplemented with inexpensive acetate.
Complete budgets for carbon and oxygen have been constructed for cyanobacterial mats dominated by Microcoleus chthonoplastes from the evaporating ponds of a salt works. We infer from the data the various sinks for O2 as well as the sources of carbon for primary production. Although seasonal variability exists, a major percentage of the O2 produced during the day did not diffuse out of the mat but was used within the mat to oxidize both organic carbon and the sulfide produced by sulfate reduction. At night, most of the O2 that diffused into the mat was used to oxidize sulfide, with O2 respiration of minor importance. During the day, the internal mat processes of sulfate reduction and O2 respiration generated as much or more inorganic carbon (DIC) for primary production as diffusion into the mat. Oxygenic photosynthesis was the most important process of carbon fixation. At night, the DIC lost from the mat was mostly from sulfate reduction. Elemental fluxes across the mat/brine interface indicated that carbon with an oxidation state of greater than zero was taken up by the mat during the day and liberated from the mat at night. Overall, carbon with an average oxidation state of near zero accumulated in the mat. Both carbon fixation and carbon oxidation rates varied with temperature by a similar amount.
Organic sulfur containing compounds represent greater than 80% of the total sulfur in sediments of eutrophic freshwater lakes. Although sedimentary sulfur is predominantly in the form of organic compounds, more sulfur is transformed by sulfate reduction than by any other process. Rates of sulfate reduction in these sediments average 7 mmol/sq m/day. This rate is 19 times greater than the net rate of production of inorganic sulfur from organic compounds on an annual basis.
Methylotrophic methanogenesis in the sulfate-rich zone of coastal and marine sediments couples with anaerobic oxidation of methane (AOM), forming the cryptic methane cycle. This study provides evidence of cryptic methane cycling in the sulfate-rich zone across a land–ocean transect of four stations–two brackish, one marine, and one hypersaline–within the Carpinteria Salt Marsh Reserve (CSMR), southern California, USA. Samples from the top 20 cm of sediment from the transect were analyzed through geochemical and molecular (16S rRNA) techniques, in-vitro methanogenesis incubations, and radiotracer incubations utilizing 35 S-SO 4 , 14 C-mono-methylamine, and 14 C-CH 4 . Sediment methane concentrations were consistently low (3 to 28 µM) at all stations, except for the marine station, where methane increased with depth reaching 665 µM. Methanogenesis from mono-methylamine was detected throughout the sediment at all stations with estimated CH 4 production rates in the sub-nanomolar to nanomolar range per cm 3 sediment and day. 16S rRNA analysis identified methanogenic archaea (Methanosarcinaceae, Methanomassiliicoccales, and Methanonatronarchaeacea) capable of producing methane from methylamines in sediment where methylotrophic methanogenesis was found to be active. Metabolomic analysis of porewater showed mono-methylamine was mostly undetectable (<3 µM) or present in trace amounts (<10 µM) suggesting rapid metabolic turnover. In-vitro methanogenesis incubations of natural sediment showed no linear methane buildup, suggesting a process limiting methane emissions. AOM activity, measured with 14 C-CH 4 , overlapped with methanogenesis from mono-methylamine activity at all stations, with rates ranging from 0.03 to 19.4 nmol cm −3 d −1 . Geochemical porewater analysis showed the CSMR sediments are rich in sulfate and iron. Porewater sulfate concentrations (9–91 mM) were non-limiting across the transect, supporting sulfate reduction activity (1.5–2,506 nmol cm −3 d −1 ). Porewater sulfide and iron (II) profiles indicated that the sediment transitioned from a predominantly iron-reducing environment at the two brackish stations to a predominantly sulfate-reducing environment at the marine and hypersaline stations, which coincided with the presence of phyla (Desulfobacterota) involved in these processes. AOM activity overlapped with sulfate reduction and porewater iron (II) concentrations suggesting that AOM is likely coupled to sulfate and possibly iron reduction at all stations. However, 16S rRNA analysis identified anaerobic methanotrophs (ANME-2) only at the marine and hypersaline stations while putative methanogens were found in sediment across all stations. In one sediment horizon at the marine station, methanogen families (Methanosarcinaceae, Methanosaetaceae, Methanomassiliicoccales, and Methanoregulaceae) and ANME 2a,2b, and 2c groups were found together. Collectively, our data suggest that at the brackish stations methanogens alone may be involved in cryptic methane cycling, while at the marine and hypersaline stations both groups may be involved in the process. Differences in rate constants from incubations with 14 C-labeled methane and mono-methylamine suggest a non-methanogenic process oxidizing mono-methylamine to inorganic carbon, likely mediated by sulfate-reducing bacteria. Understanding the potential competition of sulfate reducers with methanogens for mono-methylamine needs further investigation as it might be another important process responsible for low methane emissions in salt marshes.
Sulfur (S) is an essential macronutrient and important component of the earth’s crust, and its cycling has critical impacts on trace metal mobility, water quality, and human health. Pyrite weathering is the primary pathway by which sulfur enters surface waters. However, biogeochemical cycling of sulfur in soils and the river corridor mediates sulfate exports. In this study, we identified the major forms of sulfur across multiple compartments and scales in a pristine mountainous watershed, including shale bedrock weathering profiles, hillslope soils, and alluvial floodplain sediments, in order to provide insight into biogeochemical sulfur cycling in a hydrologically variable alpine system. X-ray absorption near-edge spectroscopy (XANES) analysis of shale weathering profiles showed clear evidence of pyrite oxidation to sulfate, with large accumulations of intermediate S(0) (20%–53%). Micro-scale XANES showed evidence of reprecipitation of pyrite at fracture surfaces within the permanently saturated zone. Organic sulfur dominated S speciation in shallow hillslope soil and floodplain sediment, with little evidence of reduced inorganic S. However, mackinawite formation, representing active sulfate reduction, was observed in saturated oxbow sediments and saturated weathered shale underlying floodplain sediments. Further evidence of sulfate reduction from aqueous sulfur isotopic analysis was observed in shallow groundwater transects across an Fe-reducing meander, whereas increases in pore water sulfate concentrations implied sulfur oxidation at other locations. Furthermore, the data present an integrated picture of sulfur cycling in a shale-dominated watershed, where riverine sulfate exports are mediated by biological cycling, particularly in redox-stratified and temporally dynamic hyporheic zone sediments.
Stromatolites offer an unparalleled geologic record of early life, because they constitute the oldest and most abundant recognizable remains of microbial ecosystems. Microbial mats are living homologs of stromatolites; thus, the physiology of the microbiota as well as the processes which create those features of mats (e.g., biomarker organic compounds, elemental and stable isotopic compositions) which are preserved in the ancient record. Observations of the carbon isotopic composition (delta C-13) of stromatolites and microbial mats were made and are consistent with the hypothesis that atmospheric CO2 concentrations have declined by at least one to two orders of magnitude during the past 2.5 Ga. Whereas delta C-13 values of carbonate carbon average about 0 permil during both the early and mid-Proterozoic, the delta C-13 values of stromatolitic organic matter increase from an average of -35 between 2.0 and 2.6 Ga ago to an average of about -28 about 1.0 Ga ago. Modern microbial mats in hypersaline environments have delta C-13 values typically in the range of -5 to -9, relative to an inorganic bicarbonate source at 0 permil. Both microbial mats and pur cultures of cyanobacteria grown in waters in near equilibrium with current atmospheric CO2 levels exhibit minimal discrimination against C-13. In contrast, hot spring cyanobacterial mats or cyanobacterial cultures grown under higher CO2 levels exhibit substantially greater discrimination. If care is taken to compare modern mats with stromatolites from comparable environments, it might be possible to estimate ancient levels of atmospheric CO2. In modern microbial mats, a tight coupling exists between photosynthetic organic carbon production and subsequent carbon oxidation, mostly by sulfate reduction. The rate of one process fuels a high rate of the other, with much of the sulfate reduction occurring within the same depth interval as oxygenic photosynthesis. Other aspects of this study are presented.
Summer sulfate reduction and sulfide and iron chemistry were monitored in marsh soils along a gradient from a creekside Spartina alterniflora region to an inland area dominated by S. patens. Sulfate reduction rates measured using S-35 increased from 0.4 in June to as much as 4.5 micro-mole/ml(sup -1)d(sup -1) in July in S. alterniflora soils with most rapid rates occurring in the upper few cm. Rates in S. patens soils were approximately equal to 5-8 fold slower with the most rapid rates occurring generally in soils deeper than 10 cm. The recovery percentage of reduced S-35 sulfur varied with depths at both locations; dissolved and acid-volatile sulfides dominating S. alterniflora regions while chromium-reducible solid phases were abundant in S. patens soils. Dissolved sulfide in S. alterniflora soils increased throughout the summer to approximately equal to 2.8 mM while sulfide in S. patens soils was abundant only in soils deeper than 15 cm. Dissolved sulfide covaried inversely with iron. Diel studies demonstrated that sulfide and iron varied approximately equal to two-fold in response to a semi-diurnal tide. Dissolved sulfide turnover was most rapid during periods of active plant growth.
Different methods for measuring the rates of processes mediated by bacteria in sediments and the rates of bacterial cell production have been compared. In addition, net production of the seagrass Zostera capricorni and bacterial production have been compared and some interrelationships with the nitrogen cycle discussed. Seagrass productivity was estimated by measuring the plastochrone interval using a leaf stapling technique. The average productivity over four seasons was 1.28 +/- 0.28 g C m-2 day-1 (mean +/- standard deviation, n = 4). Bacterial productivity was measured five times throughout a year using the rate of tritiated thymidine incorporated into DNA. Average values were 33 +/- 12 mg C m-2 day-1 for sediment and 23 +/- 4 for water column (n = 5). Spatial variability between samples was greater than seasonal variation for both seagrass productivity and bacterial productivity. On one occasion, bacterial productivity was measured using the rate of 32P incorporated into phospholipid. The values were comparable to those obtained with tritiated thymidine. The rate of sulfate reduction was 10 mmol SO4(-2) m-2 day-1. The rate of methanogenesis was low, being 5.6 mg CH4 produced m-2 day-1. A comparison of C flux measured using rates of sulfate reduction and DNA synthesis indicated that anaerobic processes were predominant in these sediments. An analysis of microbial biomass and community structure, using techniques of phospholipid analysis, showed that bacteria were predominant members of the microbial biomass and that of these, strictly anaerobic bacteria were the main components. Ammonia concentration in interstitial water varied from 23 to 71 micromoles. Estimates of the amount of ammonia required by seagrass showed that the ammonia would turn over about once per day. Rapid recycling of nitrogen by bacteria and bacterial grazers is probably important.
Bacterial sulfate reduction occurs in all marine sediments that contain organic matter. Aqueous sulfide (HS-, H2S), one of the initial products of bacterial sulfide reduction, is extremely reactive with iron bearing minerals: sulfur is fixed into sediments as iron sulfide (first FeS and then Fe2S2). A working definition is given of sulfur diagenesis in marine sediments. Controls and consequences of sulfate reduction rates in marine sediments are examined.
Long-term metal remediation in wetland treatment systems (WTSs) involves facilitating dissimilatory sulfate reduction to produce sulfide and mineralize metals in deep sediments. We evaluated seasonal sulfur cycling in two constructed wetlands (Maintained WTS constructed in 2007, and the Unmaintained WTS constructed in 2000) on the Savannah River Site in Aiken, South Carolina, USA. Significant interactions in sulfide concentration were observed between sediment depth, season, and wetland (F = 4.64, df = 11, P = 3.28 × 10 − 5). In the Maintained WTS, dissimilatory sulfate reduction dominated the surface sediments during the warm season (0–2 cm depth, t=-2.66, P = 9.70 × 10 − 3), unlike the Unmaintained system. Sulfate concentrations in pore waters increased in the warm season (F = 7.84, df = 1, P = 6.50 × 10 − 3), contrary to expectations. Sulfur limitation in the Unmaintained WTS during the warm season correlated with increased sulfur assimilation in giant bulrush. Lower sulfide concentrations in surface sediments of the Unmaintained WTS illustrated aging effects. The Maintained WTS shows potential for managing erosion, pH reduction, and sulfur limitation observed in the older Unmaintained WTS.
Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.