Anne Glerum

Anne Glerum is a postdoctoral fellow in the Geodynamic Modelling section at GFZ Potsdam. Through numerical modeling, she investigates continental rift dynamics, focusing on the effect of along-rift variations as well as magma-tectonic feedback processes on rift evolution. Anne is part of the GD blog team as an Editor. You can reach Anne via email.

Remarkable Regions – The Kenya Rift

Remarkable Regions – The Kenya Rift

Every 8 weeks we turn our attention to a Remarkable Region that deserves a spot in the scientific limelight. After looking at several convergent plate boundaries, this week the focus lies on part of a nascent divergent plate boundary: the Kenya Rift. The post is by postdoctoral researcher Anne Glerum of GFZ Potsdam.

Of course an active continental rift is worthy of the title “Remarkable Region”. And naturally I consider my own research area highly interesting. But after seeing it up-close and personal on a recent 10-day trip organized by the University of Potsdam, Roma Tre and the University of Nairobi (stay tuned for the travel log, or read that of the University of Potsdam), I must say, the Kenya Rift is a truly beautiful and fascinating region.

Figure 1. Topography (Amante and Eakins 2009) and kinematic plate boundaries (Sarah D. Stamps based on Bird 2003) of the East African Rift System (EARS). Plate boundary colors schematically indicate the western and eastern branches of the EARS.

Constituting one segment of the 5000 km long East African Rift System (EARS, Fig. 1), the Kenya Rift is host to an amazing landscape, wildlife and people, all of which somehow tie back to continental rifting processes. Although the youngest rifting phase in Kenya commenced in the Miocene, the east African region as a whole has been shaped by rifting episodes since Permian times (Bosworth and Morley 1994). The present active rift system runs from the Afar region in the north all the way south to Mozambique and is split into a western and an eastern branch that run around the Archean Tanzanian Craton (Chorowitz 2005, see Fig. 1). Generally speaking, the western branch is more seismically active, but deprived of magmatism, compared to the eastern branch, of which the Kenya Rift is part (Chorowitz 2005). Three processes characterize the EARS (Burke 1996) as well as the Kenya Rift specifically: normal faulting, volcanism and uplift.


The Tanzanian Craton together with the enveloping western and eastern EARS branches constitutes the broad, uplifted area coined the East African Plateau (~1200 m elevation, Strecker 1991; Simiyu and Keller 1997, Fig. 2). The onset of uplift of this plateau can be constrained to the Early Miocene with the help of one of the longest phonolitic lava flows on Earth (> 300 km, Wichura et al. 2010; 2011) and a whale that stranded inland 17 Ma (and was only recently found again after going missing for 30 years, Wichura et al. 2015). Plume-lithosphere interaction is thought responsible for the uplift (e.g. Wichura et al. 2010), although there is disagreement about the continuity of the low seismic velocity anomalies seen in the east African upper mantle and whether they are connected to the lower mantle. For example Ebinger and Sleep (1998), Hansen et al. (2012), Sun et al. (2017) and Torres Acosta et al. (2015) advocate for one East African superplume, while Pik et al. (2006) distinguish separate lower and upper mantle plumes and Davis and Sack (2002) and Halldórsson et al. (2014) consider a lower mantle plume splitting in the upper mantle.

Figure 2. Topography (Amante and Eakins 2009) and fault traces (GEM) of the central EARS. Triangles indicate off-rift volcanoes, dotted grey lines the three segments of the Kenya Rift.

Magmatism and volcanism

The northward motion of Africa over this hot mantle anomaly has been thought the cause of a north-to-south younging trend in the age of the ensuing EARS volcanism and rifting (e.g. Ebinger and Sleep 1998; George et al. 1998; Nyblade and Brazier 2002), although more recent studies arrive at a more spatially disparate and diachronous rifting evolution (Torres Acosta et al. 2015 and references therein). In general, massive emplacement of flood-phonolites preceded the onset of rifting in Kenya around 15 Ma (Torres Acosta et al. 2015). With ongoing rifting, and localization of faulting towards the rift axis, volcanism also migrated towards the center of the rift. Since the Miocene, massive amounts of volcanics have thus been emplaced (144,000-230,000 km3, MacDonald 1994; Wichura et al. 2011). Moreover, dyking also accommodated a significant part of the extension, with 22 to 26 % of the crust in the rift valley being composed of dykes (MacDonald 2012). Not surprisingly, the highlands directly around the rift valley, the Kenya Dome (Fig. 2) formed through a combination of volcanism and uplift (Davis and Slack 2002) with elevations of up to 1900 m.

