Abstract
An understanding of the modern relationship between diatom species and elevation is a prerequisite for using fossil diatoms to reconstruct relative sea level (RSL). We described modern diatom distributions from seven transects covering unvegetated subtidal environments to forested uplands from four tidal wetland sites (Smith Creek, Bone River, Niawiakum River, and Naselle River) of Willapa Bay, Washington, USA. We compared our diatom dataset (320 species from 104 samples) to a series of environmental variables (elevation, grain-size, total organic carbon (TOCSOM), and porewater salinity) using hierarchical clustering and ordination. While no single variable consistently explains variations in diatom assemblages at every site, elevation, salinity, and substrate (mud fraction and TOCSOM) are variables affecting diatom distribution along our transects. Elevation was the major environmental control of diatom variability (explained 27–39% variance) at four transects (Bone River Transect 1, Niawiakum River Transect 2, Naselle River Transect 1 and 2). Salinity and substrate were the major environmental controls (explained 12–32% variance) of diatom variability at three transects (Niawiakum River Transect 1: salinity; Smith Creek Transect and Bone River Transect 2: TOCSOM). Analyses of a combined regional dataset of all transects suggest that elevation is the driver of regional diatom variability in Willapa Bay, with salinity and substrate co-varying along an elevation gradient. We identify species with narrow elevation tolerances that can serve as indicator species of specific environments. Despite the site-specific variability of our modern diatom distribution, the regional dataset provides an important dataset that can be used to reconstruct RSL in Willapa Bay.
Keywords
Cascadia, Earthquakes, Relative sea level, subduction zone, Tidal wetland