Microphysics of evolving rock viscosity in the seismic and glacial cycles

Lead Research Organisation: University of Cambridge
Department Name: Earth Sciences

Abstract

Despite being the epitome of strength, the solid rocks below Earth's surface can flow surprisingly rapidly over human timescales, impacting processes of societal relevance. This project aims to deliver new equations describing this flow based on the underlying processes operating in the rocks.

Major earthquakes and the melting of ice sheets cause deflections of Earth's surface that are facilitated by viscous flow of the hot rocks below. This deformation creates important feedbacks. During the seismic cycle, earthquakes induce viscous flow of rocks beneath the fault zone that impacts the spatial and temporal distributions of future earthquakes. During the glacial cycle, viscous flow of rocks beneath melting ice sheets causes ground uplift that impacts sea-level change. Therefore, modelling these systems requires knowledge of the viscosity of rocks in Earth's lower crust and upper mantle.

Unfortunately, the viscosities of these rocks are not constant but instead undergo a transient evolution whenever there is a change in the applied forces. Whilst we know that this viscosity evolution occurs, we do not know why. Without knowing the microscale processes that control the viscosity evolution, we cannot formulate equations that reliably predict flow of rocks in the Earth over the seismic and glacial cycles. At a time when populations exposed to seismic risk are rapidly expanding and when the modelling of ice-sheet dynamics is of unprecedented importance, it is critical to develop new models for the viscosity evolution of the rocks that underpin these systems.

Deciphering the microphysical processes that control the viscosity evolution of rocks requires an ambitious multidisciplinary approach. Each element of the research will be centred on the novel adaptation of techniques from the forefront of the materials sciences to analyse key geological minerals. Experiments will be conducted at temperatures up to 1500 degrees Celcius and will induce viscosity evolution by imposing instantaneous changes in the applied forces, analogous to those imposed by earthquakes. For the first time, a subset of the experiments will be performed inside a scanning electron microscope allowing the samples to be directly imaged during the tests. The microstructures of the samples will be analysed using state-of-the-art microscopy techniques, pioneered by our group, to measure distortions of the crystal lattices and the forces trapped within them. The combined mechanical data and microstructural observations will provide the new insights necessary to determine the key processes operating in the crystal lattices of the minerals that cause their viscosities to evolve.

The interpretations from the laboratory will be subject to two critical tests. To check the consistency and robustness of the interpretations, we will employ the latest models of deforming crystalline materials. We will adapt these models, developed to simulate metals, to analyse geological materials. The relevance of the laboratory experiments to natural rocks will be tested by comparing the microstructures of minerals from both settings. We will utilise samples from the deep portions of major fault zones that provide direct records of viscous flow in the lower crust and upper mantle.

The critical information gained from experiments, microstructural analyses, and modelling will be used to construct and calibrate new equations describing viscosity evolution. For the first time, the equations will be based on rigorous analyses of the specific underlying processes. These equations will unlock the next generation of large-scale models that incorporate the impacts of viscosity evolution in the seismic and glacial cycles. The mechanical and microstructural data generated in this project will be made freely available, providing a new and unique digital resource. Similarly, the techniques pioneered in this project will open new frontiers in the dynamics of geological materials.
 
Description Although this research is at an early stage, we have made significant progress in understanding how the strength of hot rocks that flow in a ductile manner deep within Earth is controlled by microscale processes operating within the constituent crystals. The flow of these deep rocks impacts shallower processes, including the behaviour of fault zones that generate major earthquakes and also how the ground rebounds as ice sheets melt. Being able to model the flow of these deep rocks is therefore important in analysing the behaviour of faults over the earthquake cycle and how the rebound of the solid earth affects the rates of ice loss during global warming. However, these modelling efforts have been hindered by our lack of understanding of the microscale processes by which hot rocks can flow, particularly in contexts where the force applied to the rock changes rapidly, such as in the aftermath of an earthquake or during the melting of an ice sheet. Our experiments have provided new insights into how exactly the flow of hot rocks depends on the behaviour of defects in the crystal lattices of their constituent minerals. Our work on the mineral olivine has revealed that elastic forces acting among the defects control the rate at which rocks in Earth's mantle can flow and has revealed that the multiplication of these defects causes the strength of the rock to evolve dramatically during flow. A parallel set of experiments on rock salt have revealed even more complex interactions among the lattice defects demonstrating how the flow behaviour differs between minerals with different crystal structures. These results on rock salt are relevant to applications involving storage of waste or energy in underground salt caverns. For both minerals, we are formulating equations based and calibrated on our results that can be used to predict the behaviour of rocks deep in Earth's interior.
Exploitation Route Our work on the mineral olivine is providing equations that can be applied to predict the behaviour of fault zones over the earthquake cycle and rebound of the solid earth in response to melting ice sheets. These are areas of intense academic research interest and we are already working with collaborators directly to incorporate our equations into models of these large-scale processes. The improved fault-zone and ice-sheet models that we will generate could have impacts the earthquake-hazard and environmental sectors respectively. Our work on the mineral halite is providing equations that can be applied to predict the behaviour of underground salt caverns that are a key target for storage of fuel or waste from the energy sector.
Sectors Energy,Environment

 
Description Conditions for earthquake nucleation in the lower crust
Amount kr 14,415,000 (NOK)
Funding ID 334965 
Organisation Research Council of Norway 
Sector Public
Country Norway
Start 07/2023 
End 06/2027
 
Description Transient deformation in the upper mantle due to present-day deglaciation and earthquakes
Amount € 295,137 (EUR)
Funding ID ENW.GO.001.005 
Organisation Netherlands Space Office 
Sector Public
Country Netherlands
Start 01/2023 
End 12/2026
 
Title Decorated dislocations and (HR-)EBSD data from olivine of the Oman-UAE ophiolite 
Description This dataset is supplemental to the paper Wallis et al. (2021) and contains data on dislocations and their stress fields in olivine from the Oman-UAE ophiolite measured by oxidation decoration, electron backscatter diffraction (EBSD) and high-angular resolution electron backscatter diffraction (HR-EBSD). The datasets include images of decorated dislocations, measurements of lattice orientation and misorientations, densities of geometrically necessary dislocations, and heterogeneity in residual stress. Data are provided as 6 TIF files, 8 CTF files, and 37 tab-delimited TXT files. Files are organised by the figure in which the data are presented in the main paper. Data types or sample numbers are also indicated in the file names. 
Type Of Material Database/Collection of data 
Year Produced 2021 
Provided To Others? Yes  
Impact This is the first dataset on stress heterogeneity within individual grains of the mineral olivine deformed in Earth's mantle. As the analysis of such stress fields is a key research direction for my group and others, this dataset provides an important benchmark against which to assess the quality and significance of future measurements. 
URL https://dataservices.gfz-potsdam.de/panmetaworks/showshort.php?id=75f19561-5dc9-11ec-958d-ed2b0fcbcc...