Advective Heat Transport Research Articles

Relative importance of advective heat transport and boundary layer decoupling in the melt dynamics of a patchy snow cover

For a patchy snow cover the advective heat transport and the near-surface boundary layer decoupling, which have an opposite effect on sensible heat transport onto the snow surface, are both expected to increase in magnitude. The main aim of this study is to investigate the effects of locally developing atmospheric stratification over a discontinuous snow cover which can result in a decoupling from the warm atmosphere. We are particularly interested in the effect of boundary layer decoupling on the net sensible heat flux into the snow. We therefore applied local eddy flux measurements over snow patches at three different heights above the snow surface. We identified wind velocity, turbulence intensity, fetch distance and topographical curvature as the main factors driving the boundary layer depth and the efficiency of advective heat transport to contribute to snow ablation. The atmospheric decoupling is thus shown to be a key mechanism in snow patch survival.

Micrometeorological processes driving snow ablation in an Alpine catchment

Abstract. Mountain snow covers typically become patchy over the course of a melting season. The snow pattern during melt is mainly governed by the end of winter snow depth distribution and the local energy balance. The objective of this study is to investigate micro-meteorological processes driving snow ablation in an Alpine catchment. For this purpose we combine a meteorological boundary-layer model (Advanced Regional Prediction System) with a fully distributed energy balance model (Alpine3D). Turbulent fluxes above melting snow are further investigated by using data from eddy-correlation systems. We compare modeled snow ablation to measured ablation rates as obtained from a series of Terrestrial Laser Scanning campaigns covering a complete ablation season. The measured ablation rates indicate that the advection of sensible heat causes locally increased ablation rates at the upwind edges of the snow patches. The effect, however, appears to be active over rather short distances of about 4–6 m. Measurements suggest that mean wind velocities of about 5 m s−1 are required for advective heat transport to increase snow ablation over a long fetch distance of about 20 m. Neglecting this effect, the model is able to capture the mean ablation rates for early ablation periods but strongly overestimates snow ablation once the fraction of snow coverage is below a critical value of approximately 0.6. While radiation dominates snow ablation early in the season, the turbulent flux contribution becomes important late in the season. Simulation results indicate that the air temperatures appear to overestimate the local air temperature above snow patches once the snow coverage is low. Measured turbulent fluxes support these findings by suggesting a stable internal boundary layer close to the snow surface causing a strong decrease of the sensible heat flux towards the snow cover. Thus, the existence of a stable internal boundary layer above a patchy snow cover exerts a dominant control on the timing and magnitude of snow ablation for patchy snow covers.

Advective Heat Transport in Frozen Rock Clefts: Conceptual Model, Laboratory Experiments and Numerical Simulation

ABSTRACTAdvective heat transported by water percolating into discontinuities in frozen ground can rapidly increase temperatures at depth because it provides a thermal shortcut between the atmosphere and the subsurface. Here, we develop a conceptual model that incorporates the main heat‐exchange processes in a rock cleft. Laboratory experiments and numerical simulations based on the model indicate that latent heat release due to initial ice aggradation can rapidly warm cold bedrock and precondition it for later thermal erosion of cleft ice by advected sensible heat. The timing and duration of water percolation both affect the ice‐level change if initial aggradation and subsequent erosion are of the same order of magnitude. The surplus advected heat is absorbed by cleft ice loss and runoff from the cleft so that this energy is not directly detectable in ground temperature records. Our findings suggest that thawing‐related rockfall is possible even in cold permafrost if meltwater production and flow characteristics change significantly. Advective warming could rapidly affect failure planes beneath large rock masses and failure events could therefore differ greatly from common magnitude reaction‐time relations. Copyright © 2011 John Wiley & Sons, Ltd.

The role of advective heat transport in talik development beneath lakes and ponds in discontinuous permafrost

