Refine
Year of publication
- 2017 (2) (remove)
Document Type
- Article (1)
- Doctoral Thesis (1)
Language
- English (2)
Is part of the Bibliography
- yes (2)
Keywords
- numerical modeling (2) (remove)
Institute
- Institut für Geowissenschaften (2) (remove)
Inherited rheological structures in the lithosphere are expected to have large impact on the architecture of continental rifts. The Turkana depression in the East African Rift connects the Main Ethiopian Rift to the north with the Kenya rift in the south. This region is characterized by a NW-SE trending band of thinned crust inherited from a Mesozoic rifting event, which is cutting the present-day N-S rift trend at high angle. In striking contrast to the narrow rifts in Ethiopia and Kenya, extension in the Turkana region is accommodated in subparallel deformation domains that are laterally distributed over several hundred kilometers. We present both analog experiments and numerical models that reproduce the along-axis transition from narrow rifting in Ethiopia and Kenya to a distributed deformation within the Turkana depression. Similarly to natural observations, our models show that the Ethiopian and Kenyan rifts bend away from each other within the Turkana region, thus forming a right-lateral step over and avoiding a direct link to form a continuous N-S depression. The models reveal five potential types of rift linkage across the preexisting basin: three types where rifts bend away from the inherited structure connecting via a (1) wide or (2) narrow rift or by (3) forming a rotating microplate, (4) a type where rifts bend towards it, and (5) straight rift linkage. The fact that linkage type 1 is realized in the Turkana region provides new insights on the rheological configuration of the Mesozoic rift system at the onset of the recent rift episode. Plain Language Summary The Turkana depression in the Kenya/Ethiopia borderland is most famous for its several million years old human fossils, but it also holds a rich geological history of continental separation. The Turkana region is a lowland located between the East African and Ethiopian domes because its crust and mantle have been stretched in a continent-wide rift event during Cretaceous times about 140-120 Ma ago. This thin lithosphere exerted paramount control on the dynamics of East African rifting in this area, which commenced around 15 Ma ago and persists until today. Combining analog "sandbox" experiments with numerical geodynamic modeling, we find that inherited rift structures explain the dramatic change in rift style from deep, narrow rift basins north and south of the Turkana area to wide, distributed deformation within the Turkana depression. The failed Cretaceous rift is also responsible for the eastward bend of the Ethiopian rift and the westward bend of the Kenyan rift when entering the Turkana depression, which generated the characteristic hook-shaped form of present-day Lake Turkana. Combing two independent modeling techniques-analog and numerical experiments-is a very promising approach allowing to draw robust conclusions about the processes that shape the surface of our planet.
The timing and location of the two largest earthquakes of the 21st century (Sumatra, 2004 and Tohoku 2011, events) greatly surprised the scientific community, indicating that the deformation processes that precede and follow great megathrust earthquakes remain enigmatic. During these phases before and after the earthquake a combination of multi-scale complex processes are acting simultaneously: Stresses built up by long-term tectonic motions are modified by sudden jerky deformations during earthquakes, before being restored by multiple ensuing relaxation processes.
This thesis details a cross-scale thermomechanical model developed with the aim of simulating the entire subduction process from earthquake (1 minute) to million years’ time scale, excluding only rupture propagation. The model employs elasticity, non-linear transient viscous rheology, and rate-and-state friction. It generates spontaneous earthquake sequences, and, by using an adaptive time-step algorithm, recreates the deformation process as observed naturally over single and multiple seismic cycles. The model is thoroughly tested by comparing results to those from known high- resolution solutions of generic modeling setups widely used in modeling of rupture propagation. It is demonstrated, that while not modeling rupture propagation explicitly, the modeling procedure correctly recognizes the appearance of instability (earthquake) and correctly simulates the cumulative slip at a fault during great earthquake by means of a quasi-dynamic approximation.
A set of 2D models is used to study the effects of non-linear transient rheology on the postseismic processes following great earthquakes. Our models predict that the viscosity in the mantle wedge drops by 3 to 4 orders of magnitude during a great earthquake with magnitude above 9. This drop in viscosity results in spatial scales and timings of the relaxation processes following the earthquakes that are significantly different to previous estimates. These models replicate centuries long seismic cycles exhibited by the greatest earthquakes (like the Great Chile 1960 Earthquake) and are consistent with the major features of postseismic surface displacements recorded after the Great Tohoku Earthquake.
The 2D models are also applied to study key factors controlling maximum magnitudes of earthquakes in subduction zones. Even though methods of instrumentally observing earthquakes at subduction zones have rapidly improved in recent decades, the characteristic recurrence interval of giant earthquakes (Mw>8.5) is much larger than the currently available observational record and therefore the necessary conditions for giant earthquakes are not clear. Statistical studies have recognized the importance of the slab shape and its surface roughness, state of the strain of the upper plate and thickness of sediments filling the trenches. In this thesis we attempt to explain these observations and to identify key controlling parameters. We test a set of 2D models representing great earthquake seismic cycles at known subduction zones with various known geometries, megathrust friction coefficients, and convergence rates implemented. We found that low-angle subduction (large effect) and thick sediments in the subduction channel (smaller effect) are the fundamental necessary conditions for generating giant earthquakes, while the change of subduction velocity from 10 to 3.5 cm/yr has a lower effect. Modeling results also suggest that having thick sediments in the subduction channel causes low static friction, resulting in neutral or slightly compressive deformation in the overriding plate for low-angle subduction zones. These modeling results agree well with observations for the largest earthquakes. The model predicts the largest possible earthquakes for subduction zones of given dipping angles. The predicted maximum magnitudes exactly threshold magnitudes of all known giant earthquakes of 20th and 21st centuries.
The clear limitation of most of the models developed in the thesis is their 2D nature. Development of 3D models with comparable resolution and complexity will require significant advances in numerical techniques. Nevertheless, we conducted a series of low-resolution 3D models to study the interaction between two large asperities at a subduction interface separated by an aseismic gap of varying width. The novelty of the model is that it considers behavior of the asperities during multiple seismic cycles. As expected, models show that an aseismic gap with a narrow width could not prevent rupture propagation from one asperity to another, and that rupture always crosses the entire model. When the gap becomes too wide, asperities do not interact anymore and rupture independently. However, an interesting mode of interaction was observed in the model with an intermediate width of the aseismic gap: In this model the asperities began to stably rupture in anti-phase following multiple seismic cycles. These 3D modeling results, while insightful, must be considered preliminary because of the limitations in resolution.
The technique developed in this thesis for cross-scale modeling of seismic cycles can be used to study the effects of multiple seismic cycles on the long-term deformation of the upper plate. The technique can be also extended to the case of continental transform faults and for the advanced 3D modeling of specific subduction zones. This will require further development of numerical techniques and adaptation of the existing advanced highly scalable parallel codes like LAMEM and ASPECT.