An Integrated Study on Physical Properties of a KTB Gneiss Sample and Marble from Portugal: Pressure Dependence of the Permeability and Frequency Dependence of the Complex Electrical Impedance

— Pressure-induced variations in pore geometry were studied on dry- and fluid- saturated samples by means of electrical impedance spectroscopy and permeability measurements. Hydrostatic pressures (up to 120 MPa) and uniaxial pressures (up to failure) were applied. Hydrostatic pressures reduce the aspect ratio of cracks and thus cause a decrease of permeability and electrical bulk conductivity. The opposite was observed in uniaxial pressure experiments where new cracks were formed and consequently permeability and electrical conductivity were increased. More specific informations of these generated observations were derived from the interpretation of the frequency dispersion of the complex electrical conductivity. This least-squares-refinement considers electrochemical interactions between the fluid pore electrolyte and the inner surface of the sample, thus providing informations on the pore geometry and pressure-induced variations. Consequently changes in aspect ratio, size and geometry of the pore system can be detected by means of impedance spectroscopy.     

In-Situ Electrical Conductivity and Permeability of Mid-Crustal Rocks from the Ktb Drilling: Consequences for High Conductive Layers in the Earth Crust

Freshly cored samples from a microprofile (7011–7013m in depth) of the German Continental Deep Drilling Project (KTB) were taken to measure the complex electrical conductivity (1 kHz up to 1 MHz), porosity, BET-surface, permeability and density. The porosity ranged about 1 vol%, while the permeability k varied from 16.05 µD to > 0.01 µD for in-situ pressure conditions. The permeability decreased about 2 orders in magnitude up to pressures of 200 MPa. Conductivity was measured in the same pressure range on 1 M NaCl saturated samples. Thin sections and SEM analysis revealed an enrichment of carbon and ilmenite (about 1 vol%) on inner cleavage cracks of mica, thus causing an unusual high (ranging from 4.2 × 10^-3 S/m to 67 × 10^-3 S/m) being orders of magnitude higher than normally measured on such types of rocks (about 300 × 10^-6 S/m). An inverse pressure dependence of was detected on some of the samples. Electronic conduction was confirmed by least-squares-fits of model data to the frequency dispersion of the conductivity and by measuring the time dependence of the volume conductivity and its frequency dispersion. Thus the dominating role of the reconnected network of carbon and ilmenite on the enhanced volume conductivity was proved. An increase of the conductivity due to hydrofracturing by high pore fluid pressures plays a less important role.     

Changes in geophysical properties caused by fluid injection into porous rocks: analytical models

Analytical models are provided that describe how the elastic compliance, electrical conductivity, and fluid‐flow permeability of rocks depend on stress and fluid pressure. In order to explain published laboratory data on how seismic velocities and electrical conductivity vary in sandstones and granites, the models require a population of cracks to be present in a possibly porous host phase. The central objective is to obtain a consistent mean‐field analytical model that shows how each modeled rock property depends on the nature of the crack population. The crack populations are described by a crack density, a probability distribution for the crack apertures and radii, and the averaged orientation of the cracks. The possibly anisotropic nature of the elasticity, conductivity, and permeability tensors is allowed for; however, only the isotropic limit is used when comparing to laboratory data. For the transport properties of conductivity and permeability, the percolation effect of the crack population linking up to form a connected path across a sample is modeled. However, this effect is important only in crystalline rock where the host phase has very small conductivity and permeability. In general, the importance of the crack population to the transport properties increases as the host phase becomes less conductive and less permeable.     

Electrical Properties of Crustal and Mantle Rocks – A Review of Laboratory Measurements and their Explanation

