Associate Professor Rhodri Davies

Associate Professor

Positions Held

2018 - Present: Associate Professor & Associate Director (Honours & Masters), RSES, ANU.

2014 - 2018: ARC Future Fellow, RSES, ANU.

2013 - 2014: Research Fellow, RSES, ANU.

2011 - 2013: NERC Research Fellow, Imperial College London.

2009 - 2010: 1851 Research Fellow, Imperial College London.

2008 - 2009: Research Assistant, Funded by Shell, Based at the School of Earth & Ocean Sciences, Cardiff University.


2018: Anton Hales Medal, The Australian Academy of Sciences. Recognises distinguished research in the Earth sciences,

2015: Media & Outreach Award, The Australian National University.

2014: Outstanding Young Scientist Award: European Geosciences Union, Geodynamics Division. Recognises substantial, original and significant contributions to global geodynamical research.

2014: Strategic Communications & Public Affairs Award, The Australian National University.

Funding and Fellowships

2021: ARDC Platform Grant - G-Adopt - Geodynamic ADjoint OPTimization Platform.

2020: ARC Discovery Project (Valentine, Sambridge and Davies) - Seeing the Unseeable: A New Generation of Geophysical Imaging.

2018: Geoscience Australia Research Partnership - Exploring for the Future Program - Probing the Lithosphere using Mafic Geochemistry.

2017: ARC Discovery Project (Davies, Campbell and Kennett) - Earth's Intra-Plate Volcanic Engine.

2014: ARC Future Fellowship (Davies) - From Plume Source to Hotspot: Quantifying Mixing in Mantle Plumes and its Implications for the Nature of Deep Mantle Heterogeneity.

2012: NERC Standard Grant (UK) - The Mantle Transition Zone `Valve'.

2011: NERC Post-doctoral Research Fellowship (UK)  - Pulsing Mantle Plumes: Causes and Geological Consequences.

2009: Research Fellowship - Royal Commission for the Great Exhibition of 1851 (UK) - Modelling Earth's Engine: Innovative Techniques for Simulating Thermo-chemical Mantle Convection.

2009: Research Grant - Shell International Exploration and Production (UK) - Integrating Convection Models, Tomographic Observations and Tectonic Reconstructions.

Research interests

I am recognised for developing and integrating state-of-the-art tools that simulate mantle and lithosphere dynamics, with a diverse range of observational datasets, to significantly enhance understanding of mantle and lithosphere dynamics and their signature at the surface, across a range of spatial and temporal scales.

Research overview

My research focusses upon the dynamics of our planetary heat engine – the mantle – and its geophysical, geological and geochemical expression, across a range of spatial and temporal scales. I have developed some of the most advanced tools available for simulating geodynamical processes and have used these tools, alongside a variety of novel observational datasets, to significantly enhance our understanding of mantle and lithosphere dynamics and its signature at the surface.

My long-term vision is to develop a new paradigm for solid Earth evolution. The overarching goal is a reconstruction of the thermo-chemical structure and flow history of Earth’s mantle and its evolution through space and time (4-D). To achieve this, I am working towards developing ‘inverse geodynamics’, to fuse wide-ranging constraints from geophysics, geology and geochemistry with multi-resolution computational modelling. The key outcome will be a unified quantitative understanding of plate tectonics and the thermo-chemical state of Earth’s mantle in 4-D. This will reveal the interplay between the deep Earth, the lithosphere and its surface, and the physical and chemical processes that govern the tempo and multi-scale nature of Earth’s evolution, addressing a widely recognised grand challenge for the Earth Sciences. Success will require a collaborative network of scientists from several fields, including numerical modellers, inverse theorists and observational scientists. It will engage a broad cross-section of academic staff and graduate students. I am always looking for motivated and enthusiastic students to join the group.

Key focus areas of my current and future research, include:

  1. Isolating the mechanisms underpinning intra-plate volcanism on Earth, by combining cutting-edge geodynamical models with seismological, geochemical and petrological constraints – this research has already produced game-changing results, published in Geology and Nature (2015, 2017), the significance of which was formally recognised through the 2018 Anton Hales Medal from the Australian Academy of Sciences.
  2. Revealing how Earth’s surface topography responds to underlying mantle flow (dynamic topography), by reconciling observational constraints on residual topography with theoretical predictions – our recent research in this space (Nature Geosciences, 2019), provides fundamental insight into present-day dynamic topography that will underpin our ongoing efforts to reconstruct its spatio-temporal evolution.
  3. Quantifying the force-balance governing surface tectonics and rapid plate motion changes, utilising our unique coupled models of mantle and lithosphere dynamics – although it is often overlooked, our ground-breaking research in this area, which includes a 2018 publication in Nature Geosciences, highlights the key role of mantle flow in governing surface tectonics.
  4. Reconciling key geophysical and geochemical constraints on mantle structure, by validating synthetic models of mantle structure, dynamics and evolution, against diverse geophysical and geochemical datasets – this builds on almost a decade of theoretical, methodological and applied research, the novelty and significance of which was emphasized in the citation for my 2014 outstanding young scientist award from the EGU.

Fundamental progress on each of these goals relies on the development and application of `inverse geodynamics’, to fuse wide-ranging constraints from geophysics, geology and geochemistry with multi-resolution computational modelling.

