Diversion of heat by Archean cratons: a model for southern Africa

The surface heat flow in the interior of Archean cratons is typically about 40 mW m⁻² while that in Proterozoic and younger terrains surrounding them is generally considerably higher. The eighty-four heat flow observations from southern Africa provide an excellent example of this contrast in surface heat flow, showing a difference of some 25 mW m⁻² between the Archean craton and younger peripheral units. We investigate two possible contributions to this contrast: (1) a shallow mechanism, essentially geochemical, comprising a difference in crustal heat production between the two terrains, and (2) a deeper mechanism, essentially geodynamical, arising from the existence of a lithospheric root beneath the Archean craton which diverts heat away from the craton into the thinner surrounding lithosphere. A finite element numerical model which explores the interplay between these two mechanisms suggests that a range of combinations of differences in crustal heat production and lithospheric thickness can lead to the contrast in surface heat flow observed in southern Africa. Additional constraints derived from seismological observations of cratonic roots, the correlation of surface heat flow and surface heat production, petrological estimates of the mean heat production in continental crust and constraints on upper mantle temperatures help narrow the range of acceptable models. Successful models suggest that a cratonic root beneath southern Africa extends to depths of 200–400 km. A root in this thickness range can divert enough heat to account for 50–100% of the observed contrast in surface heat flow, the remainder being due to a difference in crustal heat production between the craton and the surrounding mobile belts in the range of zero to 0.35 μW m⁻³.     

Mantle heat flow and geotherms for major tectonic units in central Europe

The results of seismic measurements along the deep seismic sounding profile VII and terrestrial heat flow measurements used for construction of heat generation models for the crust in the Paleozoic Platform region, the Sudetic Mountains (Variscan Internides) and the European Precambrian Platform show considerable differences in mantle heat flow and temperatures. At the base of the crust variations from 440–510°C in the models of Precambrian Platform to 700–820°C for the Paleozoic Platform and the Variscan Internides (Sudets) are found. These differences are associated with considerable mantle heat flow variations.The calculated models show mantle heat flow of about 8.4–12.6 mW m^–2 for the Precambrian Platform and 31 mW m^–2 to 40.2 mW m^–2 for Paleozoic orogenic areas. The heat flow contribution originating from crustal radioactivity is almost the same for the different tectonic units (from 33.5 mW m^–2 to 37.6 mW m^–2). Considerable physical differences in the lower crust and upper mantle between the Precambrian Platform and the adjacent areas, produced by lateral temperature variations, could be expected. On the basis of carbon ratio data it can be concluded that the Carboniferous paleogeothermal gradient was much lower in the Precambrian Platform area than in the Paleozoic Platform region.     

On the geothermal regime of some tectonic units in Romania

Heat flow values of 33–58 mW m^–2 were found for the Transylvanian Depression, 45–57 mW m^–2 for the crystalline nucleus of the Eastern Carpathians, and 70–120 mW m^–2 for the Neogene volcanic area. Temperature-depth profile and some geophysical implications of the low values for the Transylvanian Depression are discussed, rendering evident clear-cut differences between this tectonic unit and other Noegene depressions. The heat flow values for the other two investigated tectonic units are usual ones for areas of their age.A preliminary map of the heat flow distribution over the Romanian territory is presented and its relation to other geophysical fields is discussed. A positive correlation was found between gravity and heat flow, and a negative one between crustal thickness and heat flow. A general conclusion could be drawn that the heat flow distribution over the Romanian territory seems to be governed by processes taking place in the upper mantle, rather than by the radioactive decay within the crust.     

Possible geothermal resources in the Coast Plutonic Complex of southern British Columbia,Canada