The electrical properties of rocks and minerals are controlled by thermodynamic parameters like pressure and temperature and by the chemistry of the medium in which the charge carriers move. Four different charge transport processes can be distinguished. Electrolytic conduction in fluid saturated porous rocks depends on petrophysical properties, such as porosity, permeability and connectivity of the pore system, and on chemical parameters of the pore fluid like ion species, its concentration in the pore fluid and temperature. Additionally, electrochemical interactions between water dipoles or ions and the negatively charged mineral surface must be considered. In special geological settings electronic conduction can increase rock conductivities by several orders of magnitude if the highly conducting phases (graphite or ores) form an interconnected network. Electronic and electrolytic conduction depend moderately on pressure and temperature changes, while semiconduction in mineral phases forming the Earth’s mantle strongly depends on temperature and responds less significantly to pressure changes. Olivine exhibits thermally induced semiconduction under upper mantle conditions; if pressure and temperature exceed ~ 14 GPa and 1400 °C, the phase transition olivine into spinel will further enhance the conductivity due to structural changes from orthorhombic into cubic symmetry. The thermodynamic parameters (temperature, pressure) and oxygen fugacity control the formation, number and mobility of charge carriers. The conductivity temperature relation follows an Arrhenius behaviour, while oxygen fugacity controls the oxidation state of iron and thus the number of electrons acting as additional charge carriers. In volcanic areas rock conductivities may be enhanced by the formation of partial melts under the restriction that the molten phase is interconnected. These four charge transport mechanisms must be considered for the interpretation of geophysical field and borehole data. Laboratory data provide a reproducible and reliable database of electrical properties of homogenous mineral phases and heterogenous rock samples. The outcome of geoelectric models can thus be enhanced significantly. This review focuses on a compilation of fairly new advances in experimental laboratory work together with their explanation.     

Comparison of stress-dependent geophysical,hydraulic and mechanical properties of synthetic and natural sandstones for reservoir characterization and monitoring studies

Synthetic rock samples can offer advantages over natural rock samples when used for laboratory rock physical properties studies, provided their success as natural analogues is well understood. The ability of synthetic rocks to mimic the natural stress dependency of elastic wave, electrical and fluid transport properties is of primary interest. Hence, we compare a consistent set of laboratory multi-physics measurements obtained on four quartz sandstone samples (porosity range 20–25%) comprising two synthetic and two natural (Berea and Corvio) samples, the latter used extensively as standards in rock physics research. We measured simultaneously ultrasonic (P- and S-wave) velocity and attenuation, electrical resistivity, permeability and axial and radial strains over a wide range of differential pressure (confining stress 15–50 MPa; pore pressure 5–10 MPa) on the four brine saturated samples. Despite some obvious physical discrepancies caused by the synthetic manufacturing process, such as silica cementation and anisotropy, the results show only small differences in stress dependency between the synthetic and natural sandstones for all measured parameters. Stress dependency analysis of the dry samples using an isotropic effective medium model of spheroidal pores and penny-shaped cracks, together with a granular cohesion model, provide evidence of crack closure mechanisms in the natural sandstones, seen to a much lesser extent in the synthetic sandstones. The smaller grain size, greater cement content, and cementation under oedometric conditions particularly affect the fluid transport properties of the synthetic sandstones, resulting in lower permeability and higher electrical resistivity for a similar porosity. The effective stress coefficients, determined for each parameter, are in agreement with data reported in the literature. Our results for the particular synthetic materials that were tested suggest that synthetic sandstones can serve as good proxies for natural sandstones for studies of elastic and mechanical properties, but should be used with care for transport properties studies.     

Ground penetrating radar inversion in 1-D: an approach for the estimation of electrical conductivity, dielectric permittivity and magnetic permeability