[1] Regions of warm, thin, discontinuous permafrost have been observed to be experiencing rapid changes in lake and pond dynamics in recent decades. Even though surface water and groundwater interactions are thought to play a significant role in heat transport in these regions, the effect of these interactions on permafrost remains largely unquantified. In order to examine the influence of groundwater flow on permafrost dynamics, we modeled the development of a sub-lake talik under permafrost conditions similar to those observed in the southern-central Seward Peninsula region of Alaska using a numerical solution that couples heat transport and groundwater flow, including the effect of water phase changes on soil permeability and latent heat content. A comparison of model simulations, with and without near surface subpermafrost groundwater flow, indicates that stable permafrost thicknesses are 2 to 5 times greater in the absence of groundwater flow. Simulations examining the thermal influence of lakes on underlying permafrost suggest that a through-going talik can develop in a matter of decades and that the incorporation of advective heat transport reduces the time to complete loss of ice beneath the lake by half, relative to heat transport by conduction alone. This work presents the first quantitative assessment of the rates of sub-lake permafrost response to thermal disturbances, such as talik development, in systems with near-surface groundwater flow. The results highlight the importance of coupled thermal and hydrologic processes on discontinuous permafrost dynamics.

Modeling heat fluxes in macroporous soil under sparse young forest of temperate humid climate

This research is focused on improved quantitative understanding of processes controlling the movement of water and heat below the land surface under natural vegetation. Special attention is paid to the advective transport of soil heat during rainfall events and the role of preferential flow in the soil heat transfer system. Transient flow of water and transport of heat were simulated using dual-continuum model with semi-separate flow domains representing the soil matrix and the network of preferential pathways. The model was able to reproduce quite well the measured variations of soil temperatures, and, at least in principle, also the measured temperatures of shallow subsurface discharge, closely associated with the preferential flow.

Mechanisms of mass and heat transport during Barrovian metamorphism: A discussion based on field evidence from the Central Alps (Switzerland/northern Italy)

Tectonic and metamorphic data for the Central Alps (Switzerland/Italy) are used to discuss this classic example of a Barrovian metamorphic terrain, notably the evolution of its thermal structure in space and time. Available P‐T‐t data indicate variable contributions of advective and conductive heat transport during collision and subsequent cooling and exhumation. Some areas experienced a prolonged period of partial melting while other areas, at the same time, show but moderate heating. The Barrow‐type metamorphic field gradient observed in the final orogen is the result of two distinct tectonic processes, with their related advective and conductive heat transport processes. The two tectonic processes are (1) accretion of material within a subduction channel related to decompression and emplacement of high‐pressure units in the middle crust and (2) wedging and related nappe formation in the continental lower plate. The second process postdates the first one. Wedging and underthrusting of continental lower plate material produces heat input into lower crustal levels, and this process is responsible for predominantly conductive heat transport in the overlying units. The interacting processes lead to different maximum temperatures at different times, producing the final Barrovian metamorphic field gradient. The south experienced rapid cooling, whereas the north shows moderate cooling rates. This discrepancy principally reflects differences in the temperature distribution in the deeper crust prior to cooling. Differences in the local thermal gradient that prevailed before the cooling also determined the relationships between cooling rate and exhumation rate in the different areas.

The thermal regime in the resurgent dome of Long Valley Caldera, California: Inferences from precision temperature logs in deep wells

