TS
Tectonics and Structural Geology

minds over methods

Minds over Methods: Massively dilatant faults in Iceland – from surface to subsurface structures

Minds over Methods: Massively dilatant faults in Iceland – from surface to subsurface structures
In this Minds over Methods we don’t have one, but two scientists talking about their research! Michael Kettermann and Christopher Weismüller, both from Aachen University, explain us about the multidisciplinary approach they use to understand more about massively dilatant faults. How do they form and what do they look like at depth?

Massively dilatant faults in Iceland – from surface to subsurface structures

Michael Kettermann & Christopher Weismüller, RWTH Aachen University

Michael (left) and Christopher (right) in the field. Credit: Michael Kettermann and Marianne Sophie Hollinetz.

Iceland is a volcanic island in a unique setting on the Mid-Atlantic Ridge, separating the Eurasian and North American plates. A deep mantle plume lies beneath Iceland, and the combination of rift and plume leads to very active basaltic volcanism. Ubiquitous features along the rift zone are normal faults, often exquisitely exposed at the surface. Normal faults in basalts are also common in many volcanic provinces, like Hawaii, the East African Rift, and along mid ocean ridges. These faults often form as massively dilatant faults (MDF), which show apertures up to tens of meters at the surface and supposedly have large volumes of open voids in the subsurface.

Figure 1. View along-strike a massively dilatant fault in layered basalt. The geometry of the vertical fault faces is prescribed by the cooling joints. Several basalt columns have been loosened and dropped into the fault, now being stuck in the fault (top) or filling the cavity (bottom). Opening width < 3 m. Credit: Michael Kettermann.

These openings form pathways for fluids like magma or hydrothermal waters and consequently are of importance for volcanic plumbing systems, mineralization and geothermal energy supply.

Iceland provides a perfect natural laboratory to study MDF. Due to its position on the Mid-Atlantic Ridge, Iceland is cut by extensional fault systems roughly from southwest to north. A wide range of oblique extensional to pure extensional faults can be observed mostly in flood basalts, but also in sub-glacially formed hyaloclastites (weaker volcanic sediments), pillow lavas and occasionally sediment layers formed in warmer times. Outcropping rocks in Iceland are younger than 20 Ma distal from the rift (eastern and western Iceland), while tectonic and corresponding volcanic activity at the ridge (central Iceland) constantly causes the formation of new rocks. The rough climate hinders soil formation and vegetation to overgrow faults, providing unique outcrop conditions.

While it is relatively easy to access and study the faults at surface level, investigations into the subsurface are much more challenging. Direct observations are only possible down to depths of some tens of meters by climbing into the fractures. Cavities are often filled with rubble, sediments, water or snow (Fig. 1). Steep, open fractures with meter-scale aperture are hard to detect with geophysical methods (seismic reflection/refraction, ground penetrating radar, electrical resistivity tomography) at depths greater than some meters.

We therefore started the massively dilatant fault project, a multidisciplinary, integrated project bringing together remote sensing, fieldwork, analogue modelling and numerical simulations. In essence, we utilize a modelling approach to recreate the structure and evolution of MDF at depth, using real 3D surface data as input and comparison data set.

 

Drone mapping and photogrammetry

In a first step, we capture and analyse the surface expressions of MDF at a number of representative fault areas in Iceland. To this end, we flew 27 drone surveys during our five weeks long field season in summer 2017, covering a total length of more than 42 km of faults. Luckily, in the Icelandic summer the days are very long, so the National Park Service allowed us to fly the drones early in the morning and late in the evening outside of tourist hours. For each area, we took several hundred to thousands of overlapping photographs (e.g. Fig. 2). We processed these sets with photogrammetry software, applying the Structure from Motion (SfM) technique. SfM is an increasingly popular, fast and cheap technique to reconstruct high-resolution 3D information from 2D images. This allows us to recreate digital elevation models and ortho-rectified photo-mosaics of the faults in resolutions better than 15 cm per pixel (Fig. 3). The largest area at the famous Thingvellir fissure swarm covers a length of almost 7 km with an average resolution of 11 cm per pixel.

We use these digital elevation models and ortho-photos to retrieve a wide range of structural data. Mapping the fault traces in a GIS software allows for the measurement of fault opening width, throw, orientation and length. From throw and aperture, we can then estimate the fault dip at depth. Digital elevation models further provide surface dip data that we then compare with observations from analogue models.

Figure 2. A drone photograph facing South of the Almannagjá fault in Thingvellir, where the Thingvallavegur road crosses the fault. The Almannagjá fault resembles the western shoulder of the Thingvellir graben system with locally > 50 m opening width and 40 m vertical offset. Credit: Christopher Weismüller. .

Figure 3. Digital elevation model created from drone photographs using photogrammetry software. It contains the faults at Sandvik on the Reykjanes Peninsula (SW Iceland). The detail panes (red square) show the DEM (right) at a higher zoom level and the corresponding ortho-rectified photograph (left). The bridge crossing the fault depicted in the detail panes is a famous touristic spot, known as „The bridge between the continents“, since the fault symbolically divides the North American and Eurasion plates. Credit: Michael Kettermann.

 

 

 

 

 

 

 

Figure 4. Sideview of an analogue model showing three timesteps of the development of a massivley dilatant fault and associated fractures in hemihydrate (Bücken, 2017). Note the tilted block developing at the surface of the model and the dilatant jogs and voids in the subsurface. The opening at the surface is not directly linked to the fault at depth, but caused by the rotation of the tilted block. Credit: Daniel Bücken..

Modelling approach

For the analogue modelling approach on the hundreds to thousand meter scale, we use cohesive powders as modelling material (Bücken, 2017). Especially hemihydrate powder has been proven suited to model dilatant fractures (Holland et al., 2006; Kettermann et al., 2016; van Gent et al., 2010) as it has a well characterized true cohesion and tensile strength. Faults in Iceland transform from opening mode fractures to shear mode faults at depth when overburden stress is high enough. As we are interested in the upper dilatant parts of the faults, i.e. above the shear mode faulting, we chose a basement-fault controlled approach, where a rigid basement represents the shear mode fault. It moves down-dip along a predefined surface, deforming the powder sieved on top. The basement fault dip follows the data we derived from the field and is set to 60° – 65°. The scale of the models calculates from strength and weight of the natural prototype and the modelling material. 1 cm of powder equals about 50 m of basalt.

 

Comparison of models and nature

Results show a close similarity between field and experiment at the surface structures. Open fractures form with large apertures at the surface and often we observe the formation of tilted blocks (Fig. 4). The existence and scaled dimensional similarity of fractures and tilted blocks in the field and using a scaled material suggest a validity of other observations in the models. Glass sidewalls in the analogue models provide the opportunity to examine how the faults evolve at depth. We observe that large caves form underneath these blocks and we predict that these must exist in the field as well, albeit potentially filled with rubble. Our models corroborate earlier predictions that extensional faults are open down to 800 – 1000 m (Gudmundsson and Bäckström, 1991). We also learned that below that, a hybrid failure zone exists where dilational jogs, open extensional fractures between shear mode faults, provide lateral pathways for magma or water, even at depths where the overburden stress prevents the formation of purely extensional faults.

 

Outlook

The previously shown experiments investigated MDF at a larger scale in purely dip-slip kinematics. However, faults at rifts often have strike-slip components, forming normal faults with oblique kinematics. In further experiments, we therefore explored the effect of varying basement fault obliquities, i.e. the range between dip-slip normal faults and strike-slip faults (Bitsch, 2017). As expected, early phases of faulting are dominated by Riedel shears. Surprisingly, the surface structure of mature faults, however, does not change distinctly up to obliquities of 60°, but the subsurface connectivity decreases with increasing obliquity.

Figure 5. Analogue model resembling successive layers of lava flows with cooling joints created by carefully stacking several layers of dried corn starch slurry (Winhausen, 2018). The dip of the basement fault is prescribed by the apparatus. The fault geometry generated in the model is very similar to the ones observed in Iceland. Large cavities develop and are partially refilled by loosened columns, as shown in figure 1. Tensile fractures develop on the surface of the footwall, similar to the hemihydrate model and the field. Credit: Lisa Winhausen.