This paper presents a method for inverting ground penetrating radargrams in terms of one-dimensional profiles. We resort to a special type of linearization of the damped E-field wave equation to solve the inverse problem. The numerical algorithm for the inversion is iterative and requires the solution of several forward problems, which we evaluate using the matrix propagation approach. Analytical expressions for the derivatives with respect to physical properties are obtained using the self-adjoint Green's function method. We consider three physical properties of materials; namely dielectrical permittivity, magnetic permeability and electrical conductivity. The inverse problem is solved minimizing the quadratic norm of the residuals using quadratic programming optimization. In the iterative process to speed up convergence we use the Levenberg–Mardquardt method. The special type of linearization is based on an integral equation that involves derivatives of the electric field with respect to magnetic permeability, electrical conductivity and dielectric permittivity; this equation is the result of analyzing the implication of the scaling properties of the electromagnetic field. The ground is modeled using thin horizontal layers to approximate general variations of the physical properties. We show that standard synthetic radargrams due to dielectric permittivity contrasts can be matched using electrical conductivity or magnetic permeability variations. The results indicate that it is impossible to differentiate one property from the other using GPR data.     

Permeability of Triaxially Compressed Sandstone: Influence of Deformation and Strain-rate on Permeability

— The influence of differential stress on the permeability of a Lower Permian sandstone was investigated. Rock cylinders of 50 mm in diameter and 100 mm length of a fine-grained (mean grain size 0.2 mm), low-porosity (6–9%) sandstone were used to study the relation between differential stress, rock deformation, rock failure and hydraulic properties, with a focus on the changes of hydraulic properties in the pre-failure and failure region of triaxial rock deformation. The experiments were conducted at confining pressures up to 20 MPa, and axial force was controlled by lateral strain with a rate ranging from 10^?6 to 10^?7 sec^?1. While deforming the samples, permeability was determined by steady-state technique with a pressure gradient of 1 MPa over the specimen length and a fluid pressure level between 40 and 90% of the confining pressure. The results show that permeability of low-porosity sandstones under increasing triaxial stress firstly decreases due to compaction and starts to increase after the onset of dilatancy. This kind of permeability evolution is similar to that of crystalline rocks. A significant dependence of permeability evolution on strain rate was found. Comparison of permeability to volumetric strain demonstrates that the permeability increase after the onset of dilatancy is not sufficient to regain the initial permeability up to failure of the specimen. The initial permeability, which was determined in advance of the experiments, usually was regained in the post-failure region. After the onset of dilatancy, the permeability increase displays a linear dependence on volumetric strain.     

Modelling the stress-strain behaviour of saturated rocks undergoing triaxial deformation using complex electrical conductivity measurements

Measurement of complex electrical conductivity as a function of frequency is an extremely sensitive probe for changes in pore and crack volume, crack connectivity, and crack surface topography. Such measurements have been made as a function of pore fluid chemistry, hydrostatic confining pressure, as well as uniaxial and triaxial deformation. This paper will; (1) describe the effects of triaxial deformation on the complex electrical conductivity of saturated porous rocks, (2) use the electrical data to model the mechanical stress-strain behaviour, and (3) compare the modelled behaviour with the stress-strain behaviour measured during the deformation. Experimental conductivity data tracks how the rock undergoes compaction with progressive loss of crack volume, followed by dilatation due to new crack formation, growth of existing cracks, crack interlinkage, and finally failure, as axial strain is increased. We have used the complex electrical data to produce a direction-sensitive (anisotropic) crack damage parameter, and used it to calculate the effective Young's modulus by employing the models of Walsh and Bruner. Comparison of the synthetic stress-strain curves so produced, with the experimentally derived stress-strain curves shows good agreement, particularly for undrained tests. This modelling is an improvement on similar curves produced using isotropic crack damage parameters derived from acoustic emission data. The improvement is likely to be due to the directional sensitivity of the electrical conductivity measurement, and its ability to discriminate between the formation of isolated cracks, and those cracks that contribute to the inter-connected crack space i.e. those cracks upon which transport properties of the rock such as electrical conductivity, and mechanical properties depend most critically during triaxial deformation.     

Analysis and interpretation of anomalous conductivity and magnetic permeability effects in time domain electromagnetic data: Part II: Sμ-inversion

