Current projects

Crustal magnetic source depths on Mars:


The dawn of high-resolution observations with the James Webb Space Telescope will enable spatially resolved observations of ultrashort-period rocky exoplanets. Some of these planets orbit so closely to their host star that they lack an atmosphere, which gives direct access to their surfaces and opens a window to infer their geodynamics. GJ 367 b has been observationally constrained to a planetary radius of about 0.72 to 0.75 Earth radii and a mass between 0.48 and 0.55 Earth masses, which puts this planets in a Mercury-like interior regime with a thin mantle overlying a large core. The dayside temperature ranges between 1500 and 1800 K, suggesting a permanent magma ocean or dayside magma pond on the surface, induced by stellar irradiation. We perform global 2D spherical annulus StagYY simulations of solid-state mantle convection and surface melting with the goal to constrain the geometric and compositonal properties of the planet. The combination of observations and numerical models offers a unique opportunity to constrain the compositional fractionation during magma ocean epochs and provides avenues to constrain the delivery and loss cycle of atmophile elements on strongly irradiated exoplanets.

Collaborators: T. Lichtenberg / Univ. Groningen, T. G. Meier / Univ. Oxford, D. J. Bower, P. J. Tackley / ETH Zurich


Explaining the peculiar separation of metal and troilite pockets in acapulcoites:

Acapulcoites are primitive achondrites that are derived from a partially differentiated parent body that experienced peak temperatures above the Fe-Ni-S solidus. Their textures reflect annealing driven by the minimization of surface energy, yielding smoothly curved silicate grain boundaries, common 120 degrees triple junctions and an equiaxial shape of Fe-Ni and FeS pockets. However, one peculiar obervation is not easily explained: Fe-Ni and FeS pockets are often spatially separated and seldom form common boundaries in acapulcoites. Static annealing experiments were performed in a Walker-type multianvil press using a model system of olivine with 5 wt% gold. The mixed powders were compressed to 2.5 GPa and equilibrated at 1300 degrees C before rapidly decreasing to 1000 degrees C, below the Au and FeS solidi. At this temperature the samples were statically annealed for durations between three hours and one month to mimic slow cooling of planetesimals. Our experiments successfully simulate the phase separation and the decrease of Fe-Ni-FeS boundary length observed in H3-H6 chondrites and acapulcoites.

Collaborators: A. Neri / Univ. Lille, E. Kubik, A. Bouvier / BGI Bayreuth, N. P. Walte / TU Munich

Completed projects

Probing the viscosity of Venus' mantle:

The Baltis Vallis channel on Venus preserved a record of long-wavelength deformation generated by a convecting mantle, offering a unique window into the planet's geodynamics. We statistically compare the observed dynamic topography of Baltis Vallis with dynamic topography generated by a suite of 3D stagnant-lid mantle convection models. Baltis Vallis' relatively young age and low topography indicate vigorous convection in Venus' mantle, implying a mantle viscosity 1-2 order of magnitude lower than Earth's. This difference may arise from either a water-rich, less degassed interior or a higher-temperature mantle beneath the insulating lid.

N. J. McGregor, F. Nimmo, C. Gillmann, G. J. Golabek , A. M. Plattner and J. W. Conrad (2025). Probing the viscosity of Venus's mantle from dynamic topography at Baltis Vallis. J. Geophys. Res. Planets 130, e2024JE008581.

Coverage:  EGU Blog


Tungsten isotope evolution during Earth's formation:

The Hf-W isotopic system is the reference chronometer for determining the chronology of Earth's accretion and differentiation. However, its results depend strongly on uncertain parameters, including the extent of metal-silicate equilibration and the siderophility of tungsten. Here we show that a multistage core-formation model based on N-body accretion simulations, element mass balance and metal-silicate partitioning, largely eliminates these uncertainties. We considered smoothed particle hydrodynamics estimates of the depth of melting caused by giant impacts and the isotopic evolution of ¹⁸²W. We applied two metal-silicate fractionation mechanisms: (i) the metal delivered by the cores of large impactors equilibrates with only a small fraction of the impact-induced magma pond and (ii) metal delivered by small impactors emulsifies in global magma oceans before undergoing progressive segregation. The latter is crucial for fitting the W abundance and ¹⁸²W anomaly of Earth's mantle. In addition, we show that the duration of magma ocean solidification has a major effect on Earth's tungsten isotope anomaly. Depending on the characteristics of the giant impacts, results predict that the Moon formed either 143-183 Myr or 53-62 Myr after the start of the solar system. Thus, independent evaluations of the Moon's age provide an additional constraint on the validity of accretion simulations.

D. C. Rubie, K. I. Dale, G. Nathan, M. Nakajima. E. S. Jennings, G. J. Golabek, S. A. Jacobson and A. Morbidelli (2025). Tungsten isotope evolution during Earth's formation and new constraints on the viability of accretion simulation. Earth Planet. Sci. Lett. 651, 119139.

Grain size evolution in Earth's lower mantle:

