Sediment oxygen demand (SOD) contributes immensely to hypolimnetic oxygen depletion. SOD rates thus play a key role in aquatic ecosystems' health predictions. These rates, however, can be very expensive to sample. Moreover, determination of SOD rates by sediment diagenesis modeling may require very large datasets, or may not be easily adapted to complex aquatic systems. Water quality modeling for northern aquatic systems is emerging and little is known about the seasonal trends of SOD rates for complex aquatic systems. In this study, the seasonal trend of SOD rates for a northern chained river-lake system has been assessed through the calibration of a water quality model. Model calibration and validation showed good agreement with field measurements. Results of the study show that, in the riverine section, SOD 20 rates decreased from 1.9 to 0.79 g/m 2 /day as urban effluent traveled along the river while a SOD 20 rate of 2.2 g/m 2 /day was observed in the lakes. Seasonally, the SOD 20 rates in summer were three times higher than those in winter for both river and lakes. The results of the study provide insights to the seasonal trend of SOD rates especially for northern rivers and lakes and can, thus, be useful for more complex water quality modeling studies in the region.