Exploring the mismatch between the theory and application of photosynthetic quotients in aquatic ecosystems

Abstract Estimates of primary productivity in aquatic ecosystems are commonly based on variation in O2, rather than CO2. The photosynthetic quotient (PQ) is used to convert primary production estimates from units of O2 to C. However, there is a mismatch between the theory and application of the PQ. Aquatic ecologists use PQ = 1–1.4. Meanwhile, PQ estimates from the literature support PQ = 0.1–4.2. Here, we describe the theory on why PQ may vary in aquatic ecosystems. We synthesize the current understanding of how processes such as NO3− assimilation and photorespiration can affect the PQ. We test these ideas with a case study of the Clark Fork River, Montana, where theory predicts that PQ could vary in space and time due to variation in environmental conditions. Finally, we highlight research needs to improve our understanding of the PQ. We suggest departing from fixed PQ values and instead use literature‐based sensitivity analyses to infer C dynamics from primary production estimated using O2

Methods of approximation influence aquatic ecosystem metabolism estimates

Aquatic ecologists have recently employed dynamic models to estimate aquatic ecosystem metabolism. All approaches involve numerically solving a differential equation describing dissolved oxygen (DO) dynamics. Although the DO differential equation can be solved accurately with linear multistep or Runge–Kutta methods, less accurate methods, such as the Euler method, have been applied. The methods also differ in how discrete temperature and light measurements are used to drive DO dynamics. Here, we used a representative stream DO data set to compare the metabolism estimates generated by multiple Euler based methods and an accurate numerical method. We also compared metabolism estimates using linear, piecewise constant and smoothing spline interpolation of light and temperature. Using observed DO to calculate DO saturation deficit in the Euler method results in a substantial difference in metabolism estimates compared to all other methods. If modeled DO is used to calculate DO saturation deficit, the Euler method introduces smaller error in metabolism estimates, which diminishes as logging interval decreases. Linear and smoothing spline interpolation result in similar metabolism estimates, but differ from estimates based on piecewise constant interpolation. We demonstrate how different computational methods imply distinct assumptions about process and observation error, and conclude that under the assumption of observation error, the best practice is to use the accurate numerical method of solving differential equation with a continuous interpolation of light and temperature. The Euler method will introduce minimal error if it is paired with frequently logged data and DO saturation deficit is computed using modeled DO

Variation in Detrital Resource Stoichiometry Signals Differential Carbon to Nutrient Limitation for Stream Consumers Across Biomes

Stoichiometric ratios of resources and consumers have been used to predict nutrient limitation across diverse terrestrial and aquatic ecosystems. In forested headwater streams, coarse and fine benthic organic matter (CBOM, FBOM) are primary basal resources for the food web, and the distribution and quality of these organic matter resources may therefore influence patterns of secondary production and nutrient cycling within stream networks or among biomes. We measured carbon (C), nitrogen (N), and phosphorus (P) content of CBOM and FBOM and calculated their stoichiometric ratios (C/N, C/P, N/P) from first- to fourth-order streams from tropical montane, temperate deciduous, and boreal forests, and tallgrass prairie, to compare the magnitude and variability of these resource types among biomes. We then used the ratios to predict nutritional limitations for consumers of each resource type. Across biomes, CBOM had consistently higher %C and %N, and higher and more variable C/N and C/P than FBOM, suggesting that microbial processing results in more tightly constrained elemental composition in FBOM than in CBOM. Biome-specific differences were observed in %P and N/P between the two resource pools; CBOM was lower in %P but higher in N/P than FBOM in the tropical montane and temperate deciduous forest biomes, while CBOM was higher in %P but similar in N/P than FBOM in the grassland and boreal forest biomes. Stable C-13 isotopes suggest that FBOM likely derives from CBOM in tropical and temperate deciduous forest, but that additional non-detrital components may contribute to FBOM in boreal forests and grasslands. Comparisons of stoichiometric ratios of CBOM and FBOM to estimated needs of aquatic detritivores suggest that shredders feeding on CBOM are more likely to experience nutrient (N and/or P) than C limitation, whereas collector-gatherers consuming FBOM are more likely to experience C than N and/or P limitation. Our results suggest that differences in basal resource elemental content and stoichiometric ratios have the potential to affect consumer production and ecosystem rates of C, N, and P cycling in relatively consistent ways across diverse biomes

Baseflow physical characteristics differ at multiple spatial scales in stream networks across diverse biomes

Context: Spatial scaling of ecological processes is facilitated by quantifying underlying habitat attributes. Physical and ecological patterns are often measured at disparate spatial scales limiting our ability to quantify ecological processes at broader spatial scales using physical attributes. Objective: We characterized variation of physical stream attributes during periods of high biological activity (i.e., baseflow) to match physical and ecological measurements and to identify the spatial scales exhibiting and predicting heterogeneity. Methods: We measured canopy cover, wetted width, water depth, and sediment size along transects of 1st–5th order reaches in five stream networks located in biomes from tropical forest to arctic tundra. We used hierarchical analysis of variance with three nested scales (watersheds, stream orders, reaches) to identify scales exhibiting significant heterogeneity in attributes and regression analyses to characterize gradients within and across stream networks. Results: Heterogeneity was evident at one or multiple spatial scales: canopy cover and water depth varied significantly at all three spatial scales while wetted width varied at two scales (stream order and reach) and sediment size remained largely unexplained. Similarly, prediction by drainage area depended on the attribute considered: depending on the watershed, increases in wetted width and water depth with drainage area were best fit with a linear, logarithmic, or power function. Variation in sediment size was independent of drainage area. Conclusions: The scaling of ecologically relevant baseflow physical characteristics will require study beyond the traditional bankfull geomorphology since predictions of baseflow physical attributes by drainage area were not always best explained by geomorphic power laws

Continental-scale decrease in net primary productivity in streams due to climate warming

Streams play a key role in the global carbon cycle. The balance between carbon intake through photosynthesis and carbon release via respiration influences carbon emissions from streams and depends on temperature. However, the lack of a comprehensive analysis of the temperature sensitivity of the metabolic balance in inland waters across latitudes and local climate conditions hinders an accurate projection of carbon emissions in a warmer future. Here, we use a model of diel dissolved oxygen dynamics, combined with high-frequency measurements of dissolved oxygen, light and temperature, to estimate the temperature sensitivities of gross primary production and ecosystem respiration in streams across six biomes, from the tropics to the arctic tundra. We find that the change in metabolic balance, that is, the ratio of gross primary production to ecosystem respiration, is a function of stream temperature and current metabolic balance. Applying this relationship to the global compilation of stream metabolism data, we find that a 1 °C increase in stream temperature leads to a convergence of metabolic balance and to a 23.6% overall decline in net ecosystem productivity across the streams studied. We suggest that if the relationship holds for similarly sized streams around the globe, the warming-induced shifts in metabolic balance will result in an increase of 0.0194 Pg carbon emitted from such streams every year