Page 1
Patrick Smith Academic History
PhD Project Title: Combining magma flow modelling and seismicity on Montserrat
Supervisor: Prof. Jürgen Neuberg (aka Locko)
Funded by NERC Quota Award
B.Sc. Natural Sciences from The University of Durham (2001-04)
M.Res. Physics of the Earth and Atmosphere from The University of Leeds (2005/6)
Page 2
Low frequency seismicity
Weak high frequency onset
Coda:• harmonic, slowly decaying• low frequencies (1-5 Hz)
→ Are a result of interface waves originating at the boundary between solid
rock and fluid magma
What are low-frequency earthquakes?
Specific to volcanic environments
Page 3
Source
Propagation of seismic energyConduit Resonance • Energy travels as interface waves along conduit walls at velocity controlled by magma properties
• Top and bottom of the conduit act as reflectors and secondary sources of seismic waves
• Fundamentally different process from harmonic standing waves in the conduit
Trigger Mechanism = Brittle Failure of Melt
Page 4
Propagation of seismic energy
Page 5
P-wave
S-wave
Propagation of seismic energy
Page 6
Interface waves
P-wave
S-wave
Propagation of seismic energy
Page 7
Interface waves
Propagation of seismic energy
Page 8
Interface waves
Propagation of seismic energy
Page 9
Interface waves
Propagation of seismic energy
Page 10
Interface waves
Propagation of seismic energy
Page 11
Propagation of seismic energy
Page 12
reflections
Propagation of seismic energy
Page 13
reflections
Propagation of seismic energy
Page 14
Propagation of seismic energy
Page 15
Low frequencies
High frequencies
FAST MODE: I1NORMALDISPERSION
SLOW MODE: I2INVERSEDISPERSION
Low frequencies
High frequencies
Acoustic velocity of fluid
Propagation of seismic energy
Page 16
I1
I2
Propagation of seismic energy
Page 17
I1
I2
S
Propagation of seismic energy
Page 18
S
I1
I2
Propagation of seismic energy
Page 19
S
I1
I2
Propagation of seismic energy
Page 20
‘Secondary source’
I2
Propagation of seismic energy
Page 21
Surface-wave
‘Secondary source’
Propagation of seismic energy
Page 22
Surface-wave
Propagation of seismic energy
Page 23
I1R1
Propagation of seismic energy
Page 24
I1R1
Propagation of seismic energy
Page 25
I2
I1R1
Propagation of seismic energy
Page 26
I2
‘Secondary source’
Propagation of seismic energy
Page 27
‘Secondary source’
Propagation of seismic energy
Page 28
Propagation of seismic energy
Page 29
Propagation of seismic energy
Page 30
Propagation of seismic energy
Page 31
Most of energystayswithin the conduit
Propagation of seismic energy
Page 32
Most of energystayswithin the conduit
Propagation of seismic energy
Page 33
Most of energystayswithin the conduit
Propagation of seismic energy
Page 34
Most of energystayswithin the conduit
Propagation of seismic energy
Page 35
Propagation of seismic energy
Page 36
R2
Propagation of seismic energy
Page 37
R2
Events are recorded by
seismometers as surface
waves
Propagation of seismic energy
Page 38
Why are low frequency earthquakes important?
• Have preceded several major eruptions in the past
• Provide direct link between surface observations and internal magma processes
• Correlated with the deformation and tilt - implies a close relationship with pressurisation processes (Green & Neuberg, 2006)
Page 39
Incorporating flow model data
Motivation for PhD Project
Properties of the
magma
seismic signals(surface)
Magma properties(internal)
Seismic parameters
Signal characteristics
Incorporate flow model data into wavefield models
Combining magma flow modelling and seismicity on Montserrat
Page 40
Incorporating flow model dataAims & Methodology
Derive seismic
parameters
Flow model data
2-D Finite-Element
Use in wavefield models
2-D Finite-Difference
Domain Boundary
Solidmedium
Liquid magma
Damped Zone
Free surface
Seismometers
Source Signal:
1Hz Küpperwavelet
100m below top of conduit
ρ = 2600 kgm-3
α = 3000 ms-1
β = 1725 ms-1
Domain BoundaryDomain Boundary
Solidmedium
Liquid magma
Damped Zone
Free surface
Seismometers
Source Signal:
1Hz Küpperwavelet
100m below top of conduit
ρ = 2600 kgm-3
α = 3000 ms-1
β = 1725 ms-1
Solidmedium
Liquid magma
Damped Zone
Free surface
Seismometers
Source Signal:
1Hz Küpperwavelet
100m below top of conduit
ρ = 2600 kgm-3
α = 3000 ms-1
β = 1725 ms-1