Background and Aims Nitrate is the major nitrogen source for many plants. The first step of the nitrate assimilation pathway is the reduction of nitrate to nitrite, catalysed by nitrate reductase (NR). Circadian oscillations in expression and activity of NR have been demonstrated in many plant species. The pathway by which this circadian behaviour is regulated remains to be elucidated. In this study, based on recent experimental observations, a mathematical model is proposed to explain the origin of diurnal and circadian oscillations in NR gene expression and enzyme activity. Methods The dynamic model is based on the feedback interconnections between NR and its substrate, nitrate. In the model, NR activity is regulated at the transcriptional level, in response to the balance between nitrate influx and reduction, and at the post-translational levels in response to signals from carbon assimilation. Conditions for the model system to generate self-sustained circadian oscillations are investigated by numerical simulations. Key Results and Conclusions Under light/dark cycles, the simulation results are in agreement with the observed diurnal pattern of changes in leaf nitrate concentration, NR transcript level and NR activity. Within a range of kinetic parameter values, circadian oscillation behaviour persists even under constant light, with periods of approx. 24 h. These simulation results suggest that sustained circadian oscillations can originate from the feedback interactions between NR and its substrate, nitrate, without the need to postulate the existence of an endogenous ‘circadian clock’.