The effect of flow in the benthic boundary layer on sediments in the northern Gulf of Mexico was investigated by developing a numerical flow model based on high-resolution bathymetry data and ADCP time series data in a 500 500 100 m^3 domain. We determined flow velocities from the Navier-Stokes equations and modeled turbulence using a k- model. Flow velocities were used to assess sediment resuspension by comparing simulated near-bottom flow velocities with experimentally determined critical shear velocities. The percentage of time when simulated flow velocities exceed this threshold was mapped to show the spatial distribution of resuspension probability. In-situ imagery data was compared to resuspension probability through cluster analysis. Results showed that seafloor topography leads to variations in flow velocities across the domain and with predominant flow direction. Simulated resuspension is episodic; some but not all predicted resuspension events coincide with storm events recorded with buoys.