lecture_2_handout

Comprehensive Lecture Notes on Geophysics and Tectonics

1. Introduction

The study of seismic waves, which are elastic waves propagating through the Earth, is fundamental in the field of geophysics. Seismic waves are essential for understanding various geological phenomena, including earthquakes, volcanic activity, and even anthropogenic events such as explosions or tunneling. The analysis of these waves provides critical insights into the Earth's inner structure and the dynamics of tectonic processes.

Aim:

• This resource aims to explore in detail the mechanisms of seismic wave propagation and their applications in real-world scenarios, enhancing our understanding of geophysical processes.

Importance:

• Understanding seismic waves is crucial for several reasons:

  • Seismic Imaging: Improved methods for imaging subsurface structures help locate natural resources like oil, gas, and minerals.

  • Earthquake Analysis: Insights into the mechanisms of earthquakes aid in developing predictive models, helping to mitigate risks to human life and infrastructure.

  • Civil Engineering: Knowledge of seismic wave behavior is vital in the design of buildings and bridges in earthquake-prone regions, contributing to disaster resilience.

2. Seismic Waves

Seismic waves are categorized based on their propagation characteristics, which influence their behavior as they travel through different geological media.

Definition:

Seismic waves are defined as elastic waves generated by the release of energy during geological processes, propagating through the Earth's crust and interior. They require a material medium with a linear stress-strain relationship to travel effectively.

Elastic Behavior:

• Elasticity refers to the ability of materials, such as rocks, to deform under stress and return to their original shape once the stress is removed. In geophysical contexts, small deformations that occur quickly often exhibit elastic behavior akin to that of springs.

  • Key Principle: The elastic behavior of rocks is described by Hooke’s Law, which expresses the relationship between stress (force applied per unit area) and strain (deformation) mathematically:

    • Hooke’s Law:[ F = kx ]where:

      • F = Applied force (Newtons)

      • k = Spring constant (a measure of stiffness)

      • x = Displacement (the amount of stretch or compression)

Deformation Threshold:

Under continuous stress, rocks will exhibit elastic behavior up to a certain limit known as the elastic limit or yield point. Beyond this threshold, if the applied stress exceeds the cohesive strength of the rock, permanent deformation occurs, leading to structural changes such as fracturing or faulting.

Energy Dissipation:

• As a fault ruptures during an earthquake, it emits seismic energy, which propagates as seismic waves. This energy release happens when elastic deformation is restored as the energy moves outward from the rupture point through the rock.

• Understanding energy dissipation is crucial, as this affects the amplitude and intensity of seismic waves experienced on the surface.

Typical Strain Values:

• Near-source strain: Typically measured as approximately (10^{-3}) to (10^{-4}), this indicates significant deformations close to the earthquake source.

• Far-field waves: In contrast, strain values for waves traveling further away often drop below (10^{-6}), indicating the diminishing effect of the initial rupture with distance from its source, characterized by amplitude decay.

3. Fundamental Principles

3-D Stress and Strain:

In three-dimensional space, stress can be represented mathematically by a 3x3 stress tensor, which accounts for the state of stress at a point in a material related to the forces acting on it.

Stress Tensor:

The stress tensor is expressed as: [ \sigma = \begin{bmatrix} \sigma_{11} & \sigma_{12} & \sigma_{13} \ \sigma_{21} & \sigma_{22} & \sigma_{23} \ \sigma_{31} & \sigma_{32} & \sigma_{33} \end{bmatrix} ]

• where each component represents the stress acting on a specific plane in three-dimensional space.

Stress and Elasticity in a 3D Medium:

Linear Stress-Strain Relationship

The relationship between stress and strain in a 3D context is defined using the elastic tensor, which encompasses various material properties: [ \sigma_{ij} = \sum_{k=1}^{3} \sum_{l=1}^{3} C_{ijkl} e_{ij} ]

• Here, (C_{ijkl}) is the elastic tensor with up to 81 components that define how multiple stress components contribute to overall strain in a material, providing a comprehensive model of rock behavior under stress.

P-wave and S-wave Velocities:

Two primary types of seismic waves are studied:

• P-waves (Primary waves): }}Compressional waves that travel by alternately compressing and expanding the medium, allowing them to propagate through solids, liquids, and gases.

• S-waves (Secondary waves): Shear waves that move the ground perpendicular to their direction of travel, capable of traveling only through solids.

Velocity Formulations:

• P-wave Velocity:[ V_p = \sqrt{\frac{\lambda + 2\mu}{\rho}} ]

• S-wave Velocity:[ V_s = \sqrt{\frac{\mu}{\rho}} ]

• Where:

  • (\lambda) and (\mu) are Lamé parameters, defining rheological properties of the material.

  • (\rho) is the material density.

Implication:

The derived equations indicate that P-waves are faster than S-waves (usually 1.5 to 2 times faster), highlighting the differential nature of energy transfer in geological formations. This knowledge is pivotal in earthquake engineering and wave propagation studies.

4. Anelasticity in Seismic Waves

Observation of Anelasticity:

Anelasticity describes the phenomenon where the elastic modulus of a material changes due to various conditions, such as temperature, pressure, and material composition. This variability can lead to different velocity distributions for seismic waves in heterogeneous materials.

Significance in Wave Propagation:

Anelasticity significantly impacts the energy dissipation as seismic waves propagate through various geological formations. Recognizing how these changes affect wave behavior is crucial in improving seismic models and predicting wave amplitudes and velocities under different conditions. For instance, rocks at high temperatures may exhibit softer elastic responses leading to greater energy loss, whereas colder, more rigid rocks may allow faster seismic wave propagation.