The composition of rift magmatism is bimodal, showing phonolites and trachytes on the one side and nephelinites and basalt on the other, predominantly resulting from fractional crystallization of a basaltic source. The low viscosity of these magmas allows the young volcanoes in the volcano-tectonic axis to reach significant heights (see Fig. 3; MacDonald 2012). The most impressive volcanoes are to be found outside of the rift however (Fig. 2), with Mnt. Elgon reaching 4321 m and Africa’s highest mountains Mnt. Kenya and Mnt. Kilimanjaro reaching up to 5200 m and 5964 m, respectively (Chorowitz 2005).

Figure 3. View on the crater rim of the 400 ky old Mnt. Longonot volcano in the tectono-magmatic rift axis, at 2560 m asl. Courtesy of Corinna Kallich, GFZ Potsdam.

Normal faulting

The Kenya rift itself is composed of 3 asymmetric segments, distinguished by sharp changes in their orientation (Chorowitz 2005, Fig. 2). The 2300-3000 m high Elgeyo, Mau and Nguruman escarpments result from the steep Miocene east-dipping border faults in the west, while the antithetic border faults on the eastern side formed later during the Pliocene (Strecker et al. 1990). The older border faults formed along preexisting foliation generated by the Mozambique Belt orogeny in the late Proterozoic (Shackleton 1993; Hetzel and Strecker 1994). A change in strike of this foliation from NNE in the northern and southern Kenya rifts to NW determined the change in orientation in the central Kenya rift (Strecker et al. 1990). Consequently, different generations of faults in the northern and southern rift segments run parallel, while in the central segment, the Pleistocene change in extension direction from ENE-WSW/E-W to the present-day WNW-ESE/NW-SE directed extension results in obliquely reactivated border faults and younger, en echelon arranged left-stepping NNE-striking fault zones along the rift axis (Strecker et al. 1990). Extension is transferred between the different zones by coeval normal and strike-slip faulting or dense sets of normal faults.

Figure 4. View of lake Magadi and the Nguruman escarpment. Lake Magadi is a saline, alkaline lake, commercially mined for trona. Courtesy of Corinna Kallich, GFZ Potsdam.

Human evolution

The uplift, volcanism and normal faulting together have set the stage for human and animal evolution. For example, the shift in hoofed mammals from eating predominantly woods to grazing species evidences that the large-scale uplift modified air circulation patterns resulting in aridification and savannah-expansion at the expense of forested areas (Sepulchre et al. 2006; Wichura et al. 2015). The rift basins enabled the formation of large lakes, which were subsequently compartmentalized by tectonic and volcanic morphological barriers (Fig. 4). On the short-term, lake coverage varied due to tectonically induced changes in catchment areas, drainage networks and outlets. Maslin et al. (2014) actually found a correlation between this ephemeral lake coverage and hominin diversity and dispersal. Lake highstands link with the emergence of new species and allowed the spread of hominins north and southward out of east Africa. Remarkable, or what!

Amante, C. and Eakins B. W., 2009. NOAA Technical Memorandum NESDIS NGDC-24. National Geophysical Data Center, NOAA.
Bosworth, W. and Morley, C.K., 1994.  Tectonophysics 236, 93–115.
Burke, K., 1996. S. Afr. J. Geol. 99 (4), 339–409.
Chorowitz, J., 2005. J. Afr. Earth Sci. 43, 379-410.
Davis, P. M. and Slack, P. D. 2002. Geophys. Res. Lett. 29 (7), 1117.
Ebinger, C.J. and Sleep, N.H., 1998. Nature 395, 788-791.
George, R. et al., 1998.  Geology 26, 923–926.
Halldórsson, S. A. et al., 2014. Geophys. Res. Lett. 41, 2304–2311,
Hansen, S. E. et al., 2012.  Earth Planet. Sc. Lett. 319-320, 23-34.
Hetzel, R., Strecker, M.R., 1994. J. Struct. Geol. 16, 189–201.
Macdonald, R. et al., 1994a. J. Volcanol. Geoth. Res. 60, 301–325.
Macdonald, R., et al., 1994b. J. Geol. Soc. London 151, 879–888.
MacDonald, R., 2012. Lithos 152, 11-22.
Maslin, M. A. et al., 2014. Quaternary Sci. Rev. 101, 1-17.
Nyblade, A. A. and Brazier, R. A., 2002. Geology 30 (8), 755-758.
Pik, R. et al., 2006. Chem. Geol. 266, 100-114.
Sepulchre, P. et al., 2006. Science, 1419-1423.
Shackleton, R.M., 1993. Geological Society, London, Special Publications 76, 345–362.
Simiyu, S.M., Keller, G.R., 1997. Tectonophysics 278, 291–313.
Strecker, M., 1991. Das zentrale und südliche Kenia-rift unter besonderer berücksichtigung der neotektonischen entwicklung, habilitation, Universität Fridericiana.
Sun, M. et al., 2017.  Geophys. Res. Lett. 44, 12,116–12,124.
Torres Acosta, V. et al., 2015. Tectonics 34, 2367–2386.
Wichura, H. et al., 2010. Geology 38 (6), 543–546.
Wichura, H. et al , 2011. The Formation and Evolution of Africa: A Synopsis of 3.8 Ga of Earth History, eds. D. J. J. Van Hinsbergen, S. J. H. Buiter, T. H. Torsvik, C. and Gaina, S. J.
Wichura, H. et al., 2015. P. Natl. Acad. Sci. USA 112 (13), 3910-3915.

