How and why is deformation localised in continental crust?

This Project has been filled

How and why is deformation localised in continental crust?

Sandra Piazolo, Tim Wright, Tim Craig

In this project the student will make and use geological field investigations of the exhumed lower crust to constrain models of strain localisation, and compare these against modern geodetic data from the India-Asia collision.

Over the last few years, geodetic observations of strain from the continents have revealed two modes of behaviour (Figure 1). At fault zones such as the Himalayan Frontal Thrust, Altyn Tagh and Xianshuihe Faults in the India-Eurasia collision zone high geodetic strain rates are observed at the faults themselves. In other regions, including Western Tibet and Mongolia, present-day deformation is diffuse, with little evidence for localisation on the major structures (e.g. Zheng et al., 2017). The underlying reasons for this bi-modal behaviour remain unclear and are highly debated. Understanding processes responsible for the presence or absence of strain localisation in the continents is important for diverse problems ranging from seismic hazard to the distribution of resources.

To address the question of strain localisation, the student will make new observations from exhumed lower-crustal rocks in Greenland, where bi-modal behaviour also appears to be present (Fig. 2), and use the geological and geodetic observations to assess proposed mechanisms for strain (de-)localisation in continental crust by running numerical models.

Figure 1: (a) Map of the present-day strain rate in the India-Eurasia collusion zone derived from 2576 GPS velocities (Zheng et al., 2017); (b) Annotated version of strain map, showing areas of strain localisation (red lines), undeforming regions (green blocks), and areas of diffuse strain (pink regions).


The Nagssugtoqidian Orogen in SW Greenland is a deeply eroded Paleoproterozoic, Himalaya sized orogen exposing mid to lower crustal rocks (e.g. Kolb 2014). It is characterised by geological “blocks” that differ in whether strain localization is present or absent. These exhibit striking similarities to areas of different behaviour seen in the present day India-Eurasia collision zone (Fig. 1 & 2).  Compilation of structural and metamorphic data from over 20 years of mapping in the area allow identification of areas of interest that represent the areas of strain localization (red lines), largely undeformed regions (green blocks) and areas of diffuse strain (pink regions) (Fig. 2). Using this unique dataset, freely available to us, as a guide, the student will conduct additional targeted field work to investigate the rheological, chemical and pre-deformation properties of the rocks outcropping today. They will use the results to build models that test whether the metamorphic, deformation and igneous history of areas involved in present day collisions is a deciding factor in the rheological behaviour of the whole orogen.








Figure 2 Characteristics of the Nagssugtoqidian orogeny, W Greenland. (a) General map of the Nagssugtoqidian orogen in Greenland (inset shows position within Greenland), (b) structural trends from compilation of structural data, (c) image of steep high strain zone (modified after van Gool & Piazolo, 2002).






