Waveform tomography with very large datasets reveals the upper-mantle structure of the Arctic in unprecedented detail. Using tomography jointly with computational petrology, we estimate temperature in the lithosphere–asthenosphere depth range and infer lithospheric structure and evolution. Most of the boundaries of the mantle roots of cratons in the Arctic are coincident with their geological boundaries at the surface. The thick lithospheres of the Greenland and North American cratons are separated by a corridor of thin lithosphere beneath Baffin Bay and through the middle of the Canadian Arctic Archipelago; the southern archipelago is part of the North American Craton. The mantle root of the cratonic block beneath northern Greenland may extend westwards as far as central Ellesmere Island. The Barents and Kara seas show high velocities indicative of thick lithosphere, similar to cratons. The locations of intraplate basaltic volcanism attributed to the High Arctic Large Igneous Province are all on thin, non-cratonic lithosphere. The lithosphere beneath the central part of the Siberian Traps is warmer than elsewhere beneath the Siberian Craton. This observation is consistent with lithospheric erosion associated with the large igneous province volcanism. A corridor of relatively low seismic velocities cuts east–west across central Greenland. This indicates lithospheric thinning, which appears to delineate the track of the Iceland hotspot.Supplementary material: Figures with comparisons of different tomographic models at 50 and 200 km depths are available at https://doi.org/10.6084/m9.figshare.c.3817810

Three major problems in cave micrometeorology are analysed: the concept of the temperature of a cave and its phenomenology; the internal energy flows and consequent local entropy production; and a non-hydrostatic physical model of the underground convective air circulation. A cave’s temperature is rich in information, but it is often difficult to obtain because it requires experimental accuracies to be pushed beyond the reach of common external meteorological instruments to detect a variety of factors, such as thermal sedimentation, seasonal variations and the effects of external morphology. Energy flow has an essential role in estimating the sensitivity of a cave to external inputs. A model for evaluating the local entropy production is developed. Entropy seems to be the most basic parameter, explaining both the sensitivity of the environment and its tendency to form complex structures. Analysis of the second-order terms of convective air circulation is more complex than expected; nevertheless, only such an analysis is able to explain the behaviours (e.g. continuity in the airflow) that cannot be predicted by the first-order models.

GPUSPH is a fully three-dimensional model for the simulation of the thermal and rheological evolution of lava flows that relies on the Smoothed Particle Hydrodynamics (SPH) numerical method. Thanks to the Lagrangian, meshless nature of SPH, the model incorporates a more complete physical description of the emplacement process and rheology of lava that considers the free surface, the irregular boundaries represented by the topography, the solidification fronts and the non-Newtonian rheology with temperature-dependent parameters. GPUSPH follows the very general Herschel–Bulkley rheological model, which encompasses Newtonian, power-law and Bingham flow behaviours, with both constant and temperature-dependent parameters, and can thus be used to explore in detail the impact of rheology on the behaviour of lava flows and on their emplacement. To illustrate this possibility, we present some preliminary applications of the model for studying the rheology of lava flows with different constitutive relationships and thermal regimes using the real topography of the Mt Etna volcano.