Zooming in on the faults, an inherent mechanical anisotropy (orthotropy) of basalts gains more influence on the macroscale structure of faults. Due to the shrinking during cooling of flood basalts, polygonal to blocky columns form and present regular weak zones in the rockmass. Introducing mechanical anisotropy into a stronger modelling material (dried corn-starch slurry) beautifully illustrates how the small-scale structure of the faults is affected by the layering of flood basalts and cooling fractures therein (Fig. 5; Winhausen, 2018). Close to the surface the strong basalt does not fracture, but propagating faults rather localize at the pre-existing cooling joints. This causes a jagged structure of the fault, formation of caves, and eroded basalt columns filling the opening fractures.

We are currently working on implementing all these learning points into discrete element simulations, where we can adjust material properties in a way that allows for modelling deeper parts of the faults with better mechanical control.

 

 

References

Bitsch, N.D., 2017. Massively dilatant faults in oblique rift settings – an analogue modeling study (MSc Thesis). RWTH Aachen University, Germany, Aachen.

Bücken, D.H., 2017. Effect of mechanical stratigraphy on normal fault evolution – Insights from analogue models and natural examples in Iceland (MSc Thesis). RWTH Aachen University.

Gudmundsson, A., Bäckström, K., 1991. Structure and development of the Sveinagja graben, Northeast Iceland. Tectonophysics 200, 111–125. https://doi.org/10.1016/0040-1951(91)90009-H

Holland, M., Urai, J.L., Martel, S., 2006. The internal structure of fault zones in basaltic sequences. Earth Planet. Sci. Lett. 248, 301–315. https://doi.org/10.1016/j.epsl.2006.05.035

Kettermann, M., von Hagke, C., van Gent, H.W., Grützner, C., Urai, J.L., 2016. Dilatant normal faulting in jointed cohesive rocks: a physical model study. Solid Earth 7, 843–856. https://doi.org/10.5194/se-7-843-2016

van Gent, H.W., Holland, M., Urai, J.L., Loosveld, R., 2010. Evolution of fault zones in carbonates with mechanical stratigraphy – Insights from scale models using layered cohesive powder. J. Struct. Geol. 32, 1375–1391. https://doi.org/10.1016/j.jsg.2009.05.006

Winhausen, L., 2018. Influence of columnar joints on normal fault geometry and evolution An analog modeling study Master Thesis (MSc thesis). RWTH Aachen University.

Minds over Methods: Tectonochemistry of Melting Mud in the Mantle, evidence from the Oman/UAE ophiolite

Minds over Methods: Tectonochemistry of Melting Mud in the Mantle, evidence from the Oman/UAE ophiolite

For this first Minds over Methods of 2019, we invited Christopher Spencer, Senior Research Fellow at Curtin University in Australia, to tell us something about tectonochemistry. By applying geochemistry to tectonic processes, it is possible to get more insight into the different stages of the rock cycle. By combining fieldwork and geochemical analyses of the Oman/UAE ophiolite, Chris and his co-workers believe they found the first direct and in-situ evidence of sediment melting in the mantle.

 

Credit: Christopher Spencer

Tectonochemistry of Melting Mud in the Mantle, evidence from the Oman/UAE ophiolite

Christopher Spencer, Senior Research Fellow, Curtin University, Australia

The rock cycle is the first thing we learn in Geology 101. Magma and lava cool to form igneous rocks. Igneous rocks then erode to form sediment, which forms sedimentary rocks as it is compacted. Increasing pressure and heat then create metamorphic rocks, which eventually will melt. In each of the transitions described in the rock cycle, tectonics is usually involved. Granite batholiths form in subduction zones and are uplifted and eroded in collision zones. The sediments derived therefrom are deposited along continental margins that are often then returned to subduction zones where they contribute to new magmatic systems. There is a wide array of tools that we can use to evaluate the role of tectonics in the rock cycle, of which geochemistry is able to provide insight into each stage of the process.

Applying geochemistry to tectonics is (unsurprisingly) referred to as tectonochemistry. Similar to tectonophysics, where geophysics is applied to address large-scale tectonic questions, tectonochemistry provides a unique view into geochemical proxies of tectonic processes. The melting of sediment along convergent margins is a classic tectonochemical problem, as the unique chemical signature of sediment found in a granite provides unequivocal evidence for the melting of a sedimentary rock. In collisional systems, like the Himalaya, tectonochemistry has been used to constrain the melting of meta-sedimentary rocks as crustal thickening and decompression drives dehydration of micas which leads to melting. Collisional systems provide clear and in situ evidence for sediment melting.

Figure 1: Clockwise from top left: tourmaline-bearing leucogranite from the Himalaya in NW India, leucogranite dykes intruding meta-sedimentary rocks exposed at 5000m altitude, in situ melting of meta-pelite and formation of leucogranite, incongruent melting of muscovite + plagioclase + quartz to form leucogranite but leaving the biotite behind. Credit: Christopher Spencer.

 

Sediments are also thought to melt in subduction systems, but given the difficulty of accessing the asthenosphere directly, it is more challenging to constrain the processes occurring deep in a subduction zone. The incorporation of sediment in subduction zones is often constrained using the geochemistry of the resulting magmatic rocks. The chemical signature of sediment provides a clear indication of its incorporation in the magma, but it is often unclear whether the contamination is occurring in the asthenospheric wedge or in the upper crust. For example, many granite batholiths contain zircon grains that are foreign to the host magma and whose age spectra match the detrital zircon age spectra of the adjacent sedimentary units. This relationship is a clear indication that sedimentary contamination occurred in the upper crust. Unfortunately, the geochemical proxies used to establish the sedimentary contamination only provide indirect evidence for the subduction of sedimentary material into the asthenospheric wedge. Such indirect evidence includes seismic stratigraphy showing sedimentary units being subducted beneath the forearc and whiffs of sedimentary geochemical signals in arc volcanics. Although these evidences point towards sediment being subducted deep into the asthenospheric wedge where it melts and contaminates the magmas coming off the subducting slab, they do not preserve direct evidence of sediment melting in the mantle.

To acquire direct evidence of processes happening deep in the mantle, I set my sights on the Oman/UAE ophiolite, where a thick succession of mantle peridotite is preserved beneath a complete stratigraphic section of oceanic crust. Previous work has shown that this ophiolite not only preserves an intact record of oceanic crustal stratigraphy, but also geochemical features of a subduction zone in the oceanic crust. This implies the ophiolite formed in a supra-subduction setting, where during the earliest phase of subduction, extension in the upper plate caused rifting and formation of oceanic crust above a subduction zone.

 

Figure 2: Oceanic crustal stratigraphy of the Oman/UAE ophiolite comprised of (clockwise from top left): pillow basalts, sheeted dykes, layered gabbros, and mantle peridotite. Credit: Christopher Spencer.

 

During fieldwork in the ophiolite, while traversing the 8-15 km thickness of the mantle peridotite, I encountered a number of granitoid dykes that cross cut the peridotite, but do not cross the petrologic Moho. Many of these dykes contained tourmaline, muscovite, biotite, and even andalusite, minerals that would be expected from the melting of sedimentary material. Finding these minerals in the mantle indicates these grantoid dykes formed from the melting of sedimentary material and here they were within the mantle! Subsequent analysis of zircon grains from these granitoid dykes revealed the age of these dykes was equivalent to the age of the overlying ophiolite providing bullet-proof evidence that they intruded while the ophiolite was forming above a subduction zone. To provide the nail in the coffin for a sedimentary origin, I performed oxygen isotope analysis of the zircon and quartz. Sedimentary material has a distinct oxygen isotopic composition and igneous rocks that are thought to have experienced sediment contamination have δ18O values that lie along mixing lines between a sediment end member and the mantle. The oxygen isotopic analyses of the sub-Moho granitoids of the Oman/UAE ophiolite revealed the highest δ18O values ever measured in igneous rocks, providing unequivocal evidence that these granitoids represent pure sediment melts. In a paper published in Geology (Spencer et al., 2017), my coauthors and I argue these igneous rocks represent the first direct and in situ evidence of sediment melting in the mantle. Lucky for us, we have just scratched the surface of the exciting things left to learn about these fascinating granitoids and I look forward to the opportunity to return to the Oman/UAE ophiolite.