Research highlights

Davies et al., Nat. Geo. (2019): we combine observational constraints on residual topography beneath the world’s oceans with an innovative statistical approach to performing spherical harmonic analyses to generate a robust estimate of the power spectrum of Earth’s oceanic residual topography. We demonstrate that theoretical predictions can only be reconciled with the resulting spectrum if they incorporate the effects of lithospheric structure and its contribution towards global mantle flow: both deep (long-) and shallow (shorter-wavelength) processes are crucial to generating the observed surface response. Our study resolves the long-standing debate on the spatial distribution, wavelength and amplitude of present-day dynamic topography, reconciling key predictions of mantle flow models with inferences from the geological record and emphasising the key role that Earth’s deep interior plays in shaping the surface of our planet. This will underpin our ongoing efforts to reconstruct its spatial and temporal evolution.

Inferred Oceanic Residual Topography: spherical harmonic model (regularised using ARD) up to l = 30, of our updated compilation of oceanic (point-wise and ship-track) residual topography measurements. In panel (b) we display the l = 0-3 components of this model. Contours are at 0.2 km intervals, with dashed contours denoting negative values and grey indicating the zero contour.

Iaffaldano et al., Nat. Geo. (2018): plate motions are well-studied but we lack a fundamental understanding of their impact on surface processes, particularly in regions where their boundaries are relatively diffuse. Nowhere is this more clearly apparent than the Indian Ocean floor. In this study, we identify geologically recent plate motion changes in the Indian Ocean realm and provide an explanation for their occurrence and an understanding of their impact on the development of diffuse plate margins. Specifically, we link the observed plate motion change and the consequent lithospheric deformation to the Reunion mantle plume. Our results emphasise the key role that Earth’s deep interior plays in shaping the surface of our planet.

Jones et al., Nature (2017): an inter-disciplinary study, which demonstrates that the Pacific plate’s geographically and geochemically distinct double volcanic hotspot tracks, at Hawaii, Samoa, Society, Foundation and Marquesas, emerged concurrently at ~3 Ma. The study provides a novel dynamical explanation that links their emergence to a recent and previously unidentified plate motion change, and illustrates how this change produced the systematic geochemical variations observed.

Left: Bathymetric map of recent Hawaiian volcanism, highlighting the Loa and Kea tracks; Right: Schematic diagram of the tilted Hawaiian plume, the overlying Pacific plate and associated surface volcanism.

Jones et al., EPSL (2016): we quantitatively demonstrate how deep-mantle heterogeneities are transported into a plume conduit, and determine the conditions under which deep-mantle structure can be mapped from geochemical observations at Earth’s surface.

Davies et al., Nature (2015): we identify the full-extent of Earth’s longest continental hotspot track: the Cosgrove track in eastern Australia. This study highlights a clear correlation between lithospheric thickness, volcanic outcrop and magma composition along this track, indicating that lithospheric thickness has a dominant influence on the volume and chemical composition of plume derived magmas and placing an observational constraint on the sub continental melting depth of upwelling mantle material.