In southern British Columbia the terrestrial heat flow is low (44 mW m^–2) to the west of the Coast Plutonic Complex (CPC), average in CPC (50–60 mW m^–2),and high to the east(80–90 mW m^–2). The average heat flow in CPC and the low heat generation (less than 1 W m^–3) indicate that a relatively large amount of heat flows upwards into the crust which is generally quite cool. Until two million years ago the Explorer plate underthrust this part of the American plate, carrying crustal material into the mantle. Melted crustal rocks have produced the inland Pemberton and Garibaldi volcanic belts in the CPC.Meager Mountain, a volcanic complex in the CPC 150 km north of Vancouver, is a possible geothermal energy resource. It is the product of intermittent activity over a period of 4 My, the most recent eruption being the Bridge River Ash 2440 y B.P. The original explosive eruption produced extensive fracturing in the granitic basement, and a basal explosion breccia from the surface of a cold brittle crust. This breccia may be a geothermal reservoir. Other volcanic complexes in the CPC have a similar potential for geothermal energy.Earth Physics Contribution No. 704.     

Heat flow map of the European part of the U.S.S.R.

A heat flow isoline map is presented. Low and relatively constant heat flow has been observed in the old shield areas of the East European Platform (25–40 mW/m²). Increased heat flow (>50 mW/m²) has been found in the Dniepr-Donetz depression. The area south of the East European Platform is characterized by highly variable heat flow (55–100 mW/m²). Some geophysical implications are discussed.     

Review of heat flow data from the eastern Mediterranean region

Heat flow data from the eastern Mediterranean region indicates an extensive area of low heat flow, spreading over the whole basin of the Mediterranean east of Crete (Levantine Sea), Cyprus, and northern Egypt. The average of the marine heat flow measurements in the Levantine Sea is 25.7±8.4 mW/m², and the heat flow on Cyprus is 28.0±8.0 mW/m². The estimated values of heat flow in northern Egypt range from 38.3±7.0 to 49.9±9.3 mW/m², apparently with no consistent trend. To the east, on the coast of Israel, the heat flow values increase, ranging from 36.6±22.4 to 56.7±14.2 mW/m² along a SSE trend. The trend apparently correlates with an increase in crustal thickness, which is about 23 km at the north-west base of the Nile-Delta-cone, and close to 40 km beneath Israel.Contribution No. 157, Department of Geology, Kent State University, Kent, Ohio, USA.     

Heat flow and the geothermal potential of Egypt

Preliminary heat flow values ranging from 42 to 175 mW m^–2 have been estimated for Egypt from numerous geothermal gradient determinations with a reasonably good geographical distribution, and a limited number of thermal conductivity determinations. For northern Egypt and the Gulf of Suez, gradients were calculated from oil well bottom hole temperature data; east of the Nile, and at three sites west of the Nile, gradients were calculated from detailed temperature logs in shallow boreholes. With one exception, the heat flow west of the Nile and in northern Egypt is estimated to be low, 40–45 mW m^–2, typical of a Precambrian Platform province. A local high, 175 mW m^–2, is probably due to local oxidational heating or water movement associated with a phosphate mineralized zone. East of the Nile, however, including the Gulf of Suez, elevated heat flow is indicated at several sites, with a high of 175 mW m^–2 measured in a Precambrian granitic gneiss approximately 2 km from the Red Sea coast. These data indicate potential for development of geothermal resources along the Red Sea and Gulf of Suez coasts. Water geochemistry data confirm the high heat flow, but do not indicate any deep hot aquifers. Microearthquake monitoring and gravity data indicate that the high heat flow is associated with the opening of the Red Sea.     

Combined geothermal-geophysical models of the earth's crust and upper mantle for the European continent

An attempt is made to obtain a combined geophysical model along two regional profiles: Black Sea— White Sea and Russian Platform—French Central Massif. The process of the model construction had the following stages: 1. The relation between seismic velocity (V_p, km/s) and density (σ, g/cm³) in crustal rocks was determined from seismic profiles and observed gravity fields employing the trial and error method. 2. Relations between heat production HP (μW/m³), velocity and density were established from heat flow data and crustal models of old platforms where the mantle heat flow HF_M is supposed to be constant. The HF_M value was also determined to 11 ± 5 mW/m². 3. A petrological model of the old platform crust is proposed from the velocity-density models and the observed heat flow. It includes 10–12 km of acid rocks, 15–20 km of basic/metamorphic rocks and 7–10 km of basic ones. 4. Calculation of the crustal gravity effects; its substraction from the observed field gave the mantle gravity anomalies. Extensively negative anomalies have been found in the southern part of Eastern Europe (50–70 mgal) and in Western Europe (up to 200 mgal). They correlate with high heat flow and lower velocity in the uppermost mantle. 5. A polymorphic advection mechanism for deep tectonic processes was proposed as a thermal model of the upper mantle. Deep matter in active regions is assumed to be transported (advected) upwards under the crust and in its place the relatively cold material of the uppermost mantle descends. The resulting temperature distribution depends on the type of endogeneous regime, on the age and size of geostructure. Polymorphic transitions were also taken into account.     

