
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38561458
58246
10.1038/s41598-024-58246-w
Article
Clumped isotopes reveal relationship between mussel growth and river discharge
Brewer Melanie A. 1
Grossman Ethan L. e-grossman@tamu.edu

1
Randklev Charles R. 2
1 https://ror.org/01f5ytq51 grid.264756.4 0000 0004 4687 2082 Department of Geology and Geophysics, Texas A&M University, College Station, TX 77843 USA
2 grid.264756.4 0000 0004 4687 2082 Texas A&M Natural Resources Institute, Dallas, TX 75252 USA
1 4 2024
1 4 2024
2024
14 76534 12 2023
27 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Freshwater mussels preserve valuable information about hydrology, climate, and population dynamics, but developing seasonal chronologies can be problematic. Using clumped isotope thermometry, we produced high-resolution reconstructions of modern and historic (~ 1900) temperatures and δ18Owater from mussel shells collected from an impounded river, the Brazos in Texas, before and after damming. We also performed high-resolution growth band analyses to investigate relationships between mussel growth rate, rainfall, and seasonal temperature. Reconstructed δ18Owater and temperature vary little between the modern (3R5) and historic shell (H3R). However, a positive relationship between reconstructed δ18Owater and growth rate in H3R indicates that aside from diminished growth in winter, precipitation and flow rate are the strongest controls on mussel growth in both modern and pre-dam times. Overall, our results demonstrate (1) the impact, both positive and negative, of environmental factors such as flow alteration and temperature on mussel growth and (2) the potential for clumped isotopes in freshwater mussels as a paleohydrology and paleoclimate proxies in terrestrial environments.

Subject terms

Ecology
Environmental sciences
Hydrology
http://dx.doi.org/10.13039/100007904 Texas A and M University Michel T. Halbouty Chair Michel T. Halbouty Chair Colleges of Agriculture & Life Sciences and Arts & Sciences Brewer Melanie A. Grossman Ethan L. Randklev Charles R. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Freshwater mussels (Order: Unionida) are arguably the most threatened group of animals in the world1,2, with nearly 70% of mussels in North America alone considered extinct or imperiled, and the remaining mussels living within fragile environments in need of urgent protection and proactive conservation3. Threats from habitat destruction, over-exploitation, water quality degradation, invasive species, and climate change are adding immense pressure to mussel survival, reproduction, and dispersal. Habitat destruction can be caused by alteration of natural flow regimes resulting from river impoundment, stream channelization, increases in impervious surfaces, ground water withdrawals, and changes in land use practices4–6. These environmental and anthropogenic stressors are associated with the mussel declines seen in the past and present, with alteration to natural flow regime being one of the most important because it shapes channel morphology, regulates water temperature, and influences water quality and nutrient cycling, which taken together determine mussel habitat in rivers and streams7–10. Thus, changes in flow that affect any one of these factors can result in shifts or elimination of mussel habitat, which over time can drive changes in species occurrence and community composition4.

Unfortunately, the causal mechanisms underpinning these community level changes remain largely unknown. This is due in large part to the lack of information on how population performance (i.e., growth, survivorship, and reproduction) is affected by changes to flow and temperature. A notable exception is the study by Rypel et al.11, who looked at the effect of river impoundment on shell growth in select streams in the southeastern United States. The authors found that in regulated streams, mussel growth had become decoupled from hydrologic variation, but the long-term implication of their findings is unclear. Bolotov et al.12 evaluated the effects of climate change on margaritiferid mussels, a sister family to unionids, using shell morphology and growth rates, finding that a climate warming of ~ 6 °C in mean summer temperature over the past 100 years likely contributed to their decline. These findings demonstrate the utility of using shell growth as means to evaluate the relationship between mussel growth, flow rates, and climate change.

Using shells to estimate age, growth, and mussel-environment relationships is rooted in the idea that shell growth is slower during the winter and faster in spring and summer but can be dampened or accelerated by environmental stress and certain anthropogenic practices13,14. However, these growth patterns can vary from year to year, and generally show an inverse relationship with age where shell growth slows as the organism ages15,16. Within the shell, slower growth during the winter season has typically been recognized as dark “annual” bands that are prominent on both the shell exterior and within cross-section17,18. Yet, similar banding has also been shown to occur in other seasons as a response to disturbance (i.e., sudden change in temperature, entrainment during floods, handling by surveyors, etc.) and can lead to misinterpretation of seasonality, growth rate, and age in individuals4. Therefore, accurately determining these bands is critical for reconstructing and interpreting the ecological and environmental records within freshwater mussels that can live multiple decades to more than a century16,19,20. Annual bands in sclerochronology studies have been confirmed with the help of δ18O (18O/16O) thermometry of shells15,21–25. However, in many freshwater systems, especially in the impounded river systems commonly found throughout the U.S., decoupling of the temperature and water δ18O (δ18Owater) signals makes relating annual bands to cooler δ18O temperatures problematic26–28, rendering a second unambiguous method for determining sclerochronologies critical for reconstructing and interpreting the ecological and environmental records in the shells. One such method is clumped isotopes.

Clumped isotopes, the incorporation of two rare isotopes in a molecule, has become one of the most important developments in low-temperature geochemistry in the last two decades29. Clumped isotope analyses (Δ47) of carbonates examines the bonding, or “clumping”, between 13C and 18O isotopes in carbonate molecules and their deviation from the stochastic distribution. Clumping is enhanced by thermodynamic stability of the 13C–18O bond and decreases as temperature increases30. Unlike conventional oxygen isotope thermometry, where isotopic compositions are dependent on both temperature and δ18Owater31, clumped isotope thermometry is strictly temperature dependent, making it an ideal tool for estimating growth temperatures and establishing chronologies27,32. In addition, knowing the δ18O of the carbonate and Δ47 temperature, δ18Owater can be calculated from the rearrangement of the conventional oxygen isotope thermometry equation (e.g.,33,34).

