Passy, S. I.1,2 & Freehafer, D.2 1Department of Biology, Rensselaer Polytechnic Institute, Troy, NY 12180; 2US Geological Survey, 425 Jordan Rd., Troy, NY 12180‐8349 USAOn August 27, 1999 diatoms were sampled, and current velocity was measured at 81 locations on a regular square sampling grid in an unshaded, cobble‐bottom reach of White Creek, NY. The grid had an extent of 16 m2, interval, the distance between neighboring sampling points, of 0.5 m, and grain size, the size of the elementary sampling unit, of 0.01 m2. Six of the seven dominant benthic diatoms were colonial forms, including Diatoma vulgaris, Fragilaria capucina, F. crotonensis, Gomphoneis minuta, Melosira varians, and Synedra ulna. Their morphology and distribution were investigated from the perspective of fractal geometry and stream ecology, respectively. Fractal dimension of diatom colonies, indicative of their shape complexity, ranged from 1.06 to 1.54, demonstrating vast morphological variation from simple geometric shapes to complex outlines. The relative abundance of the six diatoms was regressed against current regime, which ranged from 0.03 to 0.66 mδ s‐1. All regression models were significant at P < 0.05 and explained between 55% and 94% of the variation in diatom distribution. The diatom niche breadth, i.e. the amount of environmental variability a species can tolerate, was defined from these models and showed substantial variation, from 5 to 14. The regression model of fractal dimension against niche breadth was significant and explained 76% of colonial shape variation, revealing a strong relationship between diatom colonial complexity and habitat variability. This finding suggests that environmental variability could create highly complex colonial morphologies in benthic diatoms as an evo‐lutionary strategy for survival in unpredictable environments.