Being both strong and weak

Being both strong and weak

The Geodynamics 101 series serves to showcase the diversity of research topics and methods in the geodynamics community in an understandable manner. We welcome all researchers – PhD students to Professors – to introduce their area of expertise in a lighthearted, entertaining manner and touch upon some of the outstanding questions and problems related to their fields. For our latest ‘Geodynamics 101’ post, Postdoc Anthony Osei Tutu of GFZ Potsdam shares the outcomes of his PhD work, showing us that, like the lithosphere, it is OK to be weak sometimes!

Strength is not everything in achieving one’s goal. The lithospheric plate acts both strong and weak at times. This dual characteristic of the outermost part of the Earth, the crustal-lithospheric shell, is thought to have sustained plate tectonics throughout Earth’s history, in the presence of other controlling mechanisms such as the weak asthenospheric layer (Bercovici et al. 2000; Karato 2012). In the world of the lithospheric plates there is the saying “I might be strong and unbreakable, but sometimes and somewhere, I am very weak, soft and brittle” and this allows the plates to accommodate each other in their relative movements.

We all sometimes need to bring out the soft part in us to accommodate others such as friends, family or colleagues. For example, my graduate school, the Helmholtz-Kolleg GEOSIM, an experiment by the Helmholtz Association, GFZ-Potsdam, University of Potsdam and Free University of Berlin, brought together two or more experts in mathematics and geosciences to collaborate on and serve as PhD supervisors for answering some of Earth Sciences’ pressing questions. The many, many benefits of this multidisciplinary PhD supervising approach also came with challenges. Sometimes, the different supervisors would make opposing/contrasting suggestions to investigate a particular problem according to the experience of some students and myself. Then it falls on you as the student to stand firm (i.e. be strong) on what you believe works for your experiments and at the same time to be receptive (i.e. flexible or soft) to the different suggestions, while keeping in mind the limited time you have as a PhD student.

Figure 1: Schematic plot of the conditions in a subduction system (left) aiding or (right) hindering global plate motions.

The both strong and weak behavior of the lithospheric plates was one of the conclusions of my PhD study. Besides the strong plate interiors (Zhong and Watts 2013), weak regions along the plate boundaries, aided by sediment and water (see Fig. 1), are required to give the low friction between the subducting and overriding plates (Moresi and Solomatov 1998; Sobolev and Babeyko 2005), combined with a less viscous sublithospheric mantle. This combination was key to match the magnitude and direction of present-day global plate motions in the numerical modeling study (Osei Tutu et al. 2018). I used the global 3D lithosphere-asthenosphere numerical code SLIM3D (Popov and Sobolev 2008) with visco-elasto-plastic rheology coupled to a mantle flow code (Hager and O’Connell 1981) for the investigation. To understand the influence of intra-plate friction (brittle/plastic yielding) and asthenospheric viscosity on present-day plate motions, I tested a range of strengths of the plate boundary. Past numerical modeling studies (Moresi and Solomatov 1998; Crameri and Tackley 2015) have suggested that small friction coefficients (μ < 0.1, yield stress ~100 MPa) can lead to plate tectonics in models of mantle convection. This study shows that in order to match present-day plate motions and net rotation, the static frictional parameter must be less than 0.05 (15 MPa yield stress). I am able to obtain a good fit with the magnitude and orientation of observed plate velocities (NUVEL-1A) in a no-net-rotation reference frame with μ < 0.04 and a minimum asthenosphere viscosity of 5•1019 Pas to 1•1020 Pas (Fig. 2). The estimates of net-rotation (NR) of the lithosphere suggest that amplitudes of ~0.1– 0.2 °/My, similar to most observation-based estimates, can be obtained with asthenosphere viscosity cutoff values of ~1•1019 Pas to 5•1019 Pas and a friction coefficient μ < 0.05.