Thermal modeling of the Southern Alps,New Zealand

Finite-element modeling of the thermal regime across the Southern Alps of New Zealand has been carried out along two profiles situated near the Franz Josef and Haast valleys. The modeling involves viscous deformation beneath the Southern Alps, including both uplift and erosion, and crustal/lithospheric thickening, as a result of crustal shortening extending to 20 mm/y of a 25-km thick crust. Published uplift rates and crustal thickness variations along the two profiles are used to constrain the modeled advection of crustal material, and results are compared with the recent heat flow determinations, 190±50 mW/m² in the Franz Josef valley and 90±25 mW/m² in the Haast valley. Comparisons of the model with published K–Ar and fission track ages, show that the observed heat flow in the Franz Josef valley is consistent with observed zircon fission track ages of around 1 Ma, if the present-day uplift rate is close to 10 mm/y. Major thermal differences between the Franz Josef and Haast profiles appear to be due to different uplift and erosion rates. There is weak evidence that frictional heating close to the Alpine fault zone is not significant. The modeling provides explanations for the distribution of seismicity beneath the Southern Alps, and predicts a low surface heat flow over the eastern foothills due to the dominant thermal effect of crustal thickening beneath this region. Predicted temperatures at mid-crustal depth beneath the zone of maximum uplift rate are 50–100°C cooler than those indicated in previously published models, which implies that thermal weakening of the crust may not be the main factor causing the aseismicity of the central Southern Alps. The results of the modeling demonstrate that the different types of reset age data in the region within 25 km of the Alpine fault are critical for constraining models of the deformation and the thermal regime beneath the Southern Alps.     

Heat flow in Eastern Egypt: The thermal signature of a continental breakup

The Red Sea is a modern example of continental fragmentation and incipient ocean formation. Heat flow data have been collected from eastern Egypt to provide information relating to the mode and mechanism of Red Sea opening. Preliminary heat flow data, including new data reported here, are now available from twenty-five sites in eastern Egypt and one site in western Sinai. A pattern of low to normal heat flow (35–55 mW m⁻²) inland with high heat flow (75–100 mW m⁻²) in a zone within 30 to 40 km of the coast is indicated.Moderately high heat flow (around 70 mW m⁻²) is indicated for the Gulf of Suez. The coastal zone thermal anomaly appears continuous with high heat flow previously reported for the Red Sea shelf. Heat production data indicate that the coastal thermal anomaly is not primarily related to crustal radiogenic heat production. The effects of rapid erosion may contribute to the anomaly, but are not thought to be the primary cause of the anomaly. If the anomaly is caused by lateral conduction from hot, extended, offshore lithosphere, the extension must have been active for the last 30 Ma or so, and a minimum of 100% extension is indicated. Alternatively, the anomaly is primarily caused by high mantle heat flow causing lithospheric thinning, centred beneath the Red Sea. The Red Sea is probably underlain by dominantly basic crust, formed either by intrusion into attenuated continental crust or sea-floor spreading, and for most purposes the crust formed in these two modes of extension may be essentially indistinguishable. Fission-track ages from eastern Egypt indicate that uplift started prior to, or at latest at the time of initial Red Sea opening, and this result, together with thermo-mechanical considerations, suggests an active asthenospheric upwelling beneath the Red Sea and high temperature in the lithosphere prior to extension.     

