Revisiting carbon cycling in the Laurentian Great Lakes following dreissenid mussel invasion

Revisiting carbon cycling in the Laurentian Great Lakes following dreissenid mussel invasion | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Revisiting carbon cycling in the Laurentian Great Lakes following dreissenid mussel invasion Erin D Smith, Leigh J McGaughey, Jerome Marty, Andrea E Kirkwood, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4436844/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Since the active role of inland waters in cycling carbon (C) has been revealed, there has been a renewed interest in calculating C budgets for inland waters to understand their role with respect to global climate change. There is a lack of knowledge regarding C cycling in the Laurentian Great Lakes, the worlds largest freshwater reservoir, with current estimates neglecting the role of invasive species. For one of the most pervasive invaders, dreissenid (zebra and quagga) mussels, research has focused on filter feeding impacts on phosphorus dynamics, but there is a lack of knowledge regarding their role in C cycling, specifically, the impact of the C stored in their slowly degrading shells. As such, we set out to estimate the mass of empty shells and C stored in those shells. We calculated an estimated 1.19 E10 tonnes of empty shell mass currently sitting at the bottom of these lakes, which store approximately 1.43 E9 tonnes of C. This scale of inorganic C storage is comparable to rates of organic C storage in nature-based climate solutions. This work demonstrates the importance of a previously unexplored pathway that dreissenid mussels are altering C cycling in the Laurentian Great Lakes and the thousands of other invaded lakes and rivers. Earth and environmental sciences/Environmental sciences Earth and environmental sciences/Limnology Earth and environmental sciences/Climate sciences/Biogeochemistry carbon cycling dreissenid mussels Laurentian Great Lakes inland carbon budgets mussel shells inorganic carbon deposition Figures Figure 1 Figure 2 Figure 3 Introduction The calculation of up-to-date, accurate carbon (C) budgets is essential for climate change mitigation and adaptation plans 1 , 2 . There has been a large amount of work developing C budgets, especially in terrestrial and oceanic environments, however, more work is required to understand the C budgets of inland waters 1 , 3 , 4 . Inland waters have previously been considered as passive ‘pipes’ transporting terrestrial C to the oceans, however, recently there has been work showing the active role of inland waters in global C cycling 3 , 5 – 7 . Globally, the majority of C in inland waters (streams, rivers, lakes, reservoirs, and ponds) is from terrestrial sources (5.8 Pg C/yr; inorganic and organic C) with a small amount contributed from aquatic primary productivity (0.3 Pg C/yr; organic C), of this C supply, 4.4 Pg C/yr is emitted to the atmosphere, 1.1 Pg C/yr is exported to the oceans, while just 0.6 Pg C/yr is buried in aquatic sediments 7 . Although the annual rate of carbon burial in lakes is a relatively small, sedimentation in inland waters presents a long-term carbon storage pathway 3 , 6 , 8 . Lake size and primary productivity have been identified as key factors determining sedimentation rates, in addition, cultural eutrophication has been shown to increase C sedimentation rates by increasing primary productivity 8 – 11 . [NO_PRINTED_FORM] 11 , estimated that European lakes had an average pre-industrial organic C burial rate of 5–10 g C/m 2 /yr, which increased to an average of 60 g C/m 2 /yr post-1950. In North America, [NO_PRINTED_FORM] 10 calculated an average organic C burial rate of 88 g C/m 2 /yr across eight lakes in the Midwestern region of the United States. In the Great Lakes, Lake Erie and Ontario have higher primary productivity compared to the upper lakes (Huron, Superior, and Michigan), and organic C burial rates mirror this trend (19.07–30 vs. 0.97–4.84 g C/m 2 /yr), despite the smaller size of the downstream lakes 12 . Additionally, in the Great Lakes, dreissenid mussel filtration has been explored as a factor influencing C sedimentation as it consumes large amounts of algae and takes up calcium for shell growth, reducing its concentrations for precipitation 13 – 16 . While organic C burial is increasingly well resolved in freshwater ecosystems, the inorganic C cycle of inland waters has not been as extensively explored, despite the important role that inorganic C plays in inland water C cycling 17 – 20 . In marine systems, a study on high Arctic fjords found inorganic C burial rates ranged between 10.7–45.7 g C/m 2 /yr, with inorganic C burial dominating organic C at one of the two fjords sampled 21 . [NO_PRINTED_FORM] 21 , attributed the dominance of inorganic C burial to a greater activity of calciferous organisms and spread of carbonates in the area. The role of inorganic C burial has also been investigated for coastal blue carbon (CO 2 stored in coastal saltwater environments) and global inorganic C burial rates were estimated at 0.8 Tg C/yr for mangrove ecosystems and 15–62 Tg C/yr in seagrass ecosystems 22 . However, when determining the impact of inorganic C in bottom sediments on C budgets, it is essential to consider the C source. Allochthonous inorganic C originates from terrestrial respiration and is transported to aquatic ecosystems (via runoff and groundwater), where it may be taken up for primary production or degassed to the atmosphere 3 , 22 . Counterintuitively, the burial of inorganic C via mineralization increases CO 2 concentrations as CO 2 is a byproduct of CaCO 3 mineralization (Eq. 1). Additionally, the reverse reaction, dissolution of CaCO 3 , results in an uptake in CO 2 . Consequently, the net inorganic C burial from mineralization offsets the CO 2 sink from organic C burial, in the case of the blue carbon ecosystems described above, a ~ 30% offset was calculated 22 . Thus, it is evident that omitting inorganic C burial in inland water C budgets can result in inaccurate carbon pool evaluations, especially in systems associated with high mineralization from calciferous organisms 3 , 23 . Although the above examples provide insights on inorganic C burial from calciferous organisms on marine C budgets, there is a gap in their consideration for freshwater environments 7 , 12 , 22 . Although the indirect impacts of dreissenid filter feeding have been briefly explored in relation to organic C sedimentation 13 , the impact of shell calcification and burial has not previously been considered for inland C cycling 14 , 15 . Equation 1. \({2HCO}_{3}^{-}+{Ca}^{2+}{\leftrightarrow CaCO}_{3}+{CO}_{2}+{H}_{2}O\) In the above stated studies inorganic C from calcification was considered a source of CO 2 to the environment, however this may not be entirely accurate as it depends on the C source for calcification 24 – 26 . Previously, dissolved inorganic carbon (DIC) was considered the sole source of mussel shell C and thus that C isotopes in shell were thought to directly reflect their environmental conditions, just as stable isotope analysis of oxygen in shells has been used to reflect ambient temperature throughout a mussel’s lifetime 27 . More recently, it has become clear that metabolic C can make up a notable proportion of mussel shell C (%C M ), ~ 12% in marine mussels but up to 25–40% in freshwater mussels 25 , 26 , 28 , 29 . Specifically, a study of the North American freshwater mussel, Elliptio complanata , found an average %C M of 20 ± 10%, the large uncertainty represents the variation due to factors including ontogenic shifts, metabolic rates, and food source and selectivity 28 . Determining the relative proportion of C sources is not only important for paleo work, but also to understand the impact of modern mineralization on C cycling. While shell mineralization from DIC results in the release of a CO 2 molecule, when organic C is the original source a molecule of CO 2 from respiration is diverted from the release to the environment. Thus, the relative proportion of C sources for shell growth in dreissenid mussels need to be determined to understand their impact on inland C cycling (Fig. 1 ). The Laurentian Great Lakes (hereafter Great Lakes) contain ~ 18% of the world’s liquid, surficial freshwater, with a watershed of approximately 1.0 million km 2 30 . Due to cultural eutrophication, most of the budget calculations in the Great Lakes focused on nutrients (nitrogen and phosphorus) rather than on carbon. Preliminary calculation indicates that the carbon budget is dominated by autochthonous production, followed by inorganic and organic C inputs from stream and groundwater and surface runoff 12 . Previous work has also identified the Great Lakes as a net source of CO 2 to the atmosphere, although CO 2 saturation decreases from the upstream to downstream lakes, and varies seasonally 13 , 31 – 33 . Additionally, [NO_PRINTED_FORM] 13 , found an increase in Great Lakes pCO 2 saturation following the establishment of large quagga mussel ( Dreissena bugensis ) populations and hypothesized that this was due to the alteration of lake water biogeochemistry via decreased primary productivity, increased water clarity, and respiration from the large mussel populations. This finding draws from a large body of work examining the physicochemical impacts of invasive dreissenid mussels 34 , 35 , especially in the Great Lakes 36 – 40 . Although the impacts of living dreissenid mussels are well understood in the Great Lakes and other invaded territories, there has been very little work examining the impacts of the accumulating mass of mussel shells on the carbon budget in these lakes 38 , 41 , 42 . There has been some work exploring the scale and role of dreissenid mussel shells in inland lakes, however, as of yet, the Great Lakes, the largest freshwater ecosystem on Earth, have not been considered in this field. In Europe, a study investigated the impact of dreissenid mussel shell accumulation on a shallow lakes’ life span, and found that the invasion of D. bugensis could shorten the life span of the lake by one- to two-thirds by filling it with spent shells 43 . Additionally, [NO_PRINTED_FORM] 44 , investigated shell decay rates for native and invasive bivalves in North American waters and found that Dreissena spp. can produce standing stocks of spent shells on the scale of > 10 kg dry mass/m 2 under appropriate conditions (low flow, high calcium) 44 . [NO_PRINTED_FORM] 45 , collected benthic samples from ten Minnesota lakes and although they found high variability in empty shell mass across lakes, they determined that dead shell material accounted for an average of 322% of lake particulate organic carbon (POC). Finally, the most comprehensive study of dreissenid shell accumulation in a large lake, was conducted in Lake Simcoe, Ontario, Canada (SA = 722 km 2 ) 46 . Based on samples collected from 0–32 m depths, a total dead shell mass of 1 141 000 tonnes was estimated 46 . Additionally, the authors found that approximately 1 000 000 tonnes of shells are produced annually in Lake Simcoe, containing ~ 115 600 tonnes of C. Since their invasion of the Great Lakes in the late 1980s, dreissenid mussel populations have had consistent growth, with an initial growth of the zebra mussel ( D. polymorpha ) population followed by the delayed establishment of quagga mussels 47 – 49 . The establishment of quagga mussels not only increased dreissenid density in previously established zebra mussel habitat, but also expanded to much greater depths (> 90 m) and softer substrates 47 , 50 , 51 . The growing dreissenid mussel population was also able to withstand predation pressure from waterfowl, lake whitefish, and their native predator, the round goby ( Neogobius melanostomus ), which established Great Lake populations in the early 2000s 40 , 50 , 52 . Predation rates vary by habitat, but previous work has estimated round goby predation of dreissenid mussels in Lake Erie to be six kT/yr (270 T C shell /yr), and lake whitefish consumption of 109 and 820 kT/yr (500 and 3800 T C shell /yr) in Lakes Michigan and Huron, respectively 40 , 53 . While dreissenid predation has increased with long-term establishment, there is minimal evidence that these predators are reducing dreissenid abundance, although a shift in size structure has been observed 52 , 54 , 55 . Predation rate and the dreissenid population impact the turnover rate of dreissenid mussels (and production of empty mussel shells), through their impact on the production: biomass (P:B) ratio. The stable dreissenid population maintains biomass (B), while the increased predation rate results in increased dreissenid production (P) (Astorg et al., 2022; Kipp & Ricciardi, 2012; Naddafi & Rudstam, 2014). Therefore, the combination of stable dreissenid populations with increased predation suggests that dreissenid mussel turnover has increased, thus increasing the rate of dreissenid shell mineralization in the Great Lakes (Fig. 2 ). The evidence of large-scale inorganic C sequestration in much smaller North American lakes and a large stable dreissenid population in the Great Lakes makes a clear case for calculating C sequestered in empty mussel shells in the Great Lakes. Altogether it is evident that dreissenids are persistent ecosystem engineers in the Great Lakes 34 , 36 , 50 , however, the impact of nutrient sequestration in their shells that are continuously accumulating at the bottom of these lakes has not previously been considered as a component of their engineering role. Here we provide an estimation of the C mass in empty dreissenid shells in the four dreissenid-invaded Great Lakes (Michigan, Huron, Erie, and Ontario). This quantification provides a starting point for considering another pathway that dreissenid mussels have altered C dynamics in the Great Lakes. Finally, we suggest future work that will better resolve the impacts of this pathway on C budgets. Results To estimate the scale of mineralized C in empty dreissenid mussels in the Great Lakes we calculated the empty shell standing stock based on average pre- (1990–2005) and post (2005–2020)-goby invasion mussel densities from 1990–2020 (Table 2 ). Estimates were split by pre- and post-goby periods to account for the increased dreissenid mussel mortality (Z) following the establishment of stable round goby populations 52 – 54 . Round goby invasion also coincided with quagga mussel dominance which occurred in most lakes by the year 2000 and expanded the mussel population to deeper depths 51 , 56 . This combination of increased mussel abundance (biomass) and predation pressure (mortality rate) resulted in much higher shell turnover and thus higher mineralized C in empty shells (2.40 X 10 8 vs 1.19 X 10 9 T C). This shift was most prominent in Lakes Ontario and Michigan, where increased mussel densities in the deeper waters and increased predation resulted in much higher C mineralization in the post-round goby period 42 , 47 , 57 . An estimated total of 1.19 E10 tonnes of empty shell mass was calculated across the four lakes, containing 1.43 E9 tonnes of C. Table 2 Estimated carbon stored in empty dreissenid shells (in tonnes of carbon) in the Great Lakes before and after the invasion of its native predator, the round goby. Based on average 15 years pre-goby invasion (1990–2005) and 15 years post-goby invasion (2005–2020) mussel densities. Lake Ontario Lake Erie Lake Huron Lake Michigan Total Pre-round goby C storage (T C) 1.15 X 10 7 1.74 X 10 8 1.64 X 10 7 3.80 X 10 7 2.40 X 10 8 Post-round goby C storage (T C) 1.11 X 10 8 9.02 X 10 7 9.83 X 10 7 8.92 X 10 8 1.19 X 10 9 Total C storage (T C) 1.23 X 10 8 2.64 X 10 8 1.15 X 10 8 9.30 X 10 9 1.43 X 10 9 Discussion Dreissenid mussel shells clearly represent a large deviation in carbon flow within the lakes and need to be considered in future carbon budget calculations 12 . This work is a key step in filling the knowledge gap regarding the role of aquatic ecosystems as contributors or sinks of carbon to the atmosphere. 3 , 12 , 58 . Dreissenid mussels have existed in the Great Lakes for just over 30 years and our calculations indicate that their accumulated spent shells account for 1.43 Gt of inorganic C. Throughout time, Canadian lakes are estimate to have stored a total of 20 Gt of C, our calculations indicate that inorganic C from dreissenid shells would make up ~ 7% of this stored C 58 , 59 . On an areal basis, there was a large range of inorganic C deposition rates across the lakes (0.08–37.5 T C/ha/yr), which were primarily dependent on mussel density. This range of C storage is on the scale of estimates of organic C storage from nature-based climate solutions (NBCSs) and depending on the C source (metabolic C or DIC), may supplement or offset these NBCS’s C storage (Fig. 3 ). Based on the current literature, it is likely that DIC is the primary source, indicating that dreissenid shell mineralization will release CO 2 to the environment, offsetting the CO 2 sink from NBCS’s 26 , 28 . Dreissenid mussel inorganic C deposition is also on a similar scale to organic C sequestration rates in established coastal blue ecosystems (seagrass, mangroves, and salt marshes) (Fig. 3 ; Alongi, 2023; Claes et al., 2022). When compared to Canadian freshwater NBCSs, inorganic C is deposited from dreissenid mussel shells in the Great Lakes at a much greater rate than organic C in wetlands and peatlands. Specifically, the inorganic C from shells is deposited at a rate greater than organic C stored in peatlands (23 g/m2/yr), freshwater marshes in the Lake Simcoe watershed (average 112 g C/m 2 /yr), and other intact (62 g C/m 2 /yr) and restored (51–89 g C/m 2 /yr) wetlands in southern Ontario 62 , 63 . The inorganic C deposition of dreissenid shells is also on the scale of inorganic carbon storage calculated in other systems. [NO_PRINTED_FORM] 21 , found inorganic C burial in Arctic fjords at a rate of 10.7–45.7 g C/m 2 /yr, which is within the range of our rate of ~ 0.2–3750 g C/m 2 /yr, with a median of 201 g C/m 2 /yr. On a global scale, seagrass ecosystems have been found to accumulate inorganic C at rates from 149–756 g C/m 2 /yr 22 , 64 – 66 . However, it is important to note that these accumulation rates are based on all inorganic C deposited in bottom sediments (allochthonous and autochthonous) and not specific to calcifying organisms 22 . When looking specifically at inorganic C accumulation from local calcification in seagrass, the range is much lower, 2.2–36 g C/m 2 /yr 64 – 67 . Thus, the inorganic C accumulation from dreissenid mussels in the Great Lakes is at a scale that should be considered in future C sequestration rate and C budget calculations. The estimates provided here are based on historical turnover and shell degradation rates, however, it is essential to consider the impact of predicted future changes to the Great Lakes on dreissenid shell C burial. An updated investigation into shell decay rates in the Great Lakes is recommended to provide a more accurate value for future calculations. 44 found that shell decay rates were correlated with pH and flow, however D. bugensis occupy much deeper habitats than D. polymorpha and other depth related variables (DO, light, nutrients etc.) should be considered for future calculations, especially as climate change is expected to change these conditions 68 . Additionally, winter ice scour can result in high over-winter mussel mortality at shallow depths 69 . Climate change is predicted to alter the ice-cover period in the Great Lakes 70 , 71 , which could reduce the occurrence of ice scour mortality. Future changes to other factors that affect mussel population dynamics (nutrient levels, substrate) should also be considered since they will influence the rate of shell deposition as well. Following their initial establishment, dreissenid populations were relatively consistent in their invaded territories, even following the establishment of their native predator, the round goby 42 , 72 . However, it is possible that new predators may be introduced, or food source availability may change, altering mussel mortality rate 70 , 71 Dreissenid mussels are considered one of the most impactful invasive species in North American waters. Serving as ecosystem engineers, they modify benthic habitats and influence phosphorus cycling and nuisance algae growth 34 , 35 , 70 . We present here another ecosystem-scale impact that expands their role in nutrient cycling to include inorganic C sequestration. On the scale of billions of tonnes of inorganic carbon deposited in slowly degrading shells at the bottom of these lakes. Although our calculations provide an estimate of the scale of their C deposition capacity, physically measuring the mass of shells in the Great Lakes and connecting channels, as [NO_PRINTED_FORM] 46 did in Lake Simcoe, will provide the best estimation C stored in empty dreissenid shells. An essential part of this work will be determining the relative proportions of each shell C source (metabolic and DIC), this will provide the insight needed to determine the role of dreissenids on C cycling and CO 2 flux in the Great Lakes and other invaded territories. The storage of other elements, such as nitrogen and phosphorus, in dreissenid mussel shells should also be calculated. Although the proportions of shell nitrogen (~ 0.25%) and phosphorus (~ 0.007%) are much smaller than carbon (~ 12%) 37 , 45 , 46 , based on the shear scale of deposited shells they may be an important sink for these nutrients as well, especially when considering that phosphorus concentrations have declined in most of these lakes 70 , 73 . Finally, the creation of shell fragments from predation needs to be considered. Previous calculations of shell decay rates are based on whole shells suspended in the water column, however, most dreissenid mussel predators break the shells into fragments to consume the mussel tissue and discard shell fragments that settle into the sediments 54 . As described in the Methods, it is not clear how these small shell fragments will impact dissolution rates, they could reduce dissolution by increasing shell compaction (reducing exposed surface area) or as smaller particles they could dissolve more quickly 44 . If dissolution does increase it would result in CO 2 uptake from the water column through the reverse of Eq. 1. Evidently, determining the rate of shell dissolution with shell fragments will provide a better understanding of dreissenid mussel impacts on C cycling in the Great Lakes. This work adds to the previously explored roles of dreissenid mussels as ecosystem engineers in the Great Lakes and other invaded waterways. In addition to their phosphorus cycling role and increased respiration rates, dreissenid mussels are altering the flow of carbon through mineralization on a massive scale 13 , 35 , 36 . As mentioned above, determining the relative role of metabolic C vs. DIC in shell mineralization is key to deciphering their exact impact on C cycling and CO 2 flux in invaded waters. Previous work indicates that ambient bicarbonate is likely the primary source, and thus until the relative proportions are determined we should consider dreissenid shell mineralization as a likely source of CO 2 to invaded waters 28 . While dreissenid mussels have been found to alter CO 2 concentrations through respiration and impacts on primary productivity, their shell mineralization has not been considered as a component of this impact 13 , 74 . Regardless of the direction (sink or source) it is evident that shell mineralization has major impacts on C cycling in the Great Lakes. On a larger scale these findings have implications for calculation of C budgets for the Great Lakes and implications for the global inland waters C balance. A preliminary C budget for the Great Lakes indicates that net primary production (NPP) and respiration are the greatest C flows for all lakes except Ontario, which receives more C from upstream inputs 12 . Our calculations indicate that dreissenid mineralization alters C flow on the same scale as the processes considered in 12 calculations and should be included in future C budgets for the invaded lakes (Michigan, Huron, Erie, and Ontario). Additionally, these C fluxes should be included in global inland water C budgets, although we focused on the Great Lakes, there are thousands of invaded lakes and rivers across North America and Europe that will have altered C flows due to dreissenid mussel mineralization 34 , 35 . Finally, inclusion in these C budgets will have policy implications in the corresponding jurisdictions. Many locations have drafted adaptive management plans to minimize greenhouse gas emissions and an improved understanding of local C dynamics will make these plans more robust as dreissenid mussel populations continue to grow. Methods Calculation of dreissenid empty shell standing stock was based on 44 , where at steady state the standing stock of spent shells (SS) is equal to the production of spent shells (P) divided by the instantaneous decay rate (I). Production of spent shells was determined from living shell mass (B) and instantaneous mortality rate (Z, or P:B depending on availability), with Eq. 1. Equation 1. $$SS=\frac{B*Z}{I}$$ The instantaneous shell decay rate of 0.058 calculated by 46 was used as it is comparable to other estimates 44 , and was calculated in Lake Simcoe, a lake in the same region as the Laurentian Great Lakes, with similar geochemical characteristics. This decay rate is a very conservative estimate as it is based on shells suspended in the water column, where turbulence is greater than on the lake bottom, which would result in an overestimation of the dissolution rate. Additionally, the settlement of shells in layers on the lake bottom would reduce exposed surface area and likely result in slower decay compared to shells in the water column. Finally, the increased shell fragmentation from round goby predation could impact shell decay rate, however without published findings on this impact it is not possible to adjust decay rate for this factor. Shell fragmentation would alter the matrix of deposited shells potentially reducing decay rate by reducing exposed surface area, alternatively 44 , showed that smaller shells result in higher decay rates. Considering this contradiction, the published value was used to avoid making an incorrect assumption due to the lack of literature on this factor. Two mortality rates were used, one for the pre-round goby period (~ 1990–2005) and one for the post-round goby period (~ 2005–2020). The pre-round goby, non-predatory mortality rate of 1.26/year from 75 . The post-round goby production: biomass ratio (P: B)/mortality rate of 2.1 from 76 was applied for the post-round goby invasion period ~ 2005–2020. An average living shell biomass was determined for each period based on published density values for each depth zone (< 30 m, 31–50 m, 51–90 m) for Lake Ontario, Lake Huron, and Lake Michigan. As Lake Erie is very shallow and has more variation across its basins (East, Central and West) than depth, shell biomass was determined from mean basin-wide biomass. Lake Superior was not included in calculations as it’s dreissenid populations are limited to select areas of the lake, and extrapolating to a lake-wide estimation based on depth zones would not be appropriate 50 , 77 . Living dreissenid density data was extracted from lake-wide surveys with at least 30 sites, across a range of bottom substrates to represent the range of conditions in each lake (Table 1 ). Additionally, where a series of similar surveys were conducted the most recent published data was extracted. All surveys collected between 2–6 replicates for each sample site. Table 1 Summary of dreissenid density data sources for each lake and period (pre- and post-round goby invasion). Lake Period Source Year(s) sampled Sample sites (#) Depth range (m) Lake Ontario Pre-goby 83 1997 75 14–209 Post-goby 42 2008, 2013, 2018 55 11–209 Lake Erie Pre-goby 78 1998 30 2.0–58.5 Post-goby 84 2009, 2012 135 2–66 Lake Huron Pre-goby 48 2000, 2003 65–85 19–173 Post-goby 48 2007, 2012 80–83 12–171 Lake Michigan Pre-goby 49 1994,1995,2000 90–157 6–208 Post-goby 49 2005, 2010 140–160 13–207 Living mussel shell density was calculated from the wet mass density (g/m 2 ) assuming shells are approximately 38.13% of whole mussel wet weight (Báldi et al., 2019; Table S4), if whole mussel wet weight was not available it was calculated based on density of individuals (ind./m 2 ) using an average wet weight of 0.156 g/ind. 16 , 78 . Empty shell standing stock was multiplied by 0.12, the approximate proportion of carbon in the shell 37 , 45 , to get the areal mass of carbon in empty shells deposited each year. Lake surface area for each depth zone described above was calculated from NOAA bathymetry maps 79 – 82 . Standing stock of empty dreissenid shell C deposited per year was calculated by combining areal production with each zone’s surface area. Finally, the total mass of C stored in empty dreissenid shells was calculated by multiplying the pre- and post-round goby estimates by 15 years, the approximate time for each period. Declarations Competing Interests The authors declare no competing interests. Author Contribution EDS: Data curation, formal analysis, writing – original draft preparation. LJM: Conceptualization, supervision, writing – review & editing. JM: Writing – review & editing. AEK: Writing – review & editing. JR: Supervision, Writing – review & editing. Acknowledgement We would like to thank all members of our team for their efforts. The conceptual models were created by Lexy Harquail. We would also like to thank our funding agencies: Mitacs and RBC Foundation. Data Availability All data generated or analysed during this study are included in this published article [and its supplementary information files]. References Tranvik, L. J. et al. Lakes and reservoirs as regulators of carbon cycling and climate. Limnol Oceanogr 54, 2298–2314 (2009). Regnier, P. et al. Anthropogenic perturbation of the carbon fluxes from land to ocean. Nat Geosci 6, 597–607 (2013). Cole, J. J. et al. Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems 10, 171–184 (2007). Regnier, P., Resplandy, L., Najjar, R. G. & Ciais, P. The land-to-ocean loops of the global carbon cycle. Nature 603, 401–410 (2022). Raymond, P. A. et al. Global carbon dioxide emissions from inland waters. Nature 503, 355–359 (2013). Webb, J. R., Santos, I. R., Maher, D. T. & Finlay, K. The Importance of Aquatic Carbon Fluxes in Net Ecosystem Carbon Budgets: A Catchment-Scale Review. Ecosystems 22, 508–527 (2019). Pilla, R. M. et al. Anthropogenically driven climate and landscape change effects on inland water carbon dynamics: What have we learned and where are we going? Global Change Biology vol. 28 5601–5629 Preprint at https://doi.org/10.1111/gcb.16324 (2022). Mulholland, P. J. & Elwood, J. W. The role of lake and reservoir sediments as sinks in the perturbed global carbon cycle. Tellus 34, 490–499 (1982). Mendonça, R. et al. Organic carbon burial in global lakes and reservoirs. Nat Commun 8, (2017). Heathcote, A. J. & Downing, J. A. Impacts of Eutrophication on Carbon Burial in Freshwater Lakes in an Intensively Agricultural Landscape. Ecosystems 15, 60–70 (2012). Anderson, N. J., Bennion, H. & Lotter, A. F. Lake eutrophication and its implications for organic carbon sequestration in Europe. Glob Chang Biol 20, 2741–2751 (2014). McKinley, G., Urban, N., Bennington, V., Pilcher, D. & Mcdonald, C. Preliminary Carbon Budgets for the Laurentian Great Lakes . vol. 4 (2011). Lin, P. & Guo, L. Do invasive quagga mussels alter CO2 dynamics in the Laurentian Great Lakes? Sci Rep 6, (2016). Lu, Y., Meyers, P. A., Eadie, B. J. & Robbins, J. A. Carbon cycling in Lake Erie during cultural eutrophication over the last century inferred from the stable carbon isotope composition of sediments. J Paleolimnol 43, 261–272 (2010). Hodell, D. A., Schelske, C. L., Fahnenstiel, G. L. & Robbins, L. L. Biologically induced calcite and its isotopic composition in Lake Ontario. Limnol Oceanogr 43, 187–199 (1998). Barbiero, R. P., Tuchman, M. L. & Millard, E. S. Post-dreissenid increases in transparency during summer stratification in the offshore waters of Lake Ontario: Is a reduction in whiting events the cause? J Great Lakes Res 32, 131–141 (2006). Martin, J. B. Carbonate minerals in the global carbon cycle. Chem Geol 449, 58–72 (2017). Butman, D. et al. Aquatic carbon cycling in the conterminous United States and implications for terrestrial carbon accounting. PNAS 113, 58–63 (2016). Liu, Z., Dreybrodt, W. & Wang, H. A new direction in effective accounting for the atmospheric CO2 budget: Considering the combined action of carbonate dissolution, the global water cycle and photosynthetic uptake of DIC by aquatic organisms. Earth-Science Reviews vol. 99 162–172 Preprint at https://doi.org/10.1016/j.earscirev.2010.03.001 (2010). Minor, E. C. & Brinkley, G. Alkalinity, pH, and pCO2 in the Laurentian Great Lakes: An initial view of seasonal and inter-annual trends. J Great Lakes Res 48, 502–511 (2022). Koziorowska, K., Kuliński, K. & Pempkowiak, J. Comparison of the burial rate estimation methods of organic and inorganic carbon and quantification of carbon burial in two high Arctic fjords. Oceanologia 60, 405–418 (2018). Saderne, V. et al. Role of carbonate burial in Blue Carbon budgets. Nat Commun 10, (2019). Chapin, F. S. et al. Reconciling carbon-cycle concepts, terminology, and methods. Ecosystems 9, 1041–1050 (2006). Thomsen, J., Haynert, K., Wegner, K. M. & Melzner, F. Impact of seawater carbonate chemistry on the calcification of marine bivalves. Biogeosciences 12, 4209–4220 (2015). Gillikin, D. P., Lorrain, A., Meng, L. & Dehairs, F. A large metabolic carbon contribution to the δ13C record in marine aragonitic bivalve shells. Geochim Cosmochim Acta 71, 2936–2946 (2007). Pérez, C. A., Lagos, N. A., Saldías, G. S., Waldbusser, G. & Vargas, C. A. Riverine discharges impact physiological traits and carbon sources for shell carbonate in the marine intertidal mussel Perumytilus purpuratus. Limnol Oceanogr 61, 969–983 (2016). McConnaughey, T. A., Bljrdeti, J., Whelan, J. F. & Paull, C. K. Carbon isotopes in biological carbonates: Respiration and photosynthesis. Geochim Cosmochim Acta 61, (1997). Graniero, L. E., Gillikin, D. P. & Surge, D. Evaluating Freshwater Mussel Shell δ13C Values as a Proxy for Dissolved Inorganic Carbon δ13C Values in a Temperate River. J Geophys Res Biogeosci 126, (2021). Poulain, C. et al. Experimental shift of diet and DIC stable carbon isotopes: Influence on shell δ13C values in the Manila clam Ruditapes philippinarum. Chem Geol 272, 75–82 (2010). Waples, J. et al. The Laurentian Great Lakes. (2008). Phillips, J. C. et al. The potential for CO2-induced acidification in freshwater: A great lakes case study. Oceanography 28, 136–145 (2015). Karim, A., Dubois, K. & Veizer, J. Carbon and oxygen dynamics in the Laurentian Great Lakes: Implications for the CO2 flux from terrestrial aquatic systems to the atmosphere. Chem Geol 281, 133–141 (2011). Urban, N. R. & Desai, A. Are the Great Lakes a significant net source or sink of CO2? Verh. Internat. Verein. Limnol 30, 1283–1288 (2009). Li, J. et al. Benthic invaders control the phosphorus cycle in the world’s largest freshwater ecosystem. PNAS 118, 1–9 (2020). Higgins, S. N. & Vander Zanden, M. J. What a difference a species makes: Meta-analysis of dreissenid mussel impacts on freshwater ecosystems. Ecol Monogr 80, 179–196 (2010). Hecky, R. E., Smith, R. E. H., Barton, D. R., Guildford, S. J. & al, et. The nearshore phosphorus shunt: a consequence of ecosystem engineering by dreissenids in the Laurentian Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 61, 1285–1293 (2004). Arnott, D. L. & Vanni, M. J. Nitrogen and phosphorus recycling by the zebra mussel (Dreissena polymorpha) in the western basin of Lake Erie. Canadian Journal of Fisheries and Aquatic Sciences 53, 646–659 (1996). DeRoy, E. M. & MacIsaac, H. J. Impacts of invasive species in the Laurentian Great Lakes. in Handbook of Environmental Chemistry vol. 101 135–156 (Springer Science and Business Media Deutschland GmbH, 2020). Chase, M. E. & Bailey, R. C. The ecology of the zebra mussel (Dreissena polymorpha) in the lower Great Lakes of North America: II. Total production, energy allocation, and reproductive effort. J Great Lakes Res 25, 122–134 (1999). Madenjian, C. P. et al. Dreissenid mussels are not a ‘dead end’ in Great Lakes food webs. J Great Lakes Res 36, 73–77 (2010). Burlakova, L. E., Karatayev, A. Y., Boltovskoy, D. & Correa, N. M. Ecosystem services provided by the exotic bivalves Dreissena polymorpha, D. rostriformis bugensis, and Limnoperna fortunei. Hydrobiologia 850, 2811–2854 (2023). Karatayev, A. Y. et al. Dreissena in Lake Ontario 30 years post-invasion. J Great Lakes Res 48, 264–273 (2022). Báldi, K. et al. Sediment contributing invasive dreissenid species in a calcareous shallow lake – Possible implications for shortening life span of lakes by filling. Elementa 7, (2019). Strayer, D. L. & Malcom, H. M. Shell decay rates of native and alien freshwater bivalves and implications for habitat engineering. Freshw Biol 52, 1611–1617 (2007). Williamson, F. & Ozersky, T. Lake Characteristics, Population Properties and Invasion History Determine Impact of Invasive Bivalves on Lake Nutrient Dynamics. Ecosystems 22, 1721–1735 (2019). Ozersky, T., Evans, D. O. & Ginn, B. K. Invasive mussels modify the cycling, storage and distribution of nutrients and carbon in a large lake. Freshw Biol 60, 827–843 (2015). Karatayev, A. Y. & Burlakova, L. E. Dreissena in the Great Lakes: What have we learned in 30 years of invasion. Hydrobiologia (2022) doi: 10.1007/s10750-022-04990-x . Nalepa, T. F., Riseng, C. M., Elgin, A. K. & Lang, G. A. Abundance and Distribution of Benthic Macroinvertebrates in the Lake Huron System: Saginaw Bay, 2006–2009, and Lake Huron, Including Georgian Bay and North Channel, 2007 and 2012 . (2018). Nalepa, T. F., Fanslow, D. L., Lang, G. A., Mabrey, K. & Rowe, M. Lake-Wide Benthic Surveys in Lake Michigan in 1994-95, 2000, 2005, and 2010: Abundances of the Amphipod Diporeia Spp. and Abundances and Biomass of the Mussels Dreissena Polymorpha and Dreissena Rostriformis Bugensis . . (2014). Vanderploeg, H. A. et al. Dispersal and emerging ecological impacts of Ponto-Caspian species in the Laurentian Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 59, 1209–1228 (2002). Eifert, R. A. et al. Could quagga mussels impact offshore benthic community and surface sediment-bound nutrients in the Laurentian Great Lakes? Hydrobiologia (2023) doi: 10.1007/s10750-023-05191-w . Kipp, R. & Ricciardi, A. Impacts of the Eurasian round goby (Neogobius melanostomus) on benthic communities in the upper St. Lawrence River. Canadian Journal of Fisheries and Aquatic Sciences 69, 469–486 (2012). Bunnell, D. B., Johnson, T. B. & Knight, C. T. The impact of introduced round gobies (Neogobius melanostomus) on phosphorus cycling in central Lake Erie. Canadian Journal of Fisheries and Aquatic Sciences 62, 15–29 (2005). Naddafi, R. & Rudstam, L. G. Predation on invasive zebra mussel, Dreissena polymorpha, by pumpkinseed sunfish, rusty crayfish, and round goby. Hydrobiologia 721, 107–115 (2014). Astorg, L., Charette, C., Windle, M. J. S. & Derry, A. M. Two decades since first invasion: Revisiting round goby impacts on nearshore aquatic communities in the Upper St. Lawrence River. J Great Lakes Res (2022) doi: 10.1016/j.jglr.2022.01.017 . Karatayev, A. Y. et al. Lake morphometry determines Dreissena invasion dynamics. Biol Invasions 23, 2489–2514 (2021). Nalepa, T. F., Fanslow, D. L. & Lang, G. A. Transformation of the offshore benthic community in Lake Michigan: Recent shift from the native amphipod Diporeia spp. to the invasive mussel Dreissena rostriformis bugensis. Freshw Biol 54, 466–479 (2009). Council of Canadian Academies (CCA). Inland Freshwater Ecosystems . Nature-Based Climate Solutions (2022). Ferland, M. E., Del Giorgio, P. A., Teodoru, C. R. & Prairie, Y. T. Long-term C accumulation and total C stocks in boreal lakes in northern Québec. Global Biogeochem Cycles 26, (2012). Claes, J., Hopman, D., Jaeger, G. & Rogers, M. Blue Carbon: The Potential of Coastal and Oceanic Climate Action Nature-Based Climate Solutions in the World’s Oceans Can Play an Important Role in Conservation and Carbon Abatement Efforts Worldwide . (2022). Alongi, D. M. Current status and emerging perspectives of coastal blue carbon ecosystems. Carbon Footprints 2, (2023). Pendea, I. F., Kanavillil, N., Kurissery, S. & Chmura, G. L. Carbon Stocks and Recent Rates of Carbon Sequestration in Nutrient-Rich Freshwater Wetlands from Lake Simcoe Watershed (Southern Canada). J Geophys Res Biogeosci 128, (2023). Creed, I. F. et al. Can Restoration of Freshwater Mineral Soil Wetlands Deliver Nature-Based Climate Solutions to Agricultural Landscapes? Front Ecol Evol 10, (2022). Saderne, V. et al. Accumulation of Carbonates Contributes to Coastal Vegetated Ecosystems Keeping Pace with Sea Level Rise in an Arid Region (Arabian Peninsula). J Geophys Res Biogeosci 123, 1498–1510 (2018). Mazarrasa, I. et al. Seagrass meadows as a globally significant carbonate reservoir. Biogeosciences 12, 4993–5003 (2015). Yates, K. K. & Halley, R. B. Diurnal variation in rates of calcification and carbonate sediment dissolution in Florida Bay. Estuaries and Coasts 29, 24–39 (2006). Barron, C., Duarte, C. M., Frankignoulle, M. & Vieira Borges, A. Organic Carbon Metabolism and Carbonate Dynamics in a Mediterranean Seagrass (Posidonia oceanica) Meadow. Estuaries and Coasts 29, 417–426 (2001). Woolway, R. I., Sharma, S. & Smol, J. P. Lakes in hot water: The impacts of a changing climate on aquatic ecosystems. Bioscience 72, 1050–1061 (2022). Chase, M. E. & Bailey, R. C. The ecology of the zebra mussel (Dreissena polymorpha) in the lower Great Lakes of North America: I. Population dynamics and growth. J Great Lakes Res 25, 107–121 (1999). Mahdiyan, O., Filazzola, A., Molot, L. A., Gray, D. & Sharma, S. Drivers of water quality changes within the Laurentian Great Lakes region over the past 40 years. Limnol Oceanogr 66, 237–254 (2021). Sharma, S. et al. Widespread loss of lake ice around the Northern Hemisphere in a warming world. Nature Climate Change vol. 9 227–231 Preprint at https://doi.org/10.1038/s41558-018-0393-5 (2019). Strayer, D. L. et al. Long-term population dynamics of dreissenid mussels (Dreissena polymorpha and D. rostriformis): a cross-system analysis. Ecosphere 10, (2019). Dove, A. & Chapra, S. C. Long-term trends of nutrients and trophic response variables for the Great Lakes. Limnol Oceanogr 60, 696–721 (2015). Blagrave, K. et al. Spatial heterogeneity in water quality across the northern nearshore regions of the Laurentian Great Lakes. J Great Lakes Res 49, (2023). Boles, L. C. & Lipcius, R. N. Potential for Population Regulation of the Zebra Mussel by Finfish and the Blue Crab in North American Estuaries (1997). J Shellfish Res 16, 179–186 (1997). Kao, Y. C., Adlerstein, S. & Rutherford, E. The relative impacts of nutrient loads and invasive species on a Great Lakes food web: An Ecopath with Ecosim analysis. J Great Lakes Res 40, 35–52 (2014). Grigorovich, I. A. et al. Lake Superior: An Invasion Coldspot? Hydrobiologia vol. 499 (2003). Jarvis, P., Dow, J., Dermott, R. & Bonnell, R. Zebra (Dreissena Polymorpha) and Quagga Mussel (Dreissena Bugensis) Distribution and Density in Lake Erie, 1992–1998 . (2000). NOAA National Geophysical Data Center. Bathymetry of Lake Erie and Lake St. Clair. NOAA National Centers for Environmental Information (1999). NOAA National Geophysical Data Center. Bathymetry of Lake Huron. NOAA National Centers for Environmental Information (1999). NOAA National Geophysical Data Center. Bathymetry of Lake Michigan. NOAA National Centers for Environmental Information (1996). NOAA National Geophysical Data Center. Bathymetry of Lake Ontario. NOAA National Centers for Environmental Information (1999). Lozano, S. J., Scharold, J. V. & Nalepa, T. F. Recent declines in benthic macroinvertebrate densities in Lake Ontario. Canadian Journal of Fisheries and Aquatic Science 58, 518–529 (2001). Karatayev, A. Y. et al. Twenty five years of changes in Dreissena spp. populations in Lake Erie. J Great Lakes Res 40, 550–559 (2014). Additional Declarations No competing interests reported. Supplementary Files CSinkdata.csv Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4436844","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":308826766,"identity":"c8182999-c815-4037-89b6-893eddf2a283","order_by":0,"name":"Erin D Smith","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA/ElEQVRIiWNgGAWjYJCCA4wNNgxsIFYCELMRpeVgQxpciwQbGwNjA0E9BxsOw9kSDIS08POfTjz8ccf5fD6J5GMfHu6oq+OT7zF/8IHBTh6XFskZuRsOHDxz27JNIi15RuKZw0CH8Rg2zmBINsRllcENXqCWttsGbDxnjBkS2w6AtTTzgMIEhxb782dBWs4BtZz/DNRSB9Hyh+GAPU5bGEAOaztgwMbewwzUwgzRAgzGRFxaJG4AtZxtSwZqaQM57LBkG1ta4cweg+RkXFr4+89u/lDZZmcg38z8mPFnWx2/fPPhDR9+VNjZ4tKCCxiQqH4UjIJRMApGAQoAAK05V7FYrWPiAAAAAElFTkSuQmCC","orcid":"","institution":"St. Lawrence River Institute","correspondingAuthor":true,"prefix":"","firstName":"Erin","middleName":"D","lastName":"Smith","suffix":""},{"id":308826767,"identity":"1aa7f4a2-8794-43e9-ac88-c1af0c213396","order_by":1,"name":"Leigh J McGaughey","email":"","orcid":"","institution":"St. Lawrence River Institute","correspondingAuthor":false,"prefix":"","firstName":"Leigh","middleName":"J","lastName":"McGaughey","suffix":""},{"id":308826768,"identity":"a4443274-95c2-4766-ad99-7578bbcb503b","order_by":2,"name":"Jerome Marty","email":"","orcid":"","institution":"International Association for Great Lakes Research","correspondingAuthor":false,"prefix":"","firstName":"Jerome","middleName":"","lastName":"Marty","suffix":""},{"id":308826769,"identity":"114295af-9bec-4adc-ad34-9a00edc35d3d","order_by":3,"name":"Andrea E Kirkwood","email":"","orcid":"","institution":"Ontario Tech University","correspondingAuthor":false,"prefix":"","firstName":"Andrea","middleName":"E","lastName":"Kirkwood","suffix":""},{"id":308826770,"identity":"bf6e5d9d-176a-4d21-8d7e-3d0aa2e98a99","order_by":4,"name":"Jeff Ridal","email":"","orcid":"","institution":"St. Lawrence River Institute","correspondingAuthor":false,"prefix":"","firstName":"Jeff","middleName":"","lastName":"Ridal","suffix":""}],"badges":[],"createdAt":"2024-05-17 12:51:03","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4436844/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4436844/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":57518154,"identity":"df91151f-f87a-4911-9be8-29dfabdb7577","added_by":"auto","created_at":"2024-05-31 20:29:29","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":348419,"visible":true,"origin":"","legend":"\u003cp\u003eConceptual model of dreissenid mussel mineralization and dissolution impacts on carbon cycling in the Great Lakes. Figure created by Lexy Harquail.\u003c/p\u003e","description":"","filename":"Figure1MusselCflowLH.png","url":"https://assets-eu.researchsquare.com/files/rs-4436844/v1/ea62fe17d6428e3f6a8df6b8.png"},{"id":57518759,"identity":"f9f40634-b8c9-43ef-a368-06d167173785","added_by":"auto","created_at":"2024-05-31 20:37:29","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":197555,"visible":true,"origin":"","legend":"\u003cp\u003eConceptual model of dreissenid mussel shell deposition A) in the pre-round goby invasion and B) post-round goby invasion period. Figure created by Lexy Harquail.\u003c/p\u003e","description":"","filename":"Figure2Pre.PostGobyLH.png","url":"https://assets-eu.researchsquare.com/files/rs-4436844/v1/1965f8b1ca4a96b9ad44cd89.png"},{"id":57518156,"identity":"003e8a52-5c24-436b-b0d7-44b54ec9da88","added_by":"auto","created_at":"2024-05-31 20:29:29","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":117860,"visible":true,"origin":"","legend":"\u003cp\u003eAverage carbon sequestration rates for various marine (blue) and freshwater (green) nature-based climate solutions and dreissenid shell carbon sequestration if sourced entirely from metabolic carbon or ambient dissolved inorganic carbon.\u003c/p\u003e","description":"","filename":"Figure3CSequestrationPlot.png","url":"https://assets-eu.researchsquare.com/files/rs-4436844/v1/78d6f3ae980d453ce5a43509.png"},{"id":57586159,"identity":"2a13c428-6dd2-4628-aa7a-847eb91448fa","added_by":"auto","created_at":"2024-06-03 03:29:59","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":983641,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4436844/v1/65bd3333-8682-4376-97ce-88e29e2b102b.pdf"},{"id":57518157,"identity":"ce79cbea-19a8-42fa-97cd-1c7517e0082e","added_by":"auto","created_at":"2024-05-31 20:29:30","extension":"csv","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":9618,"visible":true,"origin":"","legend":"","description":"","filename":"CSinkdata.csv","url":"https://assets-eu.researchsquare.com/files/rs-4436844/v1/e9a6ec01dec119836a6ce351.csv"}],"financialInterests":"No competing interests reported.","formattedTitle":"Revisiting carbon cycling in the Laurentian Great Lakes following dreissenid mussel invasion","fulltext":[{"header":"Introduction","content":"\u003cp\u003eThe calculation of up-to-date, accurate carbon (C) budgets is essential for climate change mitigation and adaptation plans \u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. There has been a large amount of work developing C budgets, especially in terrestrial and oceanic environments, however, more work is required to understand the C budgets of inland waters \u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u003c/sup\u003e. Inland waters have previously been considered as passive \u0026lsquo;pipes\u0026rsquo; transporting terrestrial C to the oceans, however, recently there has been work showing the active role of inland waters in global C cycling \u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan additionalcitationids=\"CR6\" citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e. Globally, the majority of C in inland waters (streams, rivers, lakes, reservoirs, and ponds) is from terrestrial sources (5.8 Pg C/yr; inorganic and organic C) with a small amount contributed from aquatic primary productivity (0.3 Pg C/yr; organic C), of this C supply, 4.4 Pg C/yr is emitted to the atmosphere, 1.1 Pg C/yr is exported to the oceans, while just 0.6 Pg C/yr is buried in aquatic sediments \u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eAlthough the annual rate of carbon burial in lakes is a relatively small, sedimentation in inland waters presents a long-term carbon storage pathway \u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e,\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. Lake size and primary productivity have been identified as key factors determining sedimentation rates, in addition, cultural eutrophication has been shown to increase C sedimentation rates by increasing primary productivity \u003csup\u003e\u003cspan additionalcitationids=\"CR9 CR10\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. \u003csup\u003e[NO_PRINTED_FORM] 11\u003c/sup\u003e, estimated that European lakes had an average pre-industrial organic C burial rate of 5\u0026ndash;10 g C/m\u003csup\u003e2\u003c/sup\u003e/yr, which increased to an average of 60 g C/m\u003csup\u003e2\u003c/sup\u003e/yr post-1950. In North America, \u003csup\u003e[NO_PRINTED_FORM] 10\u003c/sup\u003e calculated an average organic C burial rate of 88 g C/m\u003csup\u003e2\u003c/sup\u003e/yr across eight lakes in the Midwestern region of the United States. In the Great Lakes, Lake Erie and Ontario have higher primary productivity compared to the upper lakes (Huron, Superior, and Michigan), and organic C burial rates mirror this trend (19.07\u0026ndash;30 vs. 0.97\u0026ndash;4.84 g C/m\u003csup\u003e2\u003c/sup\u003e/yr), despite the smaller size of the downstream lakes \u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e. Additionally, in the Great Lakes, dreissenid mussel filtration has been explored as a factor influencing C sedimentation as it consumes large amounts of algae and takes up calcium for shell growth, reducing its concentrations for precipitation \u003csup\u003e\u003cspan additionalcitationids=\"CR14 CR15\" citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eWhile organic C burial is increasingly well resolved in freshwater ecosystems, the inorganic C cycle of inland waters has not been as extensively explored, despite the important role that inorganic C plays in inland water C cycling \u003csup\u003e\u003cspan additionalcitationids=\"CR18 CR19\" citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e. In marine systems, a study on high Arctic fjords found inorganic C burial rates ranged between 10.7\u0026ndash;45.7 g C/m\u003csup\u003e2\u003c/sup\u003e/yr, with inorganic C burial dominating organic C at one of the two fjords sampled \u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. \u003csup\u003e[NO_PRINTED_FORM] 21\u003c/sup\u003e, attributed the dominance of inorganic C burial to a greater activity of calciferous organisms and spread of carbonates in the area. The role of inorganic C burial has also been investigated for coastal blue carbon (CO\u003csub\u003e2\u003c/sub\u003e stored in coastal saltwater environments) and global inorganic C burial rates were estimated at 0.8 Tg C/yr for mangrove ecosystems and 15\u0026ndash;62 Tg C/yr in seagrass ecosystems \u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e. However, when determining the impact of inorganic C in bottom sediments on C budgets, it is essential to consider the C source. Allochthonous inorganic C originates from terrestrial respiration and is transported to aquatic ecosystems (via runoff and groundwater), where it may be taken up for primary production or degassed to the atmosphere \u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e. Counterintuitively, the burial of inorganic C via mineralization increases CO\u003csub\u003e2\u003c/sub\u003e concentrations as CO\u003csub\u003e2\u003c/sub\u003e is a byproduct of CaCO\u003csub\u003e3\u003c/sub\u003e mineralization (Eq.\u0026nbsp;1). Additionally, the reverse reaction, dissolution of CaCO\u003csub\u003e3\u003c/sub\u003e, results in an uptake in CO\u003csub\u003e2\u003c/sub\u003e. Consequently, the net inorganic C burial from mineralization offsets the CO\u003csub\u003e2\u003c/sub\u003e sink from organic C burial, in the case of the blue carbon ecosystems described above, a\u0026thinsp;~\u0026thinsp;30% offset was calculated \u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e. Thus, it is evident that omitting inorganic C burial in inland water C budgets can result in inaccurate carbon pool evaluations, especially in systems associated with high mineralization from calciferous organisms \u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. Although the above examples provide insights on inorganic C burial from calciferous organisms on marine C budgets, there is a gap in their consideration for freshwater environments \u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e,\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eAlthough the indirect impacts of dreissenid filter feeding have been briefly explored in relation to organic C sedimentation \u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e, the impact of shell calcification and burial has not previously been considered for inland C cycling \u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e,\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eEquation 1. \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({2HCO}_{3}^{-}+{Ca}^{2+}{\\leftrightarrow CaCO}_{3}+{CO}_{2}+{H}_{2}O\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003cp\u003eIn the above stated studies inorganic C from calcification was considered a source of CO\u003csub\u003e2\u003c/sub\u003e to the environment, however this may not be entirely accurate as it depends on the C source for calcification \u003csup\u003e\u003cspan additionalcitationids=\"CR25\" citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e. Previously, dissolved inorganic carbon (DIC) was considered the sole source of mussel shell C and thus that C isotopes in shell were thought to directly reflect their environmental conditions, just as stable isotope analysis of oxygen in shells has been used to reflect ambient temperature throughout a mussel\u0026rsquo;s lifetime \u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. More recently, it has become clear that metabolic C can make up a notable proportion of mussel shell C (%C\u003csub\u003eM\u003c/sub\u003e), ~\u0026thinsp;12% in marine mussels but up to 25\u0026ndash;40% in freshwater mussels \u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e,\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e,\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e,\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. Specifically, a study of the North American freshwater mussel, \u003cem\u003eElliptio complanata\u003c/em\u003e, found an average %C\u003csub\u003eM\u003c/sub\u003e of 20\u0026thinsp;\u0026plusmn;\u0026thinsp;10%, the large uncertainty represents the variation due to factors including ontogenic shifts, metabolic rates, and food source and selectivity \u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. Determining the relative proportion of C sources is not only important for paleo work, but also to understand the impact of modern mineralization on C cycling. While shell mineralization from DIC results in the release of a CO\u003csub\u003e2\u003c/sub\u003e molecule, when organic C is the original source a molecule of CO\u003csub\u003e2\u003c/sub\u003e from respiration is diverted from the release to the environment. Thus, the relative proportion of C sources for shell growth in dreissenid mussels need to be determined to understand their impact on inland C cycling (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe Laurentian Great Lakes (hereafter Great Lakes) contain\u0026thinsp;~\u0026thinsp;18% of the world\u0026rsquo;s liquid, surficial freshwater, with a watershed of approximately 1.0\u0026nbsp;million km\u003csup\u003e2 30\u003c/sup\u003e. Due to cultural eutrophication, most of the budget calculations in the Great Lakes focused on nutrients (nitrogen and phosphorus) rather than on carbon. Preliminary calculation indicates that the carbon budget is dominated by autochthonous production, followed by inorganic and organic C inputs from stream and groundwater and surface runoff \u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e. Previous work has also identified the Great Lakes as a net source of CO\u003csub\u003e2\u003c/sub\u003e to the atmosphere, although CO\u003csub\u003e2\u003c/sub\u003e saturation decreases from the upstream to downstream lakes, and varies seasonally \u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan additionalcitationids=\"CR32\" citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. Additionally, \u003csup\u003e[NO_PRINTED_FORM] 13\u003c/sup\u003e, found an increase in Great Lakes pCO\u003csub\u003e2\u003c/sub\u003e saturation following the establishment of large quagga mussel (\u003cem\u003eDreissena bugensis\u003c/em\u003e) populations and hypothesized that this was due to the alteration of lake water biogeochemistry via decreased primary productivity, increased water clarity, and respiration from the large mussel populations. This finding draws from a large body of work examining the physicochemical impacts of invasive dreissenid mussels \u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e,\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e, especially in the Great Lakes \u003csup\u003e\u003cspan additionalcitationids=\"CR37 CR38 CR39\" citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e. Although the impacts of living dreissenid mussels are well understood in the Great Lakes and other invaded territories, there has been very little work examining the impacts of the accumulating mass of mussel shells on the carbon budget in these lakes \u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e,\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e,\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThere has been some work exploring the scale and role of dreissenid mussel shells in inland lakes, however, as of yet, the Great Lakes, the largest freshwater ecosystem on Earth, have not been considered in this field. In Europe, a study investigated the impact of dreissenid mussel shell accumulation on a shallow lakes\u0026rsquo; life span, and found that the invasion of \u003cem\u003eD. bugensis\u003c/em\u003e could shorten the life span of the lake by one- to two-thirds by filling it with spent shells \u003csup\u003e\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e. Additionally, \u003csup\u003e[NO_PRINTED_FORM] 44\u003c/sup\u003e, investigated shell decay rates for native and invasive bivalves in North American waters and found that \u003cem\u003eDreissena\u003c/em\u003e spp. can produce standing stocks of spent shells on the scale of \u0026gt;\u0026thinsp;10 kg dry mass/m\u003csup\u003e2\u003c/sup\u003e under appropriate conditions (low flow, high calcium) \u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e. \u003csup\u003e[NO_PRINTED_FORM] 45\u003c/sup\u003e, collected benthic samples from ten Minnesota lakes and although they found high variability in empty shell mass across lakes, they determined that dead shell material accounted for an average of 322% of lake particulate organic carbon (POC). Finally, the most comprehensive study of dreissenid shell accumulation in a large lake, was conducted in Lake Simcoe, Ontario, Canada (SA\u0026thinsp;=\u0026thinsp;722 km\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e) \u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e. Based on samples collected from 0\u0026ndash;32 m depths, a total dead shell mass of 1 141 000 tonnes was estimated \u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e. Additionally, the authors found that approximately 1 000 000 tonnes of shells are produced annually in Lake Simcoe, containing\u0026thinsp;~\u0026thinsp;115 600 tonnes of C.\u003c/p\u003e \u003cp\u003eSince their invasion of the Great Lakes in the late 1980s, dreissenid mussel populations have had consistent growth, with an initial growth of the zebra mussel (\u003cem\u003eD. polymorpha\u003c/em\u003e) population followed by the delayed establishment of quagga mussels \u003csup\u003e\u003cspan additionalcitationids=\"CR48\" citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. The establishment of quagga mussels not only increased dreissenid density in previously established zebra mussel habitat, but also expanded to much greater depths (\u0026gt;\u0026thinsp;90 m) and softer substrates \u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e,\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e,\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e\u003c/sup\u003e. The growing dreissenid mussel population was also able to withstand predation pressure from waterfowl, lake whitefish, and their native predator, the round goby (\u003cem\u003eNeogobius melanostomus\u003c/em\u003e), which established Great Lake populations in the early 2000s \u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e,\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e,\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e\u003c/sup\u003e. Predation rates vary by habitat, but previous work has estimated round goby predation of dreissenid mussels in Lake Erie to be six kT/yr (270 T C\u003csub\u003eshell\u003c/sub\u003e/yr), and lake whitefish consumption of 109 and 820 kT/yr (500 and 3800 T C\u003csub\u003eshell\u003c/sub\u003e/yr) in Lakes Michigan and Huron, respectively \u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e,\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eWhile dreissenid predation has increased with long-term establishment, there is minimal evidence that these predators are reducing dreissenid abundance, although a shift in size structure has been observed \u003csup\u003e\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e,\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e,\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e\u003c/sup\u003e. Predation rate and the dreissenid population impact the turnover rate of dreissenid mussels (and production of empty mussel shells), through their impact on the production: biomass (P:B) ratio. The stable dreissenid population maintains biomass (B), while the increased predation rate results in increased dreissenid production (P) (Astorg et al., 2022; Kipp \u0026amp; Ricciardi, 2012; Naddafi \u0026amp; Rudstam, 2014). Therefore, the combination of stable dreissenid populations with increased predation suggests that dreissenid mussel turnover has increased, thus increasing the rate of dreissenid shell mineralization in the Great Lakes (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe evidence of large-scale inorganic C sequestration in much smaller North American lakes and a large stable dreissenid population in the Great Lakes makes a clear case for calculating C sequestered in empty mussel shells in the Great Lakes. Altogether it is evident that dreissenids are persistent ecosystem engineers in the Great Lakes \u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e,\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e,\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e\u003c/sup\u003e, however, the impact of nutrient sequestration in their shells that are continuously accumulating at the bottom of these lakes has not previously been considered as a component of their engineering role. Here we provide an estimation of the C mass in empty dreissenid shells in the four dreissenid-invaded Great Lakes (Michigan, Huron, Erie, and Ontario). This quantification provides a starting point for considering another pathway that dreissenid mussels have altered C dynamics in the Great Lakes. Finally, we suggest future work that will better resolve the impacts of this pathway on C budgets.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003eTo estimate the scale of mineralized C in empty dreissenid mussels in the Great Lakes we calculated the empty shell standing stock based on average pre- (1990\u0026ndash;2005) and post (2005\u0026ndash;2020)-goby invasion mussel densities from 1990\u0026ndash;2020 (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Estimates were split by pre- and post-goby periods to account for the increased dreissenid mussel mortality (Z) following the establishment of stable round goby populations \u003csup\u003e\u003cspan additionalcitationids=\"CR53\" citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e\u003c/sup\u003e. Round goby invasion also coincided with quagga mussel dominance which occurred in most lakes by the year 2000 and expanded the mussel population to deeper depths \u003csup\u003e\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e,\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e\u003c/sup\u003e. This combination of increased mussel abundance (biomass) and predation pressure (mortality rate) resulted in much higher shell turnover and thus higher mineralized C in empty shells (2.40 X 10\u003csup\u003e8\u003c/sup\u003e vs 1.19 X 10\u003csup\u003e9\u003c/sup\u003e T C). This shift was most prominent in Lakes Ontario and Michigan, where increased mussel densities in the deeper waters and increased predation resulted in much higher C mineralization in the post-round goby period \u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e,\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e,\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e\u003c/sup\u003e. An estimated total of 1.19 E10 tonnes of empty shell mass was calculated across the four lakes, containing 1.43 E9 tonnes of C.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eEstimated carbon stored in empty dreissenid shells (in tonnes of carbon) in the Great Lakes before and after the invasion of its native predator, the round goby. Based on average 15 years pre-goby invasion (1990\u0026ndash;2005) and 15 years post-goby invasion (2005\u0026ndash;2020) mussel densities.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eLake Ontario\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eLake Erie\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eLake Huron\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eLake Michigan\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eTotal\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePre-round goby C storage (T C)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1.15 X 10\u003csup\u003e7\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1.74 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.64 X 10\u003csup\u003e7\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3.80 X 10\u003csup\u003e7\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2.40 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePost-round goby C storage (T C)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1.11 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e9.02 X 10\u003csup\u003e7\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e9.83 X 10\u003csup\u003e7\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e8.92 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1.19 X 10\u003csup\u003e9\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTotal C storage (T C)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1.23 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2.64 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.15 X 10\u003csup\u003e8\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e9.30 X 10\u003csup\u003e9\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1.43 X 10\u003csup\u003e9\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eDreissenid mussel shells clearly represent a large deviation in carbon flow within the lakes and need to be considered in future carbon budget calculations \u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e. This work is a key step in filling the knowledge gap regarding the role of aquatic ecosystems as contributors or sinks of carbon to the atmosphere. \u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e\u003c/sup\u003e. Dreissenid mussels have existed in the Great Lakes for just over 30 years and our calculations indicate that their accumulated spent shells account for 1.43 Gt of inorganic C. Throughout time, Canadian lakes are estimate to have stored a total of 20 Gt of C, our calculations indicate that inorganic C from dreissenid shells would make up ~\u0026thinsp;7% of this stored C \u003csup\u003e\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e,\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e\u003c/sup\u003e. On an areal basis, there was a large range of inorganic C deposition rates across the lakes (0.08\u0026ndash;37.5 T C/ha/yr), which were primarily dependent on mussel density. This range of C storage is on the scale of estimates of organic C storage from nature-based climate solutions (NBCSs) and depending on the C source (metabolic C or DIC), may supplement or offset these NBCS\u0026rsquo;s C storage (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Based on the current literature, it is likely that DIC is the primary source, indicating that dreissenid shell mineralization will release CO\u003csub\u003e2\u003c/sub\u003e to the environment, offsetting the CO\u003csub\u003e2\u003c/sub\u003e sink from NBCS\u0026rsquo;s \u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e,\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. Dreissenid mussel inorganic C deposition is also on a similar scale to organic C sequestration rates in established coastal blue ecosystems (seagrass, mangroves, and salt marshes) (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e; Alongi, 2023; Claes et al., 2022). When compared to Canadian freshwater NBCSs, inorganic C is deposited from dreissenid mussel shells in the Great Lakes at a much greater rate than organic C in wetlands and peatlands. Specifically, the inorganic C from shells is deposited at a rate greater than organic C stored in peatlands (23 g/m2/yr), freshwater marshes in the Lake Simcoe watershed (average 112 g C/m\u003csup\u003e2\u003c/sup\u003e/yr), and other intact (62 g C/m\u003csup\u003e2\u003c/sup\u003e/yr) and restored (51\u0026ndash;89 g C/m\u003csup\u003e2\u003c/sup\u003e/yr) wetlands in southern Ontario \u003csup\u003e\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e,\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe inorganic C deposition of dreissenid shells is also on the scale of inorganic carbon storage calculated in other systems. \u003csup\u003e[NO_PRINTED_FORM] 21\u003c/sup\u003e, found inorganic C burial in Arctic fjords at a rate of 10.7\u0026ndash;45.7 g C/m\u003csup\u003e2\u003c/sup\u003e/yr, which is within the range of our rate of ~\u0026thinsp;0.2\u0026ndash;3750 g C/m\u003csup\u003e2\u003c/sup\u003e/yr, with a median of 201 g C/m\u003csup\u003e2\u003c/sup\u003e/yr. On a global scale, seagrass ecosystems have been found to accumulate inorganic C at rates from 149\u0026ndash;756 g C/m\u003csup\u003e2\u003c/sup\u003e/yr \u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan additionalcitationids=\"CR65\" citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e66\u003c/span\u003e\u003c/sup\u003e. However, it is important to note that these accumulation rates are based on all inorganic C deposited in bottom sediments (allochthonous and autochthonous) and not specific to calcifying organisms \u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e. When looking specifically at inorganic C accumulation from local calcification in seagrass, the range is much lower, 2.2\u0026ndash;36 g C/m\u003csup\u003e2\u003c/sup\u003e/yr \u003csup\u003e\u003cspan additionalcitationids=\"CR65 CR66\" citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e67\u003c/span\u003e\u003c/sup\u003e. Thus, the inorganic C accumulation from dreissenid mussels in the Great Lakes is at a scale that should be considered in future C sequestration rate and C budget calculations.\u003c/p\u003e \u003cp\u003eThe estimates provided here are based on historical turnover and shell degradation rates, however, it is essential to consider the impact of predicted future changes to the Great Lakes on dreissenid shell C burial. An updated investigation into shell decay rates in the Great Lakes is recommended to provide a more accurate value for future calculations. \u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e found that shell decay rates were correlated with pH and flow, however \u003cem\u003eD. bugensis\u003c/em\u003e occupy much deeper habitats than \u003cem\u003eD. polymorpha\u003c/em\u003e and other depth related variables (DO, light, nutrients etc.) should be considered for future calculations, especially as climate change is expected to change these conditions \u003csup\u003e\u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e68\u003c/span\u003e\u003c/sup\u003e. Additionally, winter ice scour can result in high over-winter mussel mortality at shallow depths \u003csup\u003e\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e\u003c/sup\u003e. Climate change is predicted to alter the ice-cover period in the Great Lakes \u003csup\u003e\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e,\u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e\u003c/sup\u003e, which could reduce the occurrence of ice scour mortality. Future changes to other factors that affect mussel population dynamics (nutrient levels, substrate) should also be considered since they will influence the rate of shell deposition as well. Following their initial establishment, dreissenid populations were relatively consistent in their invaded territories, even following the establishment of their native predator, the round goby \u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e,\u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e72\u003c/span\u003e\u003c/sup\u003e. However, it is possible that new predators may be introduced, or food source availability may change, altering mussel mortality rate \u003csup\u003e\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e,\u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003cp\u003eDreissenid mussels are considered one of the most impactful invasive species in North American waters. Serving as ecosystem engineers, they modify benthic habitats and influence phosphorus cycling and nuisance algae growth \u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e,\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e,\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e\u003c/sup\u003e. We present here another ecosystem-scale impact that expands their role in nutrient cycling to include inorganic C sequestration. On the scale of billions of tonnes of inorganic carbon deposited in slowly degrading shells at the bottom of these lakes. Although our calculations provide an estimate of the scale of their C deposition capacity, physically measuring the mass of shells in the Great Lakes and connecting channels, as \u003csup\u003e[NO_PRINTED_FORM] 46\u003c/sup\u003e did in Lake Simcoe, will provide the best estimation C stored in empty dreissenid shells. An essential part of this work will be determining the relative proportions of each shell C source (metabolic and DIC), this will provide the insight needed to determine the role of dreissenids on C cycling and CO\u003csub\u003e2\u003c/sub\u003e flux in the Great Lakes and other invaded territories.\u003c/p\u003e \u003cp\u003eThe storage of other elements, such as nitrogen and phosphorus, in dreissenid mussel shells should also be calculated. Although the proportions of shell nitrogen (~\u0026thinsp;0.25%) and phosphorus (~\u0026thinsp;0.007%) are much smaller than carbon (~\u0026thinsp;12%) \u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e,\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e,\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e, based on the shear scale of deposited shells they may be an important sink for these nutrients as well, especially when considering that phosphorus concentrations have declined in most of these lakes \u003csup\u003e\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e,\u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e73\u003c/span\u003e\u003c/sup\u003e. Finally, the creation of shell fragments from predation needs to be considered. Previous calculations of shell decay rates are based on whole shells suspended in the water column, however, most dreissenid mussel predators break the shells into fragments to consume the mussel tissue and discard shell fragments that settle into the sediments \u003csup\u003e\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e\u003c/sup\u003e. As described in the Methods, it is not clear how these small shell fragments will impact dissolution rates, they could reduce dissolution by increasing shell compaction (reducing exposed surface area) or as smaller particles they could dissolve more quickly \u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e. If dissolution does increase it would result in CO\u003csub\u003e2\u003c/sub\u003e uptake from the water column through the reverse of Eq.\u0026nbsp;1. Evidently, determining the rate of shell dissolution with shell fragments will provide a better understanding of dreissenid mussel impacts on C cycling in the Great Lakes.\u003c/p\u003e \u003cp\u003eThis work adds to the previously explored roles of dreissenid mussels as ecosystem engineers in the Great Lakes and other invaded waterways. In addition to their phosphorus cycling role and increased respiration rates, dreissenid mussels are altering the flow of carbon through mineralization on a massive scale \u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e,\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e. As mentioned above, determining the relative role of metabolic C vs. DIC in shell mineralization is key to deciphering their exact impact on C cycling and CO\u003csub\u003e2\u003c/sub\u003e flux in invaded waters. Previous work indicates that ambient bicarbonate is likely the primary source, and thus until the relative proportions are determined we should consider dreissenid shell mineralization as a likely source of CO\u003csub\u003e2\u003c/sub\u003e to invaded waters \u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. While dreissenid mussels have been found to alter CO\u003csub\u003e2\u003c/sub\u003e concentrations through respiration and impacts on primary productivity, their shell mineralization has not been considered as a component of this impact \u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e74\u003c/span\u003e\u003c/sup\u003e. Regardless of the direction (sink or source) it is evident that shell mineralization has major impacts on C cycling in the Great Lakes.\u003c/p\u003e \u003cp\u003eOn a larger scale these findings have implications for calculation of C budgets for the Great Lakes and implications for the global inland waters C balance. A preliminary C budget for the Great Lakes indicates that net primary production (NPP) and respiration are the greatest C flows for all lakes except Ontario, which receives more C from upstream inputs \u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e. Our calculations indicate that dreissenid mineralization alters C flow on the same scale as the processes considered in \u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e calculations and should be included in future C budgets for the invaded lakes (Michigan, Huron, Erie, and Ontario). Additionally, these C fluxes should be included in global inland water C budgets, although we focused on the Great Lakes, there are thousands of invaded lakes and rivers across North America and Europe that will have altered C flows due to dreissenid mussel mineralization \u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e,\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e. Finally, inclusion in these C budgets will have policy implications in the corresponding jurisdictions. Many locations have drafted adaptive management plans to minimize greenhouse gas emissions and an improved understanding of local C dynamics will make these plans more robust as dreissenid mussel populations continue to grow.\u003c/p\u003e "},{"header":"Methods","content":"\u003cp\u003eCalculation of dreissenid empty shell standing stock was based on \u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e, where at steady state the standing stock of spent shells (SS) is equal to the production of spent shells (P) divided by the instantaneous decay rate (I). Production of spent shells was determined from living shell mass (B) and instantaneous mortality rate (Z, or P:B depending on availability), with Eq.\u0026nbsp;1.\u003c/p\u003e \u003cp\u003eEquation 1.\u003cdiv id=\"Equa\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equa\" name=\"EquationSource\"\u003e\n$$SS=\\frac{B*Z}{I}$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eThe instantaneous shell decay rate of 0.058 calculated by \u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e was used as it is comparable to other estimates \u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e, and was calculated in Lake Simcoe, a lake in the same region as the Laurentian Great Lakes, with similar geochemical characteristics. This decay rate is a very conservative estimate as it is based on shells suspended in the water column, where turbulence is greater than on the lake bottom, which would result in an overestimation of the dissolution rate. Additionally, the settlement of shells in layers on the lake bottom would reduce exposed surface area and likely result in slower decay compared to shells in the water column. Finally, the increased shell fragmentation from round goby predation could impact shell decay rate, however without published findings on this impact it is not possible to adjust decay rate for this factor. Shell fragmentation would alter the matrix of deposited shells potentially reducing decay rate by reducing exposed surface area, alternatively \u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e, showed that smaller shells result in higher decay rates. Considering this contradiction, the published value was used to avoid making an incorrect assumption due to the lack of literature on this factor.\u003c/p\u003e \u003cp\u003eTwo mortality rates were used, one for the pre-round goby period (~\u0026thinsp;1990\u0026ndash;2005) and one for the post-round goby period (~\u0026thinsp;2005\u0026ndash;2020). The pre-round goby, non-predatory mortality rate of 1.26/year from \u003csup\u003e\u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e\u003c/sup\u003e. The post-round goby production: biomass ratio (P: B)/mortality rate of 2.1 from \u003csup\u003e\u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e76\u003c/span\u003e\u003c/sup\u003e was applied for the post-round goby invasion period\u0026thinsp;~\u0026thinsp;2005\u0026ndash;2020.\u003c/p\u003e \u003cp\u003eAn average living shell biomass was determined for each period based on published density values for each depth zone (\u0026lt;\u0026thinsp;30 m, 31\u0026ndash;50 m, 51\u0026ndash;90 m) for Lake Ontario, Lake Huron, and Lake Michigan. As Lake Erie is very shallow and has more variation across its basins (East, Central and West) than depth, shell biomass was determined from mean basin-wide biomass. Lake Superior was not included in calculations as it\u0026rsquo;s dreissenid populations are limited to select areas of the lake, and extrapolating to a lake-wide estimation based on depth zones would not be appropriate \u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e,\u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e77\u003c/span\u003e\u003c/sup\u003e. Living dreissenid density data was extracted from lake-wide surveys with at least 30 sites, across a range of bottom substrates to represent the range of conditions in each lake (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Additionally, where a series of similar surveys were conducted the most recent published data was extracted. All surveys collected between 2\u0026ndash;6 replicates for each sample site.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eSummary of dreissenid density data sources for each lake and period (pre- and post-round goby invasion).\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLake\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePeriod\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSource\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eYear(s) sampled\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSample sites (#)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eDepth range (m)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eLake Ontario\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePre-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e83\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1997\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e75\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e14\u0026ndash;209\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePost-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2008, 2013, 2018\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e55\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e11\u0026ndash;209\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eLake Erie\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePre-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e78\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1998\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2.0\u0026ndash;58.5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePost-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR84\" class=\"CitationRef\"\u003e84\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2009, 2012\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e135\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2\u0026ndash;66\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eLake Huron\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePre-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2000, 2003\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e65\u0026ndash;85\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e19\u0026ndash;173\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePost-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2007, 2012\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e80\u0026ndash;83\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e12\u0026ndash;171\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eLake Michigan\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePre-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1994,1995,2000\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e90\u0026ndash;157\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e6\u0026ndash;208\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePost-goby\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2005, 2010\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e140\u0026ndash;160\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e13\u0026ndash;207\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eLiving mussel shell density was calculated from the wet mass density (g/m\u003csup\u003e2\u003c/sup\u003e) assuming shells are approximately 38.13% of whole mussel wet weight (B\u0026aacute;ldi et al., 2019; Table S4), if whole mussel wet weight was not available it was calculated based on density of individuals (ind./m\u003csup\u003e2\u003c/sup\u003e) using an average wet weight of 0.156 g/ind. \u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e,\u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e78\u003c/span\u003e\u003c/sup\u003e. Empty shell standing stock was multiplied by 0.12, the approximate proportion of carbon in the shell \u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e,\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e, to get the areal mass of carbon in empty shells deposited each year. Lake surface area for each depth zone described above was calculated from NOAA bathymetry maps \u003csup\u003e\u003cspan additionalcitationids=\"CR80 CR81\" citationid=\"CR79\" class=\"CitationRef\"\u003e79\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e82\u003c/span\u003e\u003c/sup\u003e. Standing stock of empty dreissenid shell C deposited per year was calculated by combining areal production with each zone\u0026rsquo;s surface area. Finally, the total mass of C stored in empty dreissenid shells was calculated by multiplying the pre- and post-round goby estimates by 15 years, the approximate time for each period.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e \u003ch2\u003eCompeting Interests\u003c/h2\u003e \u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eEDS: Data curation, formal analysis, writing \u0026ndash; original draft preparation. LJM: Conceptualization, supervision, writing \u0026ndash; review \u0026amp; editing. JM: Writing \u0026ndash; review \u0026amp; editing. AEK: Writing \u0026ndash; review \u0026amp; editing. JR: Supervision, Writing \u0026ndash; review \u0026amp; editing.\u003c/p\u003e\u003ch2\u003eAcknowledgement\u003c/h2\u003e\u003cp\u003eWe would like to thank all members of our team for their efforts. The conceptual models were created by Lexy Harquail. We would also like to thank our funding agencies: Mitacs and RBC Foundation.\u003c/p\u003e\u003ch2\u003eData Availability\u003c/h2\u003e\u003cp\u003eAll data generated or analysed during this study are included in this published article [and its supplementary information files].\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eTranvik, L. J. \u003cem\u003eet al.\u003c/em\u003e Lakes and reservoirs as regulators of carbon cycling and climate. Limnol Oceanogr 54, 2298\u0026ndash;2314 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRegnier, P. \u003cem\u003eet al.\u003c/em\u003e Anthropogenic perturbation of the carbon fluxes from land to ocean. Nat Geosci 6, 597\u0026ndash;607 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCole, J. J. \u003cem\u003eet al.\u003c/em\u003e Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems 10, 171\u0026ndash;184 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRegnier, P., Resplandy, L., Najjar, R. G. \u0026amp; Ciais, P. The land-to-ocean loops of the global carbon cycle. Nature 603, 401\u0026ndash;410 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRaymond, P. A. \u003cem\u003eet al.\u003c/em\u003e Global carbon dioxide emissions from inland waters. Nature 503, 355\u0026ndash;359 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWebb, J. R., Santos, I. R., Maher, D. T. \u0026amp; Finlay, K. The Importance of Aquatic Carbon Fluxes in Net Ecosystem Carbon Budgets: A Catchment-Scale Review. Ecosystems 22, 508\u0026ndash;527 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePilla, R. M. \u003cem\u003eet al.\u003c/em\u003e Anthropogenically driven climate and landscape change effects on inland water carbon dynamics: What have we learned and where are we going? \u003cem\u003eGlobal Change Biology\u003c/em\u003e vol. 28 5601\u0026ndash;5629 Preprint at \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/gcb.16324\u003c/span\u003e\u003cspan address=\"10.1111/gcb.16324\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMulholland, P. J. \u0026amp; Elwood, J. W. The role of lake and reservoir sediments as sinks in the perturbed global carbon cycle. Tellus 34, 490\u0026ndash;499 (1982).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMendon\u0026ccedil;a, R. \u003cem\u003eet al.\u003c/em\u003e Organic carbon burial in global lakes and reservoirs. Nat Commun 8, (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHeathcote, A. J. \u0026amp; Downing, J. A. Impacts of Eutrophication on Carbon Burial in Freshwater Lakes in an Intensively Agricultural Landscape. Ecosystems 15, 60\u0026ndash;70 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAnderson, N. J., Bennion, H. \u0026amp; Lotter, A. F. Lake eutrophication and its implications for organic carbon sequestration in Europe. Glob Chang Biol 20, 2741\u0026ndash;2751 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcKinley, G., Urban, N., Bennington, V., Pilcher, D. \u0026amp; Mcdonald, C. \u003cem\u003ePreliminary Carbon Budgets for the Laurentian Great Lakes\u003c/em\u003e. vol. 4 (2011).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLin, P. \u0026amp; Guo, L. Do invasive quagga mussels alter CO2 dynamics in the Laurentian Great Lakes? Sci Rep 6, (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLu, Y., Meyers, P. A., Eadie, B. J. \u0026amp; Robbins, J. A. Carbon cycling in Lake Erie during cultural eutrophication over the last century inferred from the stable carbon isotope composition of sediments. J Paleolimnol 43, 261\u0026ndash;272 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHodell, D. A., Schelske, C. L., Fahnenstiel, G. L. \u0026amp; Robbins, L. L. Biologically induced calcite and its isotopic composition in Lake Ontario. Limnol Oceanogr 43, 187\u0026ndash;199 (1998).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBarbiero, R. P., Tuchman, M. L. \u0026amp; Millard, E. S. Post-dreissenid increases in transparency during summer stratification in the offshore waters of Lake Ontario: Is a reduction in whiting events the cause? J Great Lakes Res 32, 131\u0026ndash;141 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMartin, J. B. Carbonate minerals in the global carbon cycle. Chem Geol 449, 58\u0026ndash;72 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eButman, D. \u003cem\u003eet al.\u003c/em\u003e Aquatic carbon cycling in the conterminous United States and implications for terrestrial carbon accounting. PNAS 113, 58\u0026ndash;63 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu, Z., Dreybrodt, W. \u0026amp; Wang, H. A new direction in effective accounting for the atmospheric CO2 budget: Considering the combined action of carbonate dissolution, the global water cycle and photosynthetic uptake of DIC by aquatic organisms. \u003cem\u003eEarth-Science Reviews\u003c/em\u003e vol. 99 162\u0026ndash;172 Preprint at \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.earscirev.2010.03.001\u003c/span\u003e\u003cspan address=\"10.1016/j.earscirev.2010.03.001\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMinor, E. C. \u0026amp; Brinkley, G. Alkalinity, pH, and pCO2 in the Laurentian Great Lakes: An initial view of seasonal and inter-annual trends. J Great Lakes Res 48, 502\u0026ndash;511 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKoziorowska, K., Kuliński, K. \u0026amp; Pempkowiak, J. Comparison of the burial rate estimation methods of organic and inorganic carbon and quantification of carbon burial in two high Arctic fjords. Oceanologia 60, 405\u0026ndash;418 (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaderne, V. \u003cem\u003eet al.\u003c/em\u003e Role of carbonate burial in Blue Carbon budgets. Nat Commun 10, (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChapin, F. S. \u003cem\u003eet al.\u003c/em\u003e Reconciling carbon-cycle concepts, terminology, and methods. Ecosystems 9, 1041\u0026ndash;1050 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eThomsen, J., Haynert, K., Wegner, K. M. \u0026amp; Melzner, F. Impact of seawater carbonate chemistry on the calcification of marine bivalves. Biogeosciences 12, 4209\u0026ndash;4220 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGillikin, D. P., Lorrain, A., Meng, L. \u0026amp; Dehairs, F. A large metabolic carbon contribution to the δ13C record in marine aragonitic bivalve shells. Geochim Cosmochim Acta 71, 2936\u0026ndash;2946 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eP\u0026eacute;rez, C. A., Lagos, N. A., Sald\u0026iacute;as, G. S., Waldbusser, G. \u0026amp; Vargas, C. A. Riverine discharges impact physiological traits and carbon sources for shell carbonate in the marine intertidal mussel Perumytilus purpuratus. Limnol Oceanogr 61, 969\u0026ndash;983 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcConnaughey, T. A., Bljrdeti, J., Whelan, J. F. \u0026amp; Paull, C. K. Carbon isotopes in biological carbonates: Respiration and photosynthesis. Geochim Cosmochim Acta 61, (1997).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGraniero, L. E., Gillikin, D. P. \u0026amp; Surge, D. Evaluating Freshwater Mussel Shell δ13C Values as a Proxy for Dissolved Inorganic Carbon δ13C Values in a Temperate River. J Geophys Res Biogeosci 126, (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePoulain, C. \u003cem\u003eet al.\u003c/em\u003e Experimental shift of diet and DIC stable carbon isotopes: Influence on shell δ13C values in the Manila clam Ruditapes philippinarum. Chem Geol 272, 75\u0026ndash;82 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWaples, J. \u003cem\u003eet al.\u003c/em\u003e The Laurentian Great Lakes. (2008).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePhillips, J. C. \u003cem\u003eet al.\u003c/em\u003e The potential for CO2-induced acidification in freshwater: A great lakes case study. Oceanography 28, 136\u0026ndash;145 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKarim, A., Dubois, K. \u0026amp; Veizer, J. Carbon and oxygen dynamics in the Laurentian Great Lakes: Implications for the CO2 flux from terrestrial aquatic systems to the atmosphere. Chem Geol 281, 133\u0026ndash;141 (2011).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eUrban, N. R. \u0026amp; Desai, A. Are the Great Lakes a significant net source or sink of CO2? Verh. Internat. Verein. Limnol 30, 1283\u0026ndash;1288 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi, J. \u003cem\u003eet al.\u003c/em\u003e Benthic invaders control the phosphorus cycle in the world\u0026rsquo;s largest freshwater ecosystem. PNAS 118, 1\u0026ndash;9 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHiggins, S. N. \u0026amp; Vander Zanden, M. J. What a difference a species makes: Meta-analysis of dreissenid mussel impacts on freshwater ecosystems. Ecol Monogr 80, 179\u0026ndash;196 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHecky, R. E., Smith, R. E. H., Barton, D. R., Guildford, S. J. \u0026amp; al, et. The nearshore phosphorus shunt: a consequence of ecosystem engineering by dreissenids in the Laurentian Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 61, 1285\u0026ndash;1293 (2004).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eArnott, D. L. \u0026amp; Vanni, M. J. Nitrogen and phosphorus recycling by the zebra mussel (Dreissena polymorpha) in the western basin of Lake Erie. Canadian Journal of Fisheries and Aquatic Sciences 53, 646\u0026ndash;659 (1996).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDeRoy, E. M. \u0026amp; MacIsaac, H. J. Impacts of invasive species in the Laurentian Great Lakes. in \u003cem\u003eHandbook of Environmental Chemistry\u003c/em\u003e vol. 101 135\u0026ndash;156 (Springer Science and Business Media Deutschland GmbH, 2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChase, M. E. \u0026amp; Bailey, R. C. The ecology of the zebra mussel (Dreissena polymorpha) in the lower Great Lakes of North America: II. Total production, energy allocation, and reproductive effort. J Great Lakes Res 25, 122\u0026ndash;134 (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMadenjian, C. P. \u003cem\u003eet al.\u003c/em\u003e Dreissenid mussels are not a \u0026lsquo;dead end\u0026rsquo; in Great Lakes food webs. J Great Lakes Res 36, 73\u0026ndash;77 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBurlakova, L. E., Karatayev, A. Y., Boltovskoy, D. \u0026amp; Correa, N. M. Ecosystem services provided by the exotic bivalves Dreissena polymorpha, D. rostriformis bugensis, and Limnoperna fortunei. Hydrobiologia 850, 2811\u0026ndash;2854 (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKaratayev, A. Y. \u003cem\u003eet al.\u003c/em\u003e Dreissena in Lake Ontario 30 years post-invasion. J Great Lakes Res 48, 264\u0026ndash;273 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eB\u0026aacute;ldi, K. \u003cem\u003eet al.\u003c/em\u003e Sediment contributing invasive dreissenid species in a calcareous shallow lake \u0026ndash; Possible implications for shortening life span of lakes by filling. Elementa 7, (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStrayer, D. L. \u0026amp; Malcom, H. M. Shell decay rates of native and alien freshwater bivalves and implications for habitat engineering. Freshw Biol 52, 1611\u0026ndash;1617 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWilliamson, F. \u0026amp; Ozersky, T. Lake Characteristics, Population Properties and Invasion History Determine Impact of Invasive Bivalves on Lake Nutrient Dynamics. Ecosystems 22, 1721\u0026ndash;1735 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOzersky, T., Evans, D. O. \u0026amp; Ginn, B. K. Invasive mussels modify the cycling, storage and distribution of nutrients and carbon in a large lake. Freshw Biol 60, 827\u0026ndash;843 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKaratayev, A. Y. \u0026amp; Burlakova, L. E. Dreissena in the Great Lakes: What have we learned in 30 years of invasion. Hydrobiologia (2022) doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s10750-022-04990-x\u003c/span\u003e\u003cspan address=\"10.1007/s10750-022-04990-x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNalepa, T. F., Riseng, C. M., Elgin, A. K. \u0026amp; Lang, G. A. \u003cem\u003eAbundance and Distribution of Benthic Macroinvertebrates in the Lake Huron System: Saginaw Bay, 2006\u0026ndash;2009, and Lake Huron, Including Georgian Bay and North Channel, 2007 and 2012\u003c/em\u003e. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e\u003c/span\u003e\u003cspan address=\"http://www.glerl.noaa.gov\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNalepa, T. F., Fanslow, D. L., Lang, G. A., Mabrey, K. \u0026amp; Rowe, M. \u003cem\u003eLake-Wide Benthic Surveys in Lake Michigan in 1994-95, 2000, 2005, and 2010: Abundances of the Amphipod Diporeia Spp. and Abundances and Biomass of the Mussels Dreissena Polymorpha and Dreissena Rostriformis Bugensis\u003c/em\u003e. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e\u003c/span\u003e\u003cspan address=\"http://www.glerl.noaa.gov\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e. (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVanderploeg, H. A. \u003cem\u003eet al.\u003c/em\u003e Dispersal and emerging ecological impacts of Ponto-Caspian species in the Laurentian Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 59, 1209\u0026ndash;1228 (2002).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEifert, R. A. \u003cem\u003eet al.\u003c/em\u003e Could quagga mussels impact offshore benthic community and surface sediment-bound nutrients in the Laurentian Great Lakes? \u003cem\u003eHydrobiologia\u003c/em\u003e (2023) doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s10750-023-05191-w\u003c/span\u003e\u003cspan address=\"10.1007/s10750-023-05191-w\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKipp, R. \u0026amp; Ricciardi, A. Impacts of the Eurasian round goby (Neogobius melanostomus) on benthic communities in the upper St. Lawrence River. Canadian Journal of Fisheries and Aquatic Sciences 69, 469\u0026ndash;486 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBunnell, D. B., Johnson, T. B. \u0026amp; Knight, C. T. The impact of introduced round gobies (Neogobius melanostomus) on phosphorus cycling in central Lake Erie. Canadian Journal of Fisheries and Aquatic Sciences 62, 15\u0026ndash;29 (2005).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNaddafi, R. \u0026amp; Rudstam, L. G. Predation on invasive zebra mussel, Dreissena polymorpha, by pumpkinseed sunfish, rusty crayfish, and round goby. Hydrobiologia 721, 107\u0026ndash;115 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAstorg, L., Charette, C., Windle, M. J. S. \u0026amp; Derry, A. M. Two decades since first invasion: Revisiting round goby impacts on nearshore aquatic communities in the Upper St. Lawrence River. J Great Lakes Res (2022) doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.jglr.2022.01.017\u003c/span\u003e\u003cspan address=\"10.1016/j.jglr.2022.01.017\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKaratayev, A. Y. \u003cem\u003eet al.\u003c/em\u003e Lake morphometry determines Dreissena invasion dynamics. Biol Invasions 23, 2489\u0026ndash;2514 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNalepa, T. F., Fanslow, D. L. \u0026amp; Lang, G. A. Transformation of the offshore benthic community in Lake Michigan: Recent shift from the native amphipod Diporeia spp. to the invasive mussel Dreissena rostriformis bugensis. Freshw Biol 54, 466\u0026ndash;479 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCouncil of Canadian Academies (CCA). \u003cem\u003eInland Freshwater Ecosystems\u003c/em\u003e. Nature-Based Climate Solutions (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFerland, M. E., Del Giorgio, P. A., Teodoru, C. R. \u0026amp; Prairie, Y. T. Long-term C accumulation and total C stocks in boreal lakes in northern Qu\u0026eacute;bec. Global Biogeochem Cycles 26, (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eClaes, J., Hopman, D., Jaeger, G. \u0026amp; Rogers, M. \u003cem\u003eBlue Carbon: The Potential of Coastal and Oceanic Climate Action Nature-Based Climate Solutions in the World\u0026rsquo;s Oceans Can Play an Important Role in Conservation and Carbon Abatement Efforts Worldwide\u003c/em\u003e. (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAlongi, D. M. Current status and emerging perspectives of coastal blue carbon ecosystems. Carbon Footprints 2, (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePendea, I. F., Kanavillil, N., Kurissery, S. \u0026amp; Chmura, G. L. Carbon Stocks and Recent Rates of Carbon Sequestration in Nutrient-Rich Freshwater Wetlands from Lake Simcoe Watershed (Southern Canada). J Geophys Res Biogeosci 128, (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCreed, I. F. \u003cem\u003eet al.\u003c/em\u003e Can Restoration of Freshwater Mineral Soil Wetlands Deliver Nature-Based Climate Solutions to Agricultural Landscapes? Front Ecol Evol 10, (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaderne, V. \u003cem\u003eet al.\u003c/em\u003e Accumulation of Carbonates Contributes to Coastal Vegetated Ecosystems Keeping Pace with Sea Level Rise in an Arid Region (Arabian Peninsula). J Geophys Res Biogeosci 123, 1498\u0026ndash;1510 (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMazarrasa, I. \u003cem\u003eet al.\u003c/em\u003e Seagrass meadows as a globally significant carbonate reservoir. Biogeosciences 12, 4993\u0026ndash;5003 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYates, K. K. \u0026amp; Halley, R. B. Diurnal variation in rates of calcification and carbonate sediment dissolution in Florida Bay. Estuaries and Coasts 29, 24\u0026ndash;39 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBarron, C., Duarte, C. M., Frankignoulle, M. \u0026amp; Vieira Borges, A. Organic Carbon Metabolism and Carbonate Dynamics in a Mediterranean Seagrass (Posidonia oceanica) Meadow. Estuaries and Coasts 29, 417\u0026ndash;426 (2001).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWoolway, R. I., Sharma, S. \u0026amp; Smol, J. P. Lakes in hot water: The impacts of a changing climate on aquatic ecosystems. Bioscience 72, 1050\u0026ndash;1061 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChase, M. E. \u0026amp; Bailey, R. C. The ecology of the zebra mussel (Dreissena polymorpha) in the lower Great Lakes of North America: I. Population dynamics and growth. J Great Lakes Res 25, 107\u0026ndash;121 (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMahdiyan, O., Filazzola, A., Molot, L. A., Gray, D. \u0026amp; Sharma, S. Drivers of water quality changes within the Laurentian Great Lakes region over the past 40 years. Limnol Oceanogr 66, 237\u0026ndash;254 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSharma, S. \u003cem\u003eet al.\u003c/em\u003e Widespread loss of lake ice around the Northern Hemisphere in a warming world. \u003cem\u003eNature Climate Change\u003c/em\u003e vol. 9 227\u0026ndash;231 Preprint at \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41558-018-0393-5\u003c/span\u003e\u003cspan address=\"10.1038/s41558-018-0393-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStrayer, D. L. \u003cem\u003eet al.\u003c/em\u003e Long-term population dynamics of dreissenid mussels (Dreissena polymorpha and D. rostriformis): a cross-system analysis. Ecosphere 10, (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDove, A. \u0026amp; Chapra, S. C. Long-term trends of nutrients and trophic response variables for the Great Lakes. Limnol Oceanogr 60, 696\u0026ndash;721 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBlagrave, K. \u003cem\u003eet al.\u003c/em\u003e Spatial heterogeneity in water quality across the northern nearshore regions of the Laurentian Great Lakes. J Great Lakes Res 49, (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBoles, L. C. \u0026amp; Lipcius, R. N. Potential for Population Regulation of the Zebra Mussel by Finfish and the Blue Crab in North American Estuaries (1997). J Shellfish Res 16, 179\u0026ndash;186 (1997).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKao, Y. C., Adlerstein, S. \u0026amp; Rutherford, E. The relative impacts of nutrient loads and invasive species on a Great Lakes food web: An Ecopath with Ecosim analysis. J Great Lakes Res 40, 35\u0026ndash;52 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGrigorovich, I. A. \u003cem\u003eet al. Lake Superior: An Invasion Coldspot? Hydrobiologia\u003c/em\u003e vol. 499 (2003).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJarvis, P., Dow, J., Dermott, R. \u0026amp; Bonnell, R. \u003cem\u003eZebra (Dreissena Polymorpha) and Quagga Mussel (Dreissena Bugensis) Distribution and Density in Lake Erie, 1992\u0026ndash;1998\u003c/em\u003e. (2000).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNOAA National Geophysical Data Center. Bathymetry of Lake Erie and Lake St. Clair. \u003cem\u003eNOAA National Centers for Environmental Information\u003c/em\u003e (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNOAA National Geophysical Data Center. Bathymetry of Lake Huron. \u003cem\u003eNOAA National Centers for Environmental Information\u003c/em\u003e (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNOAA National Geophysical Data Center. Bathymetry of Lake Michigan. NOAA National Centers for Environmental Information (1996).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNOAA National Geophysical Data Center. Bathymetry of Lake Ontario. \u003cem\u003eNOAA National Centers for Environmental Information\u003c/em\u003e (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLozano, S. J., Scharold, J. V. \u0026amp; Nalepa, T. F. Recent declines in benthic macroinvertebrate densities in Lake Ontario. Canadian Journal of Fisheries and Aquatic Science 58, 518\u0026ndash;529 (2001).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKaratayev, A. Y. \u003cem\u003eet al.\u003c/em\u003e Twenty five years of changes in Dreissena spp. populations in Lake Erie. J Great Lakes Res 40, 550\u0026ndash;559 (2014).\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"carbon cycling, dreissenid mussels, Laurentian Great Lakes, inland carbon budgets, mussel shells, inorganic carbon deposition","lastPublishedDoi":"10.21203/rs.3.rs-4436844/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4436844/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eSince the active role of inland waters in cycling carbon (C) has been revealed, there has been a renewed interest in calculating C budgets for inland waters to understand their role with respect to global climate change. There is a lack of knowledge regarding C cycling in the Laurentian Great Lakes, the worlds largest freshwater reservoir, with current estimates neglecting the role of invasive species. For one of the most pervasive invaders, dreissenid (zebra and quagga) mussels, research has focused on filter feeding impacts on phosphorus dynamics, but there is a lack of knowledge regarding their role in C cycling, specifically, the impact of the C stored in their slowly degrading shells. As such, we set out to estimate the mass of empty shells and C stored in those shells. We calculated an estimated 1.19 E10 tonnes of empty shell mass currently sitting at the bottom of these lakes, which store approximately 1.43 E9 tonnes of C. This scale of inorganic C storage is comparable to rates of organic C storage in nature-based climate solutions. This work demonstrates the importance of a previously unexplored pathway that dreissenid mussels are altering C cycling in the Laurentian Great Lakes and the thousands of other invaded lakes and rivers.\u003c/p\u003e","manuscriptTitle":"Revisiting carbon cycling in the Laurentian Great Lakes following dreissenid mussel invasion","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-05-31 20:29:25","doi":"10.21203/rs.3.rs-4436844/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"8d4775fe-adb1-4c88-9798-c38dc9059646","owner":[],"postedDate":"May 31st, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":32629189,"name":"Earth and environmental sciences/Environmental sciences"},{"id":32629190,"name":"Earth and environmental sciences/Limnology"},{"id":32629191,"name":"Earth and environmental sciences/Climate sciences/Biogeochemistry"}],"tags":[],"updatedAt":"2024-11-08T15:23:18+00:00","versionOfRecord":[],"versionCreatedAt":"2024-05-31 20:29:25","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4436844","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4436844","identity":"rs-4436844","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}