Figure 2: Set of predicted global plate motions for varying asthenosphere viscosity and plate boundary frictions, modified after Osei Tutu et al. (2018). Rectangular boxes show calculations with RMS velocities comparable to the observed RMS velocity of NUVEL-1A (DeMets et al. 2010).

The second part of my PhD study focused on the responses of the strong plate interiors to the convecting mantle below by evaluating the influence of shallow and deep mantle heterogeneities on the lithospheric stress field and topography. I explored the sensitivity of the considered surface observables to model parameters providing insights into the influence of the asthenosphere and plate boundary rheology on plate motion by testing various thermal-density structures to predict stresses and topography. Lithospheric stresses and dynamic topography were computed using the model setup and rheological parameters that gave the best fit to the observed plate motions (see rectangular boxes in Fig. 2). The modeled lithosphere stress field was compared the World Stress Map 2016 (Heidbach et al. 2016) and the modeled dynamic topography to models of observed residual topography (Hoggard et al. 2016; Steinberger 2016). I tested a number of upper mantle thermal-density structures. The thermal structure used to calculate the plate motions before is considered the reference thermal-density structure, see also Osei Tutu et al. (2017). This reference thermal-density structure is derived from a heat flow model combined with a sea floor age model. In addition I used three different thermal-density structures derived from global S-wave velocity models to show the influence of lateral density heterogeneities in the upper 300 km on model predictions. These different structures showed that a large portion of the total dynamic force generating stresses in the crust/lithosphere has its origin in the deep mantle, while topography is largely influenced by shallow heterogeneities. For example, there is hardly any difference between the stress orientation patterns predicted with and without consideration of the heterogeneities in the upper mantle density structure across North America, Australia and North Africa. However, inclusion of crustal thickness variations in the stress field simulations (as shown in Fig. 3a) resulted in crustal dominance in areas of high altitude in terms of stress orientation, for example in the Andes and Tibet, compared to the only-deep mantle contributions (as shown in Fig. 3b).

Figure 3: Modeled lithosphere stress field in the Andes considering (a) crustal thickness variations from the CRUST 1.0 model as well as lithospheric variations and (b) uniform crustal and lithospheric thicknesses.

The outer shell of the solid Earth is complex, exhibiting different behaviors on different scales. In our quest to understand its dynamics, we can learn from the lithospheric plate’s life cycle how to live our lives and preserve our existence as scientist-humans by accommodating one another. After all, they have existed for billions of years.



Bercovici, David, Yanick Ricard, and Mark A. Richards. 2000. “The Relation Between Mantle Dynamics and Plate Tectonics: A Primer.” 5–46.

Crameri, Fabio and Paul J. Tackley. 2015. “Parameters Controlling Dynamically Self-Consistent Plate Tectonics and Single-Sided Subduction in Global Models of Mantle Convection.” Journal of Geophysical Research: Solid Earth 120(5):3680–3706, 10.1002/2014JB011664.

DeMets, Charles, Richard G. Gordon, and Donald F. Argus. 2010. “Geologically Current Plate Motions.” Geophys. J. Int 181:1–80.

Hager, BH and RJ O’Connell. 1981. “A Simple Global Model of Plate Dynamics and Mantle Convection.” Journal of Geophysical Research: Solid Earth, 86(B6):4843–4867, 10.1029/JB086iB06p04843.

Heidbach, Oliver, Mojtaba Rajabi, Moritz Ziegler, Karsten Reiter, and Wsm Team. 2016. “The World Stress Map Database Release 2016 -Global Crustal Stress Pattern vs. Absolute Plate Motion.” Geophysical Research Abstracts EGU General Assembly 18:2016–4861.

Hoggard, M. J., N. White, and D. Al-Attar. 2016. “Global Dynamic Topography Observations Reveal Limited Influence of Large-Scale Mantle Flow.” Nature Geoscience 9(6):456–63, 10.1038/ngeo2709.

Karato, Shun-Ichiro. 2012. “On the Origin of the Asthenosphere.” Earth and Planetary Science Letters 321–322:95–103.

Moresi, Louis and Viatcheslav Solomatov. 1998. “Mantle Convection with a Brittle Lithosphere: Thoughts on the Global Tectonic Styles of the Earth and Venus.” Geophysical Journal International 133(3):669–82, 10.1046/j.1365-246X.1998.00521.x.

