Biogeochemistry

The purpose of this study was to quantify the variation of baseline carbon and nitrogen stable isotope signatures in the Lake Simcoe watershed and relate that variation to various physicochemical parameters. Particulate organic matter samples from 2009 and 2011 were used as representatives of baseline isotopic values. Temporal data from two offshore lake stations revealed that δ15N of POM was lowest mid-summer and highest after the fall turnover. POM δ13C was variable throughout the summer before declining after fall turnover. Spatial data from the lake and the tributaries revealed that POM stable isotope signatures were highly variable. Various physicochemical parameters indicative of phytoplankton biomass were significantly positively correlated with POM δ15N and significantly negatively correlated with POM δ13C. The correlations were mostly significant in the tributaries, not the lake. Moreover, many of the correlations involving δ15N of POM were driven by extreme values in Cook's Bay and its tributaries. In general, it's likely that different processes or combination of processes were affecting the δ15N and δ13C POM in the Lake Simcoe watershed as physicochemical parameters alone could not explain the variability. Measuring the δ15N of ammonium and nitrate, as well as the δ13C of DIC would help discern the dominant nitrogen and inorganic carbon cycling processes occurring in the Lake Simcoe watershed.

Author Keywords: δ13C, δ15N, isotopic baseline, particulate organic matter, spatial variation, stable isotopes

Benthic macroinvertebrtes (BMI) are functionally important in aquatic ecosystems; as such, knowledge of their community structure and function is critical for understanding these systems. BMI were sampled from ten lakes in each of two regions of south-central Ontario to investigate which chemical and physical variables could be shaping their community structure and function. Ten Precambrian Shield lakes in the Muskoka-Haliburton region, and ten St Lawrence lowland lakes in the Kawartha lakes region were sampled. These lakes are geologically and chemically distinct, creating natural chemical and physical gradients within and between both regions. Community function was assessed using stable isotope analysis to elucidate carbon transfer dynamics (δ13C) and food web interactions (δ15N). It was predicted that the BMI from Shield lakes would have a δ13C signature indicative of allochthonous carbon subsidies, whereas the lowland lake BMI signatures would reflect autochthonous production. Additionally, it was predicted that the food web length (measured in δ15N units) would be different in Shield and lowland lakes. Both of these predictions were supported; however, the data indicate that δ13C signatures are more likely influenced by catchment geology (represented by bicarbonate concentration) than the extent of allochthony. The best predictor of food web length was found to be region. To assess BMI community structure, taxonomic richness, %EPT (% Ephemeroptera, Plecoptera, Trichoptera; a water quality index), and distribution of functional feeding groups were examined. Based on chemistry it was expected that the Shield lakes would be more speciose, and of greater water quality (relatively lower nutrient levels). These predictions were rejected; since there were no significant regional differences in taxonomic richness or biologically inferred water quality (%EPT). However, sediment size was found to best explain the variability in both metrics, with greater richness and %EPT found at sites with medium and small substrates than those with large substrates. Significant regional differences were found in the distribution of functional feeding groups. Most notably, there were significantly greater proportions of scrapers and shredders in the lowland and Shield lakes, respectively. Based on the feeding mechanisms of these invertebrates it can be inferred that allochthonous subsidies are likely of greater importance to Shield lake BMI communities than those of the lowland lakes; supporting the carbon transfer prediction. These findings provide insight about the structure and function of BMI communities from two dominant lake types in Ontario, and could be useful when determining how future chemical and physical changes will impact these communities.

Author Keywords: benthic macroinvertebrates, community function, community structure, Precambrian Shield, stable isotopes, St. Lawrence lowlands

Soil is an essential natural resource that supports provisioning services such as agriculture, silviculture, and mining. However, there is limited knowledge on forest soil properties across Canada. Digital soil mapping may be used to fill these data gaps, as it can predict soil properties in areas with limited observations. The focus of this study was to develop predictive maps of select soil physicochemical properties for the Kitimat Valley, British Columbia, and apply these maps to assess the potential impacts of sulphur dioxide emissions from an aluminum smelter, on soil properties in the Valley. Exchangeable [Ex.] magnesium, organic matter, pH, coarse fragment, Ex. potassium, bulk density, Ex. calcium, Ex. acidity, and Ex. sodium were all mapped with acceptable confidence. Time to depletion of base cation pools showed that ~240 km2 of the study area had a depletion time of 50 years or less. However, sources of base cations such as atmospheric deposition and mineral weathering were not considered.

