Under certain flow conditions, fluid flow through porous media starts to deviate from the linear relationship between flow rate and hydraulic gradient. At such flow conditions, Darcy's law for laminar flow can no longer be assumed and nonlinear relationships are required to predict flow in the Forchheimer regime. To date, most of the nonlinear flow behavior data is obtained from flow experiments on packed beds of uniformly graded granular materials (C u = d 60 /d 10 < 2) with various average grain sizes, ranging from sands to cobbles. However, natural deposits of sand and gravel in the subsurface could have a wide variety of grain size distributions. Therefore, in the present study we investigated the impact of variable grain size distributions on the extent of nonlinear flow behavior through 18 different packed beds of natural sand and gravel deposits, as well as composite filter sand and gravel mixtures within the investigated range of uniformity (2.0 < C u < 17.35) and porosity values (0.23 < n < 0.36). Increased flow resistance is observed for the sand and gravel with high C u values and low porosity values. The present study shows that for granular material with wider grain size distributions (C u > 2), the d 10 instead of the average grain size (d 50) as characteristic pore length should be used. Ergun constants A and B with values of 63.1 and 1.72, respectively, resulted in a reasonable prediction of the Forchheimer coefficients for the investigated granular materials.