The latest geodetic strain rate models for the India-Eurasia collision zone are built from more than 2500 individual GPS velocities (Zheng et al., 2017). Work by COMET scientists at Leeds ( is using InSAR observations to increase the spatial density of observations even further (e.g. Wang and Wright, 2012; Wang et al., 2019), but even from the GPS data alone it is clear that there are different regimes in operation within the India-Eurasia collision zone (Figure 1). The student would have the opportunity to refine and improve geodetic models for the India-Eurasia collision zone using the latest observations. We would expect the student to interrogate the data to critically assess the degree of strain localisation that can be observed and how it varies in space (laterally and vertically), potentially also utilizing data derived from ambient noise topography (e.g. Jiang et al. 2014)

For the geological part of the project, the student will interrogate the available data from geological mapping campaigns over the last 20 years available from the Geological Survey of Denmark and Greenland to identify field areas of interest. These areas will then be investigated in detail using latest techniques to infer rheological behaviour of different rock units, field analysis and microstructural work (e.g. Svahnberg & Piazolo, 2010) to derive physico-chemical properties through time.

This data will then be used to build numerical models to forward model expected geodetic strain rate distributions arising from such areas, using methods similar to those described in Yamasaki, Wright and Houseman (2014) and Hussain et al. (2018). Comparison between the two datasets will allow for model refinement and definition of the parameters governing strain localization or strain distribution.

Specifically, the student will:

  1. Analyse existing geological data from the Nagssugtoqidian belt of West Greenland to characterise strain localisation. (Yr1)
  2. Geologically map, collect and analyse new samples from key sites to estimate viscosity distributions (in space) at locations with/without strain localisation and to establish plausible hypotheses for strain localisation. (Yr1&Yr2)
  3. Quantify the degree of strain localisation in the India-Eurasia collision zone using the latest geodetic data from InSAR and GPS, and assess how it varies in space. (Yr1&Yr2)
  4. Build geodynamic models to test the geological hypotheses against the geodetic observations and predictions built from the Greenland case study area. (Yr3)

The student will be supported by the supervisory team whose expertise covers all aspects of the projects. Depending of progression of the project, a visit to a second field area in Central Australia representing an ancient intercontinental orogen (Raimondo et al. 2014) may be organised.


This project is highly innovative as it brings together fields of research that are traditionally separated, even though their main aim is the same: Understanding of the rheological behaviour of the Earth’s continental crust. At Leeds we can offer this unique project as expertise exists in all necessary areas.  This project promises to produce cutting-edge research markedly advancing our knowledge and understanding of Earth’s deformation behaviour and the validity of classic “rigid” plate tectonic theory in the continents. This knowledge is urgently needed to be able to refine current models used to predict geohazards such as earthquake frequency. We expect research to be presented at international conferences and published in international journals.

Student profile:
The project would suit a numerate student with a background in earth sciences, geology, or geophysics who is enthusiastic about geological fieldwork in remote and challenging areas. The student will be provided with training in state of the art geological and geodetic methods, including microstructural and large geological data analysis. The student will be part of the UK Natural Environmental Research Council’s Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET) and would have the opportunity to interact with COMET students with different skills and backgrounds from across the UK.

In addition, the student will have access to a broad spectrum of training workshops offered in house e.g. image analysis,  presentation skills, through to ‘managing your degree’ and ‘preparing for your viva’ (

References / Further Reading

Hussain, E., Wright, T.J., Walters, R.J., Bekaert, D.P., Lloyd, R. and Hooper, A., 2018. Constant strain accumulation rate between major earthquakes on the North Anatolian Fault. Nature communications, 9(1), p.1392.

Jiang, C., Yang, Y, and Zheng, Y, 2014. Penetration of mid-crustal low velocity zone across the Kunlun Fault in the NE Tibetan Plateau revealed by ambient noise tomography, Earth and Planetary Science 406, 81-92.

Kolb, J. (2014). Structure of the Palaeoproterozoic Nagssugtoqidian Orogen, South-East Greenland: model for the tectonic evolution.Precambrian Research,255, 809-822.

Raimondo, T., Hand, M., & Collins, W. J. (2014). Compressional intracontinental orogens: Ancient and modern perspectives.Earth-Science Reviews,130, 128-153.

Svahnberg H; Piazolo S (2010) The initiation of strain localisation in plagioclase-rich rocks: Insights from detailed microstructural analyses,Journal of Structural Geology,32, pp.1404-1416

van Gool JAM; Piazolo S (2006) Presentation and interpretation of structural data from the Nagssugtoqidian orogen using a GIS platform: General trends and features,Geological Survey of Denmark and Greenland Bulletin, pp.125-144.

Wang H; Wright TJ (2012) Satellite geodetic imaging reveals internal deformation of western Tibet, Geophysical Research Letters, 39, . doi: 10.1029/2012GL051222

Wang, H., Wright, T.J., Liu‐Zeng, J. and Peng, L., 2019. Strain Rate Distribution in South‐Central Tibet From Two Decades of InSAR and GPS. Geophysical Research Letters, 46(10), pp.5170-5179.

Yamasaki, T., Wright, T.J. and Houseman, G.A. (2014). Weak ductile shear zone beneath a major strike‐slip fault: Inferences from earthquake cycle model constrained by geodetic observations of the western North Anatolian Fault Zone. Journal of Geophysical Research: Solid Earth119(4), pp.3678-3699.

Zheng G; Wang H; Wright TJ; Lou Y; Zhang R; Zhang W; Shi C; Huang J; Wei N (2017) Crustal deformation in the India-Eurasia collision zone from 25 years of GPS measurements, Journal of Geophysical Research: Solid Earth, . doi: 10.1002/2017JB014465