5. Ray Theory

Elastic Wave Equation:

The propagation of elastic waves is quantifiably described by the elastic wave equation:[
abla , \bullet , \sigma = \rho \frac{\partial^2 u}{\partial t^2} ]

• In this equation,

  • (u) refers to the displacement field of the wave,

  • (t) represents time,

  • (\sigma) is the stress tensor, and

  • (\rho) is the density of the medium.

• This equation is fundamental for modeling wave propagation across different rock layers within the Earth.

Challenges in Numerical Analysis:

The computational modeling of seismic wave propagation using the elastic wave equation presents challenges due to the complexity involved in resolving wavefronts across heterogeneous terrain. As grid spacing decreases to capture higher frequency content in the data, the computational resources required increase exponentially, often limiting the ability to conduct real-time analysis or to accurately analyze complex geological settings.

Applications of Ray Theory:

Ray theory provides a simplified approach to analyzing seismic wavefields under high-frequency assumptions, emphasizing the characteristics of wave propagation through various media. Key uses include:

• Predicting travel times through different geological layers.

• Understanding how seismic waves reflect or refract when encountering boundaries between materials.

• Assisting in creating seismic models that aid in assessing subsurface structures and characterizing faults.

6. Eikonal Equation and Transport Equation

Eikonal Equation: Concepts and Formulation:

The Eikonal equation is essential for understanding wave front propagation and how seismic waves travel through materials. This equation can be stated as: [ |
abla T| = U ]

• Where:

  • T = traveltime of the wavefront,

  • U = slowness, which is the inverse of wave velocity.

• This relationship is fundamental for determining the travel paths of seismic waves, especially through varying geological materials.

Transport Equation: Understanding Amplitude:

The transport equation describes changes in the amplitude of seismic waves as they propagate: [ 2
abla A \bullet
abla T + A
abla^2T = 0 ]

• In this equation,

  • A represents the amplitude of the seismic wave,

  • The equation reflects how energy spreads and decreases over distance, crucial for accurately modeling wave behavior in seismic studies.

Raypath Calculation:

Calculating the traveltime along a specified raypath involves integral calculus: [ t = \int_{S_R} U , ds ]

• By evaluating this integral along a defined path, scientists can determine the duration seismic waves take to travel through various media, providing vital data for earthquake monitoring and analysis.

7. Fermat’s Principle and Snell’s Law

Explanation of Fermat’s Principle:

Fermat’s Principle states that the actual path taken by a wave between two points is the one that minimizes the travel time. This principle holds true across multiple fields, including optics, acoustics, and seismic wave analysis, underlining the behavior of all types of wave propagation.

Derivation and Application of Snell’s Law:

• Snell’s Law provides a mathematical representation of the relationship between incident and refracted angles when seismic waves transition from one medium to another: [ n_1 \sin(\theta_1) = n_2 \sin(\theta_2) ]

• Where:

  • n₁ and n₂ refer to the indices of refraction of the respective media, which correlate to the wave velocities in those materials.

• Understanding Snell’s Law assists scientists in predicting how seismic waves behave at geological interfaces, enabling accurate modeling of subsurface structures and aiding in resource exploration.

8. Locating Paths in Heterogeneous Media

Ray Tracing Method: Fundamentals and Applications:

The ray tracing method is a tactical approach for simulating the trajectory of seismic waves as they navigate through various layers of geological media. The algorithm works iteratively, updating the wave’s trajectory until it reaches a designated receiver:

• This method is instrumental for understanding complex subsurface formations and is widely used in the fields of oil and gas exploration, mineral prospecting, and earthquake studies.

Computational Challenges:

Computational difficulties arise in heterogeneous media due to its non-linear paths and a multitude of potential solution scenarios, demanding innovative algorithms for verification. The need for advanced computational techniques is paramount for ensuring accurate seismic modeling and interpretation of geophysical data.

9. Surface Waves vs. Body Waves

Characteristics of Surface Waves:

Surface waves are seismic waves that travel along the Earth's surface, where geological interactions can significantly enhance their amplitude. Two key types of surface waves are:

• Rayleigh waves: Such waves exhibit a rolling motion, moving both vertically and horizontally as they travel, leading to significant ground motion during seismic events.

• Love waves: These waves cause horizontal shearing of the ground without vertical displacement, contributing to severe lateral movements during earthquakes.

Dispersion of Waves:

• Surface and body waves exhibit dispersion, where different frequencies travel at different velocities. This effect is driven by the properties of the medium they traverse. Generally, lower frequency (longer wavelength) waves can penetrate deeper and travel faster than higher frequency waves, affecting the arrival times recorded by seismometers and complicating the interpretation of seismic data.

10. Seismogram Structure

Components of Seismic Sensors:

3-Component Seismic Sensors are critical in capturing ground motion data. These sensors record:

• Vertical motion (up and down)

• Horizontal motion (north-south and east-west)

• The acquisition of data from multiple axes provides comprehensive insights into the nature of seismic events and facilitates detailed interpretations of the data.

Identification of Epicenter:

The epicenter of an earthquake is typically determined through the analysis of P-wave and S-wave arrival times at multiple seismic stations. By employing triangulation methods, researchers can accurately pinpoint the epicenter, essential for assessing earthquake impacts and configuring timely response strategies.

11. References

• Rawlinson, N., Hauser, J., & Sambridge, M. (2008). Seismic ray tracing and wavefront tracking. Advances in Geophysics.

• Stein, S., & Wysession, M. (2003). An Introduction to Seismology, Earthquakes, and Earth Structure. Blackwell Publishing.