Osei Tutu, A., S. V Sobolev, B. Steinberger, A. A. Popov, and I. Rogozhina. 2018. “Evaluating the Influence of Plate Boundary Friction and Mantle Viscosity on Plate Velocities.” Geochemistry, Geophysics, Geosystems n/a-n/a, 10.1002/2017GC007112.

Popov, A. A. and S. V. Sobolev. 2008. “SLIM3D: A Tool for Three-Dimensional Thermomechanical Modeling of Lithospheric Deformation with Elasto-Visco-Plastic Rheology.” Physics of the Earth and Planetary Interiors 171(1–4):55–75.

Sobolev, S. V. and A. Y. Babeyko. 2005. “What Drives Orogeny in the Andes?” Geology 33(8).

Steinberger, Bernhard. 2016. “Topography Caused by Mantle Density Variations: Observation-Based Estimates and Models Derived from Tomography and Lithosphere Thickness.” Geophysical Journal International 205(1):604–21, 10.1093/gji/ggw040.

Osei Tutu, A., B. Steinberger, S. V Sobolev, I. Rogozhina, and A. Popov. 2017. "Effects of Upper Mantle Heterogeneities on Lithospheric Stress Field and Dynamic Topography." Solid Earth Discuss.,, in review, 2017

Zhong, Shijie and A. B. Watts. 2013. “Lithospheric Deformation Induced by Loading of the Hawaiian Islands and Its Implications for Mantle Rheology.” Journal of Geophysical Research: Solid Earth 118(11):6025–48, 10.1002/2013JB010408.

Remarkable Regions – The India-Asia collision zone

Remarkable Regions – The India-Asia collision zone

Every 8 weeks we turn our attention to a Remarkable Region that deserves a spot in the scientific limelight. This week we zoom in on a particular part of the eastern Tethys as Adina Pusok discusses the India-Asia collision zone. She is a postdoctoral researcher at the Institute of Geophysics and Planetary Physics, Scripps Institution of Oceanography, UCSD, US.

Without doubt, one of the most striking features of plate tectonics and lithospheric deformation on Earth is the India-Asia collision zone, largely comprised of the Himalayan and Karakoram mountain belts and the Tibetan plateau. What makes this collision zone so remarkable? For one, Tibet is the largest, highest and flattest plateau on Earth with an average elevation exceeding 5 km, and it includes over 80% of the world’s land surface higher than 4 km. Then, the bordering Himalayas and the Karakoram Mountains include the only peaks on Earth reaching more than 8 km above sea level.

It makes one wonder, how can such a mountain belt and high plateau form? Most of the major mountain belts and orogenic plateaus on Earth are found within the overlying plate of subduction and/or collision zones (e.g. the Alps, the Andes, the American Cordilleras etc.). When an ocean closes and two continental plates meet at a destructive (subduction) boundary, the continents themselves collide. Such collisions result in intense deformation at the edges of the colliding plates. Neither continent can be subducted into the mantle due to the buoyancy of continental crust, so the forces that drive the plate movement prior to collision are brought to act directly on the continental lithosphere itself. At this stage, further convergence of the plates must be taken up by deforming one or both of the plates of continental lithosphere over hundreds of kilometres [Figure 1]. Mountain belts can form under these circumstances.

Figure 1 Global map of surface velocities and the second invariant of strain rate (from Moresi [2015]). The surface velocities show the location and extent of plates, and the strain rate map highlights the fact that most of the deformation is concentrated at plate boundaries (high strain rates), while the continental interiors have little or no deformation (low strain rates). In some places, deformation occurs over broader regions, especially following mountain belts. These boundaries are called diffuse plate boundaries. The white rectangle roughly indicates the extent of the India-Asia collision zone.

The Himalayas and the Tibetan plateau are no different. Following the closure of the Tethys ocean (see earlier blog post), the Indian continent collided with Eurasia around 50 million years ago (e.g., Patriat and Achache [1984]), thus giving rise to this anomalously high region. This tectonic boundary is complex and changes character along its length. The Tibetan plateau is a collage of continental blocks (terranes) that were added successively to the Eurasian plate during the Paleozoic and Mesozoic [Figure 2]. The boundaries between these terranes are marked by scattered occurrences of ophiolitic material, which are rocks characteristic of oceanic lithosphere. The Himalayas represent the traditional accretionary wedge formed by folding and thrusting of sediments scraped off the subducting slab.