Recent experimental investigations of grain size evolution in bridgmanite-ferropericlase assemblages have suggested very slow growth for these bimodal phases. We develop self-consistent 2-D spherical half-annulus geodynamic models of Earth's evolution using the finite volume code StagYY to assess the role of grain size on lower mantle viscosity. We are considering three scenarios: (i) uniform grain growth throughout the entire mantle with a composite rheology, (ii) different grain growth in the upper and lower mantle with a composite rheology, and (iii) different grain growth in the upper and lower mantle with purely diffusion creep rheology. In the case of different grain size evolution, the upper mantle's grain size evolution law is controlled by forsterite-enstatite grain growth, while the lower mantle's grain size evolution law is controlled by bridgmanite-ferropericlase grain growth. Our results suggest that mantle viscosity is primarily controlled by temperature, whereas grain size has a minor effect compared to the effect of temperature. To establish the robustness of this finding we vary several other model parameters, such as surface yield strength, phase transition grain size reset, different transitional stresses for creep mechanisms, pressure dependence on grain growth, and different grain damage parameters. For all our models, we consistently find that grain size has a very limited effect on controlling lower mantle viscosity in the present-day Earth. However, large grain size may have affected the lower mantle viscosity in the early Earth as larger grains of single phase bridgmanite could increase the viscosity of the early mantle delaying the onset of global convection.

J. Paul, G. J. Golabek, A. B. Rozel, P. J. Tackley, T. Katsura and H. Fei (2024). Effect of bridgmanite-ferropericlase grain size evolution on Earth's average mantle viscosity: Implications for mantle convection in early and present-day Earth. Prog. Earth Planet. Sc. 11, 64.

Combined impact and interior evolution models for Mars:

Here we use a coupled SPH-thermochemical approach to first simulate an impact event, and then use the result of this realistic model as the initial condition for the long-term mantle convection model. We find that both the impact scenario and the mantle properties affect the geometry of impact-induced crust and the subsequent state of the interior, and that the impact-induced crust extends beyond the initial magma pond. We show that a near head-on impact event with impactor radius of 750 km can best reproduce the southern highlands of Mars with a geometry similar to that of present-day Mars observations.

K. W. Cheng, H. A. Ballantyne, G. J. Golabek, M. Jutzi, A. B. Rozel & P. J. Tackley (2024). Combined impact and interior evolution models in three dimensions indicate a southern impact origin of the Martian Dichotomy. Icarus 420, 116137.

Mantle depletion after martian giant impact:

Previous giant impact models coupled with simulations of mantle convection have shown that the martian crustal dichotomy can be explained by post‐impact melt crystallization that emplaced a thick crust in the southern hemisphere. In this study, we show that the depleted residue left behind by the original post‐impact crustal formation can spread laterally, potentially persisting beneath the northern hemisphere to the present‐day. Such a large‐scale mantle province would concurrently explain both the prevalence of long‐term magmatism on Mars and its strong preference for localized equatorial regions.

K. W. Cheng, A. B. Rozel, G. J. Golabek, H. A. Ballantyne, M. Jutzi & P. J. Tackley (2024). Mars's Crustal and Volcanic Structure Explained by Southern Giant Impact and Resulting Mantle Depletion. Geophys. Res. Lett. 51, e2023GL105910.

Animation:  Movie

Induction heating of planets orbiting white dwarfs:

Recently a number of planets have been discovered that orbit white dwarfs. We investigate whether electromagnetic induction heating due to the extremely strong magnetic fields of some white dwarfs can drive volcanic activity even on small bodies orbiting such white dwarfs on short-period orbits. Our results show that induction heating can indeed rapidly partially melt the mantle of lunar-mass objects orbiting white dwarfs with strong magnetic fields. The related volcanism drives outgassing that could be observed in the future around white dwarfs.

K. G. Kislyakova, L. Noack, E. Sanchis, L. Fossati, G. G. Valyavin, G. J. Golabek & M. Guedel (2023). Induction heating of planetary interiors in white dwarf systems. Astron. Astrophys. 677, A109.


Investigating the feasibility of an impact-induced Martian Dichotomy:

Most studies attempting to explain the martian dichtomy focus on two theories; either the dichotomy formed solely through geodynamic processes or a giant impact occurred that imprinted the crustal cavity in the northern hemisphere we observe today. Recent work has proved the importance of coupling these hypotheses, introducing a hybrid exogenic-endogenic scenario with crust formation from the impact-generated melt. We model the impact using a suite of SPH simulations that explore a large parameter space. For the outcome of the simulation outcomes we apply a newly developed scheme to estimate the thickness and distribution of new formed or re-distributed post-impact crust.

H. A. Ballantyne, M. Jutzi, G. J. Golabek, L. Mishra, K. W. Cheng, A. Rozel & P. J. Tackley (2023). Investigating the feasibility of an impact-induced Martian dichotomy. Icarus 392, 115395.


Fluid flux determination for serpentinites in subduction zones:

We present a new method to determine fluid flux at subduction zone conditions. Orientated drill cores of antigorite starting material were placed into a MgO sleeve. Fluid released upon serpentine dehydration are fixed in the MgO sleeve by the formation of brucite, thus the amount and distribution of brucite in the run product serves as a proxy for fluid flux. The experiments show that dehydration results in an isotropic fluid flux. The timescales of fluid flux indicate that a large fluid pressure can build up during dehydration. Intermediate-depth earthquakes caused by brittle fracturing are thus likely to be caused by dehydration reactions in subducted serpentinites.

L. Eberhard, M. Thielmann, P. Eichheimer, A. Neri, A. Suzuki, M. Ohl, W. Fujita, K. Uesugi, M. Nakamura, G. J. Golabek & D. J. Frost (2022). A new method for determining fluid flux at high pressures applied to the dehydration of subduction zone serpentinites. Geochem. Geophys. Geosyst. 23, e2021GC010062.

Coverage:  Spektrum


Strain-weakening rheology in Earth's lower mantle:

The rheological properties of Earth's lower mantle are key for mantle dynamics and planetary evolution. The main rock-forming minerals in the lower mantle are bridgmanite (Br) and smaller amounts of ferropericlase (Fp). Previous work has suggested that the large differences in viscosity between these minerals greatly affect the bulk rock rheology. The resulting effective rheology becomes highly strain-dependent as weaker Fp minerals become elongated and eventually interconnected. This implies that strain localization may occur in Earth's lower mantle. So far, there have been no studies on global-scale mantle convection in the presence of such strain-weakening rheology. Here, we present 2D numerical models of thermo-chemical convection in spherical annulus geometry including a new strain-dependent rheology formulation for lower mantle materials, combining rheological weakening and healing terms. We find that strain-weakening rheology has several direct and indirect effects on mantle convection. The most notable direct effect is the changing dynamics of weakened plume channels as well as the formation of larger thermochemical piles at the base of the mantle. The weakened plume conduits act as lubrication channels in the mantle and exhibit a lower thermal anomaly. SW rheology also reduces the overall viscosity, notable in terms of increasing convective vigor and core-mantle boundary heat flux. Finally, we put our results into context with existing hypotheses on the style of mantle convection and mixing. Most importantly, we suggest that the new kind of plume dynamics may explain the discrepancy between expected and observed thermal anomalies of deep-seated mantle plumes on Earth.

A. J. P. Guelcher, G. J. Golabek, M. Thielmann, M. D. Ballmer & P. J. Tackley (2022). Narrow, fast and "cool" mantle plumes caused by strain-weakening rheology. Geochem. Geophys. Geosyst. 23, e2021GC010314.

Animations:  Movie1   Movie2   Movie3   Movie4   Movie5


Detailed analysis of main group pallasite slabs:



Modification of icy planetesimal interiors by early thermal evolution and collisions:
Two-stage formation of angular pallasites revealed by novel high strain-rate deformation experiments:

Pallasite formation is controversial, either sampling the core-mantle boundary or the shallower mantle of planetesimals that suffered an impact. In order to test this model, we performed high strain-rate deformation experiments employing the analogue system olivine -FeS/Au metal melt. All major textures of angular pallasites, including olivine-metal mush zones and olivine aggregates could be experimentally reproduced. Our results suggest that angular pallasites preserve evidence for a two-stage formation process, including inefficient core-mantle differentiation and an impact causing disruption, core melt injection, and fast cooling. Olivine aggregates are reinterpreted as mantle remnants containing primordial metallic melt pockets not stemming from the impactor.

N. P. Walte, G. F. D. Solferino, G. J. Golabek, D. Silva Souza & A. Bouvier (2020). Two-stage formation of pallasites and the evolution of their parent bodies revealed by deformation experiments. Earth Planet. Sci. Lett. 546, 116419.

Coverage:  TUM News   FRM News   Phys.org


Numerical and experimental study of microstructure and permeability in porous granular media:

Fluid flow on different scales is of interest for Earth science disciplines like petrophysics, hydrogeology and volcanology. To parameterize fluid flow in large-scale numerical simulations, flow properties on the microscale need to be considered. For this purpose experimental and numerical investigations of flow through porous media over a wide range of porosities are necessary. We sinter glass bead media with various porosities and determine permeability both experimentally and numerically. Using CT scans we investigate the microstructure of the samples. We test different parameterizations for isotropic porous media on their potential to predict permeability by comparing their estimations to computed and experimentally measured values. 

P. Eichheimer, M. Thielmann, W. Fujita, G. J. Golabek, M. Nakamura, S. Okumura, T. Nakatani & M. O. Kottwitz (2020). Combined numerical and experimental study of microstructure and permeability in porous granular media. Solid Earth 11, 1079-1095.

Animation: Movie


Effect of water on ringwoodite thermal conductivity and its influence on the thermal evolution of slabs:

Ringwoodite is the most abundant phase in the lowermost mantle transition zone and can host up to 1.5-2 wt% water. We studied the high-pressure thermal conductivity of dry and hydrous ringwoodite. We show that the incorporation of 1.73 wt% of water reduces the ringwoodite thermal conductivity by more than 40% at mantle transition zone pressures. Using simple 1D thermal diffusion models we demonstrate that this effect delays the decomposition of dense hydrous magnesium silicates, thus allowing them to reach Earth's lower mantle.

E. Marzotto, W.-P. Hsieh, T. Ishii, K.-H. Chao, G. J. Golabek, M. Thielmann & E. Ohtani (2020). Effect of water on lattice thermal conductivity of ringwoodite and its implications for the thermal evolution of descending slabs. Geophys. Res. Lett. 47, e2020GL087607.


N-body late accretion models coupled with atmosphere-interior models of Venus:

It remains contentious whether the late accretion after the end of core formation was rich or poor in water and other volatiles. Here we investigate the long-term evolution of Venus using self-consistent numerical simulations of global thermochemical mantle convection coupled with both an atmospheric evolution model and a late accretion N-body delivery model. The results suggest that a preferentially dry composition of the late accretion impactors is most consistent with measurements of atmospheric H₂O, CO₂ and N₂. We suggest that the late accreted material delivered to Venus was mostly dry enstatite chondrite, consistent with isotopic data for Earth, with less than 2.5% (by mass) wet carbonaceous chondrites. In this scenario, the majority of Venus's and Earth's water would have been delivered during the main accretion phase.

C. Gillmann, G. J. Golabek, S. N. Raymond, M. Schoenbaechler, P. J. Tackley, V. Dehant & V. Debaille (2020). Dry late accretion inferred from Venus's coupled atmosphere and internal evolution. Nat. Geosci. 13, 265-269.

Coverage:  News & Views   Nature Community Blog   Royal Observatory of Belgium News   Daily Science   PlanetS   Astrobiology   Astrobites


Microscale deformation of bridgmanite-ferropericlase mixtures:

The mineralogy of the lower mantle can be approximated as a bridgmanite-ferropericlase mixture. Previous work has suggested that the deformation of this mixture might be dramatically affected by the large differences in viscosity between bridgmanite and ferropericlase. We employ numerical models to establish a connection between ferropericlase morphology and the effective rheology of the Earth's lower mantle using a numerical-statistical approach. We link the statistical properties of the two-phase composite to its effective viscosity tensor using analytical approximations. We find that ferropericlase develops elongated structures within the bridgmanite matrix that result in significantly lowered viscosity. We show that significant rheological weakening can thus be already achieved even when ferropericlase does not form an interconnected network. Additionally, the alignment of weak ferropericlase leads to a pronounced viscous anisotropy that develops with total strain, which may have implications for understanding the viscosity structure of Earth's lower mantle as well as for modelling the behaviour of subducting slabs.

M. Thielmann, G. J. Golabek & H. Marquardt (2020). Ferropericlase control of lower mantle rheology: Impact of phase morphology. Geochem. Geophys. Geosyst. 21, e2019GC008688.


Desiccation of planetesimals due to ²⁶Al heating:

Here we demonstrate the power of ²⁶Al, a short-lived radionuclide abundant in the early Solar system, to control the water content of terrestrial exoplanets by rapid dehydration of planetesimals prior to accretion. Using numerical models of planet formation, evolution, and interior structure, we generate synthetic planet populations influenced by varying levels of ²⁶Al heating during accretion. We show that planet bulk water fraction and radius are anti-correlated with initial ²⁶Al levels in the planetesimal-based accretion framework. This yields a system-wide correlation of bulk abundances, consistent with the lack of a clear orbital trend in the water budgets of the TRAPPIST-1 planets.

T. Lichtenberg, G. J. Golabek, R. Burn, M. R. Mayer, Y. Alibert, T. V. Gerya & C. Mordasini (2019). A water budget dichotomy of rocky protoplanets from ²⁶Al heating. Nat. Astron. 3, 307-313.

Coverage:  Nature Community Blog   PlanetS   Phys.org   Astronomy Magazine   Discover Magazine   ETH News   Univ. Michigan News   Welt der Physik   Der Standard   Kleine Zeitung   Swissinfo   Tagesanzeiger


Two-phase flow models of melt evolution in early-formed planetesimals:

Rocky planetesimals in the early solar system melted internally and evolved chemically due to radiogenic heating from ²⁶Al. Here we quantify the parametric controls on magma genesis and transport using a coupled petrological and fluid mechanical model of reactive two-phase flow. We find the mean grain size of silicate minerals to be a key control on magma ascent.
According to our models, the evolution of partially molten planetesimal interiors falls into two categories. In the global magma ocean scenario, the whole interior of a planetesimal experiences nearly complete melting, resulting in turbulent convection and core-mantle differentiation by the rainfall mechanism. In the magma sill scenario, segregating melts gradually deplete the deep interior of the radiogenic heat source. In this case, magma may form melt-rich sills beneath a cool and stable lid, while core formation would proceed by percolation.

T. Lichtenberg, T. Keller, R. F. Katz, G. J. Golabek & T. V. Gerya (2019). Magma ascent in planetesimals: Control by grain size. Earth Planet. Sci. Lett. 507, 154-165.