In Davies et al. (2015), we identify Earth’s longest continental hotspot track in eastern Australia, the Cosgrove track, which extends from Cape Hillsborough (CH) to Cosgrove (C). Magma composition and volcanic outcrop correlate with lithospheric thickness along this track, thus constraining the sub-continental melting depth of mantle plumes. This figure shows the distribution of eastern Australian Cenozoic volcanic centers and their relationship to lithospheric thickness: (a) volcanic centers, colored using the classification of Wellman & McDougall (1974), with age-progressive hotspot-tracks denoted by black dashed lines; BU=Bunya; CA=Canabolas; CE=Central; CO=Comboyne; FI=Fraser Island; L=Liverpool; M=Monaro; MC=McBride; NU=Nulla; NVP=Newer Volcanics Province; ­­S=Sturgeon; (b) the same centers, plotted on an estimate of lithospheric thickness, highlighting a correlation between lithospheric thickness and volcanic outcrop and classific­­ation along the Cosgrove hotspot-track.

Davies & Rawlinson, Geology (2014): an inter-disciplinary study which combines high-resolution data from the transportable WOMBAT seismic array was integrated with 3-D geodynamical models, to demonstrate that several lithospheric steps combine, because of their orientation, to localize mantle melting beneath the Newer Volcanics Province of Victoria and South Australia. Although focussed on southern Australia, the study helps to solve the global puzzle of why step-changes in lithospheric thickness, which occur along craton edges and at passive margins, only produce volcanism at isolated locations: lithospheric steps are a key prerequisite for this process, but their orientation is of equal importance, due to the role of background mantle flow.

Variations in P-wave velocity at 100 km depth, beneath the Newer Volcanics region of Southeastern Australia. The horizontal limits of volcanic outcrop at the surface is denoted by a yellow dashed line. This image comes from the inter-disciplinary seismological and geodynamical study of Davies & Rawlinson (2014), which demonstrates that recent volcanism in Australia is the manifestation of edge-driven convection. It also provides a potential solution to the global puzzle of why step-changes in lithospheric thickness produce only isolated volcanism.

Garel et al., G3 (2014): competition between overriding plate strength and tensile stresses in the subducting plate controls both trench migration rates and slab sinking rates, with important consequences for slab/transition-zone interaction.

A 2-D thermo-mechanical dynamic subduction simulation from Fluidity. Temporal snapshots of: (a) temperature; (b) viscosity; (c) the dominant deformation mechanism; and (d) the underlying computational mesh, from a case where subducting/overriding plates are initially 100/20 Myr old at the trench, respectively (black squares indicate initial trench location). White lines in panel (a) mark isotherms from 600-1400 K at 200 K intervals. Black (b/c) and red (d) lines mark the location of the 1300K isotherm. In this class of subduction model, the slab’s buoyancy and distinct rheological properties arise self-consistently, through variations in temperature, pressure and strain-rate, with deformation accommodated through a composite diffusion creep (diff), dislocation creep (disl), Peierls creep (P) and yielding (YS) law. In the example shown, the slab’s excess density drives subduction and trench retreat over time. Upon interaction with the transition zone, the descending slab temporarily stalls and deforms (b/c), before slowly sinking into the lower mantle (d). Note that throughout its descent, the slab maintains a strong core (b), which is a requirement from observations of Benioff seismicity. The underlying computational mesh is adapted at fixed intervals during the simulation, with zones of high resolution analogous to regions of dynamic significance. A local resolution of 500m is required to resolve the slab’s strong core and the interface between subducting and overriding plates (Garel et al. 2014). The Garel et al. (2014) study, funded by a NERC grant, shows how competition between overriding plate strength and tensile stresses in the subducting plate controls both trench migration and slab sinking rates, with consequences for slab/transition-zone interaction.

Davies et al., EPSL (2012): synthetic seismic structures from global circulation models demonstrate that thermal heterogeneity alone can explain observed deep mantle seismic characteristics. This has important implications for the nature and distribution of mantle chemical heterogeneity.

Shear wave velocity perturbations beneath Africa from: (a) tomographic model S40RTS (Ritsema et al. 2011); (b) a purely-thermal model; and (c) a thermo-chemical model, where the dense chemical component represents only 3% of the mantle volume (Davies et al. 2012). In the isochemical model, (T,P,X) is converted into seismic velocity assuming a pyrolitic composition, whilst  the thermo-chemical model illustrated assumes a pyrolitic composition for background mantle and an iron-rich composition for the dense chemical component. The structure of a thermo-chemical model with a basaltic dense component has similar characteristics. We account for the geographic bias, smearing and damping inherent to tomographic models using the resolution operator of S40RTS, thus allowing for direct comparison between our models and S40RTS. All images include a radial surface at 2850-km depth and an isosurface at -1.0%, clipped above 1400-km to allow visualisation of lower mantle features. Continental boundaries are included for geographic reference. The Davies et al. (2012) study, fudned by a NERC post-doctoral Research Fellowship, demonstrates that thermal heterogeneity alone can explain observed deep mantle seismic characteristics. This has important implications for the nature of mantle dynamics and the likely distribution of mantle chemical heterogeneity.