Figure 2 Simplified tectonic map of Tibet and surrounding region showing approximate boundaries of the major terranes, suture zones, and strike-slip faults (from Ninomiya and Bihong [2016]). Blocks and terranes: ALT-EKL-QL: Altyn Tagh–East Kunlun–Qilian terrane; BS: Baoshan terrane; HM: Himalayan terrane; IC: Indo-China block; KA: Kohistan Arc terrane; LA: Ladakh arc terrane; LC: Lincang–Sukhothai–Chanthaburi Arc terrane; LST: Lhasa terrane; NCB: North China Block; NQT: North Qiangtang terrane; QT: Qiangtang terrane; SP: South Pamir terrane; SPGZ: Songpan-Ganze terrane; SQT: South Qiangtang terrane; TC: Tengchong terrane; TSH: Tianshuihai terrane; WB: West Burma terrane; WKL: West Kunlun terrane. Suture zones: BNS: Bangong-Nujiang Suture; EKLS: East Kunlun Suture; ITS: Indus-Tsangbo Suture; JSS: Jinsha Suture; LSS: Longmu Tso–Shuanghu–Menglian–Inthanon Suture; WKLS: West Kunlun Suture. Basins: QB: Qaidam Basin; KB: Kumkol Basin. Faults: ALT: Altyn Tagh Thrust; ALTF: Altyn Tagh Fault; KKF: Karakorum Fault; LMST: Longmen Shan Thrust; MFT: Main Frontal Thrust; NQLT: North Qilian Thrust; RRF: Red River Fault; SGF: Sagaing Fault; XXF: Xianshui River–Xiaojiang Fault.

Interestingly, the India-Asia collision orogen is not just the youngest and most spectacular active continent collision belt, it is also the most studied research area on Earth. Studies on this region span a wide range of topics and methods for over more than 100 years. I am not sure if it is the fascination with the highest mountain on Earth (Mt. Everest was actually climbed for the first time as late as 1953 by Tenzing Norgay and Edmund Hillary), similar to our fascination for exploring the Moon, Mars and the other planets in our Solar System nowadays, or the hope that studying the youngest orogeny will help us decipher the older ones (soon to realize different mountain belts evolve differently).

To understand the magnitude of the work done in the past 100 years, a simple search of the keywords “India Asia collision” on Google Scholar yielded ~90k results, and a more focused geosciences search on Web of Science (where I filtered the results to those from geophysics, geochemistry, geology, geosciences multidisciplinary only) yielded >1600 results for the same keywords (other keywords: “Himalaya” > 5600 results, “Tibet” > 6500 results, “India Asia” > 2200 results). These numbers can be intimidating to a new student taking on the topic, but it is a topic worth studying and I’ll explain why below.

From a general perspective, it is important to study the India-Asia collision zone due to the interaction between tectonics and climate and the formation of the Indian monsoon [Molnar et al., 1993], but also because it is a highly populated area (>200 million people in the Hindu Kush Himalaya region) regularly shaken by natural phenomena, such as earthquakes, floods or landslides. For example, the last large earthquake in Nepal, the Gorkha earthquake (Mw 7.8) in April 2015 caused more than 9000 deaths.

From a geophysics point of view, understanding mountain-building processes and the driving forces of plate tectonics has been one of the long-term goals of solid Earth sciences community. The India-Asia collision zone is one of the best examples in which subduction, continental collision and mountain building can be studied in a global plate tectonics perspective. Prior to plate tectonics theory, Argand [1924] and Holmes [1965] thought that the Himalayan Mountains and Tibetan Plateau had been raised due to the northern edge of the Indian craton underthrusting the entire region, causing shortening and thickening of the crust to ∼80 km. This perspective remains widely accepted, but recent ideas suggest that other processes are equally important (more below).

Today, the challenge lies in refining our understanding of the dynamics of India-Asia collision by elucidating the connections between the wealth of observations available and the underlying processes occurring at depth. Decades of study have produced data sets across various disciplines, including: active tectonics, Cenozoic geology, seismicity, global positioning system (GPS) measurements, seismic profiles, tomography, gravity anomalies, mantle-crustal anisotropy, paleomagnetism, geochemistry or magnetotelluric studies. Of these, the GPS data stands out as it clearly shows the distributed deformation across the entire collision zone and suggests that this is a highly dynamic area [Figure 3].

Figure 3 Horizontal GPS velocities of crustal motion around the Tibetan Plateau relative to stable Eurasia from Liang et al. [2013].