Author Keywords: acidification, buffering capacity, Digital soil mapping, predictive mapping, regression kriging, soil properties

Aquatic contaminant mobility and biological availability is strongly governed by the complexation of organic and inorganic ligands. Dissolved organic matter (DOM) is a complex, heterogeneous mixture of organic acids, amino acids, lipids, carbohydrates and polyphenols that vary in composition and can complex to dissolved metals thereby altering their fate in aquatic systems. The research conducted in this doctoral dissertation addresses 1) how DOM composition differs between phytoplankton taxa and 2) how DOM composition affects metal speciation and its subsequent microbial bioavailability in laboratory and field conditions. To accomplish this, a series of analytical methods were developed and applied to quantify thiols, sulphur containing DOM moieties, and the molecular composition of DOM. The works presented in this thesis represents one of the first comprehensive and multipronged analyses of the impact of phytoplankton metabolite exudates on microbial metal bioavailability. This dissertation demonstrated the analytical versatility of high-resolution mass spectrometry as a tool for compound specific information, as well as having the capabilities to obtain speciation information of organometallic complexes. The work presented in this PhD strengthens the understanding compositional differences of both autochthonous and allochthonous DOM and their effects on metal biogeochemistry.

Author Keywords: Dissolved Organic Matter, Mercury, Metal Accumulation, Phytoplankton, Spring Melts, Thiol

One possible solution to acidification and losses of base cations in central Ontario forest soils may be the application of wood ash. Wood ash is generally high in pH and contains large amounts of calcium (Ca) and other nutrients essential for ecosystem health, however it also contains trace metals. Understanding the chemistry of soils following ash application to forests is crucial for future policy recommendations and remediation efforts. In this study, soil and soil water chemistry was measured at two acidic forest sites in central Ontario. Sugar maple (Acer saccharum, Marsh.) seedling growth and chemistry, as well as understory vegetation composition, were also measured. At site one, plots (2 m x 2 m) were established with sugar maple, white pine (Pinus strobus L.) and yellow birch (Betula alleghaniensis Britt.) residential wood ash treatments and applied at rates of 0 and 6 Mg ha-1. The effects of residential wood ash on soil and understory vegetation were measured three- and 12-months following ash addition. At site two, plots (5 m x 5 m) were established with both fly and bottom industrial grade bark ash treatments of 0, 4 and 8 Mg ha-1 (n=4), and tension lysimeters were positioned in each plot at 30, 50, and 100 cm depths. The effects of industrial grade wood ash on soil, soil water and understory vegetation were measured four years following ash addition. Metal concentrations in the ashes were generally low but were higher in the fly ash and yellow birch ash types. At site one, significant increase in soil pH, and Ca and magnesium (Mg) concentrations were observed after three months, however changes varied by treatment. Some metal concentrations increased in the upper organic horizons, but metals were likely immobilized in the soil due to increases in soil pH, electrical conductivity (EC) and high organic matter content of the soil. After one year, changes to metal concentrations in soils could be seen in mineral horizons, and a few metals (aluminum (Al), zinc (Zn), copper (Cu), chromium (Cr), strontium (Sr)) increased in treatment plots. At site two, the effects of industrial-grade bark ashes on soil pH could still be seen after four years and soil water metal concentrations were not elevated relative to controls. Changes to understory vegetation composition following ash application were observed, but ash addition had no significant effect on sugar maple seedling growth (root:shoot ratio) and did not lead to significant increases in foliar metal concentrations. There were significant differences in root chemistry, suggesting metal translocation and uptake could be restricted. Mass balance estimates indicate that the organic horizon is a sink for all metals and simulated drought in this horizon led to a decrease in soil pH and increase in soil water metal concentration, but this occurred in all treatments including control. These results suggest that application of industrial and residential wood ash in moderate doses with trace metal concentrations below or near regulatory limits will increase soil pH and base cation concentrations, as well as increase seedling tissue nutrient concentrations in northern hardwood forest soils. However, depending on the parent material of the ash, increased metal availability can also occur.

Author Keywords: Acer saccharum, calcium decline, forest soil amendment, Haliburton Forest and Wildlife Reserve, heavy metal, wood ash