Long Valley Caldera in eastern California formed 0.76 Ma ago in a cataclysmic eruption that resulted in the deposition of 600 km 3 of Bishop Tuff. The total current heat flow from the caldera floor is estimated to be ~ 290 MW, and a geothermal power plant in Casa Diablo on the flanks of the resurgent dome (RD) generates ~40 MWe. The RD in the center of the caldera was uplifted by ~ 80 cm between 1980 and 1999 and was explained by most models as a response to magma intrusion into the shallow crust. This unrest has led to extensive research on geothermal resources and volcanic hazards in the caldera. Here we present results from precise, high-resolution, temperature–depth profiles in five deep boreholes (327–1,158 m) on the RD to assess its thermal state, and more specifically 1) to provide bounds on the advective heat transport as a guide for future geothermal exploration, 2) to provide constraints on the occurrence of magma at shallow crustal depths, and 3) to provide a baseline for future transient thermal phenomena in response to large earthquakes, volcanic activity, or geothermal production. The temperature profiles display substantial non-linearity within each profile and variability between the different profiles. All profiles display significant temperature reversals with depth and temperature gradients <50 °C/km at their bottom. The maximum temperature in the individual boreholes ranges between 124.7 °C and 129.5 °C and bottom hole temperatures range between 99.4 °C and 129.5 °C. The high-temperature units in the three Fumarole Valley boreholes are at the approximate same elevation as the high-temperature unit in borehole M-1 in Casa Diablo indicating lateral or sub-lateral hydrothermal flow through the resurgent dome. Small differences in temperature between measurements in consecutive years in three of the wells suggest slow cooling of the shallow hydrothermal flow system. By matching theoretical curves to segments of the measured temperature profiles, we calculate horizontal groundwater velocities in the hydrothermal flow unit under the RD that range from 1.9 to 2.8 m/yr, which corresponds to a maximum power flowing through the RD of 3–4 MW. The relatively low temperatures and large isothermal segments at the bottom of the temperature profiles are inconsistent with the presence of magma at shallow crustal levels.

Advective Heat Transport in an Unconfined Aquifer Induced by the Field Injection of an Open-Loop Groundwater Heat Pump

Problem statement: The increasing diffusion of low-enthalpy geothermal open-loop Groundwater Heat Pumps (GWHP) providing buildings air conditioning requires a careful assessment of the overall effects on groundwater system, especially in the urban areas. The impact on the groundwater temperature in the surrounding area of the re-injection well is directly linked to the aquifer properties. Physical processes affecting heat transport within an aquifer include advection (or convection) and hydrodynamic thermodispersion (diffusion and mechanical dispersion). If the groundwater flows, the advective components tend to dominate the heat transfer process within the aquifer and the diffusion can be considered negligible. This study illustrates the experimental results derived from the groundwater monitoring in the surrounding area of an injection well connected to an open-loop GWHP plant which has been installed in the "Politecnico di Torino" (NW Italy) for cooling some of the university buildings. Groundwater pumping and injection interfere only with the upper unconfined aquifer. Approach: After the description of the hydrogeological setting the authors examined the data deriving from multiparameter probes installed inside the pumping well (P2), the injection well (P4) and a downgradient piezometer (S2). Data refers to the summer 2009. To control the aquifer thermal stratification some multi-temporal temperature logs have been performed in the S2. Results: After the injection of warm water in P4 the plume arrived after 30 days in the S2. That delay is compatible with the calculated plume migration velocity (1.27 m d-1) and their respective distance (35 m). The natural temperature in the aquifer due to the switching-off of the GWHP plant has been reached after two month. The Electrical Conductivity (EC) values tend to vary out of phase with the temperature. The temperature logs in the S2 highlighted a thermal stratification in the aquifer due to a low vertical dispersion of the injected warm water. Conclusion: Experimental evidences seem to confirm the prevalence of heat advective transport component respect the dispersive phenomena. This hypothesis appears validated by the following evidences: (i) the calculated advective migration velocities are compatible with the calculated retardation factor and the temperature revealed in the S2, (ii) both the groundwater and the heat tend to flow horizontally due to the different values of horizontal and vertical hydraulic conductivity in the Unit 1 (thermal stratification) and (iii) the flowing water highlighted different geochemical characteristics during the time.

Interannual variability in North Pacific heat and freshwater budgets