Collectively, all these observation data sets stand as a different piece in the puzzle of the India-Asia collision. However, the same data sets can support a number of competing and sometimes mutually exclusive mechanisms for the uplift of the Tibetan Plateau. For example, the mantle lithosphere beneath Tibet has been proposed to be cold, hot, thickened by shortening, or thinned by viscous instability. Other controversies include the degree of mechanical coupling between the crust and deeper lithosphere and the nature of large-scale deformation. It is no surprise then, that several hypotheses emerged over time trying to explain the anomalous rise of the Himalayas and Tibetan Plateau [Figure 4]:

  1. Figure 4 Schematic cartoons of tectonic models proposed to explain the thickening and uplift of the Himalayas and the Tibetan Plateau. (Source: personal institutional web page of A. Ozacar).

    Wholescale underthrusting of the Indian plate below the Asian continent [e.g. Argand, 1924].

  2. The thin-sheet model or distributed homogeneous shortening [e.g. England and McKenzie, 1982].
  3. Homogeneous thickening of a weak, hot Asian crust, involving a large amount of magmatism [e.g. Dewey and Burke, 1973].
  4. Slip-line field model to account for the brittle deformation in and around the Tibetan Plateau and to explain extrusion of SE Tibet away from Indian continent [e.g. Molnar and Tapponnier, 1975]. The same group also proposes a time-dependent model for the growth of Tibetan plateau [e.g. Tapponnier et al., 2001], in which successive intracontinental subduction zones maintain the stepwise growth and rise of the plateau.
  5. Lower crustal flow models for the exhumation of the Himalayan units and lateral spreading of the Tibetan plateau [e.g. Royden et al., 1997, Beaumont et al., 2001].
  6. Delamination or convective removal of the lithospheric mantle that induced isostatic movement, lifting the Tibetan Plateau [e.g. Molnar, 1988].


These models were applied either to the Tibetan Plateau or the Himalayan mountain belt and were able to explain the formation of specific tectonic and geological features. However, there is no conclusive answer on which of the hypotheses works best for the entire orogen, and instead, more questions arise:

  • Which forces are at work during continental collision and mountain building?
  • What is the deformation history and evolution of this plate boundary?
  • How was the subduction accommodated in the Neo-Tethys?
  • How does subduction evolve during continental collision?
  • What drives the present-day fast convergence (~4-5 cm/yr) between India and Eurasia?
  • Which forces propagated India northwards between 70-50 million years at anomalously high speeds (up to 16 – 20cm/yr)?
  • How can you form such large elevations over such extended areas?
  • What is the effect of surface processes on uplift?
  • What is the structure at depth beneath the Himalayas and Tibetan plateau?
  • How do the Indian and Eurasia plates deform during collision?
  • How is the deformation accommodated during continental collision?
  • How do mountain belts form and why not all mountain belts look the same?
  • How did the crust beneath Himalaya and Tibet reach double-crustal thickness (normal continental crust is 35-40 km thick, whereas the crust beneath the Himalaya and Tibet is 70-100 km thick)?
  • Which mechanisms help sustain the high topographic amplitudes?
  • Why should an area as broad as the Tibetan Plateau be uplifted so high compared to other mountain belts following collision?
  • Did the Tibetan Plateau and Himalayan mountain belt rise continuously or diachronously?
  • Which the proposed models [Figure 4] can be applied, and where?
  • How do lithospheric heterogeneities and rheology affect the deformation pattern?
  • What is the degree of mechanical coupling between the crust and deeper lithosphere? Is it the “jelly sandwich” model (e.g., Burov and Watts [2006]) or the “creme-brulee” model (e.g., Jackson [2002], see earlier blog post)?
  • Why do the Himalayas have a convex curvature?
  • What about the high deformation of the prominent Himalayan syntaxes (the inflection points of the Himalayan belt): Nanga Parbat in the west and Namche Barwa in the east?
  • What is the effect of the India-Asia collision on climate? Do the Himalayas affect the Indian monsoon or is it the other way around? A chicken-and-egg question.

Seriously, can I even stop asking questions? The question that fascinated me the most during my graduate studies was “Why is the Himalayan-Tibet region so high and broad compared to other mountain belts?”. If we tune our models to Earth parameters, can we build such large elevations in computer simulations? Which factors and forces are at play? Using 3-D numerical models to address this question [Pusok and Kaus, 2015], we were able to obtain distinct topographic modes (different types of mountain belts) [Figure 5] and to show that building topography is an interplay between providing the energy to the system and the ability of that system to store it over longer periods of time. We also suggest that the reason why Himalaya-Tibet is different from the Alps, for example, is because the shape and elevation of mountain ranges can vary depending on the boundary conditions (plate driving forces that control convergence velocity and lithospheric heterogeneities such as the Tarim Basin) and internal factors (rheology), but also on the evolution stage they are in.

