The thermal performance of a laboratory-scale salt gradient solar pond has been modeled as a one-dimensional unsteady conduction heat transfer problem with heat generation. The pond is assumed to be cut into horizontal slices and finite difference heat balance equations are solved simultaneously to predict the temperature of each slice at any time. The initial conditions were the temperature profile data. The boundary conditions were determined by studying the heat balance at the bottom of the pond and by assuming the pond surface temperature to be equal to the ambient temperature. Solar radiation attenuation is calculated by the Bryant and Colbeck formula. A computer program is constructed to perform the calculations. In addition, Kooi's model was compared with our model. Similarly the salinity behavior was studied by writing the one-dimensional differential mass balance equation over a small slice with the appropriate boundary and initial conditions. The resultant set of linear equations was solved simultaneously for the unknown new concentrations. A computer program has been constructed to perform the calculations. Fair agreement between experimental and predicted profiles was obtained.