In the paper by Pavlov and Zhdanov [J. Appl. Geophys. (2001)], we demonstrated that anomalous magnetic permeability of an ore body could result in measurable anomalous effects in TDEM data, which can be used in geophysical exploration. In the current paper, we develop a new method of TDEM data interpretation, which allows simultaneous reconstruction of both electrical and magnetic properties of the rocks. The method is based on a generalization of the S-inversion technique [66th Ann. Int. Mtg., Soc. Expl. Geophys., Expanded Abstr. (1996) 1306; 68th Ann. Int. Mtg., Soc. Expl. Geophys., Expanded Abstr. (1998) 473; J. Appl. Geophys. (1999)] for models containing thin sheets with anomalous conductivity and anomalous magnetic permeability. We call this method Sμ-inversion. Numerical 3-D modeling results demonstrate that Sμ-inversion provides useful information about both subsurface conductivity and magnetic permeability distributions. It can be used as a new tool for imaging TDEM data in mineral exploration.     

基于中速-高速摩擦实验研究含碳断层带的电导率特征

野外地质调查结果显示,断层带常富集碳质.断层带中碳的分布结构是影响断层带电导率特征的一种重要参数.本文在室温、室内湿度和2MPa正应力条件下,对不同石墨含量(3,5,6和7wt%)的石英-石墨混合断层泥模拟样品开展了滑动速率介于500μm·s-1～1m·s-1的摩擦实验及相应的电导率测量,以期研究断层运动对碳分布结构的影响以及断层带电性特征对碳含量及分布的响应情况.结果显示,摩擦滑动能够显著地改变样品的电性特征(电导率大小及其各向异性).在平行滑动面方向(径向),样品电导率随着滑动位移的增加快速增加,在滑动约数十厘米之后,其电导率基本达到稳定状态;在垂直滑动面方向(轴向),样品电导率基本不随摩擦滑动速率和滑动距离而变化.SEM显微结构观测显示,摩擦滑动所引起的电导率各向异性直接反映了石墨分布结构的变化.该研究结果深化了对地震断裂带浅部电性特征的认识,为野外断层带大地电磁测深资料的解释提供了约束,同时对于了解含碳断层的力学性质和弱矿物相在剪切变形中的分布特征及其演化过程等方面也具有重要意义.     

Electrical conductivity of hydrous silicate melts and aqueous fluids: Measurement and applications

The combination of magnetotelluric survey and laboratory measurements of electrical conductivity is a powerful approach for exploring the conditions of Earth’s deep interior. Electrical conductivity of hydrous silicate melts and aqueous fluids is sensitive to composition, temperature, and pressure, making it useful for understanding partial melting and fluid activity at great depths. This study presents a review on the experimental studies of electrical conductivity of silicate melts and aqueous fluids, and introduces some important applications of experimental results. For silicate melts, electrical conductivity increases with increasing temperature but decreases with pressure. With a similar Na⁺ concentration, along the calc-alkaline series electrical conductivity generally increases from basaltic to rhyolitic melt, accompanied by a decreasing activation enthalpy. Electrical conductivity of silicate melts is strongly enhanced with the incorporation of water due to promoted cation mobility. For aqueous fluids, research is focused on dilute electrolyte solutions. Electrical conductivity typically first increases and then decreases with increasing temperature, and increases with pressure before approaching a plateau value. The dissociation constant of electrolyte can be derived from conductivity data. To develop generally applicable quantitative models of electrical conductivity of melt/fluid addressing the dependences on temperature, pressure, and composition, it requires more electrical conductivity measurements of representative systems to be implemented in an extensive P-T range using up-to-date methods.     

Synthetic radargrams from electrical conductivity and magnetic permeability variations

