Aherne, Julian

The objective of this thesis was to assess the influence of anthropogenic sulphur (S) and nitrogen (N) deposition, and harvesting on soil acidity and calcium (Ca2+), magnesium (Mg2+), potassium (K+) and N soil pools in plantation forest soils in Ireland. The response to reductions in anthropogenic S deposition was assessed using temporal trends in soil solution chemistry at two long-term monitoring plots--one on a blanket peat, the other on a peaty podzol. At the peat site, there was little evidence of a response to reductions in throughfall non marine sulphate (nmSO₄2-) and acidity; soil water acidity was determined by organic acids. In addition, temporal variation in soil water did not respond to that in throughfall. In the podzol, reductions in anthropogenic S and H+ deposition led to a significant improvement in soil water chemistry at 75 cm; pH increased and total aluminum (Al_tot) concentrations declined. The impact of harvest scenarios on exchangeable Ca2+, Mg2+ and K+ pools was assessed using input-output budgets at 40 sites (30 spruce, 10 pine). Harvest scenarios were stem-only harvest (SOH), stem plus branch harvest (SBH) and stem, branch and needle harvest (whole-tree harvesting; WTH). Average K+ and Mg2+ budgets were positive under these scenarios. However, exchangeable K+ pools were small and due to uncertainty in K+ budgets, could be depleted within one rotation. Average Ca2+ budgets for spruce were balanced under SOH, but negative under SBH and WTH. Nitrogen deposition was high, between 5 and 19 kg N ha-1 yr-1, but was balanced by N removal in SOH. However, N budgets were under SBH and WTH, indicating that these harvesting methods would lead to depletion of soil N over the long-term. Finally, monitoring of N cycling at a spruce plot indicated that N deposition was contributing to large NO₃- leaching, and as such the site was N saturated. However, N cycling did not fit the criteria of the N saturation hypothesis; instead leaching was directly related to N deposition and supported the model of kinetic N saturation.

Author Keywords: acidic deposition, base cations, input-output budgets, Ireland, nitrogen, whole-tree harvesting

Models of partitioning, uptake, and toxicity of neutral organic chemicals in fish

Alena Kathryn Davidson Celsie

A novel dynamic fugacity model is developed that simulates the uptake of chemicals in fish by respiration as applies in aquatic toxicity tests. A physiologically based toxicokinetic model was developed which calculates the time-course of chemical distribution in four tissue compartments in fish, including metabolic biotransformation in the liver. Toxic endpoints are defined by fugacity reaching a 50% mortality value. The model is tested against empirical data for the uptake of pentachloroethane in rainbow trout and from naphthalene and trichlorobenzene in fathead minnows. The model was able to predict bioconcentration and toxicity within a factor of 2 of empirical data. The sensitivity to partition coefficients of computed whole-body concentration was also investigated. In addition to this model development three methods for predicting partition coefficients were evaluated: lipid-fraction, COSMOtherm estimation, and using Abraham parameters. The lipid fraction method produced accurate tissue-water partitioning values consistently for all tissues tested and is recommended for estimating these values. Results also suggest that quantum chemical methods hold promise for predicting the aquatic toxicity of chemicals based only on molecular structure.

Author Keywords: COSMOtherm, fish model, fugacity, Partition coefficient, tissue-water, toxicokinetics

Microfibres, which are threadlike particles < 5 mm, are the most common type of microplastic reported in the environment. However, few studies have focused on their abundance in background natural environments. This study assessed the abundance of microfibres in rainfall samples (from four precipitation monitoring stations) and across three headwater lake catchments that were in remote, undeveloped areas, away from anthropogenic disturbance and anthropogenic emission sources (i.e., sites were background natural environments). Anthropogenic microfibres were observed in all samples using visual identification methods, with Raman spectroscopy confirming the presence of polyester film and synthetic pigments, e.g., indigo and hostasol green. The estimated annual average atmospheric deposition of microfibres was ~28,800 mf m-2. Meteorological variables, e.g., rain, wind direction, and relative humidity were correlated with the abundance of microfibres. The average abundance of microfibres in headwater lake catchments was 24 mf g-1 in moss, 0.70 mf m-3 in surface trawl, 9,690 mf m-3 in subsurface, 910 mf kg-1 in lake sediment and 576 mf kg-1 in lakeshore sediment.

Author Keywords: Atmospheric Deposition, Background Environments, Headwater Lake Catchments, Microfibres, Microplastics, Rainfall

Soils play a crucial role in ecosystem functioning, for example, soil minerals provide important provisioning and regulate ecosystem services. This study used major soil oxides from the North American Soil Geochemical Landscapes Project (n=560) to assess elemental associations and infer soil minerals through exploratory data analysis and to determined quantitative soil mineralogy using a normative method, Analysis to Mineralogy (n=1170). Results showed elemental variability of oxides across the provinces of Canada and strong correlations occurred between elements indicative of soil mineral composition (e.g., Silicon and Aluminium). Principal component analysis inferred soil minerals from soil oxides trends on biplots and classified minerals, generally, as carbonates, silicates, and weathered secondary oxides. Spatial variability in minerals (quartz, plagioclase, potassium feldspar, chlorite, and muscovite) was related to the underlying bedrock geology. The use of Analysis to Mineralogy led to a reliable method of quantifying soil minerals at a large scale.