Davies et al., G3 (2011): a new computational framework for geodynamical simulation, built around parallel anisotropic adaptive unstructured mesh techniques. It allows for simulations of mantle and lithospheric dynamics with an unprecedented level of realism.

(a) Anisotropic, adapted unstructured mesh (maximum element aspect ratio = 10, number of triangle vertices = 16,699), from a square convection simulation shown in Fluidity; (b, c) high‐resolution images of specific mesh sections, illustrating anisotropic elements within the thermal boundary layer; (d) a near‐isotropic unstructured mesh (maximum element aspect ratio = 1, number of triangle vertices = 94,133), with the same minimum edge length as the mesh displayed in panel a; (e, f) high‐resolution images of specific mesh sections. Placing an aspect ratio bound of 1 (isotropic) on the metric dramatically increases the number of nodes required to represent the problem.

Davies & Davies, Solid Earth (2010): a vastly improved estimate of Earth’s surface heat flux is presented, based upon GIS methodologies and an enhanced heat flow dataset. This has implications for Earth’s thermal budget and evolution.

Davies & Davies, EPSL (2009): cutting-edge global mantle convection models demonstrate that mantle plumes can reconcile a wide–range of hotspot observations, providing strong quantitative support for the mantle plume hypothesis.

A series of thermal profiles, displaying the temporal evolution of a high Rayleigh number global mantle convection simulation, with free-slip and isothermal boundaries. Snapshots are spaced 30 Myr apart and are surrounded by a latitude-longitude grid. The scale illustrates the temperature, away from the lateral average. Each snapshot shows a radial surface just above the CMB and a hot isosurface, representing regions of the mantle that are 500K hotter than average for their depth. The most prominent features are hot upwelling plumes. These display a range of characteristics. The majority are long-lived and migrate slowly across the mantle. Others are more mobile and ephemeral. New plumes form and old ones fade, some of which eventually cease. Smaller plumes merge together and coalesce over time, while long-lived plumes often pulse or vary in intensity. This snapshot comes from the study of Davies & Davies (2009), which demonstrates that mantle plumes can reconcile a wide-range of hotspot observations, providing strong quantitative support for the mantle plume hypothesis.

Davies et al., G3 (2007): first study to demonstrate applicability and benefits of adaptive mesh refinement techniques for geodynamical simulations. These schemes use intelligent algorithms to modify the computational grid automatically, enhancing resolution where required.


Link to ORCID profile

Link to Google Scholar profile

Link to Scopus profile

44. A. Beall, A. Fagerang, J. H. Davies, F. Garel and D. R. Davies. Influence of Subduction Zone Dynamics on Interface Shear Stress and Potential Relationship with Seismogenic Behavior. Geochem. Geophys. Geosys. 21, doi:10.1029/2020GC009267, In Press, 2021.

43. L. Suchoy, S. Goes, B. Maunder, F. Garel and D. R. Davies. Effects of basal drag on subduction dynamics from 2D numerical models. Solid Earth. 12, 79-93, doi:10.5194/se-12-79-2021, 2021

42. A. P. Valentine and D. R. Davies. Global models from sparse data: a robust estimate of Earth's residual topography spectrum. Geochem. Geophys. Geosys. 21, doi:10.1029/2020GC009240, 2020.

41. M. Klöcking, K. Czarnota, D.C. Champion, A.L. Jaques and D. R. Davies. Spatio-temporal evolution of Australian lithosphere-asthenosphere boundary from mafic volcanism. Exploring for the Future: extended abstracts, Geoscience Australia, Canberra, 1-4, 2020.

40. F. Garel, C. Thoraval, A. Tommasi. S. Demouchy and D. R. Davies. Using thermo-mechanical models of subduction to constrain effective mantle viscosity, Earth Planet. Sci. Lett.  539, doi: 10.1016/j.epsl.2020.116243, 2020.

39. Alsaif, M., Garel, F., Gueydan, F. & Davies, D. R. Upper plate deformation and trench-retreat modulated by subduction-driven shallow asthenospheric flows, Earth Planet. Sci. Lett. 532, doi: 10.1016/j.epsl.2019.116013, 2020.

