Abstract
Subterranean estuaries (STEs) are important coastal biogeochemical reactors facilitating unique niches for microbial communities. A common approach in determining STE greenhouse gas and nutrient fluxes is to use terrestrial endmembers, not accounting for microbially mediated transformations throughout the STE. As such, the microbial ecology and spatial distribution of specialists that cycle compounds in STEs remain largely underexplored. In this study, we applied 16S rRNA amplicon sequencing with paired biogeochemical characterisations to spatially evaluate microbial communities transforming greenhouse gases and nutrients in an STE. We show that methanogens are most prevalent at the terrestrial end (up to 2.81% relative abundance) concomitant to the highest porewater methane, carbon dioxide and dissolved organic carbon concentrations (0.41 ± 0.02 μM, 273.31 ± 6.05 μM and 0.51 ± 0.02 mM, respectively). Lower ammonium concentrations corresponded with abundant nitrifying and ammonia‐oxidising prokaryotes in the mixing zone (up to 11.65% relative abundance). Methane, ammonium and dissolved organic carbon concentrations all decreased by >50% from the terrestrial to the oceanic end of the 15 m transect. This study highlights the STE's hidden microbiome zonation, as well as the importance of accounting for microbial transformations mitigating nutrient and greenhouse gas fluxes to the coastal ecosystems.