Figure 3: Sub-Moho granitoids of the Oman/UAE ophiolite: A) Cathodoluminescence image of a zircon shown with location and result of δ18O analyses. B) Photograph of sub-Moho granitoids. C) Hand sample of granite with tourmaline and lepidolite (lithium-bearing mica). Credit: Christopher Spencer.

 

Minds over Methods: What controls the shape of oceanic ridges?

Minds over Methods: What controls the shape of oceanic ridges?

In this edition of Minds over Methods, Aurore Sibrant, postdoc at Bretagne Occidentale University (France) explains how she studies the shape of oceanic ridges, and which parameters are thought to control this shape. By using laboratory experiments combined with observations from nature, she gives new insights into how spreading rates and lithosphere thickness influence the development of oceanic ridges. 

 

Credit: Aurore Sibrant

What controls the shape of oceanic ridges? Constraints from analogue experiments

Aurore Sibrant, Post-doctoral fellow at Laboratoire Géosciences Océans, Bretagne Occidentale University, France

Mid-oceanic ridges with a total length > 70 000 km, are the locus of the most active and voluminous magmatic activity on Earth. This magmatism directly results from the passive upwelling of the mantle and decompression melting as plates separate along the ridge axis. Plate separation is taken up primarily by magmatic accretion (formation of the oceanic crust), but also by tectonic extension of the lithosphere near the mid-ocean ridge, which modifies the structure of the crust and morphology of the seafloor (Buck et al., 2005). Therefore, the morphology of the ridge is not continuous but dissected by a series of large transform faults (> 100 km) as well as smaller transform faults, overlapping spreading centres and non-transform offsets (Fig. 1). Altogether, those discontinuities form the global shape of mid-ocean ridges. While we understand many of the basic principles that govern ridges, we still lack a general framework for the governing parameters that control segmentation across all spreading rates and induce the global shape of ridges.

Geophysical (Schouten et al., 1985; Phipps Morgan and Chen, 1993; Carbotte and Macdonald, 1994) and model observations (Oldenburg and Brune, 1975, Dauteuil et al., 2002, Püthe and Gerya, 2014) suggest that segmentation of oceanic ridges reflects the effect of spreading rate on the mechanical properties and thermal structure of the lithosphere and on the melt supply to the ridge axis. To understand the conditions that control the large-scale shape of mid-ocean ridges, we perform laboratory experiments. By applying analogue results to observations made on Earth, we obtain new insight into the role of spreading velocity and the mechanical structure of the lithosphere on the shape of oceanic ridges.

 

Laboratory experiments

The analogue experiment is a lab-scale, simplified reproduction of mid-oceanic ridges system. Our set-up yields a tank filled from bottom to top by a viscous fluid (analogous to the asthenosphere) overlain by the experimental “lithosphere” that can adopt various rheologies and a thin surface layer of salted water. This analogue lithosphere is obtained using a suspension of silica nanoparticles which in contact with the salted water emplaced on the surface of the fluid causes formation of a skin or “plate” that grows by diffusion. This process is analogous to the formation of the oceanic lithosphere by cooling (Turcotte and Schubert, 1982). With increasing salinity, the rheology of the skin evolves from viscous to elastic and brittle behaviour (Di Giuseppe et al., 2012; Sibrant and Pauchard, 2016).

The plate is attached to two Plexiglas plates moving perpendicularly apart at a constant velocity. The applied extension nucleates fractures, which rapidly propagate and form a spreading axis. Underlying, less dense, fresh fluid responds by rising along the spreading axis, forming a new skin when it comes into contact with the saline solution. By separately changing the surface water salinity and the velocity of the plate separation, we independently examine the role of spreading velocity and axial lithosphere thickness on the evolution of the experimental ridges.

 

Figure 2. Close up observations of analogue mid-oceanic ridges and schematic interpretation for different spreading velocity. The grey region is a laser profile projected on the surface of the lithosphere: the laser remains straight as long as the surface is flat. Here, the large deviation from the left to centre of the image reveals the valley morphology of the axis. Credit: Aurore Sibrant.

 

Analogue mid-oceanic ridges

Over a large range of spreading rates and salinities (Sibrant et al., 2018), the morphology of the axis is different in shape. The ridge begins with a straight axis (initial condition). Then during the experiment, mechanical instabilities such as non-transform offset, overlapping spreading centres and transform faults develop (Fig. 2) and cause the spreading axis to have a non-linear geometry (Fig. 3). A key observation is the variation of the shape of the analogue ridges with the spreading rate and salinities. For similar salinity and relative slow spreading rates, each segment is offset by transform faults shaping a large tortuous ridge (i.e. non-linear geometry). In contrast, at a faster spreading rate, the ridge axis is still offset by mechanical instabilities but remains approximately linear.

Figure 3. Ridge axis morphology observed in the experiments and schematic structural interpretations of the ridge axis, transform faults (orange ellipsoids) and non-transform faults (purple ellipsoids). Measurements of lateral deviation (LD) correspond to the length of the arrows. For comparison, white squares represent the size of closeup shows in Fig 2. Credit: Aurore Sibrant.

We can quantify the ridge shape by measuring the total lateral deviation, which is the total accumulated offset of the axis, when the tortuosity amplitude becomes stable. For cases with similar salinities, the results indicate two trends. First, the lateral deviation is high at slow spreading ridges and decreases within increasing spreading rate until reaching a minimum lateral deviation value for a given critical spreading rate (Fig 4A). Then the lateral deviation remains constant despite the increasing spreading rate. Experiments with different salinities also present a transition between tortuous and linear ridges. These two trends reflect how the lithosphere deforms and fails. In the first regime, the axial lithosphere is thick and is predominantly elastic-brittle. In such cases, the plate failures occur from the surface downwards through the development of faults: it is a fault-dominated regime. In contrast, for faster spreading rate or smaller salinities, the axial lithosphere is thin and is predominantly plastic. Laboratory inspection indicates that fractures in plastic material develop from the base of the lithosphere upwards: it is a fluid-intrusion dominated regime.

 

 

Comparison with natural mid-oceanic ridge

In order to have a complete understanding of the mid-oceanic ridge system, it is essential to compare the laboratory results with natural examples. Hence, we measure the lateral deviation of nature oceanic ridges along the Atlantic, Pacific and Indian ridges. The measurements reveal the same two regimes as found in laboratory data. The remaining step consists of finding the appropriate scaling laws to superpose the natural and experiment data. This exercise requires dynamics similarity between analogue model and real-world phenomena which is demonstrated using dimensionless numbers (Sibrant et al., 2018). Particularly, the “axial failure parameter – πF” describes the predominant mechanical behaviour of the lithosphere relative to its thickness. Low-πF accretion is dominated by fractures in a predominantly elastic-brittle lithosphere: the lateral deviation of the ridges is tortuous, while at higher pF, accretion is dominated by intrusion in a predominantly plastic lithosphere: the shape of the mid oceanic ridges is mostly linear (Fig 4B).

 

Figure 4. (A) Lateral deviation values measured in the experiments in function of the spreading rate velocities and salinities. (B) Evolution of the lateral deviation of the ridge axis, normalized by the critical axial thickness (Zc) relative to the axial failure parameter. Dark grey is the laboratory experiments and the colored circles are the Earth data. Adapted from Sibrant et al., 2018.

 

Our experiments give insight into the role of axial failure mode (fault-dominated or intrusion-dominated) on the shape of mid-oceanic ridges. In the future, we want to use this experimental approach to investigate the origin of mechanical instabilities, such as transform faults or overlapping spreading centres. This experimental development and results are a collaborative work between Laboratoire FAST at Université Paris-Saclay and Department of Geological Sciences at the University of Idaho and involves E. Mittelstaedt, A. Davaille, L. Pauchard, A. Aubertin, L. Auffray and R. Pidoux.

 

 