The objective of this work is to assess the importance of electrical conductivity and magnetic permeability variations in ground penetrating radar (GPR) reflections commonly interpreted only in terms of permittivity variations. We use the matrix propagator approach to obtain the surface electric field associated with a horizontally layered model whose three electromagnetic properties vary from layer to layer. The solution is based on the plane wave boundary value problem using inverse Fourier transformation to accommodate particular GPR pulses. Our results indicate that while magnetic permeability is unimportant, reflections from electrical conductivity variations can be of the same order as those associated with electrical permittivity boundaries. In particular, we show that a realistic ground model composed of thin conductive layers can produce radargrams similar to those caused by a lossless permittivity contrast.     

Origin of High Electrical Conductivity in the Lower Continental Crust: A Review

Electromagnetic measurements have demonstrated that the lower continental crust has remarkable electrical anomalies of high conductivity and electrical anisotropy on a global scale (probably with some local exceptions), but their origin is a long-standing and controversial problem. Typical electrical properties of the lower continental crust include: (1) the electrical conductivity is usually 10⁻⁴ to 10⁻¹ S/m; (2) the overlying shallow crust and underlying upper mantle are in most cases less conductive; (3) the electrical conductivity is statistically much higher in Phanerozoic than in Precambrian areas; (4) horizontal anisotropy has been resolved in many areas; and (5) in some regions there appear to be correlations between high electrical conductivity and other physical properties such as seismic reflections. The explanation based on conduction by interconnected, highly conductive phases such as fluids, melts, or graphite films in grain boundary zones has various problems in accounting for geophysically resolved electrical conductivity and other chemical and physical properties of the lower crust. The lower continental crust is dominated by mafic granulites (in particular beneath stable regions), with nominally anhydrous clinopyroxene, orthopyroxene, and plagioclase as the main assemblages, and the prevailing temperatures are mostly 700–1,000°C as estimated from xenolith data, surface heat flow, and seismic imaging. Pyroxenes have significantly higher Fe content in the lower crust than in the upper mantle (peridotites), and plagioclase has higher Na content in the lower crust than in the shallow crust (granites). Minerals in the lower continental crust generally contain trace amounts of water as H-related point defects, from less than 100 to more than 1,000 ppm H₂O (by weight), with concentrations usually higher than those in the upper mantle. Observations of xenolith granulites captured by volcano-related eruptions indicate that the lower continental crust is characterized by alternating pyroxene-rich and plagioclase-rich layers. Experimental studies on typical lower crustal minerals have shown that their electrical conductivity can be significantly enhanced by the higher contents of Fe (for pyroxenes), Na (for plagioclase), and water (for all minerals) at thermodynamic conditions corresponding to the lower continental crust, e.g., to levels comparable to those measured by geophysical field surveys. Preferred orientation of hydrous plagioclase, e.g., due to ductile flow in the deep crust, and alternating mineral fabrics of pyroxene-rich and plagioclase-rich layers can lead to substantial anisotropy of electrical conductivity. Electrical conductivity properties in many regions of the lower continental crust, especially beneath stable areas, can mostly be accounted for by solid-state conduction due to the major constituents; other special, additional conduction mechanisms due to grain boundary phases are not strictly necessary.     

含碳结构对龙门山断层带电导率影响的实验探索

碳是影响岩石电导率大小的一个重要因素,可能是造成龙门山断层带电导率异常的重要原因之一.为了研究不同的碳含量、矿物颗粒粒径与碳晶体结构对断层带电导率的影响,在干燥、常温、0.2~300 MPa的压力条件下实验研究了人工模拟断层泥样品(石英粉末与含碳粉末混合的样品,简称模拟样品)和采自映秀-北川断层八角庙剖面的天然断层岩样品(简称天然样品)的电导率.实验结果显示,当模拟样品中的含碳粉末连通时,电导率与碳体积率的关系符合逾渗理论模型;而含碳粉末未连通时,电导率随总孔隙度降低而指数性升高.同时模拟样品的电导率也随石英颗粒粒径的变化而发生改变.相比于模拟样品中的含碳粉末主要分布于石英颗粒支撑的孔隙中,天然样品中的碳则主要以碳膜的形式赋存在颗粒边缘,导致碳体积率相同的条件下,模拟样品的电导率小于天然样品.此外,天然样品的电导率(9×10~(-4)S·m~(-1))也要小于野外大地电磁探测的结果(0.03~0.1 S·m~(-1)).在今后的实验中还需要考虑在动态摩擦条件下对含有完整含碳结构的天然样品进行电导率的实验研究.     