Terrestrial heat flow in the North Island of New Zealand

Large variations in terrestrial heat flow from 21 to 209 mW/m² have been observed over the North Island, New Zealand. This is generally in good agreement with the pattern of existing geological and geophysical observations. A high heat flow zone with a value of 92 ±3 mW/m², which corresponds to melting temperatures near the base of the crust, is delineated in the northern part of the Taranaki Basin. In the rest of the island, heat flow appears to be low to normal, but some isolated high values are also found. Observed results are interpreted in terms of crust and mantle structure in a region of plate subduction.     

Heat flow in Italy

More than fifty heat flow measurements in Italy are examined. The values, corrected only for local influences (when present), are related to the main geological features with the following results: foreland areas, 55±19 mW m^–2, foredeep areas, 45±21 mW m^–2; folded regions and intermountain depressions, 76±29 mW m^–2. In volcanic areas the heat flow rises to in excess of 600 mW m^–2. From a tectonic point of view, these values are consistent with the hypothesis that the Apennine chain is intersected by two arcuate structures: the first from Liguria to Latium is very probably a continental arc, that is an are which occurs within a continent, and the second from Campania to Calabria is very similar from geophysical evidence to the classic island arcs.     

Modern and ancient geotherms beneath southern Africa

Estimates of subsurface temperatures in the Archean craton of southern Africa during the Archean derived from diamond thermobarometry studies are remarkably similar to temperatures estimated for the same depths today, even though heat production in the earth and the mean global heat flow were probably substantially higher in the Archean. We present multi-dimensional numerical models for the thermal environment of the Archean craton in southern Africa during the Archean in which deep mantle heat is diverted away from the craton toward the surrounding oceanic lithosphere by a lithospheric root beneath the craton. Extrapolation of present-day models to thermal conditions appropriate for the Archean is inadequate to explain the similarity of present-day and Archean temperatures in the cratonic root. Reconciliation of the modern and ancient temperature estimates requires either relaxation of the constraints that the cratonic crustal heat production and/or the earth's mean mantle temperature were higher in the Archean than they are today, or that substantial “erosion” of the lithosphere comprising the cratonic root has occurred since the Archean. The latter possibility could perhaps result from revolatilization of the cratonic root in association with thermal perturbations in the mantle, for which there is evidence in southern Africa in the form of post-Archean tectonic and igneous activity.     

Crustal thickness, discontinuity depth, and upper mantle structure beneath southern Africa: constraints from body wave conversions

The technique of receiver function analysis is applied to the study of crustal and upper mantle structures beneath the Kaapvaal craton in southern Africa and its surroundings. Seismic data were recorded by the seismic array of 82 sites deployed from April 1997 to April 1999 across southern Africa, as well as a dense array of 32 sites near Kimberley, in operation from December 1998 to June 1999. Arrival times for phases converted at the Moho are used to determine crustal thickness. The Moho depth in the south–western section of the craton was found to vary between 37 and 40 km, except for one station that recorded a depth of 43 km (SA23). Farther north along the western block of the craton (into Botswana) the depth increases up to 43 km. The depth increases even further in the north–eastern section of the craton, where results vary from 40 to 52 km. Just north of the Kaapvaal craton, in the neighbouring Zimbabwe craton, the crustal thickness drops significantly. The results obtained there varied from 36 to 40 km. For the Kimberley area, using the dense array, the Moho depth was found to be 37.3 km. Arrivals of the Ps and Ppps phases were used to determine the Poisson’s ratio in the region. This was found to be 0.26±0.01. Arrivals of phases from the 410 and 660 km mantle discontinuities are used to interpret the relative positions of these discontinuities, as well as for comparison of mantle temperatures and seismic velocities in the region with global averages. In the Kimberley area the 410 and 660 km discontinuities were found at their expected depth, implying that mantle temperatures in the region are close to the global average. The seismic velocities above the ‘410’ were found up to 5% faster than the averages from the global iasp91 model, which is fast even by Precambrian standards. In other sections of the Kaapvaal craton, the velocities are also faster than global averages, but not as fast as beneath Kimberley. In these sections, the ‘410’ is also slightly elevated, while the ‘660’ is depressed, which implies a slightly lower mantle temperature relative to the global average. Beneath the Kaapvaal craton we find evidence suggesting the presence of a zone with a reduced wavespeed gradient at an upper bound of approximately 300 km, which may mark the lower chemical boundary of the craton.     