References
Buck, W.R., Lavier, L.L., Poliakov, A.N.B., 2005. Modes of faulting at mid-ocean ridges. Nature 434, 719-723.
Schouten, H., Klitgord, K.D., Whitehead, J.A., 1985. Segmentation of mid-ocean ridges. Nature 317, 225-229.
Carbotte, S.M., Macdonald, K. C., 1994. Comparison of seafloor tectonic fabric at intermediate, fast, and super fast spreading ridges: Influence of spreading rate, plate motions, and ridge segmentation on fault patterns. J. Geophys. Res. 99, 13609-13631.
Phipps Morgan, J., Chen, J., 1993. Dependence of ridge-axis morphology on magma supply and spreading rate. Nature 364, 706-708.
Oldenburg, D.W., Brune, J.N., 1975. An explanation for the orthogonality of ocean ridges and transform faults. J. Geophys. Res. 80, 2575-2585.
Dauteuil, O., Bourgeois, O., Mauduit, T., 2002. Lithosphere strength controls oceanic transform zone structure: insights from analogue models. Geophys. J. Int. 150, 706-714.
Püthe, C., Gerya, T., 2014. Dependence of mid-ocean ridge morphology on spreading rate in numerical 3-D models. Gondwana Res. 25, 270-283.
Turcotte, D., Schubert, G., Geodynamics (Cambridge Univ. Press, New York, 1982).
Di Giuseppe, E., Davaille, A., Mittelstaedt, E., Francois, M., 2012. Rheological and mechanical properties of silica colloids: from Newtonian liquid to brittle behavior. Rheologica Acta 51, 451-465.
Sibrant, A.L.R., Pauchard, L., 2016. Effect of the particle interactions on the structuration and mechanical strength of particulate materials. European Physics Lett., 116, 4, 10.1209/0295-5075/116/49002.
Sibrant, A.L.R., Mittelstaedt, E., Davaille, A., Pauchard, L., Aubertin, A., Auffray, L., Pidoux, R., 2018. Accretion mode of oceanic ridges governed by axial mechanical strength. Nature Geoscience 11, 274-279.

 

Minds over Methods: Linking microfossils to tectonics

Minds over Methods: Linking microfossils to tectonics

This edition of Minds over Methods article is written by Sarah Kachovich and discusses how tiny fossils can be used to address large scale tectonic questions. During her PhD at the University of Brisbane, Australia, she used radiolarian biostratigraphy to provide temporal constraints on the tectonic evolution of the Himalayan region – onshore and offshore on board IODP Expedition 362. Sarah explains why microfossils are so useful and how their assemblages can be used to understand the history of the Himalayas. And how are new technologies improving our understanding of microfossils, thus advancing them as a dating method?

 

                                                                          Linking microfossils to tectonics

Credit: Sarah Kachovich

Sarah Kachovich, Postdoctoral Researcher at the School of Earth and Environmental Sciences, The University of Queensland, Australia.

Radiolarians are single-celled marine organisms that have the ability to fix intricate, siliceous skeletons. This group of organism have captured the attention of artist and geologist alike due to their skeletal diversity and complexity that can be observed in rocks from the Cambrian to the present. As a virtue of their silica skeletons, small size and abundance, radiolarian skeletons can potentially exist in most fine-grained marine deposits as long as their preservation is good. This includes mudstones, hard shales, limestones and cherts. To recover radiolarians from a rock, acid digestion is commonly required. For cherts, 12-24 hours in 5 % hydrofluoric acid is needed to liberate radiolarians. Specimens are collected on a 63 µm sieve and prepared for transmitted light or scanning electron microscope analysis.

Animation of radiolarian diversity. Credit: Sarah Kachovich

Scale and diversity of modern radiolarians. Credit: Sarah Kachovich (radiolarians from IODP Expedition 362) and Adrianna Rajkumar (hair).

 

 

 

 

 

 

 

 

 

 

Improving the biostratigraphical potential of radiolarians

The radiolarian form has changed drastically through time and by figuratively “standing on the shoulders of giants”, we correlate forms from well-studied sections to determine an age of an unknown sample. A large effort of my PhD was aimed to progress, previously stagnant, research in radiolarian evolution and systematics in an effort to improve the biostratigraphical potential of spherical radiolarians, especially from the Early Palaeozoic. The end goal of this work is to improve the biostratigraphy method and its utility, thus increasing our understanding of the mountain building processes.

The main problem with older deposits is the typical states of preservation, where radiolarians partly or totally lose their transparency, which makes traditional illustration with simple transmitted light optics difficult. Micro-computed tomography (µ-CT) has been adopted in fields as diverse as the mineralogical, biological, biophysical and anatomical sciences. Although the implementation in palaeontology has been steady, µ-CT has not displaced more traditional imaging methods, despite its often superior performance.

Animation of an Ordovician radiolarian skeleton in 3D imaged through µ-CT. Credit: Sarah Kachovich

To study small complex radiolarian skeletons, you need to mount a single specimen and scan it at the highest resolution of the µ-CT. The µ-CT method is much like a CAT scan in a hospital, where X-rays are imaged at different orientations, then digitally stitched together to reconstruct a 3D model. The vital function of the internal structures provides new insights to early radiolarian morphologies and is a step towards creating a more robust biostratigraphy for radiolarians in the Early Paleozoic.

Linking radiolarian fossils to tectonics

Radiolarian chert is important to Himalayan geologists as it provides a robust tool to better document and interpret the age and consumption of oceanic lithosphere that once intervened India and Asia before their collision.The chert that directly overlies pillow basalt in the ophiolite sequence (remnant oceanic lithosphere) represents the minimum age constraint of its formation. In the Himalayas, over 2000 km of ocean has been consumed as India rifted from Western Australia and migrated north to collide with Asia. Only small slivers of ophiolite and overlying radiolarian cherts are preserved in the suture zone and it is our job to determine how these few ophiolite puzzle pieces fit together.

Another way I have been able to link microfossils to Himalayan tectonics is by studying the history and source of erosion from the Himalayas on board IODP Expedition 362. Sedimentation rates obtained from deep sea drilling can provide ages of various tectonic events related to the India-Asia collision. For example, we were able to date various events such as the collision of the Ninety East Ridge with the Sumatra subduction zone, which chocked off the sediment supply to the Nicobar basin around 2 Ma as the ridge collided with the subduction zone.

Left: Results from the McNeill et al. (2017) of the sedimentation history of Bengal Fan (green dots) and Nicobar Fan (red dots). Middle/right: Reconstruction of India and Asia for the past 9 million years showing the sediment source from the Himalayas to both basins on either side of the Ninety East Ridge.

 

 

 

 

 

 

 

 

 

 

 

 

Lastly, on board Expedition 362 we were able to use microfossils to understand how and why big earthquakes happen. We targeted the incoming sediments to the Sumatra subduction zone that were partly responsible for the globally 3rd largest recorded earthquake (Mw≈9.2). This earthquake occurred in 2004 and produced a tsunami that killed more than 250,000 people.

From the seismic profiles (see example below), we found that the seismic horizon where the pre-decollement formed coincided with a thick layer of biogenetically rich sediment (e.g. radiolarians, sponge spicules, etc.) found whilst drilling. Under the weight of the overlying Nicobar Fan sediments, this critical layer of biogenic silica is undergoing diagenesis and fresh water is being chemically released into the sediments. The fresh water within these sediments is moving into the subduction zone where it has implications to the physical properties of the sediment and the morphology of the forearc region.

The Sumatra subduction zone. The dark orange zone represents the rupture area of the 2004 earthquake. Also shown are the drill sites of IODP Expedition 362 and the location of seismic lines across the plate boundary.

Seismic profile: The fault that develops between the two tectonic plates (the plate boundary fault) forms at the red dotted line. Note the location of the drill site.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

References