Abstract Transports of volume, heat, and freshwater in the North Pacific Ocean from 1992 to 2004 are analyzed using a long-term high resolution expendable bathythermography (XBT) dataset and output from a data-assimilating model. Estimates of geostrophic transport from the data are compared with the model transport to close the volume budget north of the trans-Pacific XBT track. Advective transport from both model output and data are combined with surface fluxes to determine budgets of heat and freshwater in the closed region. The northward heat transport across the XBT track is estimated to be 0.74 ± 0.1 pW , and has variability of almost 0.5 pW on 3–4 year time scales, while freshwater transport is estimated to be − 0.1 ± 0.06 Sv . The balance between northward advective heat transport and surface heat flux gives a time-varying estimate of heat storage that compares well with observations. A similar balance is found between model estimates of advective freshwater transport and surface freshwater flux. Despite a scarcity of observations and uncertainties in all components, this analysis results in nearly closed budgets of volume, heat, and freshwater. Mean estimates of advective transport of both heat and freshwater agree with previous estimates. An analysis of each component of the heat budget with latitude indicates that a relative lack of time-variability of the surface component is consistent throughout the North Pacific. The dominant advective component is driven by changes in the wind stress curl field. For both heat and freshwater storage, strong signals occur concurrently at all latitudes. This behavior could indicate that these signals are controlled by large-scale dynamics, rather than small-scale disturbances from which signals would need to propagate to be widely felt. The analysis demonstrates the value of bringing models and data together, resulting in budgets that are consistent with observations, yet provide a comprehensive look at the variability of North Pacific heat and freshwater storage that would be unavailable from data alone.

Alpine metamorphism in the central segment of the Western Greywacke Zone (Eastern Alps)

Alpine metamorphism in the central segment of the Western Greywacke Zone (Eastern Alps)The metamorphic pattern of the central Western Greywacke Zone (Austroalpine, Eastern Alps) was investigated by organic matter reflectance, Raman spectroscopy on organic matter and clay mineralogical methods. Raman data map a 10 km wide thermal aureole along the contact zone of the Greywacke Zone to the Penninic Tauern Window. The estimated maximum temperatures of 400 °C to 200 °C decrease from South to North, that is from the contact to the uppermost parts of the Greywacke Zone. This pattern is explained by an Oligocene to Miocene thermal pulse, related to the rapid exhumation of formerly deeply buried rocks of the Penninic unit. During this event, advective heat transport and circulating fluids overprinted the Cretaceous higher anchi- to lower epizonal metamorphic pattern of the central Western Greywacke Zone.

Spring temperatures in the Sagehen Basin, Sierra Nevada, CA: implications for heat flow and groundwater circulation

AbstractHeat flow in the Sierra Nevada, CA, is low despite its young geologic age. We investigate the possibility that advective heat transport by groundwater flow leads to an underestimate of heat flow in the Sierras based purely on borehole measurements. Using temperature and discharge measurements at springs in Sagehen Basin, we find that groundwater removes the equivalent of approximately 20–40 mW m−2 of geothermal heat from the basin. This is comparable with other heat flow measurements in the region and indicates that, in this basin, at least, groundwater does transport a significant amount of geothermal heat within the basin. Additionally, we use estimates of the mean residence time of water discharged at the springs along with hourly temperature records in springs to provide constraints on groundwater flow depths within the basin. An analytical model based on these constraints indicates that the heat removed by groundwater may represent 20% to &gt;90% of the total heat flow in the basin. Without better constraints on the regional hydrogeology and the depth of circulation, we cannot determine whether the heat discharged at the springs represents a change in the mode of heat transfer, i.e. from conduction to advection at shallow depths (&lt;100 m) or whether this is a component of heat transfer that should be added to measured conductive values. If the latter is true, and Sagehen Basin is representative of the Sierras, basal heat flow in the Sierra Nevada may be higher than previously thought.

Enhanced convection and fast plumes in the lower mantle induced by the spin transition in ferropericlase