Crustal thermal models for South Siberia

Summary The results of geothermal research carried out in South Siberia (West-Siberian and Siberian Platforms, Altai-Sayan folded area and the Baikalian arched uplift zone) were employed to calculate the Earth's crustal temperatures. It is shown that temperature is a function of the heat flow value. The maps of the surface heat flow and the temperatures at the Moho discontinuity are presented and characteristic temperature-depth profiles were compiled for specific tectonic areas. The West-Siberian and Siberian Platforms, as well as the Altai-Sayan folded region are characterized by rather low heat flow, 42–50mW m ^–2 (1.0–1.2 µcal/cm ² s), and low temperatures at the Moho boundary.400–500 °C. The thermal conditions in the area of Lake Baikal are extreme, the heat flow reaches here100–142 mW m ^–2 (2.4–3.4 µcal/cm ² s), and the temperature at the crust's bottom may exceed1000 °C. The Moho-surface according to our calculations is nonisothermal.     

On the evolution of the geothermal regime of the North China Basin

Recent heat flow and regional geothermal studies indicate that the North China Basin is characterized by relatively high heat flow compared with most stable areas in other parts of the world, but lower heat flow than most active tectonic areas. Measured heat flow values range from 61 to 74 mW m⁻². The temperature at a depth of 2000 m is generally in the range 75 to 85°C, but sometimes is 90°C or higher. The geothermal gradient in Cenozoic sediments is in the range 30 to 40°C/km for most of the area. The calculated temperature at the Moho is 560 and 640°C for surface heat flow values of 63 and 71 mW m⁻², respectively. These thermal data are consistent with other geophysical observations for the North China Basin. Relatively high heat flow in this area is related to Late Cretaceous-Paleogene rifting as described in this paper.     

Thermal-rheological structure of lithosphere beneath the northern flank of Tarim Basin, western China: Implications for geodynamics

Based on the data of geo-temperature and thermophysical parameters of rocks in the Kuqa Depression and the Tabei Uplift, northern flank of the Tarim Basin, in terms of the analytical solution of 1-D heat transfer equation, the thermal structure of the lithosphere under this region is determined. Our results show that the average surface heat flow of the northern flank of the Tarim Basin is 45 mW/m2, and the mantle heat flow is between 20 and 23 mW/m2; the temperature at crust-mantle boundary (Moho) ranges from 514℃ to 603℃ and the thermal lithosphere where the heat conduction dominates is 138-182 km thick. Furthermore, in combination with the P wave velocity structure resulting from the deep seismic sounding profile across this region and rheological modeling, we have studied the local composition of the lithosphere and its rheological profile, as well as the strength distribution. We find that the rheological stratification of the lithosphere in this region is apparent. The lowermost of the lower crust is ductile; however,the uppermost of the mantle and the upper and middle parts of the crust are both brittle layers,which is typically the so-called sandwich-like structure. Lithospheric strength is also characterized by the lateral variation, and the uplift region is stronger than the depression region. The lithospheric strength of the northem flank of the Tarim Basin decreases gradually from south to north; the Kuqa Depression has the lowest strength and the south of the Tabei Uplift is strongest.The total lithospheric strength of this region is 4.77× 1012-5.03 × 1013 N/m under extension, and 6.5 × 1012-9.4× 1013 N/m under compression. The lithospheric brittle-ductile transition depth is between 20 km and 33 km. In conclusion, the lithosphere of the northern flank of the Tarim Basin is relatively cold with higher strength, so it behaves rigidly and deforms as a whole, which is also supported by the seismic activity in this region. This rigidity of the Tarim lithosphere makes it little deform interior, but only into flexure under the sedimentation and tectonic loading associated with the rapid uplift of the Tianshan at its northern margin during the Indian-Eurasian continental collision following the Late Eocene. Finally, the influences of factors, such as heat flow, temperature,crustal thickness, and especially basin sediment thickness, on the lithospheric strength are discussed here.     