The paper focuses on the Romanche transform margin on the African side of the Equatorial Atlantic, and draws from thermomechanical numerical modelling and subsequent GIS-based thermal history grid processing.Our modelling indicates that the early post-break-up thermal history of the transform margin is controlled by the cooling patterns of the adjacent pull-apart terrains, the pre-rift heat-flow regime, laterally passing seafloor spreading centres, and cooling of the newly accreted oceanic crust in two corridors located in front of and behind the spreading centres.The pre-rift thermal regime controls the background heat flow on top of which the thermal transients develop. If it is cold, the transient anomalies are very distinct on this background. If it is warm, the transient anomalies tend to blend in with background heat flow much better. The most prominent anomalies are related to thinning in the pull-apart terrains located between transforms. They are up to three times wider than anomalies related to the laterally passing spreading centre, the anomaly of which becomes no wider than 20–25 km. Together with oceanic crust accreted in corridors in front of and behind the passing centre, its heat transfer into the transform margin tends to slow down the cooling of the syn-break-up anomalies developed in the pull-apart terrains by rifting. AbstractThe paper focuses on the Romanche transform margin on the African side of the Equatorial Atlantic, and draws from thermomechanical numerical modelling and subsequent GIS-based thermal history grid processing.Our modelling indicates that the early post-break-up thermal history of the transform margin is controlled by the cooling patterns of the adjacent pull-apart terrains, the pre-rift heat-flow regime, laterally passing seafloor spreading centres, and cooling of the newly accreted oceanic crust in two corridors located in front of and behind the spreading centres.The pre-rift thermal regime controls the background heat flow on top of which the thermal transients develop. If it is cold, the transient anomalies are very distinct on this background. If it is warm, the transient anomalies tend to blend in with background heat flow much better. The most prominent anomalies are related to thinning in the pull-apart terrains located between transforms. They are up to three times wider than anomalies related to the laterally passing spreading centre, the anomaly of which becomes no wider than 20–25 km. Together with oceanic crust accreted in corridors in front of and behind the passing centre, its heat transfer into the transform margin tends to slow down the cooling of the syn-break-up anomalies developed in the pull-apart terrains by rifting.

A three-dimensional (3D) thermal–kinematic modelling approach based on finite-element techniques is used to study lower-crustal viscosity at transform margins during the continent–ocean transform development stage and after the ridge has passed by. Nine modelling scenarios combining different equilibrium surface heat flows and lower-crustal rheologies are studied. Modelling results indicate that substantial parts of the lower crust at transform margins have the potential to flow at geologically appreciable strain rates, which can lead to uplift/subsidence, as well as lateral variations, in upper- and lower-crustal thicknesses and Moho depth. These low-viscosity zones (i.e. parts of the lower crust with effective viscosities of less than 1018 Pa s) make up distinct ductility distributions that vary in space and time during margin evolution. Three basic ductility patterns and related thermal processes can be identified: reduced lower-crustal viscosities originating at the continental rift and the continent–ocean boundary (COB), respectively; reduced lower-crustal viscosities along the transform caused by the migrating ridge; and the background distribution of lower-crustal ductility resulting from the equilibrium temperature field. Superposition of all three ductility patterns and the complex interaction of the underlying perturbations of the temperature field result in distinct differences in the potential of lower-crustal flow both in space (parallel and perpendicular to the transform) and with time. Thus, modelling results provide templates for understanding lower-crustal flow at transform margins in general and await further studies comparing model predictions with actual field observations. AbstractA three-dimensional (3D) thermal–kinematic modelling approach based on finite-element techniques is used to study lower-crustal viscosity at transform margins during the continent–ocean transform development stage and after the ridge has passed by. Nine modelling scenarios combining different equilibrium surface heat flows and lower-crustal rheologies are studied. Modelling results indicate that substantial parts of the lower crust at transform margins have the potential to flow at geologically appreciable strain rates, which can lead to uplift/subsidence, as well as lateral variations, in upper- and lower-crustal thicknesses and Moho depth. These low-viscosity zones (i.e. parts of the lower crust with effective viscosities of less than 1018 Pa s) make up distinct ductility distributions that vary in space and time during margin evolution. Three basic ductility patterns and related thermal processes can be identified: reduced lower-crustal viscosities originating at the continental rift and the continent–ocean boundary (COB), respectively; reduced lower-crustal viscosities along the transform caused by the migrating ridge; and the background distribution of lower-crustal ductility resulting from the equilibrium temperature field. Superposition of all three ductility patterns and the complex interaction of the underlying perturbations of the temperature field result in distinct differences in the potential of lower-crustal flow both in space (parallel and perpendicular to the transform) and with time. Thus, modelling results provide templates for understanding lower-crustal flow at transform margins in general and await further studies comparing model predictions with actual field observations.