Author Keywords: Analysis to Mineralogy, Exploratory data analysis, Normative procedures, North American Soil Geochemical Landscapes Project, Soil geochemistry, Soil mineralogy

Recent developments at an aluminum (Al) smelter in Kitimat, BC resulted in a permitted increase of 27 to 42 tonnes of sulphur dioxide (SO2) emissions per day. Gaseous SO2 is a pollutant known to contribute to acidic deposition through processes of wet and dry deposition and can additionally react in-atmosphere to form particulate sulphate (pSO42-). Between June 2017 to October 2018, an extensive network consisting of ion exchange resin (IER) column, passive-diffusive, and active filter-pack samplers was established to provide an estimate of total annual S deposition and pSO42- variation throughout the Kitimat Valley. Filter-pack sampling determined the relative concentration of pSO42- increased downwind of the smelter. Comparison of observation-based and modelled total annual deposition suggested CALPUFF was accurate in representing the spatial viability of S deposition (R2 = > 0.85). However, the model appeared to overpredict near-field deposition suggesting the potential of underestimation further downwind of the smelter.

Author Keywords: aluminum smelter, atmospheric deposition, filter-pack sampler, ion-exchange column sampler, pSO42-, SO2

This study examined dissolved organic matter (DOM) and labile Mercury (from diffusive gradient in thin films (DGT)) in the Quesnel river, British and the Goose creek tributary of the Churchill river, Manitoba. DOM properties were measured with optical measurements of absorption coefficient (a254), spectral slopes (S275-295, SR) and fluorescence indices (HIX, BIX, FI). The DOC proxy measurements (a254) were almost 10 times higher at the Churchill site (Mean a254 116.77 cm-1) compared with the Quesnel river site (Mean a254 12.06 cm-1) during the study periods. While DGT labile Hg concentrations at the Quesnel site (2.17 to 98.97ppt) were almost 10 times more than the levels reported at the Churchill site (0.03 to 9.06 ppt). Fluorescence indices and the rise of labile Hg concentrations in spring indicated mostly terrestrial sources of DOM at both the sites. Spectral slopes and fluorescence indices substantiated that Churchill site had high molecular weight, complex and more humic DOM compared with Quesnel. DOM at both the sites was prone to temporal variation and affected by environmental conditions. Correlation between DGT labile-Hg and DOM parameters suggested that DGT collected Hg-organic complexes along with inorganic labile-Hg complexes.

Author Keywords: Churchill, Diffusive gradient in thin films, Dissolved organic matter, Labile Hg, Mercury, Quesnel

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

Anthropogenic atmospheric emissions and subsequent deposition of nitrogen (N) can affect N-sensitive habitats and lead to shifts in plant species community composition. This study assessed the effects of N deposition on plant community composition for Jack pine forests across northwestern Canada and across a smaller subset of sites surrounding the Athabasca Oil Sands Region (AOSR) using 'gradient forest' analysis. Predictor influence on community composition varied depending on the scale of the study and relatively

distinct thresholds were identified for different plant groups. In the larger scale study, a total deposited nitrogen (TDN) threshold of 1.5 – 3 kg N ha-1yr-1 was well suited to protect predominantly lichen species, consistent with lichen-based critical loads from other studies. Across the smaller scale study, a TDN threshold of 5.6 kg N ha-1yr-1 was primarily associated with vascular species changepoints but did include some important N-indicator lichen and bryophyte species.

Author Keywords: critical loads, gradientForest, Jack pine, Nitrogen deposition, species composition

Water balance models calculate water storage and movement within drainage basins, a primary concern for many hydrologists. A Thornthwaite water balance model (H2OBAAS) has shown poor accuracy in predicting low flows in the Petawawa River basin in Ontario, so lake storage and winter snow processes were investigated to improve the accuracy of the model. Lake storage coefficients, represented by the slopes of lake stage vs. lake runoff relationships, were estimated for 19 lakes in the Petawawa River basin and compared on a seasonal and inter-lake basis to determine the factors controlling lake runoff behaviour. Storage coefficients varied between seasons, with spring having the highest coefficients, summer and fall having equal magnitude, and winter having the lowest coefficients. Storage coefficients showed positive correlation with lake watershed area, and negative correlation with lake surface area during summer, fall, and winter. Lake storage was integrated into the H2OBAAS and improved model accuracy, especially in late summer, with large increases in LogNSE, a statistical measure sensitive to low flows. However, varying storage coefficients with respect to seasonal lake storage, watershed area, and surface area did not improve runoff predictions in the model. Modified precipitation partitioning and snowmelt methods using monthly minimum and maximum temperatures were incorporated into the H2OBAAS and compared to the original methods, which used only average temperatures. Methods using temperature extremes greatly improved simulations of winter runoff and snow water equivalent, with the precipitation partitioning threshold being the most important model parameter. This study provides methods for improving low flow accuracy in a monthly water balance model through the incorporation of simple snow processes and lake storages.

Author Keywords: Lake Storage, Model Calibration, Monthly Water Balance, Petawawa River, Precipitation Partitioning, Snow Melt

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