Natural establishment of vegetation on mine tailings is generally limited. Understanding the processes leading to vegetation germination and the survival mechanisms that vegetation species employ in these harsh environments is critical to future remediation efforts. As metalliferous mine tailings are generally nutrient-poor, high in harmful metals, and acidic, vegetation species require distinct mechanisms to germinate and survive in such harsh environments. In this study, edaphic and biotic factors linked to vegetation establishment and diversity were studied at two nickel-copper (Ni-Cu) tailings sites near Sudbury, Ontario. One site had experienced minimal treatment, and the second site was split into partial (hand-distribution of lime) and full (lime, fertilizer, seeding) treatment areas. Tailings were generally acidic, low in organic matter and "available" nutrients, and high in metals such as Al, Cu, Fe, and Ni, but these physical and chemical properties were extremely spatially variable. At both sites, vegetation was distributed in sparse patches, with the greatest diversity in treated areas. There was no clear link between metals and vegetation establishment/diversity at the sites. The primary limiting nutrients on the tailings were phosphorous (P) and potassium (K), and while there were areas of increased soil fertility at the sites, they were not clearly associated with increased vegetation diversity. Both traditional ecological succession and nucleation succession patterns were observed on the site, and the chief species associated with nucleation were primary colonizing trees such as B. papyrifera and P. tremuloides. The relationship between B. papyrifera nutrient retranslocation and tailings restoration was assessed and while B. papyrifera at the sites were deficient in P and K, the trees efficiently retranslocated both P and K during senescence. This research can provide insight into possibilities for future revegetation of similar tailings, enabling industry to make educated decisions when choosing where and how to revegetate, mimicking natural succession patterns.

Author Keywords: Acid-mine drainage, Betula papyrifera, ecosystem health, metals, Sudbury, tailings

Dissolved organic matter (DOM) in oceans provides nutrients and ultraviolet radiation protection to microbes. Some DOM compounds can chelate with metals, including copper, controlling their transport and bioavailability in marine systems. As copper functions as both a nutrient and toxicant, studies into the chemical structures of Cu-ligands is important, however currently limited. In this thesis, the chemical composition of total and Cu-binding DOM is investigated using Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) in the Canadian Arctic and North Pacific. Chapter 2 reveals chemical differences in DOM composition between the southern and northern Canada Basin, revealing the influence of terrestrial and biological sources. Chapter 3 shows the uniqueness of Cu-binding ligands found in the Canadian Arctic and North Pacific Ocean. Studying the composition of DOM gives insight into the chemical diversity of marine DOM, helping to predict the effects of a changing climate on marine ecosystems.

Author Keywords: biological, dissolved organic matter, fluorescence, immobilized metal-ion affinity chromatography, mass spectrometry, terrestrial

Decades of acidic deposition have depleted base cation pools in soil over large parts of eastern north America, including the Muskoka-Haliburton region of central Ontario. This region has also experienced a shift in forest species composition over the past 200 years, favouring sugar maple (Acer saccharum Marsh.) at the expense of species such as white pine (Pinus strobus L.) and eastern hemlock (Tsuga canadensis (L.) Carr.). This shift in species composition may have changed soil chemistry over time due to differences in nutrient and metal inputs in litterfall. An analysis of litterfall and soil chemistry was conducted for five tree species commonly found across central Ontario. Stands were established in the Haliburton Forest & Wild Life Reserve and were dominated by one of balsam fir (Abies balsamea (L.) Mill.), eastern hemlock, white pine, sugar maple, or yellow birch (Betula alleghaniensis Britt.). Analysis of mineral soil oxides suggested that these stands were established on similar parent material. Deciduous dominated stands (maple and birch) had greater litterfall mass compared with conifer stands (fir, hemlock, and pine), generally leading to greater macronutrient inputs to the soil. Elemental cycling through the organic horizons was more rapid in deciduous stands, with base cations having the shortest residence times. This suggests that a change from greater conifer dominance to mixed hardwood forests may lead to more rapid elemental cycling and alter the distribution of elements in soil. Forests in the region are typically mixed and the resulting differences in soil chemistry may influence model predictions of soil recovery from acidification. Laboratory leaching tests indicated that both stand type and the acidity of simulated rainwater inputs influenced soil solution chemistry, with deciduous stands generally having a greater buffering capacity than sites dominated by coniferous species. Changes in soil chemistry were examined for each stand type using the Very Simple Dynamic (VSD) biogeochemical model. Simulations showed that soil base saturation began to increase following lows reached around the year 2000, and similar patterns were observed for all stands. When sulphur (S) and nitrogen (N) deposition were held constant at present rates, soil base saturation recovery (toward pre-1900 levels) was marginal by 2100. With additional deposition reductions, further increases in base saturation were minor at all sites. In conjunction with additional deposition reductions, the elimination of future forest harvesting allowed for the greatest potential for recovery in all stands. Overall, these results suggest that changes in forest cover may influence soil chemistry over time, most notably in the organic soil horizons. However, forecasted recovery from acidification is expected to follow similar patterns among stands, since differences in soil chemistry were less significant in the mineral soil horizons which compose a greater proportion of the soil profile.

