Iris C. Anderson

Sandy sediment beaches covering 70% of non‐ice‐covered coastlines are important ecosystems for nutrient cycling along the land‐ocean continuum. Subterranean estuaries (STEs), where groundwater and seawater meet, are hotspots for biogeochemical cycling within sandy beaches. The STE microbial community facilitates biogeochemical reactions, determining the fate of nutrients, including nitrogen (N), supplied by groundwater. Nitrification influences the fate of N, oxidising reduced dissolved inorganic nitrogen (DIN), making it available for N removal. We used metabarcoding of 16S rRNA genes and quantitative PCR (qPCR) of ammonia monooxygenase ( amoA ) genes to characterise spatial and temporal variation in STE microbial community structure and nitrifying organisms. We examined nitrifier diversity, distribution and abundance to determine how geochemical measurements influenced their distribution in STEs. Sediment microbial communities varied with depth ( p ‐value = 0.001) and followed geochemical gradients in dissolved oxygen (DO), salinity, pH, dissolved inorganic carbon and DIN. Genetic potential for nitrification in the STE was evidenced by qPCR quantification of amoA genes. Ammonia oxidiser abundance was best explained by DIN, DO and pH. Our results suggest that geochemical gradients are tightly linked to STE community composition and nitrifier abundance, which are important to determine the fate and transport of groundwater‐derived nutrients to coastal waters.