Abstract: The drill hole SG-3, 12261 m deep in the Pechenga-Zapolyarny area, Kola Peninsula, Russia, is currently the deepest drill hole in the world. Geothermal measurements in the hole reveal a considerable variation (30-68 mW m-2) with depth in the vertical component of heat flow density (HFD). We simulate heat and fluid flow in the bedrock structure of the Kola deep hole site. Various potential sources for the observed HFD variation are discussed, with special emphasis on advective heat transfer, paleoclimatic ground surface temperature changes and refraction of heat flow due to thermal conductivity contrasts. A two-dimensional finite difference (FD) porous medium model of the Kola structure, constructed from all available data on lithology, hydrogeology, topography, thermal conductivity and heat production rate in the deep drilling area, is the basis of all forward model calculations. A conductive, steady-state simulation indicates that heat production and refraction create a variation of about 15 mW m-2 in the uppermost 15 km, but are insufficient to reproduce the measured HFD-depth curve in the uppermost 2-4 km. However, if topography-driven groundwater flow is considered in the model the measured HFD variation is easily explained. The most sensitive parameters in fitting the model results to the observed HFD-depth curve are the permeability of the top 4 km (10-14 - 10-15 m2) and the (constant) HFD applied at the base of the model at 15 km depth (40-50 mWm-2). The paleoclimatic effect for the Kola structure was calculated with a conductive transient simulation. A simplified ground surface temperature history (GTH) of the Kola area was simulated by varying the model's surface temperatures during different intervals of the simulation. Our results indicate that the measured variation in the vertical HFD cannot be explained by the paleoclimatic effect alone, because its amplitude decreases rapidly from about 16 mW m-2 near the surface to less than 2 mW m-2 at depths in excess of 1.5 km.

 



Steady-state simulation of heat and fluid flow in a vertical cross-sectional model showing temperature (isolines),
heat flow density (color coding) and Darcy velocity (arrows).

See full paper in GJI for details.


 



Temperature vs. depth of steady-state model simulations (blue) and measured in the Kola deep hole SG-3 (red).

See full paper in GJI for details.


 



Heat flow density vs. depth of steady-state model simulations (blue) and measured in the Kola deep hole SG-3 (red).

See full paper in GJI for details.


 



Temperature vs. depth of transient model simulations (blue) and measured in the Kola deep hole SG-3 (red).

See full paper in GJI for details.


 



Heat flow density vs. depth of transient model simulations (blue) and measured in the Kola deep hole SG-3 (red).

See full paper in GJI for details.

 

Back to Project and funding status