Hüpers, A., Torres, M. E., Owari, S., McNeill, L. C., Dugan, & Expedition 362 Scientists, 2017. Release of mineral-bound water prior to subduction tied to shallow seismogenic slip. Science, 356: 841–844. doi:10.1126/science.aal3429

McNeill, L. C., Dugan, B., Backman, J., Pickering, K. T., & Expedition 362 Scientists 2017. Understanding Himalayan Erosion and the Significance of the Nicobar Fan. Earth and Planetary Science Letters, 475: 134–142. doi:10.1016/j.epsl.2017.07.019

Minds over Methods: Experimental seismotectonics

Minds over Methods: Experimental seismotectonics

For our next Minds over Methods, we go back into the laboratory, this time for modelling seismotectonics! Michael Rudolf, PhD student at GFZ in Potsdam (Germany), tells us about the different types of analogue models they perform, and how these models contribute to a better understanding of earthquakes along plate boundaries.

 

Credit: Michael Rudolf

Experimental seismotectonics – Seismic cycles and tectonic evolution of plate boundary faults

Michael Rudolf, PhD student at Helmholtz Centre Potsdam – German Research Centre for Geosciences GFZ

The recurrence time of large earthquakes that happen along lithospheric-scale fault zones such as the San Andreas Fault or Chile subduction megathrusts, is very long (≫100 yrs.) compared to human timescales. The scarcity of such events over the instrumental record of around 60 years is unfortunate for a statistically sound analysis of the earthquake time series.

So far, only few megathrust events have been monitored in detail with near-field seismic and geodetic networks. To circumvent this lack of observational data, we at Helmholtz Tectonic Laboratory use analogue modelling to understand plate boundary faulting on multiple time-scales and the implications for seismic hazard. We use models of strike-slip zones and subduction zones, to investigate several aspects of the seismic cycle. Additionally, numerical simulations accompany and complement each experimental setup using experimental parameters.

 

Seismotectonic scale models
In my project, we develop experiments that can model multiple seismic cycles in strike-slip conditions. Our study employs two types of experimental setups both are using the same materials. The first is simpler (ring shear setup) and is able to show the on-fault rupture propagation. The second is geometrically more similar to the natural system, but only the surface deformation is observable.

To model rupture propagation, we introduce deformable sliders in a ring shear apparatus. Two cylindrical shells of ballistic gelatine (Ø20 cm), representing the side walls, rotate against each other, with a thin layer (5 mm) of glass beads (Ø355-400µm) in between representing an annular fault zone. A see-through lid connected to force sensors holds the upper shell in place, whereas the machine rotates the lower shell. Through the transparent lid and upper shell, we directly observe the fault slip. We can vary the normal stress on the fault (<20 kPa) and the loading velocity (0.0005 – 0.5 mm/s).

The next stage of analogue model, features depth-dependent normal stress and a rheological layering mimicking the strike-slip setting in the uppermost 25-30 km of the lithosphere (see also Mehmet Köküm’s blog post). A gelatine block (30x30cm) compressed in uniaxial setting represents the elastic upper crust under far-field forcing. Embedded in the block is a thin fault filled with quartz glass beads. The ductile lower crust is modelled using viscoelastic silicone oil. The model floats in a tank of dense sugar solution, to guarantee free-slip, stress-free boundaries.

 

Figure 2 – Setup and monitoring technique during an experiment. Several cameras record the displacement field and the ring shear tester records the experimental results. Credit: Michael Rudolf

 

Analogue earthquakes
Both setups generate regular stick-slip cycles including minor creep. Long phases of quiescence, where no slip or very slow creep occurs, alternate with fast slip events sometimes preceded by slow slip events. The moment magnitude of analogue earthquake events is Mw -7 to -5. The cyclic recurrence of slip events is an analogue for the natural seismic cycle of a single-fault system.

 

Figure 3 – Detailed setup and results from the ring shear tester experiments. The upper right image shows a snapshot of an analogue earthquake rupture along the fault zone. The plot shows the recorded shear forces and slip velocities over one hour of experiment. Credit: Michael Rudolf

 

Optical cameras record the slip on the fault and the deformation of the sidewalls. Using digital image correlation techniques, we are able to visualize accurately deformations on the micrometre scale at high spatial and temporal resolution. Accordingly, we can verify that analogue earthquakes behave kinematically very similar to natural earthquakes. They generally nucleate where shear stress is highest, and then propagate radially until the seismogenic width is saturated. In the end, the rupture continues laterally along the fault strike. Our experiments give insight into the role of viscoelastic relaxation, interseismic creep, and slip transients on the recurrence of earthquakes, as well as the control of loading conditions on seismic coupling and rupture dynamics.

 

Figure 4 – Setup and Results for the strike-slip geometry. The surface displacement field is very similar to natural earthquakes. The plot shows that due to technical limitations of this setup, fewer events are recorded but the slip velocities are higher. Credit: Michael Rudolf

 

Future developments
Together with our partners in the Collaborative Research Centre (CRC1114 – Scaling Cascades in Complex Systems) we employ a new mathematical and numerical description of the fault system, to simulate our experiments and get a physical understanding of the empirical friction laws. In the future, we want to use this multiscale spatial and temporal approach to model complex fault networks over many seismic cycles. The experiments serve as benchmarks and cross-validation for the numerical code, which in the future will be using natural parameters to get a better geological and mathematical understanding of earthquake slip phenomena and occurrence patterns in multiscale fault networks.

Minds over Methods: Reconstructing oceans lost to subduction

Minds over Methods: Reconstructing oceans lost to subduction

Our next Minds over Methods article is written by Derya Gürer, who just finished a PhD at Utrecht University, the Netherlands. During her PhD, she used a combination of many methods to reconstruct the evolution of the Anadolu plate, which got almost entirely lost during closure of the Neotethys in Anatolia. Here, she explains how the use of these multiple methods helped her to obtain a 3D understanding of the Anatolian double subduction system and the demise of the Anadolu plate. 

Credit: Derya Gürer

Reconstructing oceans lost to subduction

Derya Gürer, Postdoctoral Researcher, Utrecht University, the Netherlands. 

Subduction represents the single biggest recycling process on Earth and takes place at convergent plate boundaries. One plate subducts underneath another into the mantle, generating volcanism, earthquakes, tsunamis and associated hazards. Subduction zones come and go, and nearly half of the subduction zones active today formed in the Cenozoic (after ~65 Ma) (Gurnis et al., 2004). The negative buoyancy of subducted lithosphere (‘slab pull’) is thought to be the major driver of plate tectonics (Turcotte and Schubert, 2014). Changes in the configuration of subduction zones thus change the driving forces of plate tectonics, making the reconstruction of the kinematic evolution of subduction key to understanding past plate motions. Such reconstructions make use of data preserved in the modern oceans (marine magnetic anomalies and fracture zone patterns). But because subduction is a destructive process, the surface record of subduction-dominated systems is naturally incomplete, and more so backwards in time. Sometimes, relicts of subducted lithosphere are preserved in active margin mountain belts, holding valuable information to restore past plate motions and the dynamic evolution of subduction zones.

But how does one recognize a plate that has been almost entirely lost to subduction? And how do we reconstruct the evolution of subduction zones through space and time?

 

Archives of plates that were (almost) lost due to subduction

Subduction occurs in a variety of geometries and leaves behind a distinct geological record that holds key elements for the analysis of the past kinematics of now-subducted plates. Where subduction occurred below oceanic lithosphere, fragments of the leading edge of this overriding lithosphere may be left behind as remnants of oceanic crust (ophiolites). Subduction of oceanic plates may also be associated with accretion of its volcano-sedimentary cover to the overriding plate as an accretionary complex (Matsuda and Isozaki, 1991). Forearc basins associated with intra-oceanic subduction zones form on top of ophiolites and accretionary complexes and may record permanent deformation (syn-kinematic) of the overriding plate in response to tectonic interaction with the down-going plate (e.g., accretion, subduction erosion, slab roll-back) (Fig. 1).

Fig. 1: The location of archives of the evolution of “lost” oceanic plates (ophiolites, accretionary complexes, forearc basins) in a subduction zone setting.Credit: Derya Gürer.

