The relationship between diatom taxa preserved in surface soils and environmental variables at 31 sites in Water Conservation Area 2A (WCA-2A) of the Florida Everglades was explored using multivariate analyses. Surface soils were collected along a phosphorus (P) gradient and analyzed for diatoms, total P, % nitrogen (N), %carbon (C), calcium (Ca), and biogenic silica (BSi). Phosphorus varied from 315-1781 μg g-1, and was not found to be correlated with the other geochemical variables. Canonical correspondence analysis (CCA) was used to examine which environmental variables correlated most closely with the distributions in diatom taxa. Canonical correspondence analysis with forward selection, constrained and partial CCA, and Monte Carlo permutation tests of significance show the most significant changes in diatom assemblages along the P gradient (p < 0.01), with additional species differences correlated with soil C, N, Ca, and BSi.