Mantle flow and melting beneath oceanic ridge–ridge–ridge triple junctions

Plate boundary geometry likely has an important influence on crustal production at mid-ocean ridges. Many studies have explored the effects of geometrical features such as transform offsets and oblique ridge segments on mantle flow and melting. This study investigates how triple junction (TJ) geometry may influence mantle dynamics. An earlier study [Georgen, J.E., Lin, J., 2002. Three-dimensional passive flow and temperature structure beneath oceanic ridge-ridge-ridge triple junctions. Earth Planet. Sci. Lett. 204, 115–132.] suggested that the effects of a ridge–ridge–ridge configuration are most pronounced under the branch with the slowest spreading rate. Thus, we create a three-dimensional, finite element, variable viscosity model that focuses on the slowest-diverging ridge of a triple junction with geometry similar to the Rodrigues TJ. This spreading axis may be considered to be analogous to the Southwest Indian Ridge. Within 100 km of the TJ, temperatures at depths within the partial melting zone and crustal thickness are predicted to increase by ~ 40 °C and 1 km, respectively. We also investigate the effects of differential motion of the TJ with respect to the underlying mantle, by imposing bottom model boundary conditions replicating (a) absolute plate motion and (b) a three-dimensional solution for plate-driven and density-driven asthenospheric flow in the African region. Neither of these basal boundary conditions significantly affects the model solutions, suggesting that the system is dominated by the divergence of the surface places. Finally, we explore how varying spreading rate magnitudes affects TJ geodynamics. When ridge divergence rates are all relatively slow (i.e., with plate kinematics similar to the Azores TJ), significant along-axis increases in mantle temperature and crustal thickness are calculated. At depths within the partial melting zone, temperatures are predicted to increase by ~ 150 °C, similar to the excess temperatures associated with mantle plumes. Likewise, crustal thickness is calculated to increase by approximately 6 km over the 200 km of ridge closest to the TJ. These results could imply that some component of the excess volcanism observed in geologic settings such as the Terceira Rift may be attributed to the effects of TJ geometry, although the important influence of features like nearby hotspots (e.g., the Azores hotspot) cannot be evaluated without additional numerical modeling.     

The thermal state and composition of the lithospheric mantle beneath the Leizhou Peninsula, South China

Two localities on the Leizhou Peninsula, southern China (Yingfengling and Tianyang basaltic volcanoes) yield a wide variety of mantle-derived xenoliths including Cr-diopside series mantle wall rocks and two distinct types of Al-augite series pyroxenites. Metapyroxenites have re-equilibrated granoblastic microstructures whereas pyroxenites with igneous microstructures have not thermally equilibrated to the mantle conditions. An abundant suite of megacrysts and megacrystic aggregates (including garnet, plagioclase, clinopyroxene, ilmenite and apatite) is interpreted as the pegmatitic equivalents of the igneous pyroxenite suite. Layered spinel lherzolite/spinel websterite xenoliths were formed by metamorphic differentiation caused by mantle deformation, inferred to be related to lithospheric thinning. Some metapyroxenites have garnet websterite assemblages that allow calculation of their mantle equilibration temperatures and pressures and the construction of the first xenolith geotherm for the southernmost China lithosphere. Heat flow data measured at the surface in this region yield model conductive geotherms (using average crustal conductivity values) that are consistent with the xenolith geotherm for the mantle. The calculated mean surface heat flux is 110 mW/m². This high heat flux and the high geotherm are consistent with young lithospheric thinning in southern China, and with recent tomography results showing shallow low-velocity zones in this region. The xenolith geotherm allows the construction of a lithospheric rock type section for the Leizhou region; it shows that the crust–mantle boundary lies at about 30 km, consistent with seismic data, and that the lithosphere–asthenosphere boundary lies at about 100 km.