Using Continuous Underway Isotope Measurements To Map Water Residence Time in Hydrodynamically Complex Tidal Environments

Exploring factors that affect Microcystis abundance in the sacramento san joaquin delta

Cyanobacteria harmful algal blooms (cHABs) are increasing in frequency, intensity and duration in estuaries worldwide. In the upper San Francisco Estuary, also known as the Sacramento San Joaquin Delta (Delta), cHABs have been a topic of concern over the past two decades. In response, managers are urgently working to understand the factors that drive cHABs and identify feasible management options to avert ecological and human health consequences. We used a six year data set to explore relationships between flow parameters, temperature, and Microcystis biovolume to determine the potential for managing large scale hydrodynamic conditions to address Delta cHABs. We also looked at the relationship between Microcystis biovolume and the low salinity zone to see if it could be used as a proxy for residence time, because residence time is positively related to cyanobacteria abundance. We found the low salinity zone is not a useful proxy for residence time in the area of the Delta that experiences the most severe cHABs. Our finding suggest that climatic conditions (i.e., temperature and water year type) have the greatest influence on Microcystis biovolume in the Delta, with higher biovolume during years with lower flow and higher temperatures. Further, there are interannual differences in Microcystis biovolume that cannot be fully explained by flow parameters or temperature, meaning other factors not included in our model may be involved. We conclude that management actions to increase flow may be ineffective at reducing Microcystis to desired levels if water temperatures remain high.

Stable isotopes provide insight into sources and cycling of N compounds in the Sacramento-San Joaquin Delta, California, USA

River deltas and their diverse array of aquatic environments are increasingly impacted by anthropogenic inputs of nitrogen (N). These inputs can alter the N biogeochemistry of these systems and promote undesirable phenomena including harmful algae blooms and invasive aquatic macrophytes. To examine N sources and biogeochemical processes in the Sacramento-San Joaquin Delta, a river delta located in central California, USA, that is fed primarily by the Sacramento River, we utilized a multi-tracer approach that measured N species concentrations and stable isotope values monthly from April 2011 to November 2012 in samples collected from the channelized mainstem of the Sacramento River, two channelized distributaries of the Sacramento River, and the Cache Slough Complex, a network of Sacramento River tributaries and shallow water wetland habitat. We found that the Sacramento River and its channelized distributaries received N primarily in the form of NH₄⁺ from treated wastewater effluent and that NH₄⁺ was lost rapidly while NO₃^- was gained more slowly during subsequent downstream transit, driven by an array of biogeochemical processes whose identities could be constrained via examination of stable isotope values. The Cache Slough Complex, which was characterized by lower net flows and higher water residence times than the Sacramento River and its distributaries, received variable inputs of low conductivity water elevated in NH₄⁺ from the Sacramento River and higher conductivity water elevated in NO₃^- from landward tributaries. Deviations from expected conservative mixing of these sources were spatially variable but broadly indicative of local inputs of treated wastewater effluent NO₃^-, conversion of Sacramento River NH₄⁺ to NO₃^- via nitrification, uptake of NH₄⁺ and NO₃^- by phytoplankton, and remineralization of organic N. These findings highlight both the diversity in N dynamics in anthropogenically impacted river delta environments and the utility of a multi-tracer approach in constraining these processes in such complex systems.

Timescale Methods for Simplifying, Understanding and Modeling Biophysical and Water Quality Processes in Coastal Aquatic Ecosystems: A Review

In this article, we describe the use of diagnostic timescales as simple tools for illuminating how aquatic ecosystems work, with a focus on coastal systems such as estuaries, lagoons, tidal rivers, reefs, deltas, gulfs, and continental shelves. Intending this as a tutorial as well as a review, we discuss relevant fundamental concepts (e.g., Lagrangian and Eulerian perspectives and methods, parcels, particles, and tracers), and describe many of the most commonly used diagnostic timescales and definitions. Citing field-based, model-based, and simple algebraic methods, we describe how physical timescales (e.g., residence time, flushing time, age, transit time) and biogeochemical timescales (e.g., for growth, decay, uptake, turnover, or consumption) are estimated and implemented (sometimes together) to illuminate coupled physical-biogeochemical systems. Multiple application examples are then provided to demonstrate how timescales have proven useful in simplifying, understanding, and modeling complex coastal aquatic systems. We discuss timescales from the perspective of “holism”, the degree of process richness incorporated into them, and the value of clarity in defining timescales used and in describing how they were estimated. Our objective is to provide context, new applications and methodological ideas and, for those new to timescale methods, a starting place for implementing them in their own work.

Spatial variability of phytoplankton in a shallow tidal freshwater system reveals complex controls on abundance and community structure

