How Cold is the Ground II

Images from Wikipedia

Last week (ok, it was a little more than a few days ago....) I used dimensional analysis to figure out how the ground's temperature changes with time. But although dimensional analysis can give us information about the length scales in the problem, it doesn't tell us what the solution looks like. From dimensional analysis, we don't even know what the solution does at large times and distances. (Although we can usually see the asymptotic behavior directly from the equation.) So let's go ahead and solve the the heat equation exactly:
\[
\frac {\partial T}{\partial t} = a \frac {\partial ^2 T}{\partial x^2} \quad (1)
\]



Well, what type of solution do we want to this equation? We want the temperature at the Earth's surface x=0 to change with the days or the seasons. So let's start out modeling this with a sinusoidal dependence -- we'll look for a solution of the form
\[
 T(x,t) = A(x)e^{i wt}
\]
for some function A(x), then we can take the real part for our solution.

Plugging this into Eq. (1) gives A'' = iw/a * A, or

\[
A(x) = e^{ \pm \sqrt{w/2a } (1+i) x}
\]
Since we have a second-order ordinary differential equation for A, we have two possible solutions, which are like exp(+x) or exp(-x). Which one do we choose? Well, we want the temperature very far away from the surface of the ground to be constant, so we need the solution that decays with distance, A~exp(-x). Taking the real part of this solution, we find [1]
\[T(x,t) = T_0 \cos (wt + \sqrt{w/2a}\times x ) e^{-\sqrt{w/2a}x} \quad (2)
\]
Well, what does this solution say?  As we expected from our scaling arguments last week, the distance scale depends on the square root of the time scale -- if we decrease our frequency by 4 (say, looking at changes over a season vs over a month), the ground gets cooler only 2x deeper. We also see that the temperature oscillation drops off quite rapidly as we go deeper into the ground, and that there is a "lag" the farther you go into the ground. In particular, we see that at distances deep into the ground, the temperature drops to its average value at the surface. You can see this all in the pretty plot below (generated with Python):

Let's recap. To model the temperature of the ground, we looked for a solution to the heat equation which had a sinusoidally oscillating temperature at x=0, and decayed to 0 at large x. We found a solution such a solution, and it shows that the temperature decays rapidly as we go far into the ground. At this point, there are two questions that pop into mind: 1) Is the solution that we found unique? Or are there other possible solutions? 2) This is all well and good, but what if our days or seasons aren't perfect sines? Can we find a solution that describes this behavior?

I'll give one (1) VirtuosiPoint to the first commenter who can prove to what extent the above solution is unique [2]. But how about the second point? Can we solve this for non-sinusoidal time variations?

Well, at this point most of the readers are rolling their eyes and shouting "Use a Fourier series and move on." So I will. Briefly, it turns out that (more or less) any periodic function can be written as a sum of sines & cosines. So we can just add a bunch of sines and cosines together and construct our final solution.

So just for fun, here is a plot of the temperature of the ground in Ithaca (data from Wikipedia) over a year. (I used a discrete Fourier transform to compute the coefficients.)

The temperature (colorbar) is in degrees C, assuming a=0.5 mm^2/s.

Looks pretty boring, but I swear that all the frequencies are in that plot. It just turns out that the seasons in Ithaca are pretty sinusoidal. So about 20 meters below Ithaca, the temperature is a pretty constant 8 C. 

While I was postponing writing this, I wondered what the temperature on Mercury's rocks would be. If we dig deep enough, can we find an area with habitable temperatures? Some quick Googlin' shows that the daytime and nighttime temperatures on Mercury are ~550-700 K and ~110 K at the "equator." While I don't think that Mercury's temperature varies symmetrically, let's assume so for lack of better data.[3] Then we'd expect that deep into the surface, the temperature would be fairly constant in time, at the average of these two extremes. Plugging in the numbers (assuming a~0.52 mm^2/s  and using a Mercurial solar day as 176 days), we get

T=94 C, at 2.75 meters into the surface.

[1] More precisely, since the heat equation is linear and real, if T(x,t) is a solution to the equation, then so are 1/2(T+T*) or 1/2i(T-T*).

[2] Hint: It's not unique. For instance, here is another solution that satisfies the constraints, with no internal heat sources or sinks (I'll call it the "freshly buried" solution):

Can you prove that all the other solutions decay to the original solution? Or is there a second or even a spectrum of steady state solutions?

[3] If someone provides me with better data of the time variation of Mercury's surface at some specific latitude, I'll update with a full plot of the temperature as a function of depth and time.