The sedimentary infill of forearc basins implicitly records the nature and stress state of the overriding plate. Forearc basins may therefore hold the most complete record of the motion of the oceanic plate relative to the trench. However, many accretionary complexes and forearcs are deeply submerged and buried below sediments, making them highly inaccessible, and therefore expensive to study. As a consequence, our understanding of such systems is primarily based on well-studied examples in the East Pacific (e.g. Franciscan Complex, California (Wakabayashi, 2015)). Other such systems exist in the Mediterranean realm – for example in the geological record of Anatolia. The unique and direct archive of past plate motion in the geological record of Central and Eastern Anatolia is independent from constraints provided by marine magnetic anomalies, and provides a key region to unravel the evolution of destructive plate boundaries.

 

How many oceans were lost in Anatolia?

Fig. 2: The multidisciplinary approach used in my PhD research consisted of structural field analysis and stratigraphy of Anatolian sedimentary basins with focus on syn-kinematic deformation (top) with time constraints provided by absolute age dating of accessory minerals and biostratigraphy (middle). Paleomagnetic analysis (bottom left) provided information about vertical axis rotations. The combined information from these methods were integrated in a kinematic reconstruction and tested against mantle tomography (bottom right). Credit: Derya Gürer

To answer this question, I studied the deformation of sedimentary basins overlying Anatolian ophiolites (remnants of oceanic crust), and the deformation record of rocks which were buried and exhumed below these ophiolites. The Cenozoic deformation of the Anatolian orogen allowed for identifying the timing of arrest of the subduction history and revealed the simultaneous activity of two subduction zones in Late Cretaceous time. These two subduction zones bound a separate oceanic plate within the Neotethys Ocean – the Anadolu Plate (Fig. 3, Gürer et al., 2016). The aim of my PhD research was to reconstruct the birth, evolution and destruction of this oceanic plate.

Tectonic problems require a multidisciplinary approach, in order to study the evolution of orogens and associated sedimentary basins. My research involved the integration of (1) structural analysis, (2) stratigraphy, (3) geochronology, (4) paleomagnetism, (5) plate reconstruction, and (6) mantle tomography (Fig. 2). The main goal was to obtain new data on the evolution of the Central and Eastern Anatolian regions through the analysis of spatial and temporal relationships of deformation archived in the geological record.

First, I collected kinematic data from sedimentary basins (Fig. 2) overlying ophiolitic relicts of the oceanic Anadolu Plate, as well as from the underlying accretionary complex (Gürer et al., 2018a). Here, it was especially useful to focus on syn-kinematic deformation recorded by sediments. To constrain the timing of this deformation, I used geochronological data coming from absolute age dating and biostratigraphy. The integrated reconstruction of the kinematic history of basins was used to develop concepts quantitatively constraining the tectonic history of the Anadolu Plate and its surrounding trenches in 2D (Gürer et al., 2016).

 

Fig. 3: The Ulukışla Basin (Central Anatolia) represents a forearc basin in Late Cretaceous to Eocene time which recorded the evolution of the Anadolu Plate. The basin has subsequently been strongly deformed during Eocene and younger collisional processes and is juxtaposed against the Aladağ range along the Ecemiş Fault. Credit: Derya Gürer.

 

There are, however, large vertical axis rotations constrained through paleomagnetic analysis within Anatolia, not taken into account in the workflow described in the previous paragraph. Therefore, paleomagnetic data from the Late Cretaceous to Miocene sedimentary basins were collected. These data identified coherently rotating domains and major tectonic structures that accommodated differential rotations between tectonic blocks (Gürer et al., 2018b).

Fig. 4: Simplified interpretation of the Late Cretaceous double subduction geometry in Anatolia and the Anadolu Plate.Credit: Derya Gürer.

Subsequently, a kinematic reconstruction of Anatolia back to the Late Cretaceous was built (Fig. 4) incorporating the timing of deformation obtained through structural analysis, stratigraphy, geochronology, and vertical axis rotations. This reconstruction provided first-order implications for the timing and geometry of subduction zones and revealed that the demise of the Anadolu Plate and collision in Anatolia was variable along the strike of the orogen, younging from the west to the east. The exact timing of collision in Eastern Anatolia will require future studies applying structural field geology, systematic analysis of the age and nature of magmatism, and thermochronology to constrain timing of regional exhumation, as well as detrital geochronology, providing information on the relative proximity of tectonic blocks through the provenance of sediments.

 

Finally, the resulting 2D kinematic reconstruction was tested against a mantle tomographic model (UU-07, Amaru, 2007; van der Meer et al., 2017) to gain insights into its 3D geometry. Mantle tomography images the present-day structure and positive seismic anomalies (blue colours in Fig. 5), which may be interpreted as subducted slabs. Comparing the convergence estimate obtained from the kinematic reconstruction with the imaged subducted lithosphere allowed to infer that the mantle structure in the Eastern Mediterranean holds record of not only the two strands of the Neotethys Ocean that existed in Anatolia, but also of the Paleotethys Ocean.

 

Fig. 5: Map view tomographic structure below the Eastern Mediterranean region at variable depths (increasing in depth from left to right). Blue colours generally represent positive, whereas red colours represent negative wave speed anomalies. Credit: Derya Gürer & Wim Spakman.

The combination of methods to unravel the geological record of Anatolia quantitatively constrained the evolution of subduction zones and of the Anadolu Plate. The reconstruction of the Anatolian double subduction system that existed in Late Cretaceous time has implications for the dynamics of multiple simultaneously active subduction zones.

 

References

Amaru, M.L., 2007. Global travel time tomography with 3-D reference models. PhD thesis, Utrecht University, The Netherlands.

Gürer, D., van Hinsbergen, D.J.J.D.J.J., Matenco, L., Corfu, F., Cascella, A., 2016. Kinematics of a former oceanic plate of the Neotethys revealed by deformation in the Ulukışla basin (Turkey). Tectonics 35, 2385–2416. https://doi.org/10.1002/2016TC004206

Gürer, D., Plunder, A., Kirst, F., Corfu, F., Schmid, S.M., van Hinsbergen, D.J.J., 2018a. A long-lived Late Cretaceous–early Eocene extensional province in Anatolia? Structural evidence from the Ivriz Detachment, southern central Turkey. Earth Planet. Sci. Lett. 481. https://doi.org/10.1016/j.epsl.2017.10.008

Gürer, D., Hinsbergen, D.J.J. van, Özkaptan, M., Creton, I., Koymans, M.R., Cascella, A., Langereis, C.G., 2018b. Paleomagnetic constraints on the timing and distribution of Cenozoic rotations in Central and Eastern Anatolia. Solid Earth 9, 1–27. https://doi.org/10.5194/se-9-1-2018

Gurnis, M., Hall, C., Lavier, L., 2004. Evolving force balance during incipient subduction. Geochemistry Geophys. Geosystems 5, Q07001. https://doi.org/10.1029/2003GC000681

Matsuda, T., Isozaki, Y., 1991. Well-documented travel history of Mesozoic pelagic chert in Japan: from remote ocean to subduction zone. Tectonics 10, 475–499.

van der Meer, D.G., van Hinsbergen, D.J.J., Spakman, W., 2018. Atlas of the Underworld: slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity. Tectonophysics 723, 309–448.

Wakabayashi, J., 2015. Anatomy of a subduction complex: architecture of the Franciscan Complex, California, at multiple length and time scales. Int. Geol. Rev. 37–41. https://doi.org/10.1080/00206814.2014.998728

 

Minds over Methods: Making ultramylonites

Minds over Methods: Making ultramylonites

“Summer break is over, which means we will continue with our Minds over Methods blogs! For this edition we invited Andrew Cross to write about his experiments with a new rock deformation device – the Large Volume Torsion (LVT) apparatus. Andrew is currently working as a Postdoctoral Research Associate in the Department of Earth and Planetary Sciences, Washington University in St. Louis, USA. He did his PhD at the University of Otago, New Zealand, although he is originally from the UK. His main research interest lies in understanding how micro-scale deformation processes influence the evolution of Earth’s lithosphere and tectonic plate boundaries. Hopefully we will be seeing more of him in the very near future” – Subhajit Ghosh.