Using a numerical model we explore the consequences of the intrinsic density change (Δρ/ρ ≈ 2–4%) caused by the Fe2+ spin transition in ferropericlase on the style and vigor of mantle convection. The effective Clapeyron slope of the transition from high to low spin is strongly positive in pressure‐temperature space and broadens with high temperature. This introduces a net spin‐state driving density difference for both upwellings and downwellings. In 2‐D cylindrical geometry spin‐buoyancy dominantly enhances the positive thermal buoyancy of plumes. Although the additional buoyancy does not fundamentally alter large‐scale dynamics, the Nusselt number increases by 5–10%, and vertical velocities by 10–40% in the lower mantle. Advective heat transport is more effective and temperatures in the core‐mantle boundary region are reduced by up to 12%. Our findings are relevant to the stability of lowermost mantle structures.

Thermal influence of urban groundwater recharge from stormwater infiltration basins

AbstractGroundwater warming below cities has become a major environmental issue; but the effect of distinct local anthropogenic sources of heat on urban groundwater temperature distributions is still poorly documented. Our study addressed the local effect of stormwater infiltration on the thermal regime of urban groundwater by examining differences in water temperature beneath stormwater infiltration basins (SIB) and reference sites fed exclusively by direct infiltration of rainwater at the land surface. Stormwater infiltration dramatically increased the thermal amplitude of groundwater at event and season scales. Temperature variation at the scale of rainfall events reached 3 °C and was controlled by the interaction between runoff amount and difference in temperature between stormwater and groundwater. The annual amplitude of groundwater temperature was on average nine times higher below SIB (range: 0·9–8·6 °C) than at reference sites (range: 0–1·2 °C) and increased with catchment area of SIB. Elevated summer temperature of infiltrating stormwater (up to 21 °C) decreased oxygen solubility and stimulated microbial respiration in the soil and vadose zone, thereby lowering dissolved oxygen (DO) concentration in groundwater. The net effect of infiltration on average groundwater temperature depended upon the seasonal distribution of rainfall: groundwater below large SIB warmed up (+0·4 °C) when rainfall occurred predominantly during warm seasons. The thermal effect of stormwater infiltration strongly attenuated with increasing depth below the groundwater table indicating advective heat transport was restricted to the uppermost layers of groundwater. Moreover, excessive groundwater temperature variation at event and season scales can be attenuated by reducing the size of catchment areas drained by SIB and by promoting source control drainage systems. Copyright © 2009 John Wiley &amp; Sons, Ltd.

Monitoring of CO 2 plumes during storage in geological formations using temperature signals: Numerical investigation

Carbon dioxide ( CO 2 ) injection into a storage formation is accompanied by non-isothermal effects. These are caused by a CO 2 injection temperature that does not correspond to the formation temperature, cooling of the carbon dioxide due to expansion (Joule–Thomson cooling) and heat of dissolution of CO 2 in brine. During flow in the subsurface, the carbon dioxide transports energy (advective heat transport) and undergoes an equilibrating process between temperature differences (heat conduction). These non-isothermal processes can be used for the purpose of monitoring the CO 2 plume propagation in the subsurface. Temperature sensors at monitoring wells at a certain distance from the injection well can detect temperature changes and give information about the CO 2 flow in the storage site. In this study, a numerical multi-phase simulation program is used to investigate the non-isothermal effects during CO 2 injection into a storage formation. The feasibility of using temperature measurements for the observation of the carbon dioxide plume in the reservoir is addressed. Various thermal processes and their dependency on the geological characterisation of the reservoir are discussed in detail.

Modeling the Merensky Reef, Bushveld Complex, Republic of South Africa