38. Kennett, B. L. N. & Davies, D. R. Intra-plate Volcanism in North Queensland and eastern New Guinea: A Cryptic Mantle Plume? Gondwana Res. 79, 209-216, doi: 10.1016/, 2020.

37. Davies, D. R., Valentine, A. P., Kramer, S. C., Rawlinson, N., Hoggard, M. J., Eakin, C. & Wilson, C. R. Earth's multi-scale topographic response to global mantle flow. Nature Geosci. 12, 845-850, doi: 10.1038/s41561-019-0441-4, 2019.

36. Nebel, O., Sossi, P. A., Benard, A., Arculus, R. J., Yaxley, G., Woodhead, J. D., Davies, D. R., Ruttor, S. Reconciling petrological and isotopic mixing mechanisms in the Pitcairn mantle plume using stable Fe isotopes. Earth Planet. Sci. Lett., 521, 60-67, doi:10.1016/j.epsl.2019.05.037, 2019.

35. Jones, T. J., Davies, D. R., Sossi, P. A. Tungsten isotopes in mantle plumes: Heads its positive, Tails its negative. Earth Planet. Sci. Lett., 506,255-267, doi:10.1016/j.epsl.2018.11.008, 2019.

34. Ernst, R. E., Davies, D. R., Jowitt, S. M. & Campbell, I. When do mantle plumes destroy diamonds? Earth Planet. Sci. Lett. 502, 244-252, doi:10.1016/j.epsl.2018.08.058, 2018.

33. Perrin, A., Goes, S., Prytulak, J., Rondenay, S. & Davies, D. R. Mantle wedge temperatures and their potential relation to volcanic arc location. Earth Planet. Sci. Lett. 501, 67-77, doi:10.1016/j.epsl.2018.08.011, 2018.

32. Koelemeijer, P., Schuberth, B. S. A., Davies, D. R., Deuss, A. & Ritsema, J. Constraints on the presence of post-perovskite in Earth's lowermost mantle from tomographic-geodynamic model comparisons. Earth Planet. Sci. Lett. 494, 226-238, doi: 10.1016/j.epsl.2018.04.056, 2018.

31. Iaffaldano, G. & Davies, D. R. & DeMets, C. Indian Ocean floor deformation induced by the Reunion plume rather than the Tibetan Plateau. Nature Geosciences. doi:10.1038/s41561-018-0110-z, 2018.

30. Stotz, I. L., Iaffaldano, G. & Davies, D. R. Pressure driven Poiseuille flow: a major component of the torque-balance governing Pacific Plate motion. Geophys. Ref. Lett. 45, 117-125, doi:10.1002/2017GL075697, 2018.

29. Rawlinson, N., Davies, D. R. & S. Pilia. The mechanisms underpinning Cenozoic intraplate volcanism in eastern Australia: insights from seismic tomography and geodynamic modelling. Geophys. Res. Lett. doi:10.1002/2017GL074911, 2017.

28. Stotz, I. L., Iaffaldano, G. & Davies, D. R. Late Miocene Pacific Plate kinematic change explained with coupled global models of mantle and lithosphere dynamics. Geophys. Res. Lett. 44, doi:10.1002/2017GL073920, 2017.

27. Rubey, M., Brune, S., Heine, C., Davies, D. R., Williams, S. E. & Muller, R. D. Global patterns of Earth's dynamic topography since the Jurassic: the role of subducted slabs. Solid Earth, 8, 899-919, doi:10.5194/se-8-899-2017, 2017.

26. Jones, T. J., Davies, D. R., Campbell, I. H., Iaffaldano, G., Yaxley, G., Kramer, S. C. & Wilson, C. R. The concurrent emergence and causes of double volcanic hotspot tracks on the Pacific plate. Nature, 545, 472-476, doi:10.1038/nature22054, 2017.

25. Campbell, I. H. & Davies, D. R. Raising the continental crust. Earth Planet. Sci. Lett., 460, 112-122, doi:10.1016/j.epsl.2016.12.011, 2017.

24. Hedjazian, N., Garel, F., Davies, D. R. & Kaminski, E. Age-independent seismic anisotropy under oceanic plates explained by strain history in the asthenosphere. Earth Planet. Sci. Lett., 460, 135-142, doi:10.1016/j.epsl.2016.12.004, 2017.