Credit: Andrew Cross

Investigating strain-localisation processes in high-strain laboratory deformation experiments

Andrew Cross, Postdoctoral Research Associate at the Department of Earth and Planetary Sciences, Washington University in St. Louis, USA.

Below the upper few kilometres of the Earth’s surface – where rocks break and fracture under stress – elevated temperatures and pressures enable solid rocks to flow and bend, like a chocolate bar left outside on a warm day. This ductile flow of rocks and minerals plays a crucial role in many large-scale geodynamic processes, including mantle convection, the motion of tectonic plates, the flow of glaciers and ice sheets, and post-seismic and post-glacial rebound.

Fig. 1: Creep deformation occurs over very long timescales in the Earth. To replicate these processes on observable timescales, we must increase the rate of deformation in the laboratory. Credit: Andrew Cross

Unlike seismogenic slip that periodically accommodates large displacements over very short timescales, ductile flow occurs continuously, and at an almost imperceptibly slow rate: for example, rocks in the Earth’s interior creep at a rate roughly 10 billion times slower than that of the long-running pitch drop experiment1. Since few researchers are willing to wait millions of years to observe creep deformation in nature, we need ways of replicating these processes on much shorter timescales. Fortunately, by increasing temperature and the rate of deformation in the laboratory, we can generate creep behaviour in small samples of rock over timescales of a few hours, days, or weeks (Fig. 1).

In the Experimental Studies of Planetary Materials (ESPM) group at Washington University in St. Louis, we have spent the last couple of years developing a new rock deformation device – the Large Volume Torsion (LVT) apparatus (Fig. 2) – for performing torsion (twisting) experiments on geologic materials. By twisting small, disk-shaped rock samples, we are able to apply much more deformation (“strain”) than by squashing cylindrical samples end-on: this enables us to replicate deformation processes that operate in high-strain regions of the Earth (along the boundaries between tectonic plates, for instance).

Fig. 2: The Large Volume Torsion (LVT) apparatus. A 100-ton hydraulic ram applies a confining pressure, while electrical current passes through a graphite tube around the sample, generating heat through its electrical resistance. A screw actuator (typically used to raise and lower drawbridges) is used to rotate the lower platen and twist the sample, held between two tungsten-carbide anvils. Credit: Andrew Cross

Using the LVT apparatus, we are starting to investigate the microstructural and mechanical processes that lead to the formation of mylonites and ultramylonites: intensely deformed rocks that comprise the high-strain interiors of ductile shear zones and tectonic plate boundaries. It is widely thought that dramatic grain size reduction during (ultra)mylonite formation causes strain localisation, since strain-weakening deformation mechanisms (i.e., diffusion creep and grain boundary sliding) dominate at small grain sizes. However, grain size reduction (and therefore strain-weakening) is counteracted by the tendency of grains to grow over time, in the same way that bubbles in soapy water merge and grow over time.

An effective way of limiting grain growth is through “Zener pinning”, whereby the intermixing of grains of different mineral phases prevents grain boundary migration (and therefore growth). However, despite its suspected importance for ultramylonite formation and the occurrence of localised deformation on Earth (and possibly other planetary bodies), the processes leading to interphase mixing remain somewhat poorly understood and quantified.

Fig. 3: A comparison between our experimentally deformed calcite-anhydrite samples2 (backscattered electron (BSE) images), and natural metagranodiorite mylonites from Gran Paradiso, Western Alps3 (quartz grains, in black, mapped using electron backscatter diffraction (EBSD). Credit: Andrew Cross and Kilian et al., 2011.

To investigate phase mixing processes, we recently performed torsion experiments on mixtures of calcite and anhydrite. By deforming these mixtures to different amounts of strain, and then analysing the deformed samples in a scanning electron microscope, we were able to observe and quantify the evolution of deformation microstructures and mechanisms leading to ultramylonite formation. Backscattered electron (BSE) images show that clusters of the different minerals stretch out to form very thin, fine-grained layers, similar to foliation in natural shear zones (Fig. 3). At relatively large shear strains (17 < γ < 57) those layers disaggregated to form a fine-grained and homogeneously mixed aggregate. Electron backscatter diffraction (EBSD) analysis showed that calcite crystals became progressively more randomly oriented during phase mixing, indicative of a transition to the strain-weakening diffusion creep and grain boundary sliding regime.

The fact that a large amount of strain is required for phase mixing – and therefore strain-weakening – suggests that 1) only mature (highly-strained) shear zones are likely to maintain their weakness over long periods of geologic time, and 2) these features are therefore more likely to be reactivated after periods of quiescence. Inherited, long-lived mechanical weakness may well explain why tectonic plate boundaries are often reactivated over multiple cycles of continent accretion and rifting.

 

http://smp.uq.edu.au/content/pitch-drop-experiment

 Cross, A. J., & Skemer, P. (2017). Ultramylonite generation via phase mixing in high‐strain experimentsJournal of Geophysical Research: Solid Earth122(3), 1744-1759.

 3 Kilian, R., Heilbronner, R., & Stünitz, H. (2011). Quartz grain size reduction in a granitoid rock and the transition from dislocation to diffusion creepJournal of Structural Geology33(8), 1265-1284.

Minds over Methods: Block modeling of Anatolia

 

How can we use GPS velocities to learn more about present-day plate motions and regional deformation? In this edition of Minds over Methods, one of our own blogmasters Mehmet Köküm shares his former work with you! For his master thesis at Indiana University, he used block modeling to better understand the plate motion and slip rates of Anatolia and surrounding plates.

 

credit: Mehmet Köküm

Using block modeling to constrain present-day deformation of Anatolia and slip rates along the North Anatolian Fault

Mehmet Köküm, researcher at Firat University, Turkey

Until the late 1980’s, geological features such as offset of geomorphological markers were mainly used to determine historical slip rates along faults. Since the mid 1990’s, however, GPS has been widely used since it gives more accurate estimates of present-day slip rates by calculating strain accumulation at the crust. In this work, I use a GPS derived velocity field of Anatolia including data from 1988 to 2005 by Reilinger et al. (2006).

Turkey (Anatolian Plate) is located in the center of the Alpine fold and thrust belt. Due to the closure of different branches of the Neo-Tethys Ocean, main tectonic features of the Anatolian Plate are complicated by interactions between several tectonic plates.  The Arabian plate collides with the African plate in the south and the Eurasian plate in the north while the African plate subducts beneath the Anatolian plate along the Hellenic-Cyprus trench. As a result of these complex tectonic structures, the Anatolian plate displays various tectonic styles simultaneously.

Modeling and Data
Kinematic block modeling of interseismic surface motions has been used in different formats by several authors (e.g., McClusky et al. 2000; Westaway 2000; Barka and Reilingier 1997, 2006). The block modeling approach used here is described by Johnson and Fukuda (2010). In this study we used an elastic block model, which is a traditional block model that assumes no long-term deformation of the blocks. For simplicity, all faults are vertical, plates are considered as blocks and are assumed to be rigid. Block boundaries are defined from historic earthquakes, mapped faults and seismicity. Many of the major structures in Anatolia are well known except for a few submarine structures.

 

Map showing selected block model including of 14 blocks (or plates). Credit: Mehmet Köküm

 

Locking Depth
Locking depths indicate the depth for which a fault is completely locked above and creeping below. Estimates of these locking depths are output of the modeling studies and should correlate with the depth of major earthquakes along related faults. Meade and Hager (2005) suggest that there is a relation between locking depth and fault slip rates. Shallower locking depths correlate with slower slip rate estimates; therefore, GPS velocities near locked faults have slower velocities (Reilinger et al., 2006).

 

Elastic-half-space model showing fault creep at surface, locked (nonslipping) fault at depth, and freely sliding zone at great depth. (source: SFSU CREEP Project)

 

