A two-component soil water flow model was used in conjunction with an equilibrium speciation model WATEQF to study aluminum mobility in soils of a forested watershed, White Oak Run, in the Shenandoah National Park, Virginia. Soil solution samples, taken from the O, E, B, C1, and C2horizons, were collected from zero-tension lysimeters designed to collect faster gravitational macropore flow and tens