Overview of Research Interests
We live on a dynamic planet, whose upper layers respond in a brittle fashion to slow forcing from the underlying mantle. This takes the form of localised earthquake faulting on a large scale, and fracturing on a much smaller scale. Understanding the processes that lead to the mechanical, structural, hydraulic and chemical properties of rocks undergoing low-temperature deformation is a complex and fascinating challenge. The Earth is a complex material to start with, as a consequence of previous geological events. It is also opaque – we cannot seeinto its sub-surface. The acceleration from stable crack growth to catastrophic dynamic failure is extremely rapid, non-linear, and difficult to predict for individual events. There are strong feedbacks between the elements involved. Nevertheless, complex non-linear systems often ‘self-organise’ spontaneously to produce order and pattern in the population dynamics, for example the scale-invariant or fractal geometry of fault and earthquake populations.
A quantitative understanding of such processes helps: to build our cities in a way that mitigates earthquake risk; to assess the predictability of volcanoes and earthquakes; to evaluate risk of induced seismicity associated with subsurface engineering projects; and to predict fluid flow and transport in subsurface reservoirs and aquifers, including applications to geothermal and carbon storage reservoirs to mitigate climate change.
I am interested in attacking these problems with a variety of analytical, numerical, laboratory, and field techniques, usually in collaboration with colleagues from a range of disciplines. The aim is to describe patterns of earthquakes and fractures in space and time, to work out their underlying mechanisms, and where possible to provide results in a form that can be used in a practical way.
1. Earthquake hazard.
Earthquake hazard depends on the probability of occurrence of events, the nature and strength of the seismic source, and the propagation of the disturbance to the site of a particular building or facility. During my PhD I developed a method, based on statistical mechanics and information theory, to quantify earthquake hazard by synthesising short-term earthquake recurrence statistics with longer-term geological and tectonic information on deformation rates. More recent work has concentrated on improving earthquake hazard estimates by quantifying the effect of systematic and random uncertainties on estimates of earthquake occurrence rates, including issues of model selection and convergence of model parameters as more data become available.
2. Earthquake precursors, predictability and forecasting.
The notion of of deterministic earthquake “prediction” – i.e. “specifying in advance the precise location, size and time of occurrence of an individual event, above chance” (where every word counts!) – has long been the ‘holy grail’ of Earth science. However, the extreme non-linearity of the underlying physics, the complexity of the Earth, and the lack of good data at large space and time scales make this problem extremely (perhaps inherently) intractable. Along with Philip Meredith and colleagues in UCL, I am interested in developing analytical and numerical models, based on time-dependent fracture and damage mechanics, to predict and test the precise nature and form of precursors to catastrophic failure observed under controlled conditions in laboratory rock samples, and in examining the extent to which such behaviour scales in space and time. This includes a formal statistical evaluation of the degree of predictability in the population dynamics for natural events, mainly associated with the tendency of earthquakes to cluster in space and time. At present the main applications are in assessing conditional probabilities for operational forecasting of time-dependent seismic and volcanic hazard, including testing forecasting hypotheses in real time.
3. Scaling of faults and fractures. In order to underpin the models described above a direct investigation of the structure and scaling properties of fracture systems, and how they might come about, is necessary. This includes the observed scale-free or ‘fractal’ scaling of such patterns, and how they can be correlated with mechanical (e.g. the fracture toughness) and hydraulic (e.g. fluid permeability of fault gouge) properties of fault and fracture systems. Along with Ferenz Kun, and colleagues in Debrecen University, we have recently provided a physical basis for the observed scaling of rupture parameters in discrete-element models for the catastrophic failure of porous media in discrete element models, involving the competition between structural disorder and the physics of grain-scale processes.
4. Fluid-Rock Interactions. Fluid-rock interactions play a crucial role in many important geological process. The rock physics group at Edinburgh has developed a significant experimental capability to investigate the geochemistry and geophysics of this interaction, demonstrating systematic changes in fluid permeability as a function of stress state and pore-fluid chemistry, even on relatively short time scales (hours to weeks), and developing models for environmentally-assisted crack growth. We have also contributed to several European-funded projects to quantify fluid-rock interactions in boreholes drilled into the active Aegion fault in the Gulf of Corinth, Greece and the experimental evaluation of salt transport in porous fractured media. Currently the main focus is on developing a capability of directly imaging deformation and fluid flow using high-resolution CT imaging of the process, in combination with acoustic and hydraulic monitoring, and on assessing the probability of inducing earthquakes by injecting or extraction fluid into/out of the Earth.
5. Earth Structure. Earthquakes may be used as natural sources to illuminate the structure of the Earth at different scales, using techniques similar to medical imaging by tomography, the main difference being the lack of control in the source location. We have used this method to investigate earth structure in the form of subtle variations in seismic velocity and anelastic attenuation in the Aegean and the Pacific, in good agreement with tectonic models which have been suggested for these areas, and, in collaboration with the BGS Anisotropy project, using the scattered wave field to track changes in the pore pressure field using time-lapse seismic imaging in subsurface reservoirs.