The Merensky pegmatoid (normal reef) in the western Bushveld Complex is commonly characterized as a pyroxene-rich pegmatoidal unit with a base that is enriched in chromite and platinum-group element-bearing sulfides overlying a leuconorite footwall. Models for its formation have ranged from those that view it as entirely a magmatic cumulate succession to those that have suggested that it is a zone of volatile-induced remelting. The consequences of the latter interpretation are investigated using the numerical modeling program IRIDIUM, which links diffusive and advective mass and heat transport with a phase equilibration routine based on the MELTS program. The initial system consists of a simple stratigraphic succession of a partially molten leuconorite overlain by a partially molten pyroxenite, both initially at 1,190°C and 2 kbar. 2 wt% of a volatile fluid composed of 75 mol% H2O, 20 mol% CO2 and 5 mol% H2S is then added to the lower 20 cm of the pyroxenite. The system is then allowed to evolve under conditions of chemical diffusion in the liquid. The addition of the volatile components results in a modest increase in the amount of melt in the pyroxenite. However, chemical diffusion across the leuconorite–pyroxenite boundary leads to more extensive melting at and below the boundary with preferential loss of opx from the underlying leuconorite, preferential re-precipitation of sulfide and chromite and concentration of the PGE at this boundary. These results mimic actual mineral and compositional profiles across the Merensky pegmatoid and illustrate that long-term diffusion process can effectively produce mineralogical and compositional layering not present in the original assemblage.

Predictive targeting in Australian orogenic-gold systems at the deposit to district scale using numerical modelling

3D numerical models of coupled deformation and fluid flow provide a useful tool for exploration in orogenic-gold systems. Numerical modelling of ore-forming processes can lead to a reduction in targeting and detection risk, thus improving the value proposition of mineral exploration. Hydrothermal mineralisation arises from a complex interplay of deformation, fluid flow, conductive and advective heat transport, solute transport and chemical reactions. Coupled simulation of all of these processes represents a significant computational challenge that cannot be solved within the time-scale of a mineral exploration program. However, the problem can be simplified by identifying a subset of processes representing the first-order controls on mineralisation at the scale of interest. For most orogenic-gold systems, it is argued that the first-order controls on mineralisation at the camp to deposit scale are deformation-induced dilation, fluid flow and fluid focusing. Hence, numerical models of coupled deformation and fluid flow can provide a quantitative insight into the localisation of ore-forming fluids in this type of system. In two case studies, known deposits were modelled in order to determine the critical deformation and fluid-flow-related factors controlling the localisation of mineralisation in these systems. The quantitative results from the forward models were then used as a basis for constructing predictive models that were applied to regional targeting, prospect ranking and selecting the choice of detection methods. Both case studies show that numerical modelling is capable of reproducing the distribution of known anomalism, and that it can predict anomalies that were not expected or accounted for by purely empirical analysis.

Hydraulic properties of the crystalline basement

Hydraulic tests in boreholes, up to 4.5 km deep, drilled into continental crystalline basement revealed hydraulic conductivity (K) values that range over nine log-units from 10−13−10−4 m s−1. However, K values for fractured basement to about 1 km depth are typically restricted to the range from 10−8 to 10−6 m s−1. New data from an extended injection test at the KTB research site (part of the Continental Deep Drilling Program in Germany) at 4 km depth provide K=5 10−8 m s−1. The summarized K-data show a very strong dependence on lithology and on the local deformation history of a particular area. In highly fractured regions, granite tends to be more pervious than gneiss. The fracture porosity is generally saturated with Na–Cl or Ca–Na–Cl type waters with salinities ranging from 100 g L−1. The basement permeability is well within the conditions for advective fluid and heat transport. Consequently, fluid pressure is hydrostatic and a Darcy flow mechanism is possible to a great depth. Topography-related hydraulic gradients in moderately conductive basement may result in characteristic advective flow rates of up to 100 L a−1 m−2 and lead to significant advective heat and solute transfer in the upper brittle crust.

Is advective heat transport significant at the Dead Sea basin?

AbstractAn understanding of heat flux is a necessary component in reconstructing tectonic, seismic, and hydrologic models of the Dead Sea basin. Heat may be transferred by both conduction and advection by groundwater. Although the conductive heat flux in Israel has been extensively measured to be approximately 40 mW m−2, there is still a debate about the total heat flux. Recently, the discharge of hot springs along the western Dead Sea shore has been determined to be 107 m3 year−1. Simple calculations show that the heat discharged by groundwater at these hot springs is of the same order of magnitude as the measured conductive heat flux in deep boreholes. Therefore the total heat flux could be significantly higher than 40 mW m−2. However, results of numerical modeling show that the current hot‐spring heat discharge is two orders of magnitude greater than that predicted for steady‐state conditions and can be explained by the rapid recession of the Dead Sea.