Animation:  Movie


Formation of mixed-type pallasite meteorites:

The origin of pallasite meteorites remains elusive. Here we test the hypothesis of mixing of olivine fragments with Fe-Ni-S after an impact followed by annealing employing both experimental analogues and numerical models. The experiments show that the addition of sulfur to olivine + Fe-Ni accelerates olivine grain growth, though the growth rate is reduced when Fe-Ni-S is not fully molten. Numerical models satisfying available formation constraints from natural samples indicate that planetesimals with radii ≥200 km are favorable for the formation of rounded olivine-bearing pallasites by annealing of fragments in partially molten Fe-Ni-S. Moreover, early mixing in the planetesimal can form regions containing olivine grains with different grain sizes that could explain the formation of mixed-type pallasites.

G. F. D. Solferino & G. J. Golabek (2018). Olivine grain growth in partially molten Fe-Ni-S: A proxy for the genesis of pallasite meteorites. Earth Planet. Sci. Lett. 504, 38-52.

Animation:  Movie


Constraints on the metal-silicate separation on the IAB parent body:

The short-lived Hf-W decay system is a powerful chronometer for constraining the timing of metal-silicate separation and core formation in planetesimals and planets. Neutron capture effects on W isotopes, however, significantly hamper the application of this tool. In order to correct for neutron capture effects, Pt isotopes have emerged as a reliable in-situ neutron dosimeter. This study applies this method to IAB iron meteorites, in order to constrain the timing of metal segregation on the IAB parent body. Thermal models of the interior evolution that are consistent with the obtained estimates suggest that the IAB parent body underwent metal-silicate separation as a result of internal heating by short-lived radionuclides and accreted at around 1.4 Myr after CAIs with a radius of greater than 60 km.

A. C. Hunt, D. Cook, T. Lichtenberg, P. M. Reger, M. Ek, G. J. Golabek & M. Schoenbaechler (2018). Late metal-silicate separation on the IAB parent asteroid: Constraints from combined W and Pt isotopes and thermal modelling. Earth Planet. Sci. Lett. 482, 490-500.


Formation of chondrules via impact splashes:

Chondrules, mm-sized igneous-textured spherules, are the dominant bulk silicate constituent of chondritic meteorites and originate from highly energetic, local processes during the first million years after the birth of the Sun. So far, an astrophysically consistent chondrule formation scenario explaining major chemical, isotopic and textural features, in particular Fe,Ni metal abundances, bulk Fe/Mg ratios and intra-chondrite chemical and isotopic diversity, remains elusive. Here, we examine the prospect of forming chondrules from impact splashes among planetesimals heated by radioactive decay of short-lived radionuclides using thermomechanical models of their interior evolution. We examine how the observed chemical and isotopic features of chondrules constrain the dynamical environment of accreting chondrite parent bodies by interpreting the meteoritic record as an impact-generated proxy of early solar system planetesimals that underwent repeated collision and reaccretion cycles. Using a coupled evolution-collision model we demonstrate that the vast majority of collisional debris feeding the asteroid main belt must be derived from planetesimals which were partially molten at maximum. Therefore, the precursors of chondrite parent bodies either formed primarily small, from sub-canonical aluminium-26 reservoirs, or collisional destruction mechanisms were efficient enough to shatter planetesimals before they reached the magma ocean phase.

T. Lichtenberg, G. J. Golabek, C. P. Dullemond, M. Schoenbaechler, T. V. Gerya & M. R. Meyer (2018). Impact splash chondrule formation during planetesimal recycling. Icarus 302, 27-43.


Coupled giant impact and longer-term interior evolution models:

Giant impacts have been suggested to explain various characteristics of terrestrial planets and their moons. However, so far in most models only the immediate effects of the collisions have been considered, while the long-term interior evolution of the impacted planets was not studied. We combine 3-D shock physics collision calculations with 3-D thermochemical interior evolution models.
We apply the combined methods to a demonstration example of a giant impact on a Mars-sized body, using typical collisional parameters from previous studies. While the material parameters (equation of state, rheology model) used in the impact simulations can have some effect on the long-term evolution, we find that the impact angle is the most crucial parameter for the resulting spatial distribution of the newly formed crust. The results indicate that a dichotomous crustal pattern can form after a head-on collision, while this is not the case when considering a more likely grazing collision. The results underline that end-to-end 3-D calculations of the entire process are required to study in the future the effects of large-scale impacts on the evolution of planetary interiors.

G. J. Golabek, A. Emsenhuber, M. Jutzi, E. I. Asphaug & T. V. Gerya (2018). Coupling SPH and thermochemical models of planets: Methodology and example of a Mars-sized body. Icarus 301, 235-246.

Animations:  Movie1   Movie2

Formation of the Earth's oldest continental crust:

The global geodynamic regime of early Earth, which operated before the onset of plate tectonics, remains disputed. Geological and geochemical data suggest hotter Archean mantle temperature and more intense juvenile magmatism. Two alternative crust-mantle interaction modes differing in melt eruption efficiency have been proposed: the pipe tectonics regime dominated by volcanism and the squishy lid tectonics regime governed by intrusive magmatism. Both regimes are capable of producing primordial tonalite-trondhjemite-granodiorite (TTG) continental crust, however lithospheric geotherms and crust production rates as well as proportions of various TTG compositions differ significantly. This provides an opportunity to test the heat-pipe and plutonic squishy lid tectonics hypotheses on the basis of natural data. We show that the heat pipe tectonics model results is not able to produce Earth-like primordial continental crust. In contrast, the plutonic squishy lid tectonics regime dominated by intrusive magmatism results is capable of reproducing the observed proportions of various TTG-rocks. The results show that the typical modern eruption efficiency leads to the production of the expected amount of the three main primordial crustal compositions previously reported from field data.

A. B. Rozel, G. J. Golabek, C. Jain, P. J. Tackley & T. Gerya (2017). Continental crust formation of early Earth controlled by intrusive magmatism. Nature 545, 332-335.

Animation:  Movie

Coverage:  ETH News   EGU Blog   Phys.org


Effect of porosity on the thermo-mechanical evolution of planetesimals:

The thermal history and internal structure of chondritic planetesimals, assembled before the giant impact phase of chaotic growth, potentially yield important implications for the final composition and evolution of terrestrial planets. These parameters critically depend on the internal balance of heating versus cooling, which is mostly determined by the presence of short-lived radionuclides (SLRs) as well as the heat conductivity of the material. The heating by SLRs depends on their initial abundances, the formation time of the planetesimal and its size. It has been argued that the cooling history is determined by the porosity of the granular material, which undergoes dramatic changes via compaction processes and tends to decrease with time. In this study we assess the influence of these parameters on the thermo-mechanical evolution of young planetesimals with both 2D and 3D simulations. Our results indicate that powdery materials lower the threshold for melting and convection in planetesimals, depending on the amount of SLRs present. A subset of planetesimals retains a powdery surface layer which lowers the thermal conductivity and hinders cooling. However, the effect of initial porosity is small compared to those of planetesimal size and formation time, which dominate the thermo-mechanical evolution and are the primary factors for the onset of melting and differentiation.

T. Lichtenberg, G. J. Golabek, T. V. Gerya & M. R. Meyer (2016). The effects of short-lived radionuclides and porosity on the early thermo-mechanical evolution of planetesimals. Icarus 274, 350-365.

Animation:  Movie


Influence of single large impacts on the coupled mantle-atmosphere evolution of Venus:

We investigate the effect of a single large impact during the Late Veneer and Late Heavy Bombardment on the evolution of the mantle and atmosphere of Venus. We use a coupled interior-exterior model. Single vertical impacts are simulated as instant events affecting both the atmosphere and mantle of the planet by (i) eroding the atmosphere, causing atmospheric escape, and (ii) depositing energy in the crust and mantle of the planet. Main impactor parameters include timing, size/mass, velocity and efficiency of energy deposition. We observe that volatile delivery by the impactor and impact erosion of the atmosphere are both minor effects compared to melting and degassing triggered by the energy deposited in the mantle and crust. Small collisions (under 100 km radius) have only local and time-limited effects. Medium-sized impactors (100-400 km) will not have much more consequences unless the energy deposition is enhanced, for example by a fast collision. In that case, it will have comparable effects to the larger category of impacts (400-800 km): a strong thermal anomaly affecting both crust and mantle and triggering melting and a change in mantle dynamics patterns.

C. Gillmann, G. J. Golabek & P. J. Tackley (2016). Effect of a single giant impact on the coupled atmosphere-interior evolution of Venus. Icarus 268, 295-312.

Animation:  Movie

Coverage:  Smithsonian    L'Obs


Appearance of the 'ridge only' regime in mantle convection models with a visco-plastic rheology:

In this study, we report a new convection regime in which a ridge can form without destabilizing the surrounding lithosphere or forming subduction zones. Using simulations in 2D spherical annulus geometry, we show that a depth-dependent yield stress is sufficient to reach this ridge only regime. This regime occurs when the friction coffecient is close to the critical value between mobile lid and stagnant lid regime. The ridge only regime appears for both pure basal heating and mixed heating mode. For basal heating, this regime can occur for all vertical viscosity contrasts, while for mixed heating, a highly viscous deep mantle is required.

A. Rozel, G. J. Golabek, R. Naef & P. J. Tackley (2015). Formation of ridges in a stable lithosphere in mantle convection simulations with visco-plastic rheology. Geophys. Res. Lett. 42, 4770-4777.

Animation:  Movie


Formation of rounded olivine grains in pallasites:

Despite their relatively simple mineralogical composition, the origin of pallasite meteorites remains debated. It has been suggested that catastrophic mixing of olivine fragments with Fe-(Ni)-S followed by various degrees of annealing could explain pallasites bearing solely or prevalently fragmented or rounded olivines. In order to verify this hypothesis, and to quantify the grain growth rate of olivine in a metal matrix, we performed a series of annealing experiments on natural olivine plus synthetic Fe-S mixtures. Olivine grain growth in molten Fe-S is significantly faster than in solid-, sulphur-free metal. We used the experimentally determined grain growth law to model the coarsening of olivine surrounded by Fe-S melt in a 100 to 600 km radius planetesimal. Employing a 1D finite-difference model, annealing depths of up to 50 km allow for (i) average grain sizes consistent with the observed rounded olivine-bearing pallasites, (ii) a remnant magnetization of FeNi olivine inclusions as measured in natural pallasites and (iii) for the metallographic cooling rates derived from FeNi in pallasites. This is valid even if the impact occurs several millions of years after the differentiation of the target body was completed.

G. F. D. Solferino, G. J. Golabek, F. Nimmo & M. W. Schmidt (2015). Fast grain growth of olivine in liquid Fe-S and the formation of pallasites with rounded olivine grains. Geochim. Cosmochim. Acta 162, 259-275.


Is Vesta an intact and pristine protoplanet?

It is difficult to find a Vesta model of iron core, pyroxene and olivine-rich mantle, and HED crust that can match the joint constraints of (a) Vesta's density and core size as reported by the Dawn spacecraft team; (b) the chemical trends of the HED meteorites, including the depletion of sodium, the FeO abundance, and the trace element enrichments; and (c) the absence of exposed mantle material on Vesta's surface, among Vestoid asteroids, or in our collection of basaltic meteorites. These conclusions are based entirely on mass-balance and density arguments, independent of any particular formation scenario for the HED meteorites themselves. We suggest that Vesta either formed from source material with non-chondritic composition or underwent after its formation a radical physical alteration, possibly caused by collisional processes, that affected its global composition and interior structure.

G. J. Consolmagno, G. J. Golabek, D. Turrini, M. Jutzi, S. Sirono, V. Svetsov & K. Tsiganis (2015). Is Vesta an intact and pristine protoplanet? Icarus 254, 190-201.

Press coverage:  Arizona Daily Star


Size and formation time of the acapulcoite-lodranite parent body:

The acapulcoite-lodranite meteorites are members of the primitive achondrite class. The observation of partial melting and resulting partial removal of Fe-FeS indicates that this meteorite group could be an important link between achondrite and iron meteorites on the one hand and chondrite meteorites on the other. Thus a better understanding of the thermomechanical evolution of the parent body of this meteorite group can help to improve our understanding of the evolution of early planetesimals. Here we use 2D and 3D finite-difference numerical models to determine the formation time, initial radius of the parent body of the acapulcoite-lodranite meteorites and their formation depth inside the body by applying available geochronological, thermal and textural constraints to our numerical data. Our results indicate that the best fit to the data can be obtained for a parent body with 25-65 km radius, which formed around 1.3 Myr after CAI. The 2D and 3D results considering various initial temperatures and the effect of porosity indicate possible formation depths of the acapulcoite-lodranite meteorites of 9-19 and 14-25 km, respectively. Our data also suggest that other meteorite classes could form at different depths inside the same parent body.

G. J. Golabek, B. Bourdon & T. V. Gerya (2014). Numerical models of the thermomechanical evolution of planetesimals: Application to the acapulcoite-lodranite parent body. Meteorit. Planet. Sci. 49, 1083-1099.

Animation:  Movie


Single-plume state on Enceladus:

The thermal dichotomy of Enceladus suggests an asymmetrical structure in its global heat transfer. So far, most of the models proposed that obtained such a distribution have prescribed an a priori asymmetry, i.e. some anomaly in or below the south polar ice shell. We present here the first set of numerical models of convection that yield a stable single-plume state for Enceladus without prescribed mechanical asymmetry. Using the convection code StagYY in a 2D-spherical annulus geometry, we show that a non-Newtonian ice rheology is sufficient to create a localized, single hot plume surrounded by a conductive ice mantle. We obtain a self-sustained state in which a region of small angular extent has a sufficiently low viscosity to allow subcritical to weak convection to occur due to the stress-dependent part of the rheological law. We find that the single-plume state is very unlikely to remain stable if the rheology is Newtonian, confirming what has been found by previous studies. In a second set of numerical simulations, we also investigate the first-order effect of tidal heating on the stability of the single-plume state. Tidal heating reinforces the stability of the single-plume state if it is generated in the plume itself. Also we show that the likelihood of a stable single-plume state does not depend on the thickness of the ice shell.

A. Rozel, J. Besserer, G. J. Golabek, M. Kaplan and P. J. Tackley (2014). Self-consistent generation of single-plume state for Enceladus using non-Newtonian rheology. J. Geophys. Res. Planets 119, 416-439.


Oligarchic growth of Mars and its implications for Hf-W chronology:

