Simulation of wind-driven circulation in the Salton Sea: implications for indigenous ecosystems

Abstract

The Salton Sea, the largest ‘man-made’ water body wholly within California, was formed in 1904 as the result of a levee failure along the Colorado River. Initially, flow into the Salton ‘Sink’ created a fresh water lake about 24 m deep with a water surface about 85 m below the level of the ocean. Salinity of the water body, at first roughly the same as the river, rose rapidly due to solution of previously accumulated residual salt, then following levee repair, adjusted to the combined influence of agricultural drainage accretions and evaporative losses. Water levels adjusted accordingly, at first declining then rising slowly until the mid 1930s when a level about 75 m below ocean level was reached. Thereafter, both water levels and salinity gradually rose, so that at present the surface elevation of the Sea stands near −69.5 m and salinity is approaching 45 g 1−1, about 30% above ocean salinity. The Salton Se: Authority is seeking practical methods for reducing water levels and controlling salinity within ranges that will protect beneficial uses of the Sea, its adjacent lands, and its indigenous ecosystems, both aquatic and avian. Pro posed solutions include various physical changes in the bathymetry and configuration of the Sea, especially its southern basin. Because circulation in the Sea is driven primarily by wind stresses imposed on the water surface, and circulation changes are likely to affect the Sea’s quality and ecology, a methodology for quantifying the effects of specific alternatives is required. For this purpose a mathematical model for simulation of the hydrodynamic behavior of the Sea has been developed, calibrated to data gathered by a field investigation conducted in 1997, and applied to alternative schemes that will isolate sections of the southern basin, thus changing the natural wind induced circulation in areas that are ecologically sensitive. The Salton Sea Hydrodynamic/Water Quality Model is constructed using the finite element method to represent the bathymetry of the Sea as it currently exists, or may subsequently be modified, in a three-dimensional grid. Given certain boundary conditions, for example wind stresses imposed on the surface, the model solves the three dimensional equations of motion and continuity, the advection—dispersion equation, and an equation of state dependent upon temperature and salinity, to obtain temporal and spatial descriptions of velocities and temperatures over a specified period of time. Using data derived from 1997 field measurements of velocities using acoustic Doppler current profilers and temperature sensors, the model was calibrated to reproduce mathematically the historic experience of field observation. The model successfully replicated principal features of the Sea’s behavior, especially the persistence of a counterclockwise gyre in the southern basin and seasonal stratification. Once calibrated, the model was applied to evaluate the possible effects of changing water surface elevations in the Sea and altering its configuration to isolate sections for evaporative concentration of salts. These effects, evident in changes in velocity, were quantified with regard to their possible impacts on the aquatic habitat and the health of the Salton Sea ecology. A comparative evaluation of alternatives is presented.