Mechanisms Governing Interannual Variability of Upper-Ocean Temperature in a Global Ocean Hindcast Simulation

Abstract The interannual variability in upper-ocean (0–400 m) temperature and governing mechanisms for the period 1968–97 are quantified from a global ocean hindcast simulation driven by atmospheric reanalysis and satellite data products. The unconstrained simulation exhibits considerable skill in replicating the observed interannual variability in vertically integrated heat content estimated from hydrographic data and monthly satellite sea surface temperature and sea surface height data. Globally, the most significant interannual variability modes arise from El Niño–Southern Oscillation and the Indian Ocean zonal mode, with substantial extension beyond the Tropics into the midlatitudes. In the well-stratified Tropics and subtropics, net annual heat storage variability is driven predominately by the convergence of the advective heat transport, mostly reflecting velocity anomalies times the mean temperature field. Vertical velocity variability is caused by remote wind forcing, and subsurface temperature anomalies are governed mostly by isopycnal displacements (heave). The dynamics at mid- to high latitudes are qualitatively different and vary regionally. Interannual temperature variability is more coherent with depth because of deep winter mixing and variations in western boundary currents and the Antarctic Circumpolar Current that span the upper thermocline. Net annual heat storage variability is forced by a mixture of local air–sea heat fluxes and the convergence of the advective heat transport, the latter resulting from both velocity and temperature anomalies. Also, density-compensated temperature changes on isopycnal surfaces (spice) are quantitatively significant.

Composition and thermal evolution of the lithospheric mantle beneath the Harrat Uwayrid, eastern flank of the Red Sea rift (Saudi Arabia)

Harrat Uwayrid is one of several volcanic fields that developed along the western margin of the Arabian plate during opening of the Red Sea rift. Volcanic activity started in the late Miocene and Pliocene with the eruption of alkaline flood basalts. After a period of extensional faulting, numerous cinder cones formed in the central part of the volcanic field during the Quaternary. Xenoliths of both mantle and crustal origin are abundant. The mantle-derived suite is dominated by spinel lherzolite with X Mg of olivine and X Cr of spinel ranging from 0.868 to 0.911 and from 0.074 to 0.636, respectively. In addition to lherzolite, various types of pyroxenites and wehrlites occur. Most xenoliths contain primary and/or secondary amphibole (pargasite/kaersutite), indicating at least two enrichment events, both related to the formation of the Red Sea. Chemical zonation profiles of mineral grains constrain the thermal evolution of the xenoliths and their source regions in the shallow upper mantle. Ca contents in olivine are constant in the inner parts of the grains, but increase in the outermost zones (< 35 μm). X Mg values of spinel grains adjacent to olivine often show a rimward increase in their outer zones (< 150 μm). Both types of zonation patterns are interpreted to result from heating by the host magma during transport. Zonation patterns of Al and Ca/(1 − Na) in pyroxene grains record different thermal evolutions, such as (i) simple cooling, whereby a ‘near-equilibrium’ state was rarely achieved, (ii) cooling followed by heating, and (iii) simple heating. Diffusion modelling on measured zoning profiles suggests that cooling lasted for no more than ∼ 10 Myr. Temperature estimates for pyroxene core and rim compositions, calculated by using various pyroxene geothermometers range from ∼ 800 to 1140 °C (at an assumed pressure of 1.5 GPa). In conjunction with upper pressure limits derived from the composition of spinel, these values imply relatively high temperatures for the shallow upper mantle. This is corroborated by P– T estimates for spinel lherzolites obtained by combination of the two-pyroxene thermometer based on the solvus between enstatite and diopside and the Ca-in-olivine geothermobarometer. While the high- T state that preceded cooling was probably a regional phenomenon, renewed heating was of local nature, most probably related to advective heat transport by ascending and stagnating magmas.