Dauphas and Pourmand (2011) estimated the accretion timescale of Mars to be 0.8 - 2.7 Myr from the W isotopes of Martian meteorites. This timescale was derived assuming perfect metal-silicate equilibration between the impactor and the target's mantle. However, in the case of a small impactor most likely only a fraction of the target's mantle is involved in the equilibration, while only a small part of the impactor's core equilibrates in the case of a giant impact. We examined here the effects of imperfect equilibration using results of high-resolution N-body simulations for the oligarchic growth stage.

Morishima, R., G. J. Golabek and H. Samuel (2013). N-body simulations of oligarchic growth of Mars: Implications for Hf-W chronology. Earth Planet. Sci. Lett. 366, 6-16.


Early thermochemical evolution of asteroid 4 Vesta:

Diogenites are thought to represent mantle rocks formed as cumulates in magma chambers on 4 Vesta. We studied Northwest Africa (NWA) 5480, a harzburgitic diogenite, composed mainly of heterogeneously distributed olivine and orthopyroxene. Electron backscattered diffraction on the olivine grains of NWA 5480 shows that the preferred orientation can be explained neither by cumulate formation, nor by impact reprocessing near the asteroid's surface. Rather, they represent high-temperature solid-state plastic deformation by the pencil-glide slip system, which is well known from dry ultramafic rocks on Earth, typically formed by mantle shear at subsolidus temperatures. Numerical 2D models indicate that these observations may be explained by large-scale downwellings occurring in the asteroid's mantle within the first 50 Myr of Vesta's evolution. Implications include long-lasting enhanced mass exchange occurring in the dynamic interiors of differentiated asteroids.

Tkalcec, B. J, G. J. Golabek and F. E. Brenker (2013). Solid-state plastic deformation in the dynamic interior of a differentiated asteroid. Nature Geosci. 6, 93-97.

Animation:  Movie

Press coverage:  Los Angeles Times


Comparison of numerical surface topography calculations:

Topography is a direct observable of the interaction between the Earth's internal and external dynamics. Therefore, it is important for numerical models of lithospheric deformation to compute topography accurately. Earth's surface is a so-called free surface, which means that both normal and shear stress should vanish at this interface. It has been shown that correct treatment of the Earth's surface as a free surface can have a significant effect in models of lithospheric and mantle dynamics. However, a true free surface is computationally expensive which is why most mantle convection simulations until now treat the surface as a free-slip boundary. For these models, topography is computed directly from normal stresses. A free surface approximation, the so-called 'sticky air', has recently gained interest in the geodynamic community. This method requires the addition of a fluid layer in the model domain while retaining the computational advantage of a free-slip top boundary. The fluid layer is a proxy for air and should, therefore, have a near-zero density and a viscosity, which is several orders of magnitude lower than the lithosphere viscosity. Sufficiently small normal stress at the surface is ensured by the physical properties of the weak layer (i.e., the low values for density and viscosity). We present a theoretical background that provides the physical conditions under which the sticky-air approach is a valid approximation of a true free surface.

Crameri, F., H. Schmeling, G. J. Golabek, T. Duretz, R. Orendt, S. Buiter, D. A. May, B. J. P. Kaus, T. V. Gerya and P. J. Tackley (2012). A comparison of numerical surface topography calculations in geodynamic modelling: an evaluation of the 'sticky air' method. Geophys. J. Int. 189, 38-54.


Formation of the martian crustal dichotomy and long-lived volcanism:                                  

In this project we study the influence of a giant impactor on the evolution of Mars. The crustal dichotomy is the oldest and most striking surface feature on Mars. It is a large difference in elevation and crustal thickness between the southern highlands. It was formed more than 4.1 Ga ago owing to either exogenic or endogenic processes. Based on the geochemical analysis of SNC meteorites it was suggested that a primordial crust with up to 45 km thickness can be formed already during the martian core formation. As the final accretion stage of terrestrial planets is based on stochastically distributed impacts, we suggest that the sinking of the iron core, delivered by a late pre-differentiated giant impactor, induced shear heating-related temperature anomalies in the mantle, which fostered the formation of early asymmetrical primordial crust. In this study, we examine parameter sets that will likely cause an onset of hemispherical low-degree mantle convection directly after, and coupled to, an already asymmetrical core formation. For this purpose we couple 2D cylindrical simulations using the codes I2ELVIS for the core formation and the code StagYY for the long-term mantle convection in series. We apply a temperature, stress- and composition-dependent viscoplastic rheology inside a Mars-sized planet and include radioactive- and shear heating. Results show that low-degree convection can be induced by a giant impact during core formation. Furthermore, the amplitude of shear heating anomalies generally well exceeds the solidus of primitive mantle material. Therefore a hemispherical magma ocean is formed, which can be the source of a dichtomous crust. A degree-1 convective pattern establishes during the long-term evolution with the single upwelling feeding long-term volcanism as observed on Mars.

Golabek, G. J., T. Keller, T. V. Gerya, G. Zhu, P. J. Tackley & J. A. D. Connolly (2011). Origin of the martian dichotomy and Tharsis from a giant impact causing massive magmatism. Icarus 215, 346-357.

Animations:  Movie1   Movie2   Movie3

Influence of silicate rheology on core formation: