Simulation of the dissolved oxygen at the bottom of deep reservoirs using an existing one-dimensional mathematical model is discussed. A modification to the BOD and detrital components is presented and implemented in the mathematical model. A simulation of Boundary Reservoir in Saskatchewan is included to demonstrate the simulative capacity of the modified model.