Evaluation of transport and storage properties in the soil and groundwater zone from induced polarization measurements1

Spectral induced polarization as well as complex electrical measurements are used to estimate, on a non-invasive basis, hydraulic permeability in aquifers. Basic laboratory measurements on a variety of shaly sands, silts and clays showed that the main feature of their conductivity spectra in the frequency range from 10^-3 to 10³ Hertz is a nearly constant phase angle. Thus, a constant-phase-angle model of electrical conductivity is applied to interpret quantitatively surface and borehole spectral induced polarization measurements. The model allows for the calculation of two independent electrical parameters from only one frequency scan and a simple separation of electrical volume and interface effects. The proposed interpretation algorithm yields the true formation factor, the cation exchange capacity and the surface-area-to-porosity ratio, which corresponds to the inverse hydraulic radius. Using a Kozeny–Carman-like equation, the estimation of hydraulic permeability is possible.     

Multifractal characterization of airborne geophysical data at the Oak Ridge facility

Scaling analyses on geophysical measurements of electrical conductivity, gamma radiation, and magnetic fields, at the Oak Ridge Reservation were conducted. The electrical conductivity and magnetic data exhibited multifractality in the north-south and east-west directions. The radiation data were observed to be non-scaling; a variogram with a sill was found to be more appropriate. The scaling of the EC and magnetic was generally within a range smaller than the maximum distance selected, as periodicity dominated at the larger distances. The electrical conductivity had anisotropy in the scaling of their variograms. But the magnetic data appear to have an isotropic scaling. The underlying statistics of the fields were near Gaussian for the electrical conductivity, but essentially Gaussian for the magnetic data. In environmental hydrogeology, knowledge of the spatial distribution of the intrinsic permeability, K, is very helpful in understanding the transport and spreading of contaminant plumes. Our previous studies have shown that the subsurface permeability, K, is multifractal. Detailed measurement of K is costly. Hence, large data sets of value collected both on a fine scale and over large distances are rare. In this study, we hypothesize that geophysical data could be used indirectly as a surrogate measurement for K, for obtaining statistical information on scale limited K data, and perhaps, directly at sites where K and electrical conductivity are correlated.     

储层渗透性与地层因素关系的实验研究与分析

本文对渤海湾盆地不同孔隙结构样品的孔、渗、核磁、岩电、压汞、X衍射及铸体薄片等配套岩石物理实验数据进行了综合分析,通过逐一考察同一套岩芯样品的地层因素与渗透率、压汞喉径均值、储层品质指数之间的实验关系,并分别与地层因素-孔隙度交绘图进行对比分析,发现储层渗透性与地层因素之间并非简单的单调函数关系,孔隙度相近但孔隙结构类型不同、渗透率差异明显的岩芯可以具有相近的地层因素,导电能力接近.在实验数据分析的基础上通过理论分析证明了这一实验关系的合理性,并指出孔隙度及导电能力相近的岩芯,其渗透率差异与喉径均值的平方比、孔隙曲折度及几何形态相关.     

Transport properties and microstructural characteristics of a thermally cracked mylonite

An experimental study was carried out on a granitic mylonite (La Bresse, France) to analyze the influence of pore microstructure on transport properties. Different crack networks were obtained by a controlled thermal treatment. Microstructures were analyzed by means of gas adsorption and mercury porosimetry. Transport properties have been investigated by measuring gas permeability and electrical conductivity. The dependence of permeability on confining pressure shows an exponential decrease, characteristic of a porosity made of cracks. Correlations between measured parameters have been analyzed by comparing them with relations deduced from theoretical models. Linking the formation factor to the porosity leads to a rather low tortuosity value (about 2.4), characterizing a medium with a well connected porosity. Correlation between permeabilityk and formation factorF leads to a power-law relationk F ^–n wheren2.9, which is consistent with a crack model describing the behavior of the thermally treated rock.     