The application of clumped isotopes has been utilized in numerous marine studies for reconstructing paleoclimate and δ18Owater in ancient oceans (e.g.,35–39, but only recently applied to fluvial environments and freshwater mussel shells, confirming their utility as paleohydrologic proxies27. In river systems, the reconstructed δ18Owater and Δ47 temperature provide critical information for understanding the cycling of water, precipitation, and evaporation. Compared with marine environments, rivers are generally more variable due to interruptions in hydrologic connectivity, variability of the input source, and greater seasonal temperature variation. Van Plantinga and Grossman27 found that Δ47 values produced reliable temperatures and reconstructed δ18Owater values that were within measured river values. However, the resolution was low because of the small number of analyses and large sample size requirements at the time.

In our study, we combine high-resolution clumped isotope analyses with near daily-resolution growth band analyses to produce temperature, δ18Owater, and growth rate records of monthly to bi-monthly resolution in Amblema plicata, common name Threeridge (hence “3R” in specimen ID). The threeridge mussel is a common mussel distributed in a broad range of rivers, streams, and lakes within central and eastern North America3. We compare our reconstructed δ18Owater values to measured river data from Van Plantinga and Grossman40 and local precipitation data (1) to demonstrate that clumped isotopes can serve as a high-resolution proxy for reconstructing local rainfall and drought events in terrestrial environments, and (2) to evaluate the role of seasonality and rainfall in controlling mollusk growth rates. This information is critical to better understanding river hydrology, climate, and mussel ecology in prehistoric, historic, and anthropogenic times.

Study area and climate setting

The study area is the main stem of the middle Brazos River in Texas (Fig. 1), where mussel shells were collected outside of College Station in ~ 1900 and 2013. The modern shell was collected live from the shoal waters of the Brazos River. The historic shell lacks detailed information about the collection site. However, both shells were collected at almost the same location, which is represented by the marker labelled “Shell Collection Site” in Fig. 1A.Figure 1 (A) Map of the study area showing the Bryan discharge and elevation monitoring gauge (purple), water sampling site (gray), and shell sampling site (orange). (B) Inset map showing the study area in the Brazos River, TX drainage basin in green shading, river path in solid blue line, and study area in red shading. Yellow points represent the 3 major dams along the Brazos: Possum Kingdom (PK); Lake Granbury (Gb); and Lake Whitney (LW) dams. Black “x” represents the sampling site for δ18O in rainwater during 1962–65 and 1972–76. Note: Maps were created in R v4.3.1 and Rstudio v2023.− 9.1 + 494 using the ggmap and tidyverse packages with the roadmap satellite image sourced from Google Maps https://www.google.com/maps/@30.5656804,-96.4597714,11.57z?entry=ttu and the map_data [in ggplot2] function with in-built datasets. Points and shading on maps were added using ggplot2 and Microsoft PowerPoint v16.81, respectively.

The Brazos River basin is the second largest watershed in Texas and spans the longest distance for a river located entirely in Texas41. The river itself is often referred to as the boundary between arid west Texas, where evaporation often exceeds precipitation, and the more humid east Texas that receives more rain. The region’s water is supplied by precipitation, groundwater, numerous smaller tributaries and creeks, and impoundments of upstream dams.

Climate within the study area is temperate with average monthly air temperatures of 13 °C in winter to 29 °C in summer, and an annual rainfall of 89–102 cm with peak precipitation in May and September42,43. Air temperature and water temperature within the middle Brazos generally follow each other closely with 1–2 °C variation and a delay on the scale of hours to days depending on the water depth (Fig. S7;44,45). Groundwater recharge in region tends to be confined to localized portions of the river and is primarily sourced from recent event waters (i.e., flooding and precipitation) stored in the stream banks. Previous works have found that the river gains groundwater discharge from the unconfined Brazos River Alluvium during high river discharge and can receive deeper groundwater discharge from the Yegua-Jackson aquifer during low river discharge as river flow increases46–48. Both ground water discharge and release from impounded reservoirs, primarily from the upstream Lake Whitney40,47, can buffer the river temperature by providing sources of cool water in the summer and warm water during the winter49. However, shells were collected away from potential influences (i.e., impoundments and water treatment release) that may contribute to short term instability in the air and water temperature relationship, so for the purpose of this study, we can assume that monthly air temperature serves as a reasonable representation of mean monthly water temperatures.

Results

Microscopic examination of thin sections allows the discrimination of dark and light bands in the shell cross-sections, the first step in developing sclerochronologies. Dark bands are visible under transmitted light as grey bands relative to other growth periods that appear light brown to tan. Generally, dark bands continue throughout the shell as a single band. However, in the juvenile portion of our specimens, some dark bands start as a single band and bifurcate into multiple dark bands near the ventral margin of the shell (Fig. 2).Figure 2 Transmitted light image of specimen H3R (Amblema plicata) with sampling area in shaded in blue and green to represent growth for year 1 and 2, respectively. Red lines show bifurcation in opaque “annual band” towards the ventral margin. Yellow dashed lines show additional annual bands not sampled.

Our reconstructed sclerochronologies have nearly monthly resolution in the first year’s growth and are extended to bi-monthly (two samples per month) in the second year as the shell became longer and sampling area extended. Viewed from bottom to top to best show the sequence of analyses, Figs. 3 and 4 show the progression of measured δ18Oshell and calculated Δ47 temperatures (T(Δ47); Figs. 3D and 4C), comparison with available discharge records (Fig. 3C), and reconstructed δ18Owater (Figs. 3B and 4B) and monthly growth rate (Figs. 3A and 4A). Measured δ13C and δ18O in modern shell 3R5 show no significant correlation (Fig. S3B) and range from − 10.9 ± 0.6‰–− 5.2 ± 0.2‰ to − 4.6 ± 0.1‰– − 1.3 ± 0.4‰, respectively (Table S1). Measured δ13C and δ18O in historic shell H3R show a significantly negative trend (Fig. S3B) and range from − 8.9 ± 0.1‰– − 6.7 ± 0.2‰ to − 5.5 ± 0.2‰– − 1.6 ± 0.1‰, respectively (Table S1).Figure 3 Data for modern shell 3R5. From top to bottom: (A) monthly growth rate (yellow); (B) reconstructed δ18Owater (blue); (C) measured precipitation (purple) and discharge (dark grey); and (D) shell δ18O (light grey) and △47 temperature (black). The vertical shaded bands represent dark bands in the shell, with the dashed red line representing a disturbance band and the solid red line representing an annual band in the 3B-D. Modern measured values of δ18Owater are represented by the light blue horizontal rectangle in 3B. Seasons and years are shown at the top.

Figure 4 Data for historic shell H3R. From top to bottom: (A) monthly growth rate (teal); (B) reconstructed δ18Owater (blue); (C) shell δ18O (light grey) and △47 temperature (black). The vertical shaded bands represent dark bands in the shell, with the dashed red line representing a disturbance band and the solid red line representing an annual band in the 4B-C. Modern measured values of δ18Owater are represented by the light blue horizontal rectangle in 4B. Seasons and years are shown at the top.

The Δ47 values for 3R5 and H3R range from 0.579 ± 0.015 to 0.619 ± 0.012 and 0.575 ± 0.0012 to 0.620 ± 0.013, respectively. T(Δ47) values are within monthly mean temperatures of 10.8 °C in January to 29.8 °C August (Fig. S4;35) and previous mussel tolerances of ~ 12 to 35 °C2,36–38. T (Δ47) ranges from 17 ± 3 to 30 ± 5 °C in 3R5 (Fig. 3D) and 16 ± 4 to 31 ± 5 °C in H3R (Fig. 4C). Preliminary chronologies were established following the trend for monthly mean temperature data for College Station (Fig. S6), where lowest and highest clumped isotope temperatures are assigned to the months with minimum and maximum recorded temperature, January and July, respectively. The lowest temperatures were measured in the dark bands that cut through both the interior and exterior surface of the shell and confirm previous studies of “annual banding” in freshwater mussels11,21,50. Dark bands that failed to cut throughout the outermost layer in the shell were determined to be disturbance bands and occurred at both cool and warm temperatures.

Growth rate chronologies were established using T(Δ47) and high-resolution growth band analyses from thin sections stained and etched with Mutvie’s solution (51; Fig. 5). Fine banding patterns were counted near the center portion of the shells where band width was more uniform, and time could be considered equivalent to band width. The monthly boundaries determined near the center of the shell were then traced throughout the shell to the ventral margin and growth rates were estimated by comparing the measured length of the shell from a stationary point on the umbo to the ventral margin at each monthly boundary marker. Additional description of band counting can be viewed in the method and supplemental sections. In the modern shell (3R5), growth is highest in the late fall and summer and ranges from 0.4 to 2.6 mm per month (Fig. 3A). Growth rates in the historic shell (H3R) are highest in late spring and summer and range from 0.3 to 2.2 mm per month (Fig. 4A). Reconstructed δ18Owater values are within measured mean monthly values (− 5.2‰ to 0.4‰; 33) and range from − 4.7 ± 0.8‰ to − 0.2 ± 0.7‰ and − 4.9 ± 1.1‰ to − 0.5 ± 0.7‰ for 3R5 (Fig. 3B) and H3R (Fig. 4B), respectively.Figure 5 (A) Reflected light of H3R stained with Mutvei’s solution in 2.5X zoom. (B) Zoomed grayscale image of the yellow rectangle in image (A) with depictions of finer banding patterns along the ventral margin.

Comparisons of our reconstructed water values are shown in Fig. 6, where the baseline trend for precipitation, or meteoric water, was established using precipitation data collected from 1962 to 1965 and 1972 to 1976 in Waco, TX, ~ 148 km upstream of the mussel collecting site52. Reconstructed isotope values from H3R and 3R5, and the measured water values from 2012 to 2013 were corrected for the evaporation effect on the Brazos River water by subtracting 1.7‰ based on the difference between the mean 18O-enrichment in the Brazos River and the weighted mean δ18O of precipitation determined by Van Plantinga and Grossman40. This correction makes the comparison of our reconstructed values directly comparable to the measured values of precipitation values.Figure 6 Comparison of Δ47-derived δ18Owater values after applying the correction of 1.7‰ to account for evaporation effects on the Brazos River water for modern shell 3R5 (filled yellow circles), historic shell H3R (teal open circles), measured river samples from Van Plantinga and Grossman (2017; black open triangles), and average δ18O of precipitation samples from Waco, TX during 1962–65 and 1972–76 (light grey dashed line, with minima and maxima bounded by the grey shaded area). All error bars represent the standard error. The blue and red shaded arrow on the right indicates direction of environmental conditions (wet = precipitation and dry = drought).

Discussion

We interpret the isotopic record of mussel growth within the framework of mussel life cycles. Freshwater mussels have multi-stage life cycles, transitioning through the fertilized egg, glochidia, juvenile, and adult stages, respectively. A. plicata is a short-term brooder that spawns in late spring to early summer and releases glochidia a few weeks to a month later. Juvenile mussels allocate more energy towards growing and have higher growth rates before slowing as they reach sexual maturity at the age of 3 to 9 years19. Growth generally follows a logistic trend that eventually reaches an asymptotic length where growth becomes very slow or ceases. The faster growth rates of juvenile mussels can be rather plastic depending on life history and the surrounding environment16,19, but in general the faster growth rates make juvenile shell the best for detailed environmental records utilizing isotopic analyses. Specifically, clumped isotope measurements require large sample sizes (~ 1.2 mg) compared with oxygen isotopes (0.05 mg), so slower growth in the adult proportion of the shell would lead to a higher degree of averaging and less detailed environmental records.

In temperate latitudes, 23.5 to 66.5 degrees, freshwater mussel growth has typically been thought to slow and cease in the fall when temperatures reach below 6–12 °C and begin again the following spring when temperatures reach ~ 12–15 °C4; thus, peak growth occurs during spring and summer when temperatures are optimal (26–28 °C;53). However, in subtropical latitudes like our study area, seasonal temperature range is lower, which can extend the growing season and minimize growth cessation43. Haag et al.54 found that in streams across Kentucky, the onset of growth occurred at 17.8 °C. This prediction is higher than previously determined values of ~ 12–15 °C55,56 but may have been due to variation between study environments since growth is not solely temperature dependent. Instead, growth can be dependent on multiple factors including nutrient availability, alkalinity, temperature, latitude, and life stage54,57,58.

At our study site, average monthly air temperature is considered equivalent to water temperature plus or minus 1–2 °C44. Freshwater mussel growth is typically thought to slow or cease during winter when water temperatures reach below 6–12 °C and begin again the following spring when temperatures reach ~ 12–15 °C4. Mean air temperatures near the study site reach below 12 °C during the months of December and January but can experience average daily highs up to 16.5 °C and 17.4 °C, respectively. These average highs are not much lower than the predicted value for the onset of growth (~ 17.8 °C) determined by Haag et al.54 and higher than previously determined ranges (~ 12–15 °C) that can trigger the onset of growth. Therefore, we cannot confidently assume a complete cessation of growth from December to January, and instead assume slower growth at these lower temperatures. Based on our understanding of the study site and local temperature trends, we assume that the lowest calculated Δ47 temperatures occur in January and that Δ47 temperature follow the seasonal trend for air temperature based on monthly norms from NOAA climatological reports (i.e., means for 1991 to 2020 and the air–water temperature relationship for the Brazos River44,45).

For the modern shell (3R5), we determined the years sampled by the isotope analyses by first estimating the age of the specimen (~ 4 years) based on criteria described in Haag and Commens-Carson17, then subtracted the number of dark bands characteristic of annual bands within the sampled area from known date of collection. Since 3R5 was collected in August 2013, this yielded an age for isotopic samples of summer 2009 to summer 2010. Based on comparisons with local precipitation and discharge records, we find that the dark bands that fail to extend to the exterior of the shell align with periods of high river discharge and precipitation (Fig. 3C). During elevated stream flows, increased turbidity and sediment suspension can displace or stress mussels and lead to dissolution or precipitation misalignments that cause deposition of dark “disturbance bands” as the organism’s valves close and reopen9,17,59,60.

Reconstructed δ18Owater values in 3R5 for the period of high rainfall and increased natural discharge from mid-August to late-September (2009) decrease from − 0.6 to − 3.8‰ (Fig. 3B), revealing a significant freshening signal with the increased amount of 18O-depleted rainwater discharging to the river. Measured Δ47 temperature also show a gradual decrease, most likely in response to the cooler rainfall entering the river and decreased residence time that the water has available to warm during faster streamflow61–63. After that, both reconstructed δ18Owater and T(Δ47) increase before declining in fall and winter. The second dark band, or annual band, in 3R5 also correlates with the onset of high discharge based on our chronology. However, the observed increase in discharge in the 2009–2010 winter does not correlate with an increase in measured local precipitation, implying the major increase in flow was a response to release from the upstream impoundment at Lake Whitney and more northern precipitation events (Fig. S9).

The temperature corresponding with the annual band during this time recorded the coldest temperature, 17 °C, but was still higher than the expected 11 °C and 15 °C winter minimum based on the seasonal trend for mean air temperature (Fig. S6) and available water (Fig. S7) temperature, respectively. Winter releases from impoundments are known to increase water temperatures when released from the bottom waters of reservoirs and elevate δ18Owater due to enhanced evaporation in impounded water40,46,64, which would contribute to the warmer water temperatures observed and the potential to lengthen peak growing seasons. Water temperatures in the Lake Whitney reservoir are typically thermally stratified with cooler bottom water throughout the year until the fall turnover in October–November when waters become isothermal and lasts until the following spring. During this isothermal period, the water becomes well mixed with depth and water temperatures follow the surface and air temperature65,66. However, bottom waters after overturning can become thermally insulated from cooler surface temperatures in winter and contribute to warmer downstream temperature that deviates from cooler air temperatures during Lake Whitney releases44.

For the historic shell, H3R, we could not define an exact age of the sampling period because the archive lists an imprecise collection date (~ 1880–1900), but we were able to define it as growing within an environment that predates the damming of the Brazos River. Despite the lack of impoundment, reconstructed δ18Owater values calculated from H3R were within the monthly mean of modern measured values and those predicted by 3R5. Disturbance bands in H3R also correlated with decreases in δ18Owater before and after the summer season, suggesting that increased rainfall and associated runoff inhibited growth. The timing of these events also validates our methods used to reconstruct the chronology since peak precipitation in the Brazos River basin occurs in May and September42. During H3R’s late spring rainfall event in Year 2, Δ47 temperature measured ~ 25 °C, while Δ47 temperature during the rainfall event at the summer-fall transition decreased to ~ 18 °C (Fig. 4C,D). In both H3R and 3R5, growth rates slowed during the initial increases in stream flow but rebounded afterwards, most likely due to increased availability of nutrients via runoff and increased particulate settling with decreased turbulence67–69.

To compare our reconstructed δ18Owater results to measured δ18O values of meteoric water, we corrected the reconstructed and measured δ18Owater to account for evaporation effect on the surface water, and created a baseline using precipitation data collected from 1962–1965 to 1972–1976 in Waco, TX, ~ 148 km upstream of the mussel collecting site40,52. Overall, our reconstructed δ18Owater values show the expected trend for δ18O in the river when meteoric water is the main source (Fig. 6). Mean measured δ18Owater from Van Plantinga and Grossman40 also followed the meteoric trend, except for fall and winter measurements in 2012 when drought conditions were present. Though our baseline relies on a small sample set due to limited number of measured δ18O in precipitation recorded near our study site, deviations from the meteoric baseline can be an indication for drought and flood conditions.

Growth trends for both shells agree with previous findings of high growth in spring and summer, but also confirm that rapid growth can occur well into late fall4,55. There is no significant difference between modern and historic growth rates (two tailed t-test, p > 0.05; Fig. S7). However, only the historic shell (H3R) shows a significant correlation between growth rate and both reconstructed δ18Owater and T(Δ47) (R = 0.67, p < 0.001, and R = 0.43, p = 0.046, respectively; Fig. 8C,D). The stronger correlation between growth and reconstructed δ18Owater suggests that historic growth is more sensitive to changes in δ18Owater, or precipitation, than to temperature (p ≥ 0.05 for both H3R and 3R5). Of course, these results are based on only two specimens; nevertheless, they contribute valuable insight for understanding the factors driving community level changes in mussel populations within the Brazos River and reveal the potential of clumped isotopes in freshwater mussels as a proxy in rivers that experience flow alterations.

Conclusions and Implications

This study successfully used clumped isotope analyses to reconstruct seasonal δ18Owater values for A. plicata from the Brazos River, located in Central Texas. Disturbance bands correlate with rainfall and high streamflow events based on reconstructed δ18Owater and known relationships between river water δ18O and precipitation. Calculated growth rates ranged from 0.3 to 2.5 mm per month with high growth in the spring through fall, suggesting an extended growing season. No significant relationship with growth rate was found in the modern shell (3R5), but the significant relationship between reconstructed δ18Owater and growth rates in the historic shell (H3R) suggests that precipitation and flow rate are stronger drivers than temperature in historic growth that predates river impoundment and represents natural flow conditions. The results from these two shells demonstrate that clumped isotope analyses can provide valuable understanding of local hydrologic cycles and climate events such as droughts and floods on multiyear time scales. However, broader application using museum, fossil, and live-collected specimens would make it possible to examine historic and prehistoric climate change over a wide swath of the US mid-continent, while at the same time revealing the impact of river impoundment and climate change on mussel growth and survival.

In addition to understanding the relationship between mussel growth, flow rates, and climate change, clumped isotope analyses can be helpful for validating growth and longevity estimates. This is important because these biological end points are often used to make inferences on how a species will respond to environmental disturbance or different management actions. For example, Randklev et al.9 evaluated how mussels cope with flow rates and found that species with high growth rates and reduced longevity were more tolerant of streambed disturbance during high flow events. In contrast, species with lower growth rates and longer lifespans were less resilient. Such inferences can be used to guide environmental flow recommendations70. However, it is important that age and growth estimates are validated to minimize inaccurate determinations that could lead to inappropriate management recommendations that are not protective during floods or droughts. This could be done by additional coupled isotope-growth band studies of shells from the Brazos and other rivers as an alternative to the time-consuming caged experiments that can create artificial growth hiatuses with each handling. The literature is replete with examples of the impact of inaccurate estimates for these endpoints. For example, the orange roughy (Hoplostethus atlanticus) was intensively harvested off New Zealand based on the presumption that it lived 20 to 30 years and had a fast growth rate. More recent research has shown that it is much longer lived, over 100 years, with a slow growth rate71.

For unionid mussels, examples of inaccurate determinations of age and growth and their impact on management actions are not as well documented. This is due, in part, to the fact that growth and longevity data are lacking for most mussel species. To date, these data are available for only 35 of the 300 mussel species known to occur in North America. For these species, very few if any have been validated and to the best of our knowledge clumped isotopes have only been used on two species (Amblema plicata and Crytonias tampicoensis27). Thus, expanding this knowledge base and validating age and growth estimates is critical for the long-term conservation of mussels as a group. Future efforts should focus not only on evaluating the ecology of different species but also on creating records spanning historic drought periods such as the Texas drought of record (1950–57) and deeper time (e.g., Pliocene warm period) with high temperatures and drought events that resemble predictions for future climate.

Materials and methods

Samples and sampling

Shell specimens were collected from the main stem of the middle Brazos River west of College Station, TX (30.55 latitude, − 96.42 longitude; Fig. 1) in 2013 and ~ 1900 (30.6 latitude, − 96.4 longitude; Fig. 1). The modern sample was collected live from the shoal water by the third author and the historic specimen was provided by the Singley Collection from the University of Texas at Austin Non-Vertebrate Paleontology Laboratory (collector unknown). Discharge data were compiled from the USGS gauge 08,108,700 near Bryan, TX72 and precipitation δ18O values for Waco, TX were retrieved from the Global Network of Isotopes in Precipitation (GNIP;52).

Upon arrival, shells were cleaned with double deionized water and the right valves were sectioned along the long axis of growth and then sequentially cut into thick (~ 2–3 mm) and thin sections (300–400 μm). Thin sections were then examined under the microscope for physical discoloration and damage that could indicate chemical alterations. A New Wave micromill equipped with a ~ 400 μm dental bur, operated at low speed to prevent heating and reordering of the clumped isotopes, was used to mill carbonate powders parallel to growth bands of juvenile growth (first 2–3 years) in 300–350 μm increments along the ventral margin of thick sections (Fig. 2).

Isotope measurements

Isotope analyses were performed using a Kiel IV carbonate device coupled to a Thermo Fisher Scientific 253 + isotope ratio mass spectrometer (IRMS) in the Stable Isotope Geosciences Facility (SIGF) at Texas A&M University using methodology described in Meinicke et al.73 and Barney and Grossman35. We analyzed 5–10 replicates of each individual powdered sample weighing ~ 120 µm and used a running mean of three samples to obtain 15–30 replicates for each data point and increase precision. Analytical precisions (1σ) for ETH standards during the period of sample analyses were ± 0.03‰, ± 0.05‰, and 0.040‰ for δ13C, δ18O, and Δ47, respectively. For sample analyses of unknowns, δ18Oshell, δ13Cshell, and Δ47 uncertainties were calculated based on the standard error for each sample and its 5–10 replicates and then averaged for each three-sample mean (15–30 replicates) that define each data point. Errors for T(Δ47) and δ18Owater were estimated using the propagation of the minima and maxima of the Δ47 uncertainty (Table S1). Overall, the average standard error (SE) for measured Δ47 equates to ~  ± 3.5 °C. Clumped (Δ47) temperatures were calculated using the equation of Anderson et al.74:1 Δ47=0.0391±0.0004·106T2+0.154±0.004

where T is the water temperature in kelvin and Δ47 is standardized within the “InterCarb-Carbon Dioxide Equilibrium Scale” (I-CDES). Reconstructed δ18Owater values were calculated using the rearranged δ18O paleotemperature equation from Grossman and Ku26 corrected for water δ18O in terms of VSMOW75:2 T∘C=19.7-4.34δ18OaragoniteVPDB-δ18OwaterVSMOW

Chronologies and growth rates

Monthly chronologies were assigned using a combination of the Δ47 temperature, shell banding, and monthly air temperature norms, or monthly averages over a 30-year period, for College Station based on NOAA records. Preliminary chronologies were determined using the shell length at individual milling transects and the local trend for monthly air temperature norms. The shell intervals with the coolest Δ47 temperatures were assigned to the month of January, while the interval with the maximum Δ47 temperature was assigned to the month of July, matching the local minimum and maximum of average seasonal air temperatures. This preliminary chronology was then refined using thin sections stained with Mutvei’s solution to enhance the fine banding patterns and better resolve the distance between annual bands and monthly boundaries (Figs. S1 and S2;51). Sub-annual bands were counted in the center of the shell near the beak area where bands had almost uniform band width. Assuming equal bandwidth and time, we subdivided the sub-annual growth bands between the previously determined temperature minima and maxima within a years’ time, or area between two annual bands. The monthly boundaries between sub-annual bands in the center of the shell were then traced throughout the shell and marked on scan images of each shell thin section. Monthly growth rates were then calculated by measuring the difference in shell length from the umbo to the edge of the ventral margin at the assigned monthly boundaries. See supplemental material for further discussion on chronology anchoring and sub-annual banding for the individual shells.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-58246-w.

Acknowledgements

We gratefully acknowledge funding from the Michel T. Halbouty Chair; the W.G. Mills Memorial Fellowship in Hydrology administered through the Texas Water Resources Institute, part of Texas Agrilife Research; the Texas Agrilife Extension Service; and the Colleges of Agriculture & Life Sciences and Arts & Sciences at Texas A&M University. The clumped isotope mass spectrometer system was purchased through an award from the Chancellor’s $5 million commitment to mass spectrometry. We thank the Singley Collection from the University of Texas at Austin Non-Vertebrate Paleontology Laboratory for providing the historic sample. We thank Chris Maupin, Co-Director of the Stable Isotope Geoscience Facility (SIGF), and graduate students Zeyang Sun and Bryce Barney for analytical assistance and instrument maintenance. Lastly, the manuscript benefitted from careful reviews by Niels de Winter and Gregory Henkes.

Author contributions

Conceptualization and visualization were contributed equally by M.B. and E.G. The first draft was written by M.B. with revisions, editing, and supervision from E.G. and C.R.

Data availability

Data are provided in the supplemental material and more detailed methods and raw data will be made available upon request from the authors M.A. Brewer and E.L. Grossman.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Newton, T. J., Johnson, N. A. & Hu, D. H. U.S. Geological Survey science vision for native freshwater mussel research in the United States. (2023; https://pubs.usgs.gov/publication/cir1511).
2. Lopes-Lima M Burlakova L Karatayev A Gomes-dos-Santos A Zieritz A Froufe E Bogan AE Revisiting the North American freshwater mussel genus Quadrula sensu lato (Bivalvia Unionidae): Phylogeny, taxonomy and species delineation Zool. Scr. 2019 48 313 336 10.1111/zsc.12344
Lopes-Lima, M. et al. Revisiting the North American freshwater mussel genus Quadrula sensu lato (Bivalvia Unionidae): Phylogeny, taxonomy and species delineation. Zool. Scr. 48, 313–336 (2019).10.1111/zsc.12344
3. Williams JD Warren ML Cummings KS Harris JL Neves RJ Conservation status of freshwater mussels of the United States and Canada Fisheries 1993 18 6 22 10.1577/1548-8446(1993)018<0006:CSOFMO>2.0.CO;2
Williams, J. D., Warren, M. L., Cummings, K. S., Harris, J. L. & Neves, R. J. Conservation status of freshwater mussels of the United States and Canada. Fisheries 18, 6–22 (1993).10.1577/1548-8446(1993)018<0006:CSOFMO>2.0.CO;2
4. Haag WR North American freshwater mussels: Natural history, ecology, and conservation 2012 Cambridge University Press
Haag, W. R. North American freshwater mussels: Natural history, ecology, and conservation (Cambridge University Press, 2012).
5. Ferreira-Rodríguez N Akiyama YB Aksenova OV Araujo R Christopher Barnhart M Bespalaya YV Bogan AE Bolotov IN Budha PB Clavijo C Clearwater SJ Darrigran G Do VT Douda K Froufe E Gumpinger C Henrikson L Humphrey CL Johnson NA Klishko O Klunzinger MW Kovitvadhi S Kovitvadhi U Lajtner J Lopes-Lima M Moorkens EA Nagayama S Nagel KO Nakano M Negishi JN Ondina P Oulasvirta P Prié V Riccardi N Rudzīte M Sheldon F Sousa R Strayer DL Takeuchi M Taskinen J Teixeira A Tiemann JS Urbańska M Varandas S Vinarski MV Wicklow BJ Zając T Vaughn CC Research priorities for freshwater mussel conservation assessment Biol. Conserv. 2019 231 77 87 10.1016/j.biocon.2019.01.002
Ferreira-Rodríguez, N. et al. Research priorities for freshwater mussel conservation assessment. Biol. Conserv. 231, 77–87 (2019).10.1016/j.biocon.2019.01.002
6. Strayer DL Downing JA Haag WR King TL Layzer JB Newton TJ Nichols JS Changing perspectives on pearly mussels, North America’s most imperiled animals BioScience 2004 54 429 439 10.1641/0006-3568(2004)054[0429:CPOPMN]2.0.CO;2
Strayer, D. L. et al. Changing perspectives on pearly mussels, North America’s most imperiled animals. BioScience 54, 429–439 (2004).10.1641/0006-3568(2004)054[0429:CPOPMN]2.0.CO;2
7. Brainwood M Burgin S Byrne M The impact of small and large impoundments on freshwater mussel distribution in the Hawkesbury-Nepean River, Southeastern Australia River Res. Appl. 2008 24 1325 1342 10.1002/rra.1087
Brainwood, M., Burgin, S. & Byrne, M. The impact of small and large impoundments on freshwater mussel distribution in the Hawkesbury-Nepean River, Southeastern Australia. River Res. Appl. 24, 1325–1342 (2008).10.1002/rra.1087
8. Randklev CR Ford N Wolverton S Kennedy JH Robertson C Mayes K Ford D The influence of stream discontinuity and life history strategy on mussel community structure: A case study from the Sabine River, Texas Hydrobiologia 2016 770 173 191 10.1007/s10750-015-2586-5
Randklev, C. R. et al. The influence of stream discontinuity and life history strategy on mussel community structure: A case study from the Sabine River, Texas. Hydrobiologia 770, 173–191 (2016).10.1007/s10750-015-2586-5
9. Randklev CR Hart MA Khan JM Tsakiris ET Robertson CR Hydraulic requirements of freshwater mussels (Unionidae) and a conceptual framework for how they respond to high flows Ecosphere 2019 10 e02975 10.1002/ecs2.2975
Randklev, C. R., Hart, M. A., Khan, J. M., Tsakiris, E. T. & Robertson, C. R. Hydraulic requirements of freshwater mussels (Unionidae) and a conceptual framework for how they respond to high flows. Ecosphere 10, e02975 (2019).10.1002/ecs2.2975
10. Khan JM Dudding J Hart M Tsakiris E Randklev CR Linking life history strategies and historical baseline information shows effects of altered flow regimes and impoundments on freshwater mussel assemblages Freshw. Biol. 2020 65 1950 1961 10.1111/fwb.13591
Khan, J. M., Dudding, J., Hart, M., Tsakiris, E. & Randklev, C. R. Linking life history strategies and historical baseline information shows effects of altered flow regimes and impoundments on freshwater mussel assemblages. Freshw. Biol. 65, 1950–1961 (2020).10.1111/fwb.13591
11. Rypel AL Haag WR Findlay RH Validation of annual growth rings in freshwater mussel shells using cross dating Can. J. Fish. Aquat. Sci. 2008 65 2224 2232 10.1139/F08-129
Rypel, A. L., Haag, W. R. & Findlay, R. H. Validation of annual growth rings in freshwater mussel shells using cross dating. Can. J. Fish. Aquat. Sci. 65, 2224–2232 (2008).10.1139/F08-129
12. Bolotov IN Makhrov AA Gofarov MY Aksenova OV Aspholm PE Bespalaya YV Kabakov MB Kolosova YS Kondakov AV Ofenböck T Climate warming as a possible trigger of keystone mussel population decline in oligotrophic rivers at the continental scale Sci. Rep. 2018 8 35 10.1038/s41598-017-18873-y 29311629
Bolotov, I. N. et al. Climate warming as a possible trigger of keystone mussel population decline in oligotrophic rivers at the continental scale. Sci. Rep. 8, 35 (2018).29311629 10.1038/s41598-017-18873-y
13. Nobles T Zhang Y Survival, growth and condition of freshwater mussels: Effects of municipal wastewatere effluent PLOS ONE 2015 10 e0128488 10.1371/journal.pone.0128488 26042840
Nobles, T. & Zhang, Y. Survival, growth and condition of freshwater mussels: Effects of municipal wastewatere effluent. PLOS ONE 10, e0128488 (2015).26042840 10.1371/journal.pone.0128488
14. Hornbach DJ Stutzman HN Hove MC Kozarek JL MacGregor KR Newton TJ Ries PR Influence of surrounding land-use on mussel growth and glycogen levels in the St. Croix and Minnesota river Basins Hydrobiologia 2021 848 3045 3063 10.1007/s10750-019-04016-z
Hornbach, D. J. et al. Influence of surrounding land-use on mussel growth and glycogen levels in the St. Croix and Minnesota river Basins. Hydrobiologia 848, 3045–3063 (2021).10.1007/s10750-019-04016-z
15. Goodwin DH Gillikin DP Banker R Watters GT Dettman DL Romanek CS Reconstructing intra-annual growth of freshwater mussels using oxygen isotopes Chem. Geol. 2019 526 7 22 10.1016/j.chemgeo.2018.07.030
Goodwin, D. H. et al. Reconstructing intra-annual growth of freshwater mussels using oxygen isotopes. Chem. Geol. 526, 7–22 (2019).10.1016/j.chemgeo.2018.07.030
16. Sansom BJ Atkinson CL Vaughn CC Growth and longevity estimates for mussel populations in three Ouachita Mountain rivers Freshw. Mollusk Biol. Conserv. 2016 19 19–26 18
Sansom, B. J., Atkinson, C. L. & Vaughn, C. C. Growth and longevity estimates for mussel populations in three Ouachita Mountain rivers. Freshw. Mollusk Biol. Conserv. 19(19–26), 18 (2016).
17. Haag WR Commens-Carson AM Testing the assumption of annual shell ring deposition in freshwater mussels Can. J. Fish. Aquat. Sci. 2008 65 493 508 10.1139/f07-182
Haag, W. R. & Commens-Carson, A. M. Testing the assumption of annual shell ring deposition in freshwater mussels. Can. J. Fish. Aquat. Sci. 65, 493–508 (2008).10.1139/f07-182
18. Neves RJ Moyer Steven N Evaluation of techniques for age determination of freshwater mussels (Unionidae) Am. Malacol. Bull. 1988 6 179 188
Neves, R. J., Moyer, & Steven, N. Evaluation of techniques for age determination of freshwater mussels (Unionidae). Am. Malacol. Bull. 6, 179–188 (1988).
19. Haag WR Rypel AL Growth and longevity in freshwater mussels: evolutionary and conservation implications Biol. Rev. 2011 86 225 247 10.1111/j.1469-185X.2010.00146.x 20608928
Haag, W. R. & Rypel, A. L. Growth and longevity in freshwater mussels: evolutionary and conservation implications. Biol. Rev. 86, 225–247 (2011).20608928 10.1111/j.1469-185X.2010.00146.x
20. Valery Z Eduardo San M Richard JN Angeles L Carlos F Rafaela A Victor B Ekaterina P Sviatoslav K Torbjörn J Life span variation of the freshwater pearl shell: A model species for testing longevity mechanisms in animals AMBIO J. Hum. Environ. 2000 29 102 105 10.1579/0044-7447-29.2.102
Valery, Z. et al. Life span variation of the freshwater pearl shell: A model species for testing longevity mechanisms in animals. AMBIO J. Hum. Environ. 29, 102–105 (2000).10.1579/0044-7447-29.2.102
21. Veinott GI Cornett RJ Identification of annually produced opaque bands in the shell of the freshwater mussel Elliptio complanata using the seasonal cycle of δ18Ο Can. J. Fish. Aquat. Sci. 1996 53 372 379
Veinott, G. I. & Cornett, R. J. Identification of annually produced opaque bands in the shell of the freshwater mussel Elliptio complanata using the seasonal cycle of δ18Ο. Can. J. Fish. Aquat. Sci. 53, 372–379 (1996).
22. Ricken W Steuber T Freitag H Hirschfeld M Niedenzu B Recent and historical discharge of a large European river system–oxygen isotopic composition of river water and skeletal aragonite of Unionidae in the Rhine Palaeogeogr. Palaeoclimatol. Palaeoecol. 2003 193 73 86 10.1016/S0031-0182(02)00713-7
Ricken, W., Steuber, T., Freitag, H., Hirschfeld, M. & Niedenzu, B. Recent and historical discharge of a large European river system–oxygen isotopic composition of river water and skeletal aragonite of Unionidae in the Rhine. Palaeogeogr. Palaeoclimatol. Palaeoecol. 193, 73–86 (2003).10.1016/S0031-0182(02)00713-7
23. Versteegh EAA Vonhof HB Troelstra SR Kroon D Can shells of freshwater mussels (Unionidae) be used to estimate low summer discharge of rivers and associated droughts? Int. J. Earth Sci. 2011 100 1423 1432 10.1007/s00531-010-0551-0
Versteegh, E. A. A., Vonhof, H. B., Troelstra, S. R. & Kroon, D. Can shells of freshwater mussels (Unionidae) be used to estimate low summer discharge of rivers and associated droughts?. Int. J. Earth Sci. 100, 1423–1432 (2011).10.1007/s00531-010-0551-0
24. Kelemen Z Gillikin DP Graniero LE Havel H Darchambeau F Borges AV Yambélé A Bassirou A Bouillon S Calibration of hydroclimate proxies in freshwater bivalve shells from Central and West Africa Geochim. et Cosmochim. Acta 2017 208 41 62 10.1016/j.gca.2017.03.025
Kelemen, Z. et al. Calibration of hydroclimate proxies in freshwater bivalve shells from Central and West Africa. Geochim. et Cosmochim. Acta 208, 41–62 (2017).10.1016/j.gca.2017.03.025
25. Schöne BR Meret AE Baier SM Fiebig J Esper J McDonnell J Pfister L Freshwater pearl mussels from northern Sweden serve as long-term, high-resolution stream water isotope recorders Hydrol. Earth Syst. Sci. 2020 24 673 696 10.5194/hess-24-673-2020
Schöne, B. R. et al. Freshwater pearl mussels from northern Sweden serve as long-term, high-resolution stream water isotope recorders. Hydrol. Earth Syst. Sci. 24, 673–696 (2020).10.5194/hess-24-673-2020
26. Grossman EL Ku T-L Oxygen and carbon isotope fractionation in biogenic aragonite: Temperature effects Chem. Geol. Isot. Geosci. Sect. 1986 59 59 74 10.1016/0168-9622(86)90057-6
Grossman, E. L. & Ku, T.-L. Oxygen and carbon isotope fractionation in biogenic aragonite: Temperature effects. Chem. Geol. Isot. Geosci. Sect. 59, 59–74 (1986).10.1016/0168-9622(86)90057-6
27. Van Plantinga AA Grossman EL Stable and clumped isotope sclerochronologies of mussels from the Brazos River, Texas (USA): Environmental and ecologic proxy Chem. Geol. 2018 502 55 65 10.1016/j.chemgeo.2018.10.012
Van Plantinga, A. A. & Grossman, E. L. Stable and clumped isotope sclerochronologies of mussels from the Brazos River, Texas (USA): Environmental and ecologic proxy. Chem. Geol. 502, 55–65 (2018).10.1016/j.chemgeo.2018.10.012
28. Pfister L Thielen F Deloule E Valle N Lentzen E Grave C Beisel J-N McDonnell JJ Freshwater pearl mussels as a stream water stable isotope recorder Ecohydrology 2018 11 e2007 10.1002/eco.2007
Pfister, L. et al. Freshwater pearl mussels as a stream water stable isotope recorder. Ecohydrology 11, e2007 (2018).10.1002/eco.2007
29. Eiler JM “Clumped-isotope” geochemistry—the study of naturally-occurring, multiply-substituted isotopologues Earth Planet. Sci. Lett. 2007 262 309 327 10.1016/j.epsl.2007.08.020
Eiler, J. M. “Clumped-isotope” geochemistry—the study of naturally-occurring, multiply-substituted isotopologues. Earth Planet. Sci. Lett. 262, 309–327 (2007).10.1016/j.epsl.2007.08.020
30. Ghosh P Adkins J Affek H Balta B Guo W Schauble EA Schrag D Eiler JM 13C–18O bonds in carbonate minerals: A new kind of paleothermometer Geochim. et Cosmochim. Acta 2006 70 1439 1456 10.1016/j.gca.2005.11.014
Ghosh, P. et al. 13C–18O bonds in carbonate minerals: A new kind of paleothermometer. Geochim. et Cosmochim. Acta 70, 1439–1456 (2006).10.1016/j.gca.2005.11.014
31. Epstein S Buchsbaum R Lowenstam HA Urey HC Revised carbonate-water isotopic temperature scale Geol. Soc. Am. Bull. 1953 64 1315 1326 10.1130/0016-7606(1953)64[1315:RCITS]2.0.CO;2
Epstein, S., Buchsbaum, R., Lowenstam, H. A. & Urey, H. C. Revised carbonate-water isotopic temperature scale. Geol. Soc. Am. Bull. 64, 1315–1326 (1953).10.1130/0016-7606(1953)64[1315:RCITS]2.0.CO;2
32. Eiler JM Paleoclimate reconstruction using carbonate clumped isotope thermometry Quat. Sci. Rev. 2011 30 3575 3588 10.1016/j.quascirev.2011.09.001
Eiler, J. M. Paleoclimate reconstruction using carbonate clumped isotope thermometry. Quat. Sci. Rev. 30, 3575–3588 (2011).10.1016/j.quascirev.2011.09.001
33. Affek HP Clumped isotope paleothermometry: Principles, applications, and challenges Paleontol. Soc. Papers 2012 18 101 114 10.1017/S1089332600002576
Affek, H. P. Clumped isotope paleothermometry: Principles, applications, and challenges. Paleontol. Soc. Papers 18, 101–114 (2012).10.1017/S1089332600002576
34. Huntington KW Petersen SV Frontiers of carbonate clumped isotope thermometry Ann. Rev. Earth Planet. Sci. 2023 51 611 641 10.1146/annurev-earth-031621-085949
Huntington, K. W. & Petersen, S. V. Frontiers of carbonate clumped isotope thermometry. Ann. Rev. Earth Planet. Sci. 51, 611–641 (2023).10.1146/annurev-earth-031621-085949
35. Barney BB Grossman EL Reassessment of ocean paleotemperatures during the Late Ordovician Geology 2022 50 572 576 10.1130/G49422.1
Barney, B. B. & Grossman, E. L. Reassessment of ocean paleotemperatures during the Late Ordovician. Geology 50, 572–576 (2022).10.1130/G49422.1
36. Came RE Brand U Affek HP Clumped isotope signatures in modern brachiopod carbonate Chem. Geol. 2014 377 20 30 10.1016/j.chemgeo.2014.04.004
Came, R. E., Brand, U. & Affek, H. P. Clumped isotope signatures in modern brachiopod carbonate. Chem. Geol. 377, 20–30 (2014).10.1016/j.chemgeo.2014.04.004
37. de Winter NJ Witbaard R Kocken IJ Müller IA Guo J Goudsmit B Ziegler M Temperature dependence of clumped isotopes (∆47) in aragonite Geophys. Res. Lett. 2022 49 e2022GL099479 10.1029/2022GL099479 36589267
de Winter, N. J. et al. Temperature dependence of clumped isotopes (∆47) in aragonite. Geophys. Res. Lett. 49, e2022GL099479 (2022).36589267 10.1029/2022GL099479
38. Henkes GA Passey BH Grossman EL Shenton BJ Yancey TE Pérez-Huerta A Temperature evolution and the oxygen isotope composition of Phanerozoic oceans from carbonate clumped isotope thermometry Earth Planet. Sci. Lett. 2018 490 40 50 10.1016/j.epsl.2018.02.001
Henkes, G. A. et al. Temperature evolution and the oxygen isotope composition of Phanerozoic oceans from carbonate clumped isotope thermometry. Earth Planet. Sci. Lett. 490, 40–50 (2018).10.1016/j.epsl.2018.02.001
39. Huyghe D Daëron M de Rafelis M Blamart D Sébilo M Paulet Y-M Lartaud F Clumped isotopes in modern marine bivalves Geochim. et Cosmochim. Acta 2022 316 41 58 10.1016/j.gca.2021.09.019
Huyghe, D. et al. Clumped isotopes in modern marine bivalves. Geochim. et Cosmochim. Acta 316, 41–58 (2022).10.1016/j.gca.2021.09.019
40. Van Plantinga AA Grossman EL Roark EB Chemical and isotopic tracer evaluation of water mixing and evaporation in a dammed Texas river during drought River Res. Appl. 2017 33 450 460 10.1002/rra.3080
Van Plantinga, A. A., Grossman, E. L. & Roark, E. B. Chemical and isotopic tracer evaluation of water mixing and evaporation in a dammed Texas river during drought. River Res. Appl. 33, 450–460 (2017).10.1002/rra.3080
41. K. Hendrickson Jr. (Texas State Historical Association, 1952).
42. Lee K Singh VP Analysis of uncertainty and non-stationarity in probable maximum precipitation in Brazos River Basin J. Hydrol. 2020 590 125526 10.1016/j.jhydrol.2020.125526
Lee, K. & Singh, V. P. Analysis of uncertainty and non-stationarity in probable maximum precipitation in Brazos River Basin. J. Hydrol. 590, 125526 (2020).10.1016/j.jhydrol.2020.125526
43. Nielsen-Gammon JW The changing climate of Texas 2011 University of Texas Press 39 68
Nielsen-Gammon, J. W. The changing climate of Texas 39–68 (University of Texas Press, 2011).
44. Rawson, J. Reconnaissance of water temperature of selected streams in Southeastern Texas. vol. 105, (1970).
45. Segura C Caldwell P Sun G McNulty S Zhang Y A model to predict stream water temperature across the conterminous USA Hydrol. Process. 2015 29 2178 2195 10.1002/hyp.10357
Segura, C., Caldwell, P., Sun, G., McNulty, S. & Zhang, Y. A model to predict stream water temperature across the conterminous USA. Hydrol. Process. 29, 2178–2195 (2015).10.1002/hyp.10357
46. Chowdhury, A., Osting, T., Furnans, J. & Mathews, R. Groundwater-surface water interaction in the Brazos River Basin: Evidence from lake connection history and chemical and isotopic compositions. Texas Water Dev. Rep. 375, (2010).
47. Prince BH Marcantonio F Knappett PSK Characterizing groundwater–surface water interaction across the Brazos River watershed, Texas, with uranium isotopes Appl. Geochem. 2022 147 105491 10.1016/j.apgeochem.2022.105491
Prince, B. H., Marcantonio, F. & Knappett, P. S. K. Characterizing groundwater–surface water interaction across the Brazos River watershed, Texas, with uranium isotopes. Appl. Geochem. 147, 105491 (2022).10.1016/j.apgeochem.2022.105491
48. Rhodes KA Proffitt T Rowley T Knappett PSK Montiel D Dimova N Tebo D Miller GR The importance of bank storage in supplying sbseflow to rivers flowing through compartmentalized, alluvial aquifers Water Resour. Res. 2017 53 10539 10557 10.1002/2017WR021619
Rhodes, K. A. et al. The importance of bank storage in supplying sbseflow to rivers flowing through compartmentalized, alluvial aquifers. Water Resour. Res. 53, 10539–10557 (2017).10.1002/2017WR021619
49. Kaandorp VP Doornenbal PJ Kooi H Peter Broers H de Louw PGB Temperature buffering by groundwater in ecologically valuable lowland streams under current and future climate conditions J. Hydrol. X 2019 3 100031
Kaandorp, V. P., Doornenbal, P. J., Kooi, H., Peter Broers, H. & de Louw, P. G. B. Temperature buffering by groundwater in ecologically valuable lowland streams under current and future climate conditions. J. Hydrol. X 3, 100031 (2019).
50. Versteegh, E. A. A., Vonhof, H. B., Troelstra, S. R., Kaandorp, R. J. G. & Kroon D. Seasonally resolved growth of freshwater bivalves determined by oxygen and carbon isotope shell chemistry. Geochem. Geophys. Geosyst. 11, (2010).
51. Schöne BR Dunca E Fiebig J Pfeiffer M Mutvei’s solution: An ideal agent for resolving microgrowth structures of biogenic carbonates Palaeogeogr. Palaeoclimatol. Palaeoecol. 2005 228 149 166 10.1016/j.palaeo.2005.03.054
Schöne, B. R., Dunca, E., Fiebig, J. & Pfeiffer, M. Mutvei’s solution: An ideal agent for resolving microgrowth structures of biogenic carbonates. Palaeogeogr. Palaeoclimatol. Palaeoecol. 228, 149–166 (2005).10.1016/j.palaeo.2005.03.054
52. W. IAEA. (The GNIP database, 2023).
53. Carey CS Jones JW Hallerman EM Butler RS Determining optimum temperature for growth and survival of laboratory-propagated juvenile freshwater mussels North Am. J. Aquac. 2013 75 532 542 10.1080/15222055.2013.826763
Carey, C. S., Jones, J. W., Hallerman, E. M. & Butler, R. S. Determining optimum temperature for growth and survival of laboratory-propagated juvenile freshwater mussels. North Am. J. Aquac. 75, 532–542 (2013).10.1080/15222055.2013.826763
54. Haag WR Culp JJ McGregor MA Bringolf R Stoeckel JA Growth and survival of juvenile freshwater mussels in streams: Implications for understanding enigmatic mussel declines Freshw. Sci. 2019 38 753 770 10.1086/705919
Haag, W. R., Culp, J. J., McGregor, M. A., Bringolf, R. & Stoeckel, J. A. Growth and survival of juvenile freshwater mussels in streams: Implications for understanding enigmatic mussel declines. Freshw. Sci. 38, 753–770 (2019).10.1086/705919
55. Dettman DL Reische AK Lohmann KC Controls on the stable isotope composition of seasonal growth bands in aragonitic fresh-water bivalves (Unionidae) Geochim. et Cosmochim. Acta 1999 63 1049 1057 10.1016/S0016-7037(99)00020-4
Dettman, D. L., Reische, A. K. & Lohmann, K. C. Controls on the stable isotope composition of seasonal growth bands in aragonitic fresh-water bivalves (Unionidae). Geochim. et Cosmochim. Acta 63, 1049–1057 (1999).10.1016/S0016-7037(99)00020-4
56. Hanlon, S. D. Release of juvenile mussels into a fish hatchery raceway: a comparison of techniques. Thesis, Virginia Tech (2000).
57. DuBose TP Patten MA Holt AS Vaughn CC Latitudinal variation in freshwater mussel potential maximum length in Eastern North America Freshw. Biol. 2022 67 1020 1034 10.1111/fwb.13898
DuBose, T. P., Patten, M. A., Holt, A. S. & Vaughn, C. C. Latitudinal variation in freshwater mussel potential maximum length in Eastern North America. Freshw. Biol. 67, 1020–1034 (2022).10.1111/fwb.13898
58. White DEJ Haag WR McGregor MA Price SJ Effects of food abundance on juvenile freshwater mussel survival and growth in aquaculture, and comparison with growth in streams Aquaculture 2022 560 738473 10.1016/j.aquaculture.2022.738473
White, D. E. J., Haag, W. R., McGregor, M. A. & Price, S. J. Effects of food abundance on juvenile freshwater mussel survival and growth in aquaculture, and comparison with growth in streams. Aquaculture 560, 738473 (2022).10.1016/j.aquaculture.2022.738473
59. Lopez JW Vaughn CC A review and evaluation of the effects of hydrodynamic variables on freshwater mussel communities Freshw. Biol. 2021 66 1665 1679 10.1111/fwb.13784
Lopez, J. W. & Vaughn, C. C. A review and evaluation of the effects of hydrodynamic variables on freshwater mussel communities. Freshw. Biol. 66, 1665–1679 (2021).10.1111/fwb.13784
60. Lutz RA Rhoads DC Anaerobiosis and a theory of growth line formation Science 1977 198 1222 1227 10.1126/science.198.4323.1222 17741686
Lutz, R. A. & Rhoads, D. C. Anaerobiosis and a theory of growth line formation. Science 198, 1222–1227 (1977).17741686 10.1126/science.198.4323.1222
61. Gu R Montgomery S Austin TA Quantifying the effects of stream discharge on summer river temperature Hydrol. Sci. J. 1998 43 885 904 10.1080/02626669809492185
Gu, R., Montgomery, S. & Austin, T. A. Quantifying the effects of stream discharge on summer river temperature. Hydrol. Sci. J. 43, 885–904 (1998).10.1080/02626669809492185
62. Sinokrot BA Gulliver JS In-stream flow impact on river water temperatures J. Hydraul. Res. 2000 38 339 349 10.1080/00221680009498315
Sinokrot, B. A. & Gulliver, J. S. In-stream flow impact on river water temperatures. J. Hydraul. Res. 38, 339–349 (2000).10.1080/00221680009498315
63. Yang ZF Yan Y Liu Q The relationship of streamflow-precipitation-temperature in the Yellow River Basin of China during 1961–2000 Procedia Environ. Sci. 2012 13 2336 2345 10.1016/j.proenv.2012.01.222
Yang, Z. F., Yan, Y. & Liu, Q. The relationship of streamflow-precipitation-temperature in the Yellow River Basin of China during 1961–2000. Procedia Environ. Sci. 13, 2336–2345 (2012).10.1016/j.proenv.2012.01.222
64. Sylvester, R. O. Paper presented at the Proceedings of the twelfth Pacific Northwest symposium on water pollution research, (1963).
65. Doyle RDB Bruce W Lake Whitney comprehensive water quality assessment, Phase 1B-Physical and biological assessment (USDOE) 2009 Baylor University
Doyle, R. D. B. & Bruce, W. Lake Whitney comprehensive water quality assessment, Phase 1B-Physical and biological assessment (USDOE) (Baylor University, 2009).
66. Strause JLAA Freeman L Water quality of Lake Whitney, north-central Texas 1984 U.S. Geological Survey 172
Strause, J. L. A. A. & Freeman, L. Water quality of Lake Whitney, north-central Texas 172 (U.S. Geological Survey, 1984).
67. Dycus J. C., Wisniewski, J. M. & Peterson, J. T. The effects of flow and stream characteristics on the variation in freshwater mussel growth in a Southeast US river basin. Freshw. Biol. 60, (2015).
68. Sansom BJ Hornbach DJ Hove MC Kilgore JS Effects of flow restoration on mussel growth in a wild and scenic North American river Aquat. Biosyst. 2013 9 6 10.1186/2046-9063-9-6 23452382
Sansom, B. J., Hornbach, D. J., Hove, M. C. & Kilgore, J. S. Effects of flow restoration on mussel growth in a wild and scenic North American river. Aquat. Biosyst. 9, 6 (2013).23452382 10.1186/2046-9063-9-6
69. Sotola VA Sullivan KT Littrell BM Martin NH Stich DS Bonner TH Short-term responses of freshwater mussels to floods in a southwestern U.S.A. river estimated using mark–recapture sampling Freshw. Biol. 2021 66 349 361 10.1111/fwb.13642
Sotola, V. A. et al. Short-term responses of freshwater mussels to floods in a southwestern U.S.A. river estimated using mark–recapture sampling. Freshw. Biol. 66, 349–361 (2021).10.1111/fwb.13642
70. Gates KK Vaughn CC Julian JP Developing environmental flow recommendations for freshwater mussels using the biological traits of species guilds Freshw. Biol. 2015 60 620 635 10.1111/fwb.12528
Gates, K. K., Vaughn, C. C. & Julian, J. P. Developing environmental flow recommendations for freshwater mussels using the biological traits of species guilds. Freshw. Biol. 60, 620–635 (2015).10.1111/fwb.12528
71. Campana S Accuracy, precision and quality control in age determination, including a review of the use and abuse of age validation methods J. Fish Biol. 2001 59 197 242 10.1111/j.1095-8649.2001.tb00127.x
Campana, S. Accuracy, precision and quality control in age determination, including a review of the use and abuse of age validation methods. J. Fish Biol. 59, 197–242 (2001).10.1111/j.1095-8649.2001.tb00127.x
72. USGS, Brazos RV at SH nr Bryan TX: TX USGS Water data for the nation. (2023).
73. Meinicke N Ho SL Hannisdal B Nürnberg D Tripati A Schiebel R Meckler AN A robust calibration of the clumped isotopes to temperature relationship for foraminifers Geochim. et Cosmochim. Acta 2020 270 160 183 10.1016/j.gca.2019.11.022
Meinicke, N. et al. A robust calibration of the clumped isotopes to temperature relationship for foraminifers. Geochim. et Cosmochim. Acta 270, 160–183 (2020).10.1016/j.gca.2019.11.022
74. Anderson N Kelson JR Kele S Daëron M Bonifacie M Horita J Mackey TJ John CM Kluge T Petschnig P A unified clumped isotope thermometer calibration (0.5–1,100°C) using carbonate-based standardization Geophys. Res. Lett. 2021 48 e2020GL092069 10.1029/2020GL092069
Anderson, N. et al. A unified clumped isotope thermometer calibration (0.5–1,100°C) using carbonate-based standardization. Geophys. Res. Lett. 48, e2020GL092069 (2021).10.1029/2020GL092069
75. Hudson J Anderson T Ocean temperatures and isotopic compositions through time Earth Environ. Sci. Trans. Royal Soc. Edinb. 1989 80 183 192
Hudson, J. & Anderson, T. Ocean temperatures and isotopic compositions through time. Earth Environ. Sci. Trans. Royal Soc. Edinb. 80, 183–192 (1989).