Author Keywords: base cation decline, forest harvesting, litterfall, mineral weathering, soil acidification, VSD model

Lake of the Woods (LOW) is a large international waterbody which suffers from frequent and widespread algae blooms. Previous studies have highlighted the importance of the lake's largest tributary, the Rainy River (RR) and its significance in total phosphorus (TP) delivery to the LOW. Unfortunately, little is known about TP contributions from the RR and its tributaries within the Canadian portion of the watershed. This thesis examines patterns and sources of TP from four tributaries on the Canadian side of the lower RR region, two of which are predominantly natural, and two that are predominantly agricultural. Relationships between water quality parameters, land use and geologic characteristics were observed over a complete hydrologic year (Oct 1, 2018 - Sept 31, 2019), and through an intensive sampling campaign using a nested watershed approach during the spring high flow and summer low flow periods. Results revealed that TP and total suspended sediment (TSS) concentrations (>100 µg/L and >20 mg/L respectively), and loads (>20 kg/km2 and >3500 kg/km2, respectively), were greater at agricultural sites compared with natural sites (<65 µg/L TP and <15 mg/L TSS concentration, and <20 kg/km2 TP and <4000 kg/km2 TSS export). Total P, TSS, Fe, and Al were significantly positively correlated (R2= 0.26-0.59; p<0.05) and intensive sampling revealed that these relationships were strongest during the spring and at the agricultural sites (R2= 0.73-0.98; p<0.05). In contrast, the summer intensive sampling revealed that TP and redox sensitive Fe were significantly correlated (R2= 0.72; p<0.005), whereas redox insensitive Al and TSS were not, suggesting TP may be sourced via redox processes in the summer due to favourable hydrologic conditions. This was observed not only at sites with high wetland influence, but also at sites with more agricultural presence suggesting that redox sourced TP may also originate from mineral stream bed sediment during low flow periods. This research suggested two primary TP sources in the lower RR region: erosion in the spring, and redox processes (internal release) in the summer. It is recommended that intensive monitoring continue in Canada, and further research be conducted to fully understand the significance of internal P release in the tributaries.

Author Keywords: erosion, land use, nutrients, particulates, redox, water quality

In this study we seek to better understand the potential effects of short-term (5-year) N fertilization on jack pine forest biogeochemistry, vascular plant community composition and to project a temporal endpoint of nitrogen leaching below the major rooting zone. Aqueous ammonium nitrate (NH4NO3) was applied above the forest canopy across five treatment plots (20 x 80 m) four times annually. The experimental deposition gradient followed those known for localized areas around the major open pit operations at 0, 5, 10, 15, 20 and 25 kg N ha-1 yr-1 over a five-year period (2011 – 2015). Nitrate recovery in throughfall was significantly higher than NH4+ (p < 0.05), indicating canopy NH4+ immobilization. There was a strong treatment effect (p < 0.05) of N on the epiphytic lichen thalli concentrations of Hypogymnia physodes and Evernia mesomorpha after five years. The canopy appeared to approach saturation at the highest deposition load (25 kg N ha-1 yr-1) during the fifth year of N additions as most N added above the canopy was accounted for in throughfall and stemflow. The non-vascular (lichen and moss) vegetation pool above the forest floor was the largest receptor of N as cryptogam foliar and thalli N concentrations showed a significant treatment effect (p < 0.05). Nitrogen in decomposing litter (25 kg N ha-1 yr-1) remained immobilized after five years, while treatments ≤ 20 kg N ha-1 yr-1 started to mobilize. Understory vascular plant cover expansion was muted when deposition was ≥ 10 kg N ha-1 yr-1. Finally, modeling suggests the forest may not leach N below the rooting zone until around 50 years after chronic addition begin (25 kg N ha-1 yr-1). The modeling results are consistent with empirical data from a high exposure (~20 - 25 kg N ha-1 yr-1) jack pine site approximately 12 km west of the experimental site that has not yet experienced N leaching.

Author Keywords: Biogeochemistry, Canopy, Deposition, Jack Pine, Nitrogen, Understory