Laboratory Electrical Conductivity Measurement of Mantle Minerals

Electrical conductivity structures of the Earth’s mantle estimated from the magnetotelluric and geomagnetic deep sounding methods generally show increase of conductivity from 10⁻⁴–10⁻² to 10⁰ S/m with increasing depth to the top of the lower mantle. Although conductivity does not vary significantly in the lower mantle, the possible existence of a highly conductive layer has been proposed at the base of the lower mantle from geophysical modeling. The electrical properties of mantle rocks are controlled by thermodynamic parameters such as pressure, temperature and chemistry of the main constituent minerals. Laboratory electrical conductivity measurements of mantle minerals have been conducted under high pressure and high temperature conditions using solid medium high-pressure apparatus. To distinguish several charge transport mechanisms in mantle minerals, it is necessary to measure the electrical conductivity in a wider temperature range. Although the correspondence of data has not been yet established between each laboratory, an outline tendency of electrical conductivity of the mantle minerals is almost the same. Most of mineral phases forming the Earth’s mantle exhibit semiconductive behavior. Dominant conduction mechanism is small polaron conduction (electron hole hopping between ferrous and ferric iron), if these minerals contain iron. The phase transition olivine to high-pressure phases enhances the conductivity due to structural changes. As a result, electrical conductivity increases in order of olivine, wadsleyite and ringwoodite along the adiabat geotherm. The phase transition to post-spinel at the 660 km discontinuity further can enhance the conductivity. In the lower mantle, the conductivity once might decrease in the middle of the lower mantle due to the iron spin transition and then abruptly increase at the condition of the D″ layer. The impurities in the mantle minerals strongly control the formation, number and mobility of charge carriers. Hydrogen in nominally anhydrous minerals such as olivine and high-pressure polymorphs can enhance the conductivity by the proton conduction. However, proton conduction has lower activation enthalpy compared with small polaron conduction, a contribution of proton conduction becomes smaller at high temperatures, corresponding to the mantle condition. Rather high iron content in mantle minerals largely enhances the conductivity of the mantle. This review focuses on a compilation of fairly new advances in experimental laboratory work together with their explanation.     

Temperature and Pressure Dependencies of Thermal Transport Properties of Rocks: Implications for Uncertainties in Thermal Lithosphere Models and new Laboratory Measurements of High-Grade Rocks in the Central Fennoscandian Shield

Measurements on thermal conductivity and diffusivity as functions of temperature (up to 1150 K) and pressure (up to 1000 MPa) are presented for Archaean and Proterozoic mafic high-grade rocks metamorphosed in middle and lower crustal pressures, and situated in eastern Finland, central Fennoscandian Shield. Decrease of 12–20% in conductivity and 40–55% in diffusivity was recorded between room temperature and 1150 K, which can be considered as typical of phonon conductivity. Radiative heat transfer effects were not detected in these samples. Pressure dependencies of the samples are weak if compared to crystalline rocks in general, but relatively typical for mafic rocks.The temperature and pressure dependencies of thermal transport properties (data from literature and the present study) were applied in an uncertainty analysis of lithospheric conductive thermal modellings with random (Monte Carlo) simulations using a 4-layer model representative of shield lithosphere. Model parameters were varied according to predetermined probability functions and standard deviations were calculated for lithospheric temperature and heat flow density after 1500 independent simulations. The results suggest that the variations (uncertainties) in calculated temperature and heat flow density values due to variations in the temperature and pressure dependencies of conductivity are minor in comparison to the effects produced by typical variations in the room temperature value of conductivity, heat production rate or lower boundary condition values.