23. Perrin, A., Goes, S., Prytulak, J., Davies, D. R., Kramer, S. C. & Wilson, C. R. Reconciling mantle wedge thermal structure with arc lava thermobarometric determinations in oceanic subduction zones, Geochem. Geophys. Geosys. 17, 4105-4127, doi:10.1002/2016GC006527, 2016.

22. Sossi, P. A., Eggins, S. M., Nesbitt, R. W., Nebel, O., Hergt, J. M., Campbell, I. H., O'Neill, H. C., van Kranendonk, M. & Davies, D. R. Petrogenesis and Geochemistry of Archean Komatiites. J. Petrol. doi:10.1093/petrology/egw004, 2016.

21. Jones, T. J., Davies, D. R., Campbell, I. H., Wilson, C. R. & Kramer, S. C. Do mantle plumes preserve the heterogeneous structure of their deep-mantle source? Earth. Planet. Sci. Lett. 434, 10-17, doi:10.1016/j.epsl.2015.11.016, 2016.

20. Davies, D. R., Le Voci, G., Goes, S., Kramer, S. C. & Wilson, C. R. The mantle wedge's transient 3-D flow regime and thermal structure. Geochem. Geophys. Geosys. 17, 1-23, doi:10.1002/2015GC006125, 2016.

19. Tkalcic, H., Young, M., Muir, J. B., Davies, D. R. & Mattesini, M. Strong multi-scale heterogeneity in Earth's lowermost mantle. Sci. Rep. 5, 18416, doi:10.1038/srep18416, 2015.

18. Davies, D. R., Rawlinson, N., Iaffaldano, G. & Campbell, I. H. Lithospheric controls on magma composition along Earth's longest continental hotspot track. Nature, 525, 511–514, doi:10.1038/nature14903, 2015.

17. Tosi, N., Stein, C., Noack, L., Huttig, C., Maierova, P., Samuel, H., Davies, D. R., Wilson, C. R., Kramer, S. C., Thieulot, C., Glerum, A., Fraters, M., Spakman, W., Rozel, A., Tackley, P. J. A community benchmark for viscoplastic thermal convection in a 2-D square box. Geochem. Geophys. Geosys., 16, 2175–2196, doi:10.1002/2015GC005807, 2015.

16. Davies, D. R., Goes, S. & Lau, H. C. P. Thermally dominated deep mantle LLSVPs: a review. In: The Earth's heterogeneous mantle, Eds. A. Khan, F. Deschamps & K. Kawai. Springer, 441-478, doi:10.1007/978-3-319-15627-9_14, 2015.

15. Davies, D. R., Goes, S. & Sambridge, M. On the relationship between volcanic hotspot locations, the reconstructed eruption sites of LIPs and deep mantle seismic structure. Earth Planet. Sci. Lett. 411, 121-131, doi:10.1016/j.epsl.2014.11.052, 2015.

14. Davies, D. R. & Rawlinson, N. On the origin of recent intra-plate volcanism in Australia. Geology, 42, 1031-1034, doi: 10.1130/G36093.1, 2014.

13. Garel, F., Goes, S., Davies, D. R., Davies, J. H., Kramer, S. C. & Wilson, C. R. Interaction of subducted slabs with the mantle transition zone: A regime diagram from 2-D thermo-mechanical models with a mobile trench and an overriding plate. Geochem. Geophys. Geosys. 15, 1739-1765, doi:10.1002/2014GC005257, 2014.

12. Le Voci, G., Davies, D. R., Goes, S., Kramer, S. C. & Wilson, C. R. A systematic 2-D investigation into the mantle wedge's transient flow regime and thermal structure: Complexities arising from a hydrated rheology and thermal buoyancy. Geochem. Geophys. Geosys. 15, 28-51, doi:10.1002/2013GC005022, 2014.

11. Davies, D. R., Davies, J. H., Bollada, P. C., Morgan, K., Hassan, O. & Nithiarasu, P. A hierarchical mesh refinement technique for global 3-D spherical mantle convection modelling. Geosci. Mod. Dev. 6, 1095-1107, doi:10.5194/gmd-6-1095-2013, 2013.

10. Davies, D. R., Goes, S., Davies, J. H., Schuberth, B. S. A., Bunge, H.-P. & Ritsema, J. Reconciling dynamic and seismic models of Earth's lower mantle: the dominant role of thermal heterogeneity. Earth Planet Sci. Lett. 353-354, 253-269, doi:10.1016/j.epsl.2012.08.016, 2012.