Estuaries worldwide are undergoing changes to patterns of aquatic productivity because of human activities that alter flow, impact sediment delivery and thus the light field, and contribute nutrients and contaminants like pesticides and metals. These changes can influence phytoplankton communities, which in turn can alter estuarine food webs. We used multiple approaches-including high-resolution water quality mapping, synoptic sampling, productivity and nitrogen uptake rates, Lagrangian parcel tracking, enclosure experiments and bottle incubations-over a short time period to take a "spatial snapshot" of conditions in the northern region of the San Francisco Estuary (California, USA) to examine how environmental drivers like light availability, nutrients, water residence time, and contaminants affect phytoplankton abundance and community attributes like size distribution, taxonomic structure, and nutrient uptake rates. Zones characterized by longer residence time (15-60 days) had higher chlorophyll-a concentrations (9 ± 4 µg L^-1) and were comprised primarily of small phytoplankton cells (<5 µm, 74 ± 8%), lower ammonium concentrations (1 ± 0.8 µM), higher nitrate uptake rates, and higher rates of potential carbon productivity. Conversely, zones characterized by shorter residence time (1-14 days) had higher ammonium concentration (13 ± 5 µM) and lower chlorophyll-a concentration (5 ± 1 µg L^-1) with diatoms making up a larger percent contribution. Longer residence time, however, did not result in the accumulation of large (>5 µm) cells considered important to pelagic food webs. Rather, longer residence time zones had a phytoplankton community comprised primarily of small cells, particularly picocyanobacteria that made up 38 ± 17% of the chlorophyll-a - nearly double the concentration seen in shorter residence time zones (22 ± 7% picocyanobacterial of chlorophyll-a). Our results suggest that water residence time in estuaries may have an effect as large or larger than that experimentally demonstrated for light, contaminants, or nutrients.

Complex life histories discovered in a critically endangered fish

Effective conservation of endangered species requires knowledge of the full range of life-history strategies used to maximize population resilience within a stochastic and ever-changing environment. California’s endemic Delta Smelt (Hypomesus transpacificus) is rapidly approaching extinction in the San Francisco Estuary, placing it in the crossfire between human and environmental uses of limited freshwater resources. Though managed as a semi-anadromous species, recent studies have challenged this lifecycle model for Delta Smelt, suggesting the species is an estuarine resident with several localized “hot-spots” of abundance. Using laser-ablation otolith strontium isotope microchemistry, we discovered three distinct life-history phenotypes including freshwater resident (FWR), brackish-water resident (BWR), and semi-anadromous (SA) fish. We further refined life-history phenotypes using an unsupervised algorithm and hierarchical clustering and found that in the last resilient year-class, the FWR (12%) and BWR (7%) comprised a small portion of the population, while the majority of fish were SA (81%). Furthermore, the semi-anadromous fish could be clustered into at least four additional life-history phenotypes that varied by natal origin, dispersal age and adult salinity history. These diverse life-history strategies should be incorporated into future conservation and management efforts aimed at preventing the extinction of Delta Smelt in the wild.

Characterizing macroinvertebrate community composition and abundance in freshwater tidal wetlands of the Sacramento-San Joaquin Delta

Restored tidal wetlands may provide important food web support for at-risk fish species in the Sacramento-San Joaquin Delta (Delta) of California, including Delta Smelt (Hypomesus transpacificus) and Chinook Salmon (Oncorhynchus tshawytscha). Since many tidal wetland restoration projects are planned or have recently been constructed in the Delta, understanding the diversity and variability of wetland invertebrates that are fish prey items is of increasing importance. During this study, two different invertebrate sampling techniques were tested (leaf packs and sweep nets) in four habitat types within three different wetland areas to evaluate which sampling technique provided the most reliable metric of invertebrate abundance and community composition. Sweep nets provided a better measure of fish food availability than leaf packs and were better able to differentiate between habitat types. Generalized linear models showed submerged and floating vegetation had higher abundance and taxa richness than channel habitats or emergent vegetation. Permutational multivariate analysis of variance showed significantly different communities of invertebrates in different habitat types and in different wetlands, and point-biserial correlation coefficients found a greater number of mobile taxa associated with sweep nets. There were more taxa associated with vegetated habitats than channel habitats, and one area had more taxa associated with it than the other two areas. These results suggest that restoration sites that contain multiple habitat types may enhance fish invertebrate prey diversity and resilience. However, the effect of habitat diversity must be monitored as restoration sites develop to assess actual benefits to at-risk fish species.

The Use of Stable Isotope-Based Water Age to Evaluate a Hydrodynamic Model

Transport time scales are common metrics of the strength of transport processes. Water age is the time elapsed since water from a specific source has entered a study area. An observational method to estimate water age relies on the progressive concentration of the heavier isotopes of hydrogen and oxygen in water that occurs during evaporation. The isotopic composition is used to derive the fraction of water evaporated, and then translated into a transport time scale by applying assumptions of representative water depth and evaporation rate. Water age can also be estimated by a hydrodynamic model using tracer transport equations. Water age calculated by each approach is compared in the Cache Slough Complex, located in the northern San Francisco Estuary, during summer conditions in which this region receives minimal direct freshwater inflow. The model’s representation of tidal dispersion of Sacramento River water into this backwater region is evaluated. In order to compare directly to isotopic estimates of the fraction of water evaporated (“fractional evaporation”) in addition to age, a hydrodynamic model-based property tracking approach analogous to the water age estimation approach is proposed. The age and fractional evaporation model results are analyzed to evaluate assumptions applied in the field-based age estimates. The generally good correspondence between the water age results from both approaches provides confidence in applying the modeling approach to predict age through broader spatial and temporal scales than are practical to assess using the field method, and discrepancies between the two methods suggest aspects of both approaches that may be improved. Model skill in predicting water age is compared to skill in predicting salinity. Compared to water age, salinity observations are shown to be a less useful diagnostic of transport in this low salinity region in which salt inputs are poorly constrained.

Calibrating temperature reconstructions from fish otolith oxygen isotope analysis for California's critically endangered Delta Smelt

RATIONALE

Oxygen isotope ratios (δ¹⁸ O values) of fish otoliths (ear bones) are valuable geochemical tracers of water conditions and thermal life history. Delta Smelt (Hypomesus transpacificus) are osmerid forage fish endemic to the San Francisco Estuary, California, USA, that are on the verge of extinction. These fish exhibit a complex life history that allows them to survive in a dynamic estuarine environment; however, a rapidly warming climate threatens this thermally sensitive species. Here we quantify the accuracy and precision of using δ¹⁸ O values in otoliths to reconstruct the thermal life histories of Delta Smelt.

METHODS

Delta Smelt were reared for 360 days using three different water sources with different ambient δ¹⁸ O_water values (-8.75‰, -5.28‰, and -4.06‰) and different water temperatures (16.4°C, 16.7°C, 18.7°C, and 20.5°C). Samples were collected after 170 days (n = 28) and 360 days (n = 14) post-hatch. In situ δ¹⁸ O values were measured from the core of the otolith to the dorsal edge using secondary ion mass spectrometry (SIMS) to reconstruct temporally resolved thermal life histories.

RESULTS

The δ¹⁸ O_otolith values for Delta Smelt varied as a linear inverse function of water temperature: 1000 ln α = 18.39 (±0.43, 1SE)(10³ TK^-1 ) - 34.56 (±1.49, 1SE) and δ¹⁸ O_{otolith(VPDB)}  - δ¹⁸ O_{water (VPDB)}  = 31.34(±0.09, 1SE) - 0.19(±0.01, 1SE) × T ° C. When the ambient δ¹⁸ O_water value is known, this species-specific temperature-dependent oxygen isotope fractionation model facilitated the accurate (0.25°C) and precise (±0.37°C, 2σ) reconstruction of the water temperature experienced by the fish. In contrast, the use of existing general fractionation equations resulted in inaccurate temperature reconstructions.

CONCLUSIONS

The species-specific δ¹⁸ O_otolith fractionation equation allowed for accurate and precise reconstructions of water temperatures experienced by Delta Smelt. Characterization of ambient δ¹⁸ O_water values remains a critical next step for reconstructing thermal life histories of wild Delta Smelt. This tool will provide new insights into habitat utilization, potential thermal refugia, and resilience to future warming for this critically endangered fish.

A novel high-frequency groundwater quality monitoring system

High-frequency, long-term monitoring of water quality has revolutionized the study of surface waters in recent years. However, application of these techniques to groundwater has been limited by the ability to remotely pump and analyze groundwater. This paper describes a novel autonomous groundwater quality monitoring system which samples multiple wells to evaluate temporal changes and identify trends in groundwater chemistry. The system, deployed near Fresno, California, USA, collects and transmits high-frequency data, including water temperature, specific conductance, pH, dissolved oxygen, and nitrate, from supply and monitoring wells, in real-time. The system consists of a water quality sonde and optical nitrate sensor, manifold, submersible three-phase pump, variable frequency drive, data collection platform, solar panels, and rechargeable battery bank. The manifold directs water from three wells to a single set of sensors, thereby reducing setup and operation costs associated with multi-sensor networks. Sampling multiple wells at high frequency for several years provided a means of monitoring the vertical distribution and transport of solutes in the aquifer. Initial results show short period variability of nitrate, specific conductivity, and dissolved oxygen in the shallow aquifer, while the deeper portion of the aquifer remains unchanged-observations that may be missed with traditional discrete sampling approaches. In this aquifer system, nitrate and specific conductance are increasing in the shallow aquifer, while invariant changes in deep groundwater chemistry likely reflect relatively slow groundwater flow. In contrast, systems with high groundwater velocity, such as karst aquifers, have been shown to exhibit higher-frequency groundwater chemistry changes. The stability of the deeper aquifer over the monitoring period was leveraged to develop estimates of measurement system uncertainty, which were typically lower than the manufacturer's stated specifications, enabling the identification of subtle variability in water chemistry that may have otherwise been missed.