In 2012, the Rio Tinto aluminum smelter in Kitimat, British Columbia increased sulphur dioxide (SO2) emissions from 27 to 42 tonnes/day. An initial study was conducted to investigate the effect of the increased sulphur (S) deposition on forest soils. A key uncertainty of the initial study was mineral surface area estimations that were applied to critical load calculations. The current study investigates the effect of organic matter (OM) removal techniques on mineral surface area and the ability to predict mineral surface area using pedotransfer functions (PTFs). Mineral surface area was measured on bulk soil samples using BET gas-adsorption. Organic matter was removed from soil samples prior to surface area measurements using a sodium hypochlorite treatment (NaOCl), loss on ignition (LOI) and no treatment. Removal techniques were found to affect surface area measurements; decreasing in the order of LOI> untreated> NaOCl. Particle-size based PTFs developed from other regions were not significantly correlated with measured surface area. A regionally-specific particle-size based function had stronger predictive value of surface area measurements (adjusted R2=0.82). The PTF that best reflected surface area measurements of bulk soil for the Kitimat area used particle-size data as well as kaolinite, the most abundant clay mineral in the region. Surface area values estimated using the particle-size PTF were applied to the PROFILE model to calculate weathering rates. Weathering rates were then input to critical load calculations using steady-state mass balance. These estimates predicted that none of the 24 measured sites are receiving SO2 deposition in exceedance of their critical load.
Author Keywords: acid deposition, critical loads, mineral surface area, mineral weathering, pedotransfer functions, PROFILE