Results
On the basis of the GPS velocity field, the Anatolia and Aegean blocks show counterclockwise motion with respect to the Eurasian plate and the rate of the motion increases towards the west. The locking depth variations of the work are between 20-25 km, which correlates with the focal depths of significant earthquakes. The major fault slip rates are consistent with some of the geological slip rate estimates.

 

Results of the model. Figure shows Anatolian plate motion and slip rate estimates of major faults. Credit: Mehmet Köküm

Minds over Methods: Reconstruction of salt tectonic features

Minds over Methods: Reconstruction of salt tectonic features

What is the influence of salt tectonics on the evolution of sedimentary basins and how can we reconstruct such salt features? Michael Warsitzka, PhD student at the Friedrich Schiller University of Jena, explains which complementary methods he uses to better understand salt structures and their relation to sedimentary basins. Enjoy!

 

Credit: Michael Warsitzka

Reconstruction of salt tectonic features from analogue models and geological cross-sections

Michael Warsitzka, PhD student, Institute of Geosciences, Friedrich Schiller University Jena

Salt tectonics, as a sub-discipline of structural geology, describe deformation structures developing due to the special deformation behaviour of salt (as synonym for a sequence of evaporitic rocks). Salt behaves like a viscous fluid over geological time scales and, therefore, it may flow due to lateral differences in thickness and density of the supra-salt layers. This influences the structural evolution of sedimentary basins, because salt flow can modify the amount of regional subsidence of the basin. Local sinks (“minibasins”) develop in regions from where salt is squeezed out and salt structure uplifts, e.g. diapirs or pillows evolve in regions of salt influx. Unfortunately, temporal changes of salt flow patterns are often difficult to reconstruct owing to enigmatic ductile deformation structures in salt layers. Understanding the evolution of salt-related structures requires either forward modelling techniques (e.g. physically scaled sandbox experiments) or restoration of sedimentary and tectonic structures of the supra-salt strata.

In my PhD thesis, I tried to integrate both, analogue modelling and restoration, to investigate salt structures and related minibasins developed in the realm of extensional basins. The sandbox model is a lab-scale, simplified representative of natural salt-bearing grabens, e.g. the Glückstadt Graben located in the North German Basin (Fig. 1). A viscous silicone putty and dry, granular sand were used to simulate ductile salt and brittle overburden sediments. Cross sections were cut through the model at the end of each experiment to conduct reconstruction of the final experimental structures. The material movements were monitored with a particle tracking velocimetry (PIV) technique at the sidewalls of the experimental box.

(text continues below figure)

Fig 1: 2D restoration of the supra-salt (post-Permian) strata in the Glückstadt Graben (Northern Germany). Credit: Michael Warsitzka

Using experimental and geological cross sections, structures in the overburden of the ductile layer can be reconstructed, if present-day layer geometries and lithologies of the overburden strata can be identified. From natural clastic and carbonatic sediments we know that they compact with burial, reducing the layer thickness. Therefore, the reconstruction procedure sequentially removes the uppermost layer and layers beneath are decompacted and shifted upwards to a horizontal surface (Fig. 2). The sequence of decompaction and upward shifting is then repeated until the earliest, post-salt stage is reached (Fig. 1). It intends to restore the initial position, shape and thickness of each reconstructed layer.

In analogue experiments, no decompaction is necessary, because the compressibility of the granular material is insignificant for depths of a few centimetre. Restoration can be directly applied to coloured granular layers revealing detailed layer geometries for each experimental period (Fig. 2a). The PIV technique displays coeval material movement and strain patterns occurring during the subsidence of the experimental minibasins (Fig. 2b). Based on the observation that the experimental structures resemble those reconstructed from the natural example (Glückstadt Graben during the Early Triassic, Fig. 1), it can be inferred that strain patterns observed in the experiments took place in a similar manner during the early stage of extensional basins. This demonstrates the advantage of applying both methods. First, original geometries of basin structures can be determined from the restoration and then reproduced in the model. If the restored geometries are suitably validated by the models, the kinematics observed in the model can be translated back to nature and help to understand the effect of salt flow on the regional subsidence pattern.

Fig 2: Result of an analogue model showing (a) reconstructed sand layers restored from a central cross section, and (b) monitored displacement and strain patterns in the viscous layer above the left basal normal fault. Credit: Michael Warsitzka

Minds over Methods: Sensing Earth’s gravity from space

Minds over Methods: Sensing Earth’s gravity from space

How can we learn more about the Earth’s interior by going into space? This edition of Minds over Methods discusses using satellite data to study the Earth’s lithospere. Anita Thea Saraswati, PhD student at the University of Montpellier, explains how information on the gravity of the Earth is obtained by satellites and how she uses this information to get to know more about the lithosperic structure in subduction zones.

 

Sensing Earth’s gravity from space

Anita Thea Saraswati – PhD student, Géosciences Montpellier

From the basic physics we all know that the value of the gravity is a constant 9.81 meter per second squared. This assumption would be true if the Earth were a smooth nonrotating spherical symmetric body made of uniform element and material. However, because of the Earth’s rotation, internal lateral density variation, and the diversity of the topography (including mountains, valleys, oceans and glaciers), the gravity  varies all over the surface. These tiny changes in gravity due to the mass variations could be a crucial hint for understanding the structure of the Earth, both on the surface and at depth.

The determination of Earth’s gravity field has benefited from various gravity satellite missions that have been launched recently. Among them are the Challenging Minisatellite Payload (CHAMP) (2000-2010), the Gravity Recovery and Climate Experiment (GRACE) (2002-recent), and most recently the Gravity field and steady-state Ocean Circulation Explorer (GOCE) (2009-2013). From these missions, finally a global high quality coverage of Earth’s gravity field became available. (Yay!)

GRACE observation data are very useful for the temporal analysis of changes in gravity. For example to detect the gravity signal before and after a big earthquake, like the Sumatra Mw 9.1 (2004) and Tohoku Mw 9.1 (2011) ones. By analyzing the changes of gravity signal during a certain period of time, it could also be used to detect the drought over a large scale area, which is used in several areas in Africa and Australia.

(text continues below figure)

Design of GOCE satellite observation. A geoid’s shape is showed on the bottom left. On the top right, the GOCE gravity gradients in six components. (Source : ESA)

 

Meanwhile, GOCE is very suitable for the construction of a static model of Earth’s gravity field. Since this satellite has a very low orbit, ~250 km above mean sea level, it has a better spatial resolution. Its accuracy is also better than the previous missions, up to 1 mGal. GOCE is equipped with a gradiometer, which measures the gravity acceleration in three directions (x, y, and z). Afterwards this information is processed into a gravity-gradient dataset containing six components (XX, XY, XZ, YY, YZ, ZZ).

This gravity gradient is the first derivative of the gravity acceleration, which provides us better information about the geometry of the earth’s structure than the gravity acceleration itself. For my PhD, I use this gravity gradient dataset to analyze the lithospheric structure of subduction zones. Before treating the GOCE observation data, I am developing a computational code to calculate the gravity and gravity gradient due to the effect of topography, also called the topographic reduction. The observed gravity and gravity gradient values will be reduced by this topography effect in order to get the anomaly signal. This means that only the signal due to other geodynamic phenomena over the observed area (e.g. slab, isostasy, mantle plum, etc.) is left. By doing further processing, we can obtain the lateral variations of the lithospheric structure in the study areas and then investigate the correlation with the occurrence of mega-earthquakes in these subduction zones.

Since there is still some ambiguity about the information that is produced by gravity data only, it is better to combine the use of them with others geophysical or geological measurements, e.g. seismic tomography measurements and magnetic field observations.

 

Global coverage of GOCE gravity gradient (in milliEötvös) in radial direction (ZZ) (Panet, I. et al., 2014)

 

Reference:

Panet, I., Pajot-Métivier, G., Greff-Lefftz, M., Métivier, L., Diament, M. and Mandea, M., 2014. Mapping the mass distribution of Earth/’s mantle using satellite-derived gravity gradients. Nature Geoscience7(2), pp.131-135.