Thermal Engine Slides The Tectonic Engine
Part 1: Mantle Thermodynamics
Unit 04 12th Grade Physics
What powers a planet?
To move 100-km thick slabs of solid rock, you need a mechanism of immense power. Like any machine, Earth's tectonic system requires an energy source and a heat sink .
Energy Source
Internal Thermal Energy
Heat Sink
Deep Space (via the Crust)
Mechanical Work
Plate Motion
Thermal Budget: \( Q_{total} \)
1. Primordial Heat
Residual energy from planetary accretion and core formation (the "big splash" of Earth's birth).
2. Radiogenic Heat
Energy released by the decay of unstable isotopes in the mantle and crust: \( ^{238}U \), \( ^{232}Th \), and \( ^{40}K \).
~47 TW
Global Heat Flow
Equivalent to 1.5 trillion light bulbs constantly burning.
Will it Convect?
Fluid motion (convection) occurs when buoyancy forces overcome viscous resistance and thermal diffusion.
\[ Ra = \frac{\alpha g \Delta T d^3}{\kappa \nu} \]
\( \alpha \): Thermal expansion coefficient
\( g \): Gravitational acceleration
\( \Delta T \): Temperature difference
\( d \): Depth of the layer
\( \kappa \): Thermal diffusivity
\( \nu \): Kinematic viscosity
The Cycle of Motion
1 Heating: Lower mantle absorbs heat from the core.
2 Expansion: Rock expands, density decreases (\( \rho = \frac{m}{V} \)).
3 Ascent: Buoyant material rises as a "plume" or upwelling.
4 Cooling: At the surface, heat is lost; material becomes dense and sinks.
COOL SURFACE
HOT CORE
In the mantle, rock flows like an extremely viscous liquid.
Mechanical Link
If the mantle is moving horizontally at the top of a convection cell, how does that translate to the movement of a solid lithospheric plate?
Discuss with a partner: Is the plate being pushed by the mantle, or is the plate part of the cell?
Heat Budget Worksheet Thermal Engine Analysis
Earth Systems Physics | Module 01: Thermodynamics
NAME:
DATE:
01
The Internal Thermal Budget
Total Earth heat flow is estimated at ~47 Terawatts. Analyze the two primary sources of this energy.
A. Primordial Heat Source
Describe the physical processes that generated this heat during Earth's formation.
B. Radiogenic Heat Source
Identify the key isotopes involved and the nuclear process responsible for energy release.
Critical Thinking: Heat Flux Variation
Why is heat flux significantly higher at mid-ocean ridges compared to stable continental cratons? Use the concept of lithospheric thickness and thermal gradients in your answer.
02
Convection Stability
The Rayleigh Number (\( Ra \)) determines whether a fluid layer will convect or remain stagnant.
\[ Ra = \frac{\alpha g \Delta T d^3}{\kappa \nu} \]
Proportionality Analysis:
If the viscosity (\( \nu \)) of the mantle increases significantly (e.g., due to cooling), what happens to the Rayleigh Number and the vigor of convection?
Depth Relationship:
The term \( d^3 \) implies that the depth of the convection layer is critical. How does the Ra for whole-mantle convection differ from upper-mantle-only convection?
03
The Convection Circuit Blueprint
In the space below, construct a cross-sectional diagram of a mantle convection cell. Your model must explicitly label: Thermal boundary layers (top/bottom), adiabatic core of the cell, upwelling plumes, and downwelling slabs.
04
Synthesis: Energy Transformation
Conservation of Energy
Explain how thermal energy (heat) is converted into kinetic energy (plate motion). Where in the "engine" does this energy transformation occur most efficiently?
Prediction
Billions of years from now, Earth's internal heat will diminish. Predict how the Rayleigh Number will change and what the eventual surface expression of Earth will look like (tectonically speaking).
REF_DOC: TE-L1-WS EARTH_CORE_TEMP: ~6000K MANTLE_VISCOSITY: ~10^21 Pa·s
Thermodynamics Teacher Guide Thermodynamics Facilitation Guide
Lesson 01: The Mantle Engine
Level 12th Physics
Instructional Overview
The Hook: Convection Simulation
Before slides, set up a clear tank of water with a heat source (bunsen burner or heating pad) at one end. Add potassium permanganate crystals or rheoscopic fluid. Prompt: "We often hear the mantle is 'solid'. If it's solid, how can it move like this water? What must be different about the physics of rock over millions of years?"
Key Misconception Alert
Students often believe the mantle is liquid magma. Emphasize that it is a solid that flows (ductile/plastic behavior). Use the analogy of glacier ice or high-viscosity glass. Magma only exists in tiny pockets (1-2%) or specific decompression zones.
Learning Targets
Differentiate between primordial and radiogenic heat sources.
Calculate/evaluate Rayleigh Number variables for mantle stability.
Model the energy transfer loop of a convection cell.
Physics Breakdown: Rayleigh Number
The critical Rayleigh Number (\( Ra_c \)) for the onset of convection is typically ~1700. The Earth's mantle has an estimated \( Ra \) of \( 10^6 \) to \( 10^8 \). This means the mantle is highly unstable and convective motion is turbulent and vigorous on geologic scales.
Driving Factors (Numerator)
Thermal Expansion (\( \alpha \)): More expansion = more buoyancy.Temperature Gradient (\( \Delta T \)): Higher heat difference = stronger drive.Depth (\( d^3 \)): The most influential factor. Large volumes of fluid are much more likely to convect.
Resisting Factors (Denominator)
Thermal Diffusivity (\( \kappa \)): If heat leaks away too fast, buoyancy is lost.Viscosity (\( \nu \)): Internal friction. Higher viscosity strongly resists motion.
Worksheet Answer Key
Section 01: Heat Budget
A. Primordial Heat:
Focus on accretionary energy (kinetic energy of impactors converted to heat) and core formation (gravitational potential energy of sinking iron converting to thermal energy).
B. Radiogenic Heat:
Must mention radioactive decay of Uranium-238, Thorium-232, and Potassium-40. Explain that as these atoms decay to stable daughter products, they release kinetic energy (alpha/beta particles) and gamma rays which are absorbed as heat.
Motion Mechanics Slides Gravity at Work
Part 2: Ridge Push & Slab Pull
Kinematics Gravitational Dynamics
Is Convection Enough?
Earlier models suggested plates were simply "riding" on top of mantle convection cells like a conveyor belt. Modern physics tells a different story: The plates are not just passengers; they are part of the engine itself.
"The tectonic plate is the cold, top boundary layer of a convection cell."
1. Ridge Push
At mid-ocean ridges, the lithosphere is hot and elevated. As it cools and moves away, it thickens and subsides.
The Mechanism:
Gravitational sliding off the topographic high of the ridge.
\[ F_{rp} \propto \Delta h \cdot g \cdot \rho \]
2. Slab Pull
The dominant force in plate tectonics. Cold, old oceanic lithosphere is denser than the underlying warm asthenosphere.
Negative Buoyancy
As the slab sinks, it pulls the rest of the plate behind it. It is essentially falling through the mantle in slow motion.
PULL
The Power of the Slab
We can estimate the work done by a sinking slab using Gravitational Potential Energy (\( U_g \)):
\[ \Delta U_g = m g \Delta h \]
But wait—the slab is in a fluid (the mantle). We must use buoyant mass (\( m_{eff} \)):
\[ m_{eff} = V(\rho_{slab} - \rho_{mantle}) \]
Key Data Point:
Slab Density (\( \rho_s \)) \(\approx\) 3300 kg/m³
Mantle Density (\( \rho_m \)) \(\approx\) 3250 kg/m³
The Result:
That small \( 50 \, \text{kg/m}^3 \) difference across a slab 100km thick and 1000km long generates enough force to move entire continents.
Who is the Driver?
Hypothesis A
Ridge Push
Active mid-ocean ridges "push" the plates apart. Plate velocity should be proportional to ridge length.
Hypothesis B
Slab Pull
Slabs pull the plates into trenches. Plate velocity should be proportional to the length of subduction zones.
Observation: Plates attached to large subduction zones move 5-10x faster.
Gravity Drivers Worksheet Gravity Drivers Analysis
Kinematics of Plate Motion | Module 02: Gravitational Forces
NAME:
DATE:
01
Ridge Push (Gravitational Sliding)
The Mid-Atlantic Ridge (MAR) is approximately 2,500m higher than the surrounding abyssal plains. This topographic gradient creates a "sliding" force.
Vector Diagram
LITHOSPHERE AT RIDGE
Decompose the gravitational force into parallel and perpendicular components relative to the lithosphere-asthenosphere boundary.
The Cooling Effect
As lithosphere moves away from the ridge, it cools and becomes thicker/denser. Explain how this isostatic change reinforces the ridge push force.
02
Slab Pull: Energy Transformation
The Problem:
Old oceanic lithosphere (\( \rho \approx 3300 \, \text{kg/m}^3 \)) is subducting into the asthenosphere (\( \rho \approx 3250 \, \text{kg/m}^3 \)). Calculate the excess mass (\( m_{excess} \)) of a slab that is 100km thick, 1000km wide, and extending 500km into the mantle.
STEP 1: CALCULATE VOLUME (m³)
STEP 2: CALCULATE EXCESS MASS (kg)
Work and Force:
Using the excess mass above, explain why Slab Pull is considered the "primary" driver of plate motion compared to Ridge Push or Mantle Drag.
03
Correlating Force with Velocity
Plate Name Velocity (cm/yr) % Perimeter Subducting Pacific 8.0 35% Nazca 7.5 30% Eurasian 1.5 5% African 2.0 2%
Data adapted from Forsyth and Uyeda (1975).
Analysis Prompt:
Analyze the correlation between velocity and subduction percentage. What does this empirical data suggest about the relative strength of Slab Pull vs. Ridge Push?
04
The Metamorphic Engine
As the slab subducts, minerals in the basaltic crust undergo phase transitions to denser minerals like Eclogite .
Physics Synthesis
Explain how the Basalt-to-Eclogite transition creates a "runaway" effect in slab pull. How does this change the net force (\( \sum F \)) acting on the plate over time?
REF: L2-PLATE-DYNAMICS G = 9.81 m/s² 1 cm/yr = 3.17e-10 m/s
Motion Teacher Guide Motion Mechanics Guide
Lesson 02: Gravitational Drivers
Physics Focus Mechanics
Lesson Facilitation
The Core Debate
Students often enter this lesson assuming convection "pushes" the plates. This lesson shifts the agency to the plates themselves. Key Question: "Does the engine move the car, or is the plate the engine?"
Vector Analysis Tip
For Section 01 of the worksheet, ensure students decompose gravity (\( F_g \)) correctly. The component parallel to the slope (\( F_{parallel} = F_g \sin \theta \)) is what drives "sliding." Even a very small angle (\( \approx 1-2^\circ \)) results in massive forces due to the mass of the lithosphere.
Physics Concepts
Negative Buoyancy
Topographic Gravitational Potential
Phase Transitions (Eclogitization)
Advanced Context: Slab pull is estimated to provide ~90% of the total driving force for plate motion.
Answer Key & Calculations
Worksheet Section 02: Slab Pull Excess Mass
VOLUME CALCULATION
\( V = 100\,\text{km} \times 1000\,\text{km} \times 500\,\text{km} \)
\( V = 5 \times 10^7 \, \text{km}^3 = 5 \times 10^{16} \, \text{m}^3 \)
EXCESS MASS CALCULATION
\( \Delta \rho = 3300 - 3250 = 50 \, \text{kg/m}^3 \)
\( m_{excess} = 5 \times 10^{16} \, \text{m}^3 \times 50 \, \text{kg/m}^3 \)
\( m_{excess} = 2.5 \times 10^{18} \, \text{kg} \)
Worksheet Section 03: Correlation Analysis
Students should observe a strong positive correlation between subduction length and plate velocity. The Pacific and Nazca plates (highly subducting) move significantly faster than the Eurasian and African plates (largely continent-bearing with few subduction zones). This supports Hypothesis B: Slab Pull Dominance.
Worksheet Section 04: Phase Transitions
The Basalt-to-Eclogite transition occurs at depths of ~40-60km. Eclogite is significantly denser than basalt/gabbro. This "densification" increases the negative buoyancy of the slab as it descends, creating a positive feedback loop that accelerates slab pull.
Required Classroom Materials
String and weights (to demo slab pull)
Inclined plane (to demo ridge push)
High-friction surface (to demo mantle drag)
Scientific calculators
Flowing Rock Lab Guide Lab: Flowing Rock Dynamics
Geomaterials Rheology & Viscoelastic Modeling
PHYSICS LAB 03
Abstract
While the mantle is often described as "solid" because it transmits S-waves (shear waves), it behaves as a fluid on geologic timescales (\( 10^6 \) to \( 10^8 \) years). In this lab, we use Silly Putty as an analogue for mantle material to explore the transition between elastic and plastic behavior.
Essential Question:
How can a substance act as both a solid and a fluid simultaneously depending on the timescale of the applied force?
Required Materials
Silly Putty (approx. 50g)
Stopwatch / Timer
Metric Ruler (mm precision)
500g and 1kg weights
Thermometer
Physics Framework: The Maxwell Model
The mantle is a viscoelastic material. Its behavior can be modeled as a spring (elasticity) and a dashpot (viscosity) in series.
SHORT TIMESCALE
Elastic (Spring-like)
LONG TIMESCALE
Viscous (Fluid-like)
A Experiment 1: Stress-Strain Relationship
Roll the putty into a uniform cylinder approximately 5cm long.
High-Stress Event: Quickly pull the ends of the putty apart. Record the nature of the break (brittle or ductile).Low-Stress Event: Slowly and steadily pull the ends of the putty apart over 60 seconds. Observe the deformation.Relaxation Time (\( \tau \)): Hold a stretched piece of putty fixed. Observe how long it takes for the internal stress to dissipate.
B Experiment 2: Gravitational Creep
This experiment models the "flattening" of mantle upwellings or the sinking of cold slabs under their own weight.
Roll the putty into a perfect sphere. Measure its initial height (\( H_0 \)).
Place the sphere on a flat surface. Every 2 minutes for 10 minutes, measure the height of the sphere as it flattens under gravity.
Repeat this experiment with a 500g mass placed on top of the sphere to simulate increased lithospheric pressure.
Temperature Variable: Chill one sample of putty in ice and warm another in your hands. Compare their "sag rates" over 5 minutes.
Observation Data
Experimental Qualitative Notes
Temperature Impact Hypothesis
Rheology Lab Report Rheology Lab Report
Quantifying Viscoelasticity and Creep
STUDENT:
STATION ID: ________
01
Strain Rate Visualization
Plot the Percent Strain (\( \epsilon \)) over time for your control sample and the loaded sample (\( +500\text{g} \)). Strain is calculated as: \( \epsilon = \frac{H_0 - H_t}{H_0} \times 100 \).
Graph: Strain (%) vs. Time (min)
Time (minutes)
Percent Strain (%)
02
Viscosity Estimate
Viscosity (\( \eta \)) is the ratio of stress (\( \sigma \)) to strain rate (\( \dot{\epsilon} \)). Based on the slope of your graph, which sample showed the highest viscosity? Explain using the relationship \( \sigma = \eta \dot{\epsilon} \).
03
Deborah Number
The Deborah Number (\( De \)) is defined as \( De = \frac{t_c}{t_p} \), where \( t_c \) is the relaxation time and \( t_p \) is the observation time. What was the \( De \) for your "quick pull" experiment vs. the "slow sag"?
04
Upscaling the Model
Critical Synthesis Prompt
Earthquakes release energy in seconds (short timescale), while mantle convection moves rock over millions of years (long timescale). Based on your lab results, explain why the mantle acts as a solid to earthquake waves (S-waves) but as a fluid to tectonic plates.
05
Uncertainty and Calibration
Identify two potential sources of error in your measurements (e.g., surface tension, temperature fluctuations, measurement parallax) and how they may have affected your calculated strain rates.
Error Source 1
Error Source 2
END OF EXPERIMENTAL RECORD // RHEO-REPORT-L3
Plume Physics Slides Deep Earth Plumes
Part 4: Hotspots & Plate Vectors
Geospatial Analysis Relative Motion
The Anomaly
If volcanoes are caused by plate boundaries (subduction or rifting), why are there chains of volcanoes in the middle of tectonic plates?
Hawaii
Pacific Plate (Intraplate)
Yellowstone
North American Plate (Intraplate)
Mantle Plumes
Narrow columns of hot mantle rock that rise from the Core-Mantle Boundary (CMB).
Fixed Position: Unlike moving plates, plumes are relatively stationary in the mantle.Decompression Melting: As the plume reaches the lithosphere, it melts to create a "hotspot."Tracking: As the plate moves over the plume, it leaves a trail of extinct volcanoes.
The Kinematic Tape Measure
A hotspot track is a vector record of plate motion over time.
Magnitude (\( v \)) \( \frac{\text{Distance between volcanoes}}{\text{Age difference}} \)
Direction (\( \theta \)) Orientation of the track (Azimuth)
10.5 cm/yr
Average Pacific Speed
Fun Fact: The Pacific Plate moves at roughly the same speed as your fingernails grow.
Change of Plans
The Hawaiian-Emperor chain features a dramatic "bend." This indicates a major shift in the direction of the Pacific Plate approximately 47 million years ago .
Why the bend?
Hypothesis: The subduction of another plate (the Izanagi Plate) changed the torque and forces acting on the Pacific Plate, forcing a redirection.
Pacific Plate Track Visualization
Your Turn
Using real bathymetric data from the Hawaii-Midway-Emperor chain, you will calculate the velocity vectors for the last 80 million years.
Required: Ruler, Protractor, and Geospatial Data Sheets.
Hotspot Vector Worksheet Hotspot Vector Analysis
Kinematics of the Pacific Plate
NAME:
DATE:
Experimental Context
The Hawaiian-Emperor hotspot track provides a high-resolution record of the Pacific Plate's velocity. By measuring the distance between volcanic centers and their radiometric ages, we can determine the speed (\( v \)) and direction (\( \theta \)) of tectonic motion over millions of years.
Table 4.1: Hawaiian-Emperor Chain Data
Location Age (Ma) Dist. from Kilauea (km) Interval Vel. (cm/yr) Kilauea (Hawaii) 0 0 --- Molokai 1.8 150 Kauai 5.1 520 Midway 27.7 2,430 Daikakuji (The Bend) 47.0 3,500 Suiko 64.7 4,800 Meiji 80.1 6,000
Data Source: Clague and Dalrymple (1987). Age in Ma (Millions of years).
01
Magnitude Calculation
Calculate the velocity (cm/yr) for each interval. Remember: \( 1 \, \text{km/Ma} = 0.1 \, \text{cm/yr} \). Show your work for the Midway to Daikakuji interval below:
02
Directional Analysis
Observe a map of the Hawaiian-Emperor chain. The modern chain (Hawaii to Midway) trends roughly N 70° W . The older Emperor chain (Daikakuji to Meiji) trends roughly N 10° W .
Vector Diagram: 0-47 Ma
Vector Diagram: 47-80 Ma
03
The Physics of "The Bend"
1. Velocity Comparison
Did the Pacific Plate speed up or slow down after the 47 Ma event? Support your claim with evidence from your Table 4.1 calculations.
2. Dynamic Causes
Assuming mantle plumes are fixed, a change in plate direction requires a change in the net force (\( \sum F \)) acting on the plate. Based on your knowledge of Ridge Push and Slab Pull, what geological events could cause a plate to change direction so suddenly?
REF_DOC: PLUME-VEC-L4 PLATE_ID: PACIFIC_701 1 km/Ma = 0.1 cm/yr
Plume Teacher Guide Plume Physics Guide
Lesson 04: Hotspots & Kinematic Records
Physics Focus Kinematics
Facilitation Notes
Concept: Stationary Reference Frames
Students often struggle with relative motion. Emphasize that the Mantle Plume is our "stationary" marker (fixed reference frame), while the Plate is the moving object. Analogy: A moving sheet of paper (plate) over a fixed spray-paint can (plume).
Vector Addition vs. Trail Analysis
Note that a trail pointing Northwest means the plate is moving Northwest . Some students mistakenly think the trail points in the direction the plate "came from."
Advanced Discussion
"Are plumes really fixed?"
Recent research suggests "plume wander"—that plumes can actually drift slightly (1-2 cm/yr). However, for 12th-grade physics models, we assume they are stationary to simplify the kinematic calculations.
Worksheet Answer Key
Interval Age Diff (Ma) Dist Diff (km) Velocity (cm/yr) Hawaii to Kauai 5.1 520 10.2 Kauai to Midway 22.6 1,910 8.45 Midway to Daikakuji 19.3 1,070 5.54 Daikakuji to Suiko 17.7 1,300 7.34
Analysis 3.1: Speed Shift
The Pacific Plate significantly accelerated after the bend. Pre-bend (Emperor) speeds average ~7 cm/yr; post-bend (Hawaii) speeds average ~9-10 cm/yr.
Analysis 3.2: Direction Shift
Likely cause: Collision or subduction of the Izanagi-Kula ridge system, or the initiation of new subduction zones in the Western Pacific (Izu-Bonin-Mariana), which changed the Slab Pull vector's direction.
End of Lesson 04 Facilitation Materials | Tectonic Engine Mechanics
Tectonic Project Brief The Tectonic Engine Blueprint
Culminating Synthesis Project | Earth Systems Engineering
Scenario
You have been hired as a "Planetary Architect" for a new interstellar mission to terraform a silicate-based planet. Your task is to design and present a functional Tectonic Engine that ensures long-term carbon cycling and heat dissipation.
You must synthesize everything you've learned about thermodynamics, gravitational forces, and mantle rheology to create a "Technical Blueprint" and a "Physical Defense" of your model.
Deliverables
1. Technical Blueprint (A3 size)
2. Force Vector Supplement
3. 5-Minute Technical Defense
01
The Technical Blueprint
Your blueprint must be a cross-sectional diagram (Core to Crust) that identifies and illustrates:
Thermodynamic Systems
Primordial & Radiogenic heat zones.
Convection cell geometry and boundary layers.
Adiabatic upwellings (plumes).
Mechanical Drivers
Subduction zones (Slab Pull vectors).
Divergent ridges (Ridge Push vectors).
Lithospheric thickness variations.
02
The Force Supplement
Accompanying your blueprint must be a 1-page document containing three mathematical "Proof of Concepts":
PROOF A
A Rayleigh Number calculation (estimated) that proves your mantle will actually convect (\( Ra > 2000 \)).
PROOF B
A vector decomposition showing how ridge height translates to lateral push force.
PROOF C
A density calculation explaining how the Basalt-to-Eclogite transition powers your subduction zones.
03
The Technical Defense
During your 5-minute presentation to the "Architecture Board" (your peers), you must answer one of the following "Crisis Scenarios" for your planet:
Scenario Alpha
Your planet's internal heat drops by 30% suddenly. How does your mantle rheology change, and how will you restart the engine?
Scenario Beta
The planet's ocean evaporates. With no hydration to soften the mantle minerals, how will your subduction zones remain functional?
Preparation Checklist
Review Lesson 1 (Rayleigh Number)
Review Lesson 2 (Slab Density)
Review Lesson 3 (Viscosity)
Review Lesson 4 (Hotspots)
Gather Blueprint Materials
Draft Force Supplement
Tectonic Model Rubric Blueprint Evaluation Rubric
Tectonic Engine Synthesis Project
Total Possible: 40 pts
Criterion Advanced (10 pts) Proficient (8 pts) Developing (6 pts) Thermodynamic Accuracy Blueprint perfectly illustrates heat flow from core to crust, including radiogenic sources and adiabatic cooling. Illustrates internal heat budget and convection cells with minor labeling errors. Convection cells are shown but lack heat source/sink differentiation. Mechanical Drivers Clear, accurate vector representation of Slab Pull and Ridge Push. Accurately shows density variations. Shows gravitational forces; vector magnitudes are generally proportional to model scale. Forces are listed but not visually integrated into the blueprint or lack directionality. Mathematical Proofs All three proofs (Ra, Vector, Density) are mathematically sound and use appropriate scientific notation. Two or more proofs are complete and correct. Minor calculation errors present. Calculations are incomplete or do not correctly apply the relevant physical formulas. Technical Defense Crisis response uses rheological concepts (viscosity, strain rate) to defend a logical tectonic outcome. Response is scientifically plausible but relies on qualitative descriptions rather than physics. Response lacks scientific justification or fails to address the specific crisis scenario.
Feedback & Observations
Total Score: / 40
Evaluator Signature: ____________________________________
