Between Earth's crust and core lies the mantle, a 2,900 km-thick layer of hot rock that comprises greater than 80% of Earth's volume. Carrying Earth's internal heat to the surface, the convecting mantle creeps like tar on a hot day. This overturning is the 'engine' that drives our dynamic Earth; ultimately, all large- scale geological activity is driven by mantle convection (Davies, 1992). Quantitative modelling of this process is thus essential for understanding Earth's dynamics, structure, and evolution, from earthquakes and volcanoes to the forces that build mountains and break continents apart.
While sluggish in human terms, mantle convection is actually an extremely vigorous process. Indeed, the mantle's Rayleigh number, a dimensionless parameter quantifying the vigour of convection, is of order 109, which is almost that of a pot of boiling water. Accurate modelling of mantle convection is thus extremely challenging, because kilometre-scale features, such as subducting plates/slabs and upwelling mantle plumes, must be resolved within the thousand-kilometre-scale mantle. Consequently, despite the large computing clusters available today, with current methods it is difficult to accurately simulate the whole mantle at its true dynamical regime.
To overcome this obstacle to a better understanding of mantle convection on Earth and other planets, the proposed research will further enhance and test a recently developed computational framework for mantle dynamics, Fluidity (Davies et al. 2011, Kramer et al. 2012). The code:
- Employs an unstructured mesh with non-uniform resolution, so that regions of dynamical activity, or specific interest, can be modeled at a higher resolution than is required by the large-scale flow.
- Employs state-of-the-art adaptive mesh methods, to automatically resolve fine-scale features as they develop. This minimizes computational cost by reducing resolution where possible; computational power is automatically concentrated in the most active regions of the domain (Davies et al. 2007).
- It takes advantage of load-balanced domain decomposition algorithms, in order to run on parallel machines with distributed memory.
The modificatios central to the proposed project include:
- Incorporation of 500 Myr of plate motion histories as a surface boundary condition.
- Incorporation of complex rheologies (grain-size and histroy dependent - damage), such that 'plates' can be generated and analysed self-consistently within a global framework.
The modified code will subsequently be applied in examining a range of processes relating to mantle convection and plate tectonics, focusing upon mantle/lithosphere interaction, including: (i) how surface plates modulate underlying mantle flow; (ii) how plate tectonics has evolved over Earth's history; (iii) why plate tectonics is unique to Earth, among the terrestrial planets.
The successful applicant will hold a 1st class degree (or equivalent) and have a strong record in a numerate subject such as mathematics, computer sciences, physics, geophysics or engineering. Experience in computer programming is desirable. The student will join a large and active international group, with collaborators at Imperial College London, Columbia University, New York and the University of Sydney.