Denitrification beds are being used to reduce the transport of water-soluble nitrate via subsurface drainage systems to surface water. Only recently has the non-linearity of water flow through woodchips been ascertained. To successfully design and model denitrification beds with optimum nitrate removal, a better understanding of flow in denitrification beds is needed. The main objectives of this study were to characterize the hydraulic properties of old degraded woodchips and provide a better understanding of the factors affecting flow. To achieve this goal, we conducted constant-head column experiments using old woodchips that were excavated from a four-year old denitrification bed near Willmar, Minnesota, USA. For Izbash's equation, the non-Darcy exponent (n) ranged from 0.76 to 0.87 that indicates post-linear regime, and the permeability coefficient (M-10) at 10 degrees C ranged from 0.9 to 2.6 cm s(-1). For Forchheimer's equation, the intrinsic permeability of 5.6 x 10(-5) cm(2) and omega constant of 0.40 (at drainable porosity of 0.41) closely resembled the in-situ properties found in a previous study. Forchheimer's equation was better than that of Izbash's for describing water flow through old woodchips, and the coefficients of the former provided stronger correlations with drainable porosity. The strong correlation between intrinsic permeability and drainable porosity showed that woodchip compaction is an important factor affecting water flow through woodchips. Furthermore, we demonstrated the importance of temperature effects on woodchip hydraulics. In conclusion, the hydraulic properties of old woodchips should be characterized using a non-Darcy equation to help design efficient systems with optimum nitrate removal. (C) 2016 Elsevier B.V. All rights reserved.