To sum up, it is clear that many of the above questions remain unanswered. But I think this is good news, meaning that in the future, exciting new results will shape our understanding of this remarkable region.

Figure 5 3-D Simulation results showing different modes of surface expressions in continental collision models. Modified from Pusok and Kaus [2015].

Argand, E. (1924). La tectonique de l’Asie. Proc. 13th Int. Geol. Cong., 7:171–372.

Beaumont, C., Jamieson, R. A., Nguyen, M. H., and Lee, B. (2001). Himalayan Tectonics Explained by Extrusion of a Low-Viscosity Crustal Channel Coupled to Focused Surface Denudation. Nature, 414:738–742.

Burov, E. B. and Watts, A. B. (2006). The long-term strength of continental lithosphere: “jelly sandwich” or “crème brûlée”? GSA Today, 16(1):4.

Dewey, J. F. and Burke, K. (1973). Tibetan, Variscan, and Precambrian Basement Reactivation: Products of Continental Collision. The Journal of Geology, 81(6):683–692.

England, P. and McKenzie, D. (1982). A Thin Viscous Sheet Model for Continental Deformation. Geophys. J. R. astr. Soc., 70:295–321.

Holmes, A. (1965). Principles of Physical Geology. The Ronald Press Company, New York, second edition.

Jackson, J. (2002). Strength of the Continental Lithosphere: Time to Abandon the Jelly Sandwich? GSA Today, 4–9.

Liang, S., Gan, W., Shen, C., Xiao, G., Liu, J., Chen, W., Ding, X., and Zhou, D. (2013). Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements. Journal of Geophysical Research: Solid Earth, 118:1–11.

Molnar, P. and Tapponnier, P. (1975). Cenozoic Tectonics of Asia: Effects of a Continental Collision. Science, 189:419–426.

Molnar, P. (1988). A Review of Geophysical Constraints on the Deep Structure of the Tibetan Plateau, the Himalaya and the Karakoram, and their Tectonic Implications. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 326(1589):33–88.

Molnar, P., England, P., and Martinod, J. (1993). Mantle Dynamics, Uplift of the Tibetan Plateau, and the Indian Monsoon. Reviews of Geophysics, 31:357–396.

Moresi, L. (2015). Computational Plate Tectonics and the Geological Record in the Continents. SIAM News, 48:1–6.

Ninomiya, Y. and Bihong Fu, B. (2016). Regional Lithological Mapping Using ASTER-TIR Data: Case Study for the Tibetan Plateau and the Surrounding Area. Geosciences 2016, 6(3), 39; doi:10.3390/geosciences6030039.

Patriat, P. and Achache, J. (1984). India-Eurasia collision chronology has implications for crustal shortening and driving mechanism of plates. Nature, 311:615–621.

Pusok, A. E. and Kaus, B. J. P. (2015). Development of topography in 3-D continental-collision models. Geochemistry, Geophysics, Geosystems, 16(5):1378–1400.

Royden, L. H., Burch el, B. C., King, R., Wang, E., Chen, Z., Shen, F., and Liu, Y. (1997). Surface Deformation and Lower Crustal Flow in Eastern Tibet. Science, 276(5313):788–790.

Tapponnier, P., Zhiqin, X., Roger, F., Meyer, B., Arnaud, N., Wittlinger, G., and Jingsui, Y. (2001). Oblique Stepwise Rise and Growth of the Tibet Plateau. Science, 294(5547):1671–1677.


One month to AGU!

One month to AGU!
As the leaves are falling; the sun is going down before you leave the office; and the stores are stacking up on Christmas decorations, it’s time to face the facts: it’s almost AGU! It shouldn’t come as a surprise, but just in case. Don’t worry, there is still time to reread your abstract to see what you’re supposed to be presenting, figure out how to do that in the several weeks that are left and wrap it all up in a convincing poster or talk.
But first, check out these tips on presenting your work: our EGU GD blogs on creating prize-winning posters and selecting proper color schemes, EGU resources on making posters and these Science magazine posts on creating and giving great oral presentations. And while you’re at it, research New Orleans  highlights, such as the French quarter and the excellent cuisine. There is still time to register for some of the AGU field trips too! And wait, did you book your hotel and flights? Then, get to work!