9. Kramer, S. C., Wilson, C. R. & Davies, D. R. An implicit free-surface algorithm for geodynamical simulations. Phys. Earth Planet Int. 194-195, 25-37, doi:10.1016/j.pepi.2012.01.001, 2012.

8.  Hunt, S. A., Davies, D. R., Walker, A. M., McCormack, R. J., Wills, A. S., Dobson, D. P. & Li, L. On the increase in thermal diffusivity caused by the perovskite to post-perovskite phase transition and its implications for mantle dynamics. Earth Planet. Sci. Lett., 319-320, 96-103, doi:10.1016/j.epsl.2011.12.009, 2012.

7. Davies, D. R., Wilson, C. R. & Kramer, S. C. Fluidity: A fully unstructured anisotropic adaptive mesh computational modeling framework for geodynamics. Geochem. Geophys. Geosys., 12, Q06001, 20 PP., doi:10.1029/2011GC003551, 2011.

6. Styles, E., Davies, D. R., Goes, S. Mapping spherical seismic into physical structure: biases from 3-D phase-transition and thermal boundary-layer heterogeneity. Geophys. J. Int., 184, 1371-1378, doi: 10.1111/j.1365-246X.2010.04914.x, 2011.

5. Davies, J. H & Davies, D. R. Earth's surface heat flux. Solid Earth,  1, 5 - 24, 2010, doi:10.5194/se-1-5-2010

4. Wolstencroft, M., Davies, J. H. & Davies, D. R. Nusselt-Rayleigh number scaling for spherical shell Earth mantle simulations up to a Rayleigh number of 10^9. Phys. Earth Planet. Int., 176, 132-141, doi:10.1016/j.pepi.2009.05.002, 2009.

3.  Davies, D. R. & Davies, J. H. Thermally-driven mantle plumes reconcile multiple hot-spot observations. Earth Planet. Sci. Lett., 278, 50-54, doi:10.1016/j.epsl.2008.11.027, 2009.

2.  Davies, D. R., Davies, J. H., Hassan, O., Morgan, K. & Nithiarasu, P. Adaptive finite element methods in geodynamics; Convection dominated mid-ocean ridge and subduction zone simulations. Int. J. Num. Meth. Heat Fluid Flow, 7-8, 1015-1035, doi:10.1108/09615530810899079, 2008.

1.  Davies, D. R., Davies, J. H., Hassan, O., Morgan, K., & Nithiarasu, P. Investigations into the applicability of adaptive finite element methods for infinite Prandtl number thermal and thermo-chemical convection. Geochem. Geophys. Geosyst., 8, Q05010, doi:10.1029/2006GC001470, 2007.


EMSC 3034/6034: Dynamic Earth -- Plates, Plumes and Mantle Convection 

Mantle convection is the fundamental agent driving many of the geological features observed at Earth's surface, including plate tectonics and volcanism. However, many geologists have an incomplete understanding of the process, whilst there are many misconceptions about how it relates to surface processes. A broad background to the physics and fluid dynamics of mantle convection will be provided in this course, by explaining what it is, how it works, and how to quantify it in simple terms. It assumes no specialist background: mechanisms will be explained simply and the required basic physics will be fully reviewed and explained. The distinctive forms that convection takes within Earth's mantle will be described within the context of tectonic plates and mantle plumes, whilst the implications for geochemistry and Earth's tectonic evolution will be explored. Common misconceptions and controversies will be addressed, providing a straightforward, but rigorous, explanation of this key process. Emerging insights into the fundamental links between climate, surface processes (weathering, erosion and sediment transport), plate tectonics and underlying mantle flow will also be covered, thus providing a complete overview of Earth’s dynamic engine.

EMSC 8034: Research Orientation - Big Questions in the Earth Sciences

A course for Master of Earth Sciences students that will reveal the biggest unanswered questions in the Earth Sciences. Students will be introduced to the breadth of science undertaken at the Research School of Earth Sciences and will gain an appreciation for the different research approaches utilised and facilities available. Further, this course will provide a framework for students to make an informed choice of research supervisor and project for their subsequent studies.


2016: Fellow of the Higher Education Academy (FHEA). Recognises attainment against the UK professional standards framework for teaching and learning support in higher education.

Supervised students