In this article we consider a linear programming model for assessing the aggregate impact of land‐use activities scattered over a large area on the resultant pollutant concentrations in coastal waters. The dispersion to coastal waters of the adverse environmental loads generated by the land uses is described by a set of transport coefficients that measure the attenuation of pollutants, such as industrial BOD, carried to the coast along surface drainage basins. Further dispersion in the waters caused by tidal action is then described by the simple procedure know as “pollution susceptibility.”; The model seeks to minimize the steady‐state concentrations of pollutants by establishing an optimal spatial configuration of residential, commercial, and industrial land uses. This configuration is constrained by a number of restrictions based on local and regional targets for growth and development. The methodology discussed in this article is intended to be useful to regional planners, and is based on a study conducted for the Long Island area.