A variational approach to identify the depth distribution of radiogeni
c heat generation in the crust is utilized. It is inferred that in the
crust, in which at least the upper section represents a porous medium
where fluid circulation may exist, the exponential distribution, incl
uding an additional term corresponding to the removed part of radioele
ments, is optimal.