-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathLecture9.html
More file actions
472 lines (386 loc) · 19.5 KB
/
Copy pathLecture9.html
File metadata and controls
472 lines (386 loc) · 19.5 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
<!DOCTYPE html>
<html lang="en">
<head>
<meta charset="UTF-8">
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<link rel="icon" href="favicon.png" type="image/x-icon">
<link href="https://fonts.googleapis.com/css2?family=Noto+Sans+JP&display=swap" rel="stylesheet">
<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/highlight.js/11.8.0/styles/default.min.css">
<link rel="stylesheet" href="style.css">
<script src="https://cdnjs.cloudflare.com/ajax/libs/highlight.js/11.8.0/highlight.min.js"></script>
<script>hljs.highlightAll();</script>
<script>
MathJax = {
tex: {
inlineMath: [['$', '$'], ['\\(', '\\)']],
displayMath: [['$$', '$$'], ['\\[', '\\]']]
},
options: {
renderActions: {
addMenu: []
}
},
svg: {
fontCache: 'global',
scale: 1.0
}
};
document.addEventListener("DOMContentLoaded", () => {
const sidebar = document.querySelector('.sidebar');
const main = document.querySelector('.main-content');
const wrapper = document.createElement('div');
wrapper.className = 'layout';
sidebar.parentNode.insertBefore(wrapper, sidebar);
wrapper.appendChild(sidebar);
wrapper.appendChild(main);
});
</script>
<script src="https://cdn.jsdelivr.net/npm/mathjax@3/es5/tex-svg.js" async></script>
</head>
<body>
<div class="sidebar">
<img src="Logo_NS-blackboard-Expand.png" alt="Sidebar Image">
<nav>
<ul class="nav-menu">
<li><a href="index.html">Welcome to NFYK22003U</a></li>
<ul class="sub-menu">
<li><a href="CourseInformation.html">Course Information</a></li>
</ul>
<li><a href="Welcome.html">Scope and Sequence</a></li>
<ul class="sub-menu">
<li><a href="CoursePrerequisitesExam.html">Prerequisites Exam</a></li>
<li><a href="IntroductiontoPythonProgramming.html">Introduction to Python Programming</a></li>
<li><a href="MathsSpeedrun.html">Maths Speedrun</a></li>
</ul>
<li>
<a href="Lecture1.html">Lecture 1 – Introduction and Transport Processes</a></li>
<ul class="sub-menu">
<li><a href="Project1.html">Project 1 – Oxygen Budget</a></li>
</ul>
<li><a href="Lecture2.html">Lecture 2 – The Dynamics of Rotating Planets</a></li>
<ul class="sub-menu">
<li><a href="Project2.html">Project 2 – The Sverdrup Relation</a></li>
</ul>
<li><a href="Lecture3.html">Lecture 3 – Origin of Carbon and Water on Earth</a></li>
<li><a href="Lecture4.html">Lecture 4 – The Gulfstream</a></li>
<ul class="sub-menu">
<li><a href="Project3.html">Project 3 – Western Boundary Currents</a></li>
</ul>
<li><a href="Lecture5.html">Lecture 5 – Carbon and Plankton</a></li>
<li><a href="Lecture6.html">Lecture 6 – Quasi-Geostrophy</a></li>
<li><a href="Lecture7.html">Lecture 7 – Theory of the Ventilated Thermocline</a></li>
<ul class="sub-menu">
<li><a href="Project4.html">Project 4 – Shadow Zone</a></li>
</ul>
<li><a href="Lecture8.html">Lecture 8 – Waves</a></li>
<li><a href="Lecture9.html">Lecture 9 – Atmospheric Thermodynamics</a></li>
<ul class="sub-menu">
<li><a href="Project5.html">Project 5 – Hadley Cell</a></li>
</ul>
<li><a href="Lecture10.html">Lecture 10 – Surface and Abyssal Circulation</a></li>
<ul class="sub-menu">
<li><a href="Project6.html">Project 6 – Overturning</a></li>
</ul>
<li><a href="Lecture11.html">Lecture 11 – Ocean Mixing and Overturning Circulation</a></li>
<li><a href="Lecture12.html">Lecture 12 – Stability of the Atlantic Meridional Overturning Circulation</a></li>
<ul class="sub-menu">
<li><a href="Project7.html">Project 7 – Oscillations</a></li>
<li><a href="TestExam.html">Test Exam</a></li>
</ul>
<li><a href="Lecture13.html">Lecture 13 – Applications of Machine Learning in CFD and GFD</a></li>
<ul class="sub-menu">
<li><a href="SailbyNightPhysics.html">Sail by Night Physics</a></li>
</ul>
<li><a href="Appendix.html">Appendix</a></li>
<li><a href="References.html">References</a></li>
</li>
</ul>
</nav>
</div>
<div class="main-content">
<h1>Lecture 9 – Atmospheric Thermodynamics</h1>
<h2>Climate Reanalysis</h2>
Climate reanalyses combine past observations with models to generate consistent time series of multiple climate variables.
A climate reanalysis gives a numerical description of the recent climate. It contains estimates of atmospheric parameters such as air temperature, pressure and wind at different altitudes, and surface parameters such as rainfall, soil moisture content, ocean-wave height and sea-surface temperature.
The estimates are produced for all locations on earth, and they span a long time period that can extend back several decades or more.
Newly available remotely-sensed data, the development of Earth System models, and the benefits that would derive from more consistent states across analysis components, all point to a considerable interest in integrated Earth System reanalysis systems, including some level of coupling among the various analysis components (e.g., atmosphere and ocean).
<ul>
<li>Climate reanalyses combine past observations with models to generate consistent time series of multiple
climate variables.</li>
<li>Reanalyses are among the most-used datasets in the geophysical sciences. They provide a comprehensive
description of the observed climate as it has evolved during recent decades, on 3D grids at sub-daily intervals
(ECMWF, copernicus).</li>
<li>Reanalysis products are used extensively in climate research and services, including for monitoring and
comparing current climate conditions with those of the past, identifying the causes of climate variations and
change, and preparing climate predictions. Information derived from reanalyses is also being used increasingly
in commercial and business applications in sectors such as energy, agriculture, water resources, and insurance.</li>
</ul>
<h2>Inter-Tropical Convergence Zone (ITCZ)</h2>
The ITCZ appears as a band of clouds consisting of showers and occasional thunderstorms that encircles the globe near the equator.
Scientists call this area intertropical convergence zone to reflect the influences of wind convergence on conditions near the equator.
Strong heating in the ITCZ causes surface air to expand and rise. The humid, rising, expanding air loses moisture as rainm some of which contributes to the success of tropical rain forests.
The position of the ITCZ varies seasonally because it follows the Sun; it moves north in the Northern Hemisphere summer and south in the Northern Hemisphere winter. As a result, the ITCZ is responsible for the wet and dry seasons in the tropics.
<h2>Hadley Cells</h2>
Because of the differential heating over the globe, Hadley already knew in 1735 that some kind of energy transport must occur from the equator to the poles. He described a circulation cell that would extend all the way to the poles as shown below:
<img src="Figures/Lecture9_HadleyCells.png" alt="HadleyCells" style="margin-top: 1.5em; max-width: 100%; display: block;">
<div class="figure-caption">(Hadley cells could extend all the way to the poles. Credit: Lyndon State College Atmospheric Sciences.)</div>
With very simple arguments of momentum conservation (Held and Hou, 1980) we are going to show that Hadley's vision is not physically reasonable.
The angular momentum of the air parcel is given by:
\[
M = \Omega r^2 + u r
\]
Since \( r = a \cos \varphi \), the moment arm is:
\[
r_{\bot} = a \cos(\varphi)
\]
Substituting into the equation for \( M \):
\[
M = \Omega a^2 \cos^2 \varphi + u a \cos \varphi
\]
If the parcel is not moving at the equator, then \( u = 0 \), and the angular momentum at the equator is:
\[
M_{\text{Eq}} = \Omega a^2
\]
Conservation of angular momentum implies:
\[
M = M_{\text{Eq}} = \Omega a^2 = \Omega a^2 \cos^2 \varphi + u a \cos \varphi
\]
Solving for \( u(\varphi) \), using the identity \( \sin^2 \varphi = 1 - \cos^2 \varphi \), we get:
\[
u(\varphi) = \Omega a \frac{\sin^2 \varphi }{\cos \varphi }
\]
You can see how this is related to the classical definition of angular momentum
$$ L = I \omega$$
The moment of inertia per unit mass is in units of meters$^2$
The classical angular momentum for a rotating object is:
$$
\underbrace{L}_{\text{angular momentum}} = \underbrace{I}_{\text{moment of inertia}} \cdot \underbrace{\omega}_{\text{angular velocity}}
$$
For an air parcel of unit mass at latitude $\varphi$, its absolute angular momentum about the Earth's rotation axis is:
$$
M =
\underbrace{\Omega}_{\text{rotation}} \cdot
\underbrace{r^2}_{\text{moment of inertia}} +
\underbrace{u}_{\text{zonal velocity}} \cdot
\underbrace{r}_{\substack{\text{moment arm} \\ r = \underbrace{a}_{\text{Earth's radius}} \cos \varphi}}
$$
So:
$$
M = \underbrace{\Omega}_{\text{planetary rotation}} \cdot \underbrace{a^2 \cos^2 \varphi}_{r^2} + \underbrace{u}_{\text{relative velocity}} \cdot \underbrace{a \cos \varphi}_{r}
$$
If we treat the air parcel as a point mass $m = 1$, then:
$$
I = m r^2 = \underbrace{a^2 \cos^2 \varphi}_{\text{moment of inertia per unit mass}}
$$
The total angular velocity is the sum of Earth's rotation and the parcel’s relative motion:
$$
\omega = \underbrace{\Omega}_{\text{planetary rotation}} + \underbrace{\frac{u}{r}}_{\text{relative contribution}} = \Omega + \frac{u}{a \cos \varphi}
$$
So the total angular momentum becomes:
$$
L = I \omega = \underbrace{a^2 \cos^2 \varphi}_{I} \left( \underbrace{\Omega}_{\text{rotation}} + \underbrace{\frac{u}{a \cos \varphi}}_{\text{relative motion}} \right)
$$
$$
L = \Omega a^2 \cos^2 \varphi + u a \cos \varphi = M
$$
This shows that the atmospheric expression for angular momentum,
$$
M = \Omega a^2 \cos^2 \varphi + u a \cos \varphi
$$
is equivalent to the classical expression $L = I \omega$, when the moment of inertia and angular velocity are appropriately defined for a rotating air parcel
<h2>The Radiative Equilibrium Distribution</h2>
<ul>
<li>The radiative equilibrium temperature is a hypothetical temperature distribution
if there were no atmospheric and oceanic
motion. Only based upon the radiative
properties of the atmosphere (e.g. clouds
and moisture, etc).</li>
<li>The radiative equilibrium temperature has a
much stronger pole-to-equator gradient
than does the effective emitting
temperature, indicating that there is a
poleward transport of heat.</li>
<li>Radiative forcing thus seeks to maintain a
pole-to-equator temperature gradient, and
the ensuing circulation seeks to reduce this
gradient.</li>
</ul>
A black body subject to a net incoming radiation of $S$ (watts per square metre) has a radiative-equilibrium temperature $T_{\text{rad}}$ given by
$$
\sigma T_{\text{rad}}^4 = S
$$
this being Stefan’s law with Stefan–Boltzmann constant
$$
\sigma = 5.67 \times 10^{-8} \, \text{W} \, \text{m}^{-2} \, \text{K}^{-4}
$$
Thus, for the Earth, we have, at each latitude,
$$
\sigma T_{\text{rad}}^4 = S(\vartheta)(1 - \alpha)
$$
where $\alpha$ is the albedo of the Earth and $S(\vartheta)$ is the incoming solar radiation at the top of the atmosphere, and its solution is shown in the figure below.
<img src="Figures/Lecture9_RadiativeEquilibriumDistribution.png" alt="RadiativeEquilibriumDistribution" style="margin-top: 1.5em; max-width: 100%; display: block;">
<div class="figure-caption">((a) The (approximate) observed net average incoming solar radiation and outgoing infrared radiation at the top of the atmosphere, as a function of latitude (plotted on a sine scale).
(b) The temperatures associated with these fluxes, calculated using $T = \left( \frac{R}{\sigma} \right)^{1/4}$,
where $R$ is the solar flux for the radiative equilibrium temperature and $R$ is the infrared flux for the effective emitting temperature.
Thus, the solid line is an approximate radiative equilibrium temperature. Credit: G. Vallis, 2017.)</div>
<h2>Potential Temperature</h2>
Potential temperature is defined as the temperature that a parcel of air would have if it were expanded or compressed adiabatically (without exchange of heat or matter with the surroundings) from its current pressure to a standard reference pressure, typically 1000 hPa (or 1 bar).
Thus, if air flows across a mountain range it will first cool down and then heat up, but the potential temperature will stay constant. Thus potential temperature is the dynamically relevant quantity. The equation to calculate potential temperature (Θ) from the actual temperature (T) is:
$$
\Theta = T \cdot \left( \frac{P_0}{P} \right)^{\frac{R}{c_p}}
$$
<ul>
<li>$\Theta$ is the potential temperature</li>
<li>$T$ is the actual temperature</li>
<li>$P$ is the atmospheric pressure</li>
<li>$P_0$ is the reference (surface) pressure (often 100,000 Pa)</li>
<li>$R$ is the specific gas constant for dry air</li>
<li>$c_p$ is the specific heat capacity of air at constant pressure</li>
</ul>
<h2>Adiabatic Flow</h2>
$$
\frac{D}{Dt} \log \left[ \frac{T}{T_0} \left( \frac{p_0}{p} \right)^{R/c_p} \right] = \frac{J}{c_p T}
$$
Adiabatic flow → $J = 0$
$$
\frac{D}{Dt} \left[ T \left( \frac{p_0}{p} \right)^{R/c_p} \right] = 0
$$
Potential temperature $
\theta = T \left( \frac{p_0}{p} \right)^{R/c_p}
$
is the temperature an air parcel would have if its pressure were adiabatically adjusted to the reference pressure $p_0$
$$
\Rightarrow \frac{D\theta}{Dt} = 0 \quad \text{(potential temperature is conserved under adiabatic flow)}
$$
Typically $p_0$ represents sea-level pressure ($10^5$ Pa = 1000 hPa). In this case, potential temperature is defined as the temperature an air parcel (with temperature $T$ and pressure $p$) would have if it was adiabatically brought to sea-level pressure.<br>
Potential temperature is closely associated with <i>entropy</i> (constant potential temperature is the same as constant entropy).
<ul>
<li>Adiabatic flow is the one that does not involve heat transfer, i.e., there is no heat added or taken away from the flow</li>
<li>Entropy should increase along an adiabatic flow
$$
ds \geq \frac{\delta q}{T} = 0
\Rightarrow ds \geq 0
$$</li>
</ul>
<h2>Ideal Gas Law</h2>
The ideal gas law is an equation which relates pressure, density and temperature and applies for many fluids where intermolecular attractions are negligible.
Atmosphere is composed of air, which is a mixture of gases (treated as an ideal gas). Below an altitude of 100 km, the atmosphere behaves as a fluid (under the <i>continuum hypothesis</i>).
$$
\Rightarrow \quad p = \rho R_d T
$$
<ul>
<li>$p$: Pressure (Pa)</li>
<li>$\rho$: Density (kg/m³)</li>
<li>$R_d$: Ideal gas constant for dry air</li>
<li>$T$: Temperature (Kelvin)</li>
</ul>
For dry air:
$$
R_d = 287.0 \ \text{J K}^{-1} \text{kg}^{-1}
$$
<h2>Horizontal Momentum</h2>
$$
\underbrace{\frac{D\mathbf{u}}{Dt}}_{\text{Acceleration (change in momentum)}}
=
\overbrace{
\underbrace{-\frac{1}{\rho} \nabla p}_{\text{Pressure Gradient Force: Initiates motion}}
+
\underbrace{\nu \nabla^2 \mathbf{u}}_{\text{Friction / Viscosity: Opposes motion}}
}^{\text{Surface}}
+
\overbrace{
\underbrace{-g\mathbf{k}}_{\text{Gravity: Stratification / buoyancy}}
}^{\text{Body}}
+
\overbrace{
\underbrace{f v \mathbf{i} - f u \mathbf{j}}_{\text{Coriolis: Modifies motion}}
}^{\text{Apparent}}
$$
For an atmosphere in hydrostatic balance, pressure corresponds to the total mass of the atmosphere above that point.
Consequently, pressure must decrease monotonically with height.
This makes pressure a good choice as a vertical coordinate.
The use of pressure as a vertical coordinate greatly simplifies the equations of motion, the pressure gradient force simply reduces to the gradient of the geopotential on constant pressure surfaces!<br>
Horizontal momentum equations, <strong>height coordinate</strong>, no viscosity:
$$
\frac{Du}{Dt} = -\frac{1}{\rho} \left( \frac{\partial p}{\partial x} \right)_z + fv
$$
$$
\frac{Dv}{Dt} = -\frac{1}{\rho} \left( \frac{\partial p}{\partial y} \right)_z - fu
$$
$$
\boxed{ \frac{D\mathbf{u}}{Dt} + f \mathbf{k} \times \mathbf{u} = -\frac{1}{\rho} \nabla_z p }
$$
Horizontal momentum equations, <strong>pressure coordinate</strong>, no viscosity:
$$
\frac{Du}{Dt} = - \left( \frac{\partial \Phi}{\partial x} \right)_p + fv
$$
$$
\frac{Dv}{Dt} = - \left( \frac{\partial \Phi}{\partial y} \right)_p - fu
$$
$$
\boxed{ \frac{D\mathbf{u}}{Dt} + f \mathbf{k} \times \mathbf{u} = -\nabla_p \Phi }
$$
Observe density no longer explicitly present.<br>
Air is compressible, altitude difficult to measure, and air cools as it expands.
Thus, atmospheric scientists typically describe the atmosphere in terms of
potential temperature (to remove the adiabatic effect of expansion which affects
regular temperature), and pressure instead of depth.
<h2>The Overturning Circulation</h2>
<img src="Figures/Lecture9_HadleyCellOverturning.png" alt="HadleyCellOverturning" style="margin-top: 1.5em; max-width: 100%; display: block;">
<div class="figure-caption">(A simple model of the Hadley
Cell. Rising air near the equator moves
poleward near the tropopause, descending in the subtropics and returning.
The poleward moving air conserves its angular momentum, leading to a shear of
the zonal wind that increases away from
the equator. By thermal wind the temperature of the air falls as it moves poleward,
and to satisfy the thermodynamic budget
it sinks in the subtropics.
Credit: G. Vallis, 2017.)</div>
Hadley Cell (Assumptions):
<ul>
<li>1. The circulation is steady</li>
<li>2. The polewards moving air conserves its angular momentum</li>
<li>3. The equatorwards moving flow is frictionally retarded and weak</li>
<li>4. The circulation is in thermal wind balance</li>
<li>5. There is solar heating from above, symmetric about the equator</li>
<li>6. The flow is symmetric about the equator</li>
<li>7. Potential temperature is conserved</li>
<li>8. The zonal velocity is zero at the equator (rising motion)</li>
</ul>
Food for thought: <a href="https://arxiv.org/pdf/1306.2418">Atmospheric Circulation of Terrestrial Exoplanets</a>
<h2>The Held-Hou Model of the Hadley Circulation</h2>
<ul>
<li>Uses thermodynamic constraints to predict the width and strength of the Hadley Cell</li>
<li>Angular Momentum is conserved</li>
<li>Thermal wind balance needs to be true</li>
<li>Predicts the temperature from the angular momentum balance</li>
<li>Theory predicts the poleward edge of the Hadley Cell to be around 20°</li>
</ul>
<script>
function copyToClipboard(btn) {
const code = btn.nextElementSibling.innerText;
navigator.clipboard.writeText(code).then(() => {
btn.classList.add('copied');
setTimeout(() => btn.classList.remove('copied'), 2000);
});
}
</script>
<p><strong>Next</strong>: <a href="Project5.html">Project 5 – Hadley Cell</a></p>
<div class="footer">
Qi-fan based on the course material by Markus Jochum & Guido Vettoretti <br>
It is very much a work in progress! Have you spotted a mistake or an error on this page?
Click <a href="mailto:qifan.wu@nbi.ku.dk">here</a> to tell me!<br>
© 2025 TeamOcean | NBI/KU
</div>
</div>
<script>
function copyToClipboard(btn) {
const code = btn.nextElementSibling.innerText;
navigator.clipboard.writeText(code).then(() => {
btn.classList.add('copied');
setTimeout(() => btn.classList.remove('copied'), 2000);
});
}
</script>
</body>
</html>