Temperature Profile in Lake

Lakes become stratified in temperature over the course of the year because of the way that heat transfer occurs. In the spring, the temperature of a lake as a function of depth is reasonably constant. During the summer, as heat enters the lake from solar radiation, the upper layers begin to warm. Since warm fluids tend to rise and the warmest water is already on the surface, there is not much mixing in the lake during this time. As a result, a temperature profile like that attached below tends to form.

Later in the year, the air becomes cooler than the water and the lake loses heat. Eventually, the surface waters cool enough that they sink, causing the lake to mix during the winter and setting up a new cycle for the following year.

During the time of year that the lake is stratified it is of interest to know the depth at which the separation of layers occurs – because physical processes such as transport of pollutants or nutrients are very different in the different layers. The separation of layers is considered to occur at the thermocline, which is defined as the location of the steepest slope in the temperature gradient. Mathematically, this occurs at the inflection point – so the position of the thermocline can be found from the following criterion:

(1)∂2T/∂y2=0

where y is the depth (measured from the lake surface) and T is the temperature.

A mathematical model for temperature as a function of depth y (in m) and time t (in days) is:

(2)T(y,t)-T0/Tsurf(t)-T0=exp(-y2/4αt)

where Tsurf(t) is the water temperature of the lake surface at time t, α is a property called the “eddy thermal diffusivity” and T0 is the lake temperature at time zero. Time zero must be chosen to be on a day when the lake temperature is more or less uniform.

Here are the specific tasks:

(1) Apply equation (1) to equation (2) and develop an expression for the location ytc of the thermocline as a function of time.

(2) The speed at which the thermocline moves vtc is defined as

vtc=∂ytc/∂t

Use your results from (1) to obtain an expression for vtc as a function of time.

2. Relevant equations

[STRIKE]if you need to see the picture it is attached here: [/STRIKE]
https://www.physicsforums.com/attachment.php?attachmentid=57866&d=1365975296

3. The attempt at a solution
First I expressed the function as T(y,t)=e^(-y2/4αt)(Tsurf(t)-T0)+T0....

Then I got the second partial derivative with respect to y which gave me:
ytc=1/4α2t2(Tsurf(t)-T0)e-y2/4αt=0.

Im not sure if this is the way to go so I just wanted to check to see if I was in the right path. Also for number 2 im not sure how to differentiate with respect to t when I have Tsurf(t) as a function of time. Any help would be appreciated.