Riverine organic matter functional diversity increases with catchment size

A large amount of dissolved organic matter (DOM) is transported to the ocean from terrestrial inputs each year (~0.95 Pg C per year) and undergoes a series of abiotic and biotic reactions, causing a significant release of CO2. Combined, these reactions result in variable DOM characteristics (e.g., nominal oxidation state of carbon, double-bond equivalents, chemodiversity) which have demonstrated impacts on biogeochemistry and ecosystem function. Despite this importance, however, comparatively few studies focus on the drivers for DOM chemodiversity along a riverine continuum. Here, we characterized DOM within samples collected from a stream network in the Yakima River Basin using ultrahigh-resolution mass spectrometry (i.e., FTICR-MS). To link DOM chemistry to potential function, we identified putative biochemical transformations within each sample. We also used various molecular characteristics (e.g., thermodynamic favorability, degradability) to calculate a series of functional diversity metrics. We observed that the diversity of biochemical transformations increased with increasing upstream catchment area and landcover. This increase was also connected to expanding functional diversity of the molecular formula. This pattern suggests that as molecular formulas become more diverse in thermodynamics or degradability, there is increased opportunity for biochemical transformations, potentially creating a self-reinforcing cycle where transformations in turn increase diversity and diversity increase transformations. We also observed that these patterns are, in part, connected to landcover whereby the occurrence of many landcover types (e.g., agriculture, urban, forest, shrub) could expand DOM functional diversity. For example, we observed that a novel functional diversity metric measuring similarity to common freshwater molecular formulas (i.e., carboxyl-rich alicyclic molecules) was significantly related to urban coverage. These results show that DOM diversity does not decrease along stream networks, as predicted by a common conceptual model known as the River Continuum Concept, but rather are influenced by the thermodynamic and degradation potential of molecular formula within the DOM, as well as landcover patterns.