পৃথিবী ও মহাকাশ বিজ্ঞান · বই 3 · স্নাতক প্রথম বর্ষ

University Earth and Space Sciences — Year 1

বইগুলো এখনো বাংলা ভাষায় পাওয়া যায় না — ইংরেজি সংস্করণ দেখানো হচ্ছে।

University Earth and Space Sciences — Year 1 · স্নাতক প্রথম বর্ষ

26Structure of the Atmosphere and the Ocean

Twice a day, at weather stations around the world, a balloon is released carrying a radiosonde, a small instrument that radios back the pressure, the temperature and the humidity of the air as it rises; it climbs at about 300 m300\,\mathrm{m} a minute and can pass 35 km35\,\mathrm{km} before the balloon bursts. Far out at sea, a research ship lowers a frame of instruments and water bottles on a cable to the abyss, measuring the conductivity, the temperature and the depth all the way. The two profiles, one of the air and one of the water, are the raw material of this chapter. Both envelopes are layered, and in both the layering is set by density: light fluid above heavy fluid. What makes a layer light or heavy is less obvious than it looks. Air that is colder than the air below it can still be the lighter; water that is warmer than the water above it can still be the heavier. The chapter builds the two tools that sort this out, the potential temperature and the equation of state of seawater.

A CTD and rosette being brought back on board in the Bering Strait in June 2011. The grey bottles round the frame are closed at chosen depths on the way up to sample the water; the instruments at the bottom of the frame measure conductivity (from which salinity is computed), temperature and depth (in fact pressure) continuously as the frame is lowered. (NASA Goddard Space Flight Center, public domain, via Wikimedia Commons.)
A CTD and rosette being brought back on board in the Bering Strait in June 2011. The grey bottles round the frame are closed at chosen depths on the way up to sample the water; the instruments at the bottom of the frame measure conductivity (from which salinity is computed), temperature and depth (in fact pressure) continuously as the frame is lowered. (NASA Goddard Space Flight Center, public domain, via Wikimedia Commons.)

You already know

The physics volumes give the hydrostatic relation, dp/dz=−ρg\mathrm dp/\mathrm dz = -\rho g, for a fluid at rest; for an ideal gas of molar mass MM, p=ρRT/Mp = \rho R T/M, and in an isothermal atmosphere the pressure falls as e−z/He^{-z/H} with the scale height H=RT/(Mg)H = RT/(Mg). They state the first law for a closed system and, for an ideal gas, Laplace’s relation along a reversible adiabatic transformation, pVγ=pV^{\gamma} = constant, with γ=cp/cv\gamma = c_p/c_v and cp−cv=R/Mc_p - c_v = R/M per unit mass. A body immersed in a fluid receives an upward thrust equal to the weight of fluid it displaces (Archimedes). Book 1 named the layers of the atmosphere and the environmental lapse rate of the standard atmosphere, 6.5 ∘C6.5\,{}^{\circ}\mathrm{C} per kilometre in the troposphere; Book 2 admitted the dry adiabatic rate, 9.8 ∘C9.8\,{}^{\circ}\mathrm{C} per kilometre, called air stable or unstable according to the fate of a parcel pushed upwards, and described the mixed layer, the thermocline, the halocline and the pycnocline of the ocean. The salinity of seawater is close to 35 grams of dissolved salts per kilogram.

26.1 The vertical structure of the atmosphere

No real sounding is the same as the next one. To compare them, a reference profile is fixed by convention.

Definition 26.1 (Standard atmosphere)

A standard atmosphere is a profile of temperature, pressure and density against height fixed by convention: dry air of constant molar mass, an ideal gas in hydrostatic balance, with given values at sea level and a temperature that varies linearly with height in successive layers. The 1976 standard atmosphere takes 288.15 K288.15\,\mathrm{K} and 101 325 Pa101\,325\,\mathrm{Pa} at sea level, M0=28.9644 kg kmol−1M_0 = 28.9644\,\mathrm{kg}\,\mathrm{kmol}^{-1}, R∗=8314.32 J kmol−1 K−1R^* = 8314.32\,\mathrm{J}\,\mathrm{kmol}^{-1}\,\mathrm{K}^{-1} and g0=9.806 65 m s−2g_0 = 9.806\,65\,\mathrm{m}\,\mathrm{s}^{-2}.

Below 47 km47\,\mathrm{km} the 1976 standard has four layers (Table 26.1). The lowest, the troposphere, cools at 6.5 K6.5\,\mathrm{K} per kilometre; from its top, the tropopause, the temperature stays at 216.65 K216.65\,\mathrm{K} for nine kilometres, then rises. The heights of the standard are geopotential heights, the heights that a constant gravity g0g_0 would give the same potential energy; below 32 km32\,\mathrm{km} they differ from geometric heights by less than 0.6 %0.6\,\% (11 km11\,\mathrm{km} of the standard lie at 11.02 km11.02\,\mathrm{km}).

LayerBase (km)dT/dz\mathrm dT/\mathrm dz (K/km)TT at base (K)pp at base (hPa)
1 (troposphere)0−6.5-6.5288.151013.25
2 (stratosphere)110.0216.65226.32
3 (stratosphere)20+1.0+1.0216.6554.75
4 (stratosphere)32+2.8+2.8228.658.68
Table 26.1. The four lowest layers of the 1976 standard atmosphere. The gradients and the sea-level values are defining values; the temperatures and pressures at the bases are computed from them with Proposition 26.2.

In each layer the temperature is a linear function of height, so the hydrostatic relation can be integrated in closed form.

Proposition 26.2 (Pressure in a layer of constant gradient)

In a layer of dry air where T=Tb−L (z−zb)T = T_b - L\,(z - z_b) with a constant L≠0L \neq 0, the pressure is

p=pb(TTb)gM/(RL),p = p_b \left(\frac{T}{T_b}\right)^{gM/(RL)},

where TbT_b and pbp_b are the values at the base zbz_b. When L=0L = 0 the pressure falls exponentially with the scale height H=RTb/(Mg)H = RT_b/(Mg).

Proof. The hydrostatic relation and the ideal gas give dp/p=−(gM/R) dz/T\mathrm dp/p = -(gM/R)\,\mathrm dz/T. Since dT=−L dz\mathrm dT = -L\,\mathrm dz, this is dp/p=(gM/(RL)) dT/T\mathrm dp/p = (gM/(RL))\,\mathrm dT/T, which integrates to ln⁡(p/pb)=(gM/(RL))ln⁡(T/Tb)\ln(p/p_b) = (gM/(RL)) \ln(T/T_b). For L=0L = 0, T=TbT = T_b and ln⁡(p/pb)=−(z−zb)/H\ln(p/p_b) = -(z - z_b)/H: the isothermal atmosphere of the physics volumes. ∎

In the troposphere of the standard, gM/(RL)=9.80665×28.9644/(8314.32×0.0065)=5.256gM/(RL) = 9.80665 \times 28.9644 / (8314.32 \times 0.0065) = 5.256, and at 11 km11\,\mathrm{km}, T=216.65T = 216.65 K: p=101 325×(216.65/288.15)5.256=22 632 Pap = 101\,325 \times (216.65/288.15)^{5.256} = 22\,632\,\mathrm{Pa}, a little under a quarter of the sea-level value. In the isothermal layer above, H=8314.32×216.65/(28.9644×9.80665)=6.34 kmH = 8314.32 \times 216.65/(28.9644 \times 9.80665) = 6.34\,\mathrm{km}, and nine kilometres higher the pressure has fallen by a further factor e9/6.34=4.13e^{9/6.34} = 4.13, to 5475 Pa5475\,\mathrm{Pa}. Since the pressure at a level is the weight, per unit area, of the air above it, half of the mass of the standard atmosphere lies below 5.5 km5.5\,\mathrm{km} and nine tenths below 16 km16\,\mathrm{km}.

26.2 Lapse rates and potential temperature

Definition 26.3 (Environmental and dry adiabatic lapse rates)

The environmental lapse rate is Γ=−dT/dz\Gamma = -\mathrm dT/\mathrm dz, the rate at which the temperature of the air at rest decreases with height, as a sounding measures it. The dry adiabatic lapse rate Γd\Gamma_d is the rate at which the temperature of a parcel of unsaturated air decreases as the parcel is lifted.

A lapse rate is positive when the temperature falls with height. The two rates are of different natures: the first describes the state of the atmosphere at one moment; the second, the transformation of a piece of air that moves.

Definition 26.4 (Air parcel)

An air parcel is a volume of air followed in its motion, small enough to have a single temperature and pressure and not to disturb the air around it, large enough to keep its identity: it exchanges no mass and, during a quick displacement, no heat with its surroundings, and its pressure is at every moment that of the surrounding air.

The parcel is an idealisation, like the closed system of thermodynamics: the results below hold for real air to the extent that a moving piece of it exchanges little heat and little mass with its surroundings.

Theorem 26.5 (The dry adiabatic lapse rate)

A parcel of dry air, lifted adiabatically and reversibly through air at rest of nearly the same temperature, cools at the rate

Γd=gcp,\Gamma_d = \frac{g}{c_p},

where cpc_p is the specific heat of air at constant pressure. With cp=1004.7 J kg−1 K−1c_p = 1004.7\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, Γd=9.76 K km−1\Gamma_d = 9.76\,\mathrm{K}\,\mathrm{km}^{-1}.

Proof. For a kilogram of the parcel, of temperature TpT_p and density ρp\rho_p, the first law of a reversible transformation without heat exchange, written with the enthalpy, is cp dTp−dp/ρp=0c_p\,\mathrm dT_p - \mathrm dp/\rho_p = 0. The pressure of the parcel is that of the air around it, which is hydrostatic: dp=−ρeg dz\mathrm dp = -\rho_e g\,\mathrm dz, with ρe\rho_e the density of the surrounding air. Hence

dTpdz=−gcp ρeρp=−gcp TpTe,\frac{\mathrm dT_p}{\mathrm dz} = -\frac{g}{c_p}\,\frac{\rho_e}{\rho_p} = -\frac{g}{c_p}\,\frac{T_p}{T_e},

since at equal pressure the densities of two ideal gases of the same molar mass are in the inverse ratio of their temperatures. When the parcel is at nearly the temperature of its surroundings, Tp/Te≈1T_p/T_e \approx 1 and −dTp/dz=g/cp-\mathrm dT_p/\mathrm dz = g/c_p. For dry air, an ideal diatomic gas, γ=1.4\gamma = 1.4, so cp=γ/(γ−1)×R/M=3.5×8314.32/28.9644=1004.7 J kg−1 K−1c_p = \gamma/(\gamma - 1) \times R/M = 3.5 \times 8314.32/28.9644 = 1004.7\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, and Γd=9.80665/1004.7=0.009 76 K m−1\Gamma_d = 9.80665/1004.7 = 0.009\,76\,\mathrm{K}\,\mathrm{m}^{-1}. ∎

A parcel 10 K10\,\mathrm{K} warmer than its surroundings, near 280 K280\,\mathrm{K}, cools at Γd×290/280\Gamma_d \times 290/280, 3.6 %3.6\,\% faster; the correction is neglected in what follows. The rate does not depend on the pressure or the temperature of the parcel, and a parcel that sinks warms at the same rate, as Book 2 followed air over a mountain.

Integrating the same equation without the hydrostatic step gives a quantity that the parcel carries with it.

Proposition 26.6 (Potential temperature)

Along a reversible adiabatic transformation of dry air,

θ=T(p0p)κ,κ=RMcp=γ−1γ=27,\theta = T\left(\frac{p_0}{p}\right)^{\kappa},\qquad \kappa = \frac{R}{Mc_p} = \frac{\gamma - 1}{\gamma} = \frac{2}{7},

stays constant, whatever the reference pressure p0p_0.

Proof. From cp dT=dp/ρc_p\,\mathrm dT = \mathrm dp/\rho and ρ=pM/(RT)\rho = pM/(RT), cp dT/T=(R/M) dp/pc_p\,\mathrm dT/T = (R/M)\,\mathrm dp/p, so ln⁡T−κln⁡p\ln T - \kappa \ln p is constant, and so is Tp−κT p^{-\kappa}; multiplying by the constant p0κp_0^{\kappa} gives θ\theta. This is Laplace’s relation written with TT and pp. ∎

Definition 26.7 (Potential temperature)

The potential temperature θ\theta of a parcel of dry air is the temperature it would have if brought adiabatically to the reference pressure p0=1000 hPap_0 = 1000\,\mathrm{hPa}: θ=T (1000/p)2/7\theta = T\,(1000/p)^{2/7} with pp in hectopascals and TT in kelvins.

The potential temperature is a label that a dry parcel keeps when it rises or sinks without exchanging heat, as the salinity of a water parcel is. Heating by radiation, by condensation of water or by mixing changes it; the rest does not.

Example 26.8 (The tropopause of the standard atmosphere)

At the surface of the standard, θ=288.15×(1000/1013.25)2/7=287.07 K\theta = 288.15 \times (1000/1013.25)^{2/7} = 287.07\,\mathrm{K}. At the tropopause, T=216.65 KT = 216.65\,\mathrm{K} and p=226.32 hPap = 226.32\,\mathrm{hPa}: θ=216.65×(1000/226.32)2/7=216.65×1.529=331.2 K\theta = 216.65 \times (1000/226.32)^{2/7} = 216.65 \times 1.529 = 331.2\,\mathrm{K}. Air from the tropopause brought down adiabatically to 1000 hPa1000\,\mathrm{hPa} would arrive at 331.2 K331.2\,\mathrm{K}, 44 K44\,\mathrm{K} warmer than surface air brought to the same pressure. The temperature falls by 71.5 K71.5\,\mathrm{K} from the ground to the tropopause, yet the potential temperature rises by 44 K44\,\mathrm{K}.

26.3 Static stability

Definition 26.9 (Static stability; temperature inversion)

The static stability of a layer at rest is its response to a small vertical displacement of a parcel. The layer has a stable stratification if the displaced parcel is pushed back towards its level, an unstable stratification if it is pushed further away, and a neutral stratification if no force acts on it. A temperature inversion is a layer in which the temperature increases with height.

Theorem 26.10 (The criterion of static stability)

A layer of dry air at rest, of temperature TT and environmental lapse rate Γ\Gamma, is stable if Γ<Γd\Gamma < \Gamma_d, neutral if Γ=Γd\Gamma = \Gamma_d and unstable if Γ>Γd\Gamma > \Gamma_d; equivalently, it is stable if and only if its potential temperature increases with height. A parcel displaced by a small δz\delta z in a stable layer oscillates about its level at the angular frequency NN, where

N2=gT (Γd−Γ)=gθ dθdz.N^2 = \frac{g}{T}\,(\Gamma_d - \Gamma) = \frac{g}{\theta}\,\frac{\mathrm d\theta}{\mathrm dz}.

Proof. A parcel at the temperature TT of its level is lifted by δz\delta z. It arrives at T−Γd δzT - \Gamma_d\,\delta z, while the air around it is at T−Γ δzT - \Gamma\,\delta z. Per unit mass, the parcel feels its weight and the thrust of Archimedes, a net upward acceleration g (ρe−ρp)/ρp=g (Tp−Te)/Teg\,(\rho_e - \rho_p)/\rho_p = g\,(T_p - T_e)/T_e, the pressures being equal. To first order in δz\delta z this is

d2δzdt2=−gT (Γd−Γ) δz=−N2 δz.\frac{\mathrm d^2 \delta z}{\mathrm dt^2} = -\frac{g}{T}\,(\Gamma_d - \Gamma)\,\delta z = -N^2\,\delta z .

If Γ<Γd\Gamma < \Gamma_d, N2>0N^2 > 0 and the parcel oscillates at the angular frequency NN; if Γ>Γd\Gamma > \Gamma_d the displacement grows exponentially; if Γ=Γd\Gamma = \Gamma_d no force acts. For the second form, take the logarithmic derivative of θ=T(p0/p)κ\theta = T (p_0/p)^{\kappa} and use the hydrostatic relation:

1θdθdz=1TdTdz−κpdpdz=−ΓT+κ ρgp=−ΓT+κgMRT=Γd−ΓT,\frac{1}{\theta}\frac{\mathrm d\theta}{\mathrm dz} = \frac{1}{T}\frac{\mathrm dT}{\mathrm dz} - \frac{\kappa}{p}\frac{\mathrm dp}{\mathrm dz} = -\frac{\Gamma}{T} + \frac{\kappa\,\rho g}{p} = -\frac{\Gamma}{T} + \frac{\kappa g M}{RT} = \frac{\Gamma_d - \Gamma}{T},

since κM/R=1/cp\kappa M/R = 1/c_p. ∎

The angular frequency NN is called the buoyancy frequency. In the troposphere of the standard, near the ground, N2=(9.80665/288.15)×(9.761−6.5)×10−3=1.11×10−4 s−2N^2 = (9.80665/288.15) \times (9.761 - 6.5)\times10^{-3} = 1.11 \times 10^{-4}\,\mathrm{s}^{-2}, N=0.0105 s−1N = 0.0105\,\mathrm{s}^{-1}, and the period 2π/N2\pi/N of the oscillation is about 10 min10\,\mathrm{min}. In the isothermal lower stratosphere Γ=0\Gamma = 0: N=0.0210 s−1N = 0.0210\,\mathrm{s}^{-1}, twice as large, a period of five minutes. The stratosphere is a lid: it resists vertical motion much more strongly than the troposphere below it.

The 1976 standard atmosphere to 32\, km. Left: temperature, with the dry adiabat that a parcel lifted from the surface would follow (dashed). The standard cools more slowly than the adiabat, so it is stable everywhere. Right: potential temperature, which rises with height throughout, gently in the troposphere and steeply above the tropopause (dot).
The 1976 standard atmosphere to 32 km32\,\mathrm{km}. Left: temperature, with the dry adiabat that a parcel lifted from the surface would follow (dashed). The standard cools more slowly than the adiabat, so it is stable everywhere. Right: potential temperature, which rises with height throughout, gently in the troposphere and steeply above the tropopause (dot).

On a diagram of temperature against height, the criterion compares two slopes (figure below): the parcel follows the dry adiabat through its starting point, and the question is on which side of the sounding it finds itself after a small displacement.

The parcel method. Solid blue: the sounding; dashed red: the dry adiabat through P. A parcel starts at P at the temperature of the air around it and is lifted to z + z along the dry adiabat (red dot); the sounding gives the temperature of the air there (blue dot). Left: the sounding cools more slowly than the adiabat; the parcel arrives colder, hence denser, than its surroundings and sinks back. Right: the sounding cools faster than the adiabat; the parcel arrives warmer and keeps rising.
The parcel method. Solid blue: the sounding; dashed red: the dry adiabat through PP. A parcel starts at PP at the temperature of the air around it and is lifted to z+δzz + \delta z along the dry adiabat (red dot); the sounding gives the temperature of the air there (blue dot). Left: the sounding cools more slowly than the adiabat; the parcel arrives colder, hence denser, than its surroundings and sinks back. Right: the sounding cools faster than the adiabat; the parcel arrives warmer and keeps rising.

Method 26.11 (Reading the stability of a sounding)

  1. Divide the sounding into layers between levels where the slope changes.
  2. In each layer compute Γ=−(T2−T1)/(z2−z1)\Gamma = -(T_2 - T_1)/(z_2 - z_1), in kelvins per kilometre.
  3. Compare with Γd=9.8 K/km\Gamma_d = 9.8\,\mathrm{K}/\mathrm{km}: Γ>Γd\Gamma > \Gamma_d unstable, Γ\Gamma close to Γd\Gamma_d neutral, Γ<Γd\Gamma < \Gamma_d stable; a negative Γ\Gamma is an inversion, the most stable case.
  4. Where pressures are given rather than heights, compute θ=T(1000/p)2/7\theta = T(1000/p)^{2/7} at each level: θ\theta decreasing upwards marks an unstable layer, constant a neutral one, increasing a stable one.

An unstable layer does not last: a displacement in it grows exponentially, at the rate −N2\sqrt{-N^2}, and the parcels it releases overturn it, mixing it until its potential temperature is uniform and its lapse rate is Γd\Gamma_d. A layer heated from below, by warm ground or a warm sea, is pushed towards instability; one cooled from below, towards stability. An inversion caps whatever lies below it: a parcel rising from below finds itself colder than its surroundings as soon as it enters it. Above 20 km20\,\mathrm{km} the whole stratosphere of the standard is an inversion.

All of this assumes dry air. A rising parcel that becomes saturated releases latent heat as its water condenses and cools more slowly than Γd\Gamma_d; the saturated adiabat and the conditional instability it brings are the subject of next year’s volume.

26.4 The layered ocean

The ocean exchanges heat and water with the atmosphere only at its surface, through heating and cooling, evaporation and rain. Its vertical structure has three ingredients, temperature, salinity and density.

Definition 26.12 (Mixed layer, thermocline, halocline, pycnocline)

The mixed layer is the surface layer of the ocean whose temperature, salinity and density are nearly uniform with depth; its depth is read on a profile as the depth where the temperature first differs from its surface value by a chosen threshold, often 0.5 ∘C0.5\,{}^{\circ}\mathrm{C}. The thermocline, the halocline and the pycnocline are the layers below it in which the temperature, the salinity and the density respectively change fastest with depth.

The pressure in the ocean is quoted in decibars, because one decibar, 104 Pa10^{4}\,\mathrm{Pa}, is nearly the pressure of one metre of seawater: at 30∘30{}^{\circ} of latitude, 1000 dbar1000\,\mathrm{dbar} lie at 990.8 m990.8\,\mathrm{m} and 10 000 dbar10\,000\,\mathrm{dbar} at 9712.7 m9712.7\,\mathrm{m}.

Water is nearly incompressible, but not quite, and a parcel of water brought up from the deep expands and cools a little, as a parcel of air does. The adiabatic lapse rate of seawater is small, about 0.04 ∘C0.04\,{}^{\circ}\mathrm{C} per 1000 dbar1000\,\mathrm{dbar} at the surface at 0 ∘C0\,{}^{\circ}\mathrm{C} and 0.09 ∘C0.09\,{}^{\circ}\mathrm{C} per 1000 dbar1000\,\mathrm{dbar} at 3000 dbar3000\,\mathrm{dbar}, against 9.89.8 per kilometre in air; but over five kilometres it adds up to a few tenths of a degree, more than the differences between deep waters. Oceanographers therefore compare waters by their potential temperature too: the temperature a parcel would have if raised adiabatically, without change of salinity, to the surface. A water of 2 ∘C2\,{}^{\circ}\mathrm{C} at 4000 dbar4000\,\mathrm{dbar} and salinity 35 has a potential temperature of 1.66 ∘C1.66\,{}^{\circ}\mathrm{C}. In the deepest levels of the South Atlantic at 10.5∘S10.5{}^{\circ}\mathrm{S}, the measured temperature rises from 0.730 ∘C0.730\,{}^{\circ}\mathrm{C} at 5400 m5400\,\mathrm{m} to 0.731 ∘C0.731\,{}^{\circ}\mathrm{C} at 5500 m5500\,\mathrm{m} while the potential temperature falls from 0.261 ∘C0.261\,{}^{\circ}\mathrm{C} to 0.250 ∘C0.250\,{}^{\circ}\mathrm{C}: the warming is compression, not heat.

Potential temperature, salinity and potential density anomaly () in the upper 2000\, m, in the annual mean, at three places: the tropical South Atlantic (10.5 S, 25.5 W), the subtropical South Atlantic (30.5 S, 25.5 W) and the Canada Basin of the Arctic Ocean (75.5 N, 150.5 W). Data: World Ocean Atlas 2023.
Potential temperature, salinity and potential density anomaly (Definition 26.13) in the upper 2000 m2000\,\mathrm{m}, in the annual mean, at three places: the tropical South Atlantic (10.5∘S10.5{}^{\circ}\mathrm{S}, 25.5∘W25.5{}^{\circ}\mathrm{W}), the subtropical South Atlantic (30.5∘S30.5{}^{\circ}\mathrm{S}, 25.5∘W25.5{}^{\circ}\mathrm{W}) and the Canada Basin of the Arctic Ocean (75.5∘N75.5{}^{\circ}\mathrm{N}, 150.5∘W150.5{}^{\circ}\mathrm{W}). Data: World Ocean Atlas 2023.

The three columns of the figure above show the three kinds of upper ocean. In the tropics, a warm mixed layer near 26.5 ∘C26.5\,{}^{\circ}\mathrm{C} sits on a sharp thermocline: between 100 m100\,\mathrm{m} and 200 m200\,\mathrm{m} the potential temperature drops from 24.0 ∘C24.0\,{}^{\circ}\mathrm{C} to 14.4 ∘C14.4\,{}^{\circ}\mathrm{C}, and the density follows it, so the thermocline is also the pycnocline. At 30.5∘S30.5{}^{\circ}\mathrm{S} the surface is cooler, 21.1 ∘C21.1\,{}^{\circ}\mathrm{C}, and the fall more gradual, spread over the upper 800 m800\,\mathrm{m}. In both, salinity is highest near the surface and decreases downwards; it works against the stratification, which the temperature more than makes up for. The Arctic column is the opposite case. Its surface water, at −1.4 ∘C-1.4\,{}^{\circ}\mathrm{C}, is the coldest of the three, yet the lightest, σθ=22.5\sigma_\theta = 22.5, because its salinity is only 28.128.1; salinity rises to 34.534.5 at 300 m300\,\mathrm{m}. The halocline is the pycnocline here. Below it lies a layer warmer than the water above, near 0.7 ∘C0.7\,{}^{\circ}\mathrm{C} at 475 m475\,\mathrm{m}: a temperature inversion that is perfectly stable, held down by its salt.

26.5 Seawater density and T–S diagrams

Definition 26.13 (Equation of state of seawater; haline contraction coefficient; density anomaly)

The equation of state of seawater gives its density ρ(S,T,p)\rho(S, T, p) as a function of salinity, temperature and pressure. The haline contraction coefficient is β=(1/ρ) ∂ρ/∂S\beta = (1/\rho)\,\partial\rho/\partial S at constant temperature and pressure. The density anomaly σ\sigma is a density minus 1000 kg m−31000\,\mathrm{kg}\,\mathrm{m}^{-3}: σt=ρ(S,T,0)−1000\sigma_t = \rho(S, T, 0) - 1000 uses the measured temperature at zero sea pressure, and the potential density anomaly σθ=ρ(S,θ,0)−1000\sigma_\theta = \rho(S, \theta, 0) - 1000 the potential temperature.

The equation of state is not derived from first principles but fitted to laboratory measurements. The international equation of 1980 at the pressure of the surface is a polynomial with fifteen coefficients, ρ=ρw(T)+A(T) S+B(T) S3/2+C S2\rho = \rho_w(T) + A(T)\,S + B(T)\,S^{3/2} + C\,S^2, where ρw\rho_w is the density of pure water; it gives ρ(35,0 ∘C)=1028.1063 kg m−3\rho(35, 0\,{}^{\circ}\mathrm{C}) = 1028.1063\,\mathrm{kg}\,\mathrm{m}^{-3}. A thermodynamic description of seawater, TEOS-10, was adopted in 2009 to replace it as the official one; the 1980 equation is used here because its polynomial fits on a page. Its thermal expansion coefficient α=−(1/ρ) ∂ρ/∂T\alpha = -(1/\rho)\,\partial\rho/\partial T is the coefficient of the physics volumes.

TT (∘C{}^{\circ}\mathrm{C})ρ(35,T,0)\rho(35, T, 0) (kg m−3\mathrm{kg}\,\mathrm{m}^{-3})α\alpha (K−1\mathrm{K}^{-1})β\beta
01028.1060.53×10−40.53 \times 10^{-4}7.85×10−47.85 \times 10^{-4}
101026.9521.67×10−41.67 \times 10^{-4}7.61×10−47.61 \times 10^{-4}
251023.3432.97×10−42.97 \times 10^{-4}7.39×10−47.39 \times 10^{-4}
Table 26.2. Density and expansion coefficients of seawater of salinity 35 at the surface pressure, from the international equation of state of 1980.

The expansion coefficient is six times larger in warm water than in water near freezing, while β\beta hardly changes (Table 26.2). One unit of salinity weighs as much as about 4.6 ∘C4.6\,{}^{\circ}\mathrm{C} at 10 ∘C10\,{}^{\circ}\mathrm{C}, but as much as 15 ∘C15\,{}^{\circ}\mathrm{C} near 0 ∘C0\,{}^{\circ}\mathrm{C}: in cold polar water, density is set by salinity, as the Arctic column showed.

Proposition 26.14 (Linear equation of state)

Near a reference state (S0,T0)(S_0, T_0) of density ρ0\rho_0, at fixed pressure,

ρ≈ρ0[1−α (T−T0)+β (S−S0)],\rho \approx \rho_0\bigl[1 - \alpha\,(T - T_0) + \beta\,(S - S_0)\bigr],

with α\alpha and β\beta evaluated at the reference state; the error is of second order in T−T0T - T_0 and S−S0S - S_0.

Proof. This is the first-order Taylor expansion of ρ(S,T)\rho(S, T) at (S0,T0)(S_0, T_0), since ∂ρ/∂T=−ρ0α\partial\rho/\partial T = -\rho_0\alpha and ∂ρ/∂S=ρ0β\partial\rho/\partial S = \rho_0\beta there. ∎

About (35,10 ∘C)(35, 10\,{}^{\circ}\mathrm{C}) the linear form gives the density of water at (35.2,11 ∘C)(35.2, 11\,{}^{\circ}\mathrm{C}) within 0.01 kg m−30.01\,\mathrm{kg}\,\mathrm{m}^{-3}, but that of a surface water at (35.8,22.4 ∘C)(35.8, 22.4\,{}^{\circ}\mathrm{C}) 0.74 kg m−30.74\,\mathrm{kg}\,\mathrm{m}^{-3} too high: the equation of state is curved, and its curvature has consequences.

Definition 26.15 (T–S diagram; water mass)

A T–S diagram is a plot of the potential temperature of the samples of a profile against their salinity, with the depths written along the curve and the lines of equal σθ\sigma_\theta drawn underneath. A water mass is a body of water with a common formation history, recognised by its range of potential temperature and salinity.

A water mass acquires its temperature and salinity at the surface, through heating and cooling, evaporation and rain, in the region where it forms; once it sinks below the mixed layer, out of contact with the atmosphere, it can change them only by mixing with its neighbours. Potential temperature and salinity are then tracers.

Proposition 26.16 (Mixing lines are straight)

When masses m1m_1 and m2m_2 of two waters (S1,θ1)(S_1, \theta_1) and (S2,θ2)(S_2, \theta_2) mix completely, the mixture has

S=m1S1+m2S2m1+m2,θ≈m1θ1+m2θ2m1+m2:S = \frac{m_1 S_1 + m_2 S_2}{m_1 + m_2},\qquad \theta \approx \frac{m_1 \theta_1 + m_2 \theta_2}{m_1 + m_2}:

on a T–S diagram it lies on the segment joining the two waters, at distances from them in the inverse ratio of their masses. Mixtures of three waters fill the triangle they form.

Proof. The salinity is a mass fraction of salt, and salt is conserved: the mass of salt in the mixture, m1S1+m2S2m_1 S_1 + m_2 S_2, divided by its mass, gives SS. Mixing at constant pressure conserves the enthalpy; with a specific heat that varies little, this is a weighted mean of the potential temperatures, to within a relative error of about 3.4×10−43.4\times10^{-4} per unit of salinity difference. The point (S,θ)(S, \theta) is then the barycentre of the two points with the weights m1m_1 and m2m_2. ∎

T–S diagram of the two South Atlantic columns of the profiles figure, from the surface to the sea floor, over lines of equal potential density anomaly _ (grey, labelled in kg\, m-3 at their upper ends). Dots: the cores of the water masses read on the column at 30.5 S. Blue: two waters of equal density, _ = 27.0, at 2\, C and 14\, C, and their equal-parts mixture, on the straight line between them and denser than either. Data: World Ocean Atlas 2023.
T–S diagram of the two South Atlantic columns of the profiles figure, from the surface to the sea floor, over lines of equal potential density anomaly σθ\sigma_\theta (grey, labelled in kg m−3\mathrm{kg}\,\mathrm{m}^{-3} at their upper ends). Dots: the cores of the water masses read on the column at 30.5∘S30.5{}^{\circ}\mathrm{S}. Blue: two waters of equal density, σθ=27.0\sigma_\theta = 27.0, at 2 ∘C2\,{}^{\circ}\mathrm{C} and 14 ∘C14\,{}^{\circ}\mathrm{C}, and their equal-parts mixture, on the straight line between them and denser than either. Data: World Ocean Atlas 2023.

The two South Atlantic columns draw the same figure on the T–S diagram: the similarity of curves from places far apart is what reveals water masses. From the surface down to about 500 m500\,\mathrm{m}, the curves follow a nearly straight line, steep in temperature: this is Central Water, which fills the upper kilometre of the subtropical ocean. Below it, both curves turn to a salinity minimum, 34.3034.30 at 900 m900\,\mathrm{m} and 4.3 ∘C4.3\,{}^{\circ}\mathrm{C} at 30.5∘S30.5{}^{\circ}\mathrm{S}: Antarctic Intermediate Water (AAIW), formed under the rainy westerlies of the Southern Ocean, sinking to about 1000 m1000\,\mathrm{m} and spreading north. Deeper, a salinity maximum, 34.9034.90 at 2900 m2900\,\mathrm{m}: North Atlantic Deep Water (NADW), formed in the Greenland and Labrador Seas, sinking to about 3000 m3000\,\mathrm{m} and moving south through the whole Atlantic. At the bottom the curves converge on the coldest water, 0.8 ∘C0.8\,{}^{\circ}\mathrm{C} at 4900 m4900\,\mathrm{m}: Antarctic Bottom Water (AABW), formed along the Antarctic continent, which fills the oceans below about 4000 m4000\,\mathrm{m}. At 10.5∘S10.5{}^{\circ}\mathrm{S} the salinity minimum lies at 750 m750\,\mathrm{m} and the maximum at 1800 m1800\,\mathrm{m}, with a higher maximum, 34.9534.95: the NADW is closer to its source there.

The extremes are themselves mixing products. AAIW, as it moves north, mixes with the saltier water above and below it; by Proposition 26.16 its core moves towards them on the diagram and the minimum erodes. A salinity minimum or maximum on a profile is always the trace of a water that came from elsewhere, carried sideways at its own density level.

Method 26.17 (Reading water masses on a T–S diagram)

  1. Convert the measured temperatures to potential temperatures (in the upper kilometre the difference is below 0.1 ∘C0.1\,{}^{\circ}\mathrm{C} and may be neglected).
  2. Plot θ\theta against SS, writing the depths along the curve.
  3. Read the straight segments as mixtures of the two waters at their ends, the extremes of salinity or temperature as the cores of water masses, and the end of the curve near the bottom as the densest water present.
  4. Read the density on the isopycnals: along a stable profile, σθ\sigma_\theta increases with depth.

Because the isopycnals are curved, two waters of the same density do not mix to a water of that density. Water at (33.79,2 ∘C)(33.79, 2\,{}^{\circ}\mathrm{C}) and water at (36.05,14 ∘C)(36.05, 14\,{}^{\circ}\mathrm{C}) both have σθ=27.00\sigma_\theta = 27.00; mixed in equal parts, they give (34.92,8 ∘C)(34.92, 8\,{}^{\circ}\mathrm{C}), of σθ=27.21\sigma_\theta = 27.21 (figure above). The mixture is denser than both parents and sinks below them. This densification on mixing, called cabbeling, is the cost of the curvature of the equation of state: with the linear form of Proposition 26.14, the density of a mixture would be the mean of the densities. Along the isopycnal surfaces of the ocean, where eddies stir waters of different temperature and salinity, cabbeling drives a slow sinking across them.

A last caution. Referred to the surface, the AABW at 30.5∘S30.5{}^{\circ}\mathrm{S} comes out at σθ=27.854\sigma_\theta = 27.854, slightly lighter than the NADW above it, 27.85727.857; at the bottom of the column at 10.5∘S10.5{}^{\circ}\mathrm{S}, σθ\sigma_\theta even decreases downwards. The ratio α/β\alpha/\beta of the equation of state changes with pressure, so that the order of two waters of different temperatures can depend on the pressure at which their densities are compared. Waters thousands of metres deep are compared with a density referred to a pressure nearer their own, such as σ2\sigma_2, referred to 2000 dbar2000\,\mathrm{dbar}; the method is that of next year’s volume.

26.6 Exercises

Exercise 26.1 ★

(a) Compute the scale height of air at 250 K250\,\mathrm{K} with the constants of the standard atmosphere. (b) By what factor does the pressure fall over 10 km10\,\mathrm{km} in an isothermal layer at that temperature? (c) Why is the scale height smaller at the tropopause than at the ground?

Solution

Solution of Exercise 26.1.

(a) H=R∗T/(M0g0)=8314.32×250/(28.9644×9.80665)=7.32 kmH = R^*T/(M_0 g_0) = 8314.32 \times 250/(28.9644 \times 9.80665) = 7.32\,\mathrm{km}. (b) e10/7.32=3.92e^{10/7.32} = 3.92: the pressure falls to about a quarter. (c) HH is proportional to the temperature: at 216.65 K216.65\,\mathrm{K} it is 6.34 km6.34\,\mathrm{km}, against 8.43 km8.43\,\mathrm{km} at 288.15 K288.15\,\mathrm{K}. Cold air is denser at a given pressure, so the pressure falls faster with height.

Exercise 26.2 ★

A sounding gives, in kelvins per kilometre, the lapse rates Γ=11.0\Gamma = 11.0 (0 to 0.3 km0\text{ to }0.3\,\mathrm{km}), 9.79.7 (0.3 to 1.2 km0.3\text{ to }1.2\,\mathrm{km}), −4.0-4.0 (1.2 to 1.5 km1.2\text{ to }1.5\,\mathrm{km}), 6.06.0 (1.5 to 8 km1.5\text{ to }8\,\mathrm{km}) and 0.00.0 (above 8 km8\,\mathrm{km}). Classify each layer as unstable, neutral, stable or an inversion, and say which is the most stable.

Solution

Solution of Exercise 26.2.

With Γd=9.76 K/km\Gamma_d = 9.76\,\mathrm{K}/\mathrm{km}: 0 to 0.3 km0\text{ to }0.3\,\mathrm{km}, 11.0>9.7611.0 > 9.76, unstable; 0.3 to 1.2 km0.3\text{ to }1.2\,\mathrm{km}, 9.79.7, neutral to within the precision of the data; 1.2 to 1.5 km1.2\text{ to }1.5\,\mathrm{km}, an inversion; 1.5 to 8 km1.5\text{ to }8\,\mathrm{km}, stable; above 8 km8\,\mathrm{km}, isothermal, stable. The inversion is the most stable, since Γd−Γ=13.8\Gamma_d - \Gamma = 13.8 is largest there.

Exercise 26.3 ★

(a) Seawater has a density of 1026.95 kg m−31026.95\,\mathrm{kg}\,\mathrm{m}^{-3}. What is its density anomaly? (b) Of two waters, A with σθ=27.12\sigma_\theta = 27.12 and B with σθ=26.78\sigma_\theta = 26.78, which lies above the other when both are at rest in the same column? (c) Why do oceanographers write σ\sigma rather than ρ\rho?

Solution

Solution of Exercise 26.3.

(a) σ=1026.95−1000=26.95 kg m−3\sigma = 1026.95 - 1000 = 26.95\,\mathrm{kg}\,\mathrm{m}^{-3}. (b) B, the lighter, lies above A. (c) All seawater densities begin with the same digits, 10; σ\sigma keeps only the digits that change from one water to another.

Exercise 26.4 ★

Water P (S=35.40,θ=15.0 ∘C)(S = 35.40, \theta = 15.0\,{}^{\circ}\mathrm{C}) mixes with water Q (S=34.30,θ=4.0 ∘C)(S = 34.30, \theta = 4.0\,{}^{\circ}\mathrm{C}). (a) Compute the salinity and the potential temperature of a mixture of one part of P with three parts of Q. (b) Where does it lie on the T–S diagram? (c) A third water R (34.90,2.5 ∘C)(34.90, 2.5\,{}^{\circ}\mathrm{C}) is added. Can the mixture of the three have (34.60,8.0 ∘C)(34.60, 8.0\,{}^{\circ}\mathrm{C})?

Solution

Solution of Exercise 26.4.

(a) S=0.25×35.40+0.75×34.30=34.575S = 0.25 \times 35.40 + 0.75 \times 34.30 = 34.575; θ=0.25×15.0+0.75×4.0=6.75 ∘C\theta = 0.25 \times 15.0 + 0.75 \times 4.0 = 6.75\,{}^{\circ}\mathrm{C}. (b) On the segment PQ, a quarter of the way from Q to P (Proposition 26.16). (c) The mixtures of P, Q and R fill the triangle PQR. Writing (34.60,8.0)=a P+b Q+c R(34.60, 8.0) = a\,\mathrm P + b\,\mathrm Q + c\,\mathrm R with a+b+c=1a + b + c = 1 gives a=0.345a = 0.345, b=0.788b = 0.788, c=−0.133c = -0.133: a negative proportion of R, so the point lies outside the triangle and no mixture of the three can produce it.

Exercise 26.5 ★★

(a) At 500 hPa500\,\mathrm{hPa} the air is at −20 ∘C-20\,{}^{\circ}\mathrm{C}. Compute its potential temperature. (b) At the ground, at 1000 hPa1000\,\mathrm{hPa}, the air is at 25 ∘C25\,{}^{\circ}\mathrm{C}. Is the layer between the two stable for dry air? (c) What would be the temperature of the 500 hPa500\,\mathrm{hPa} air brought down adiabatically to the ground?

Solution

Solution of Exercise 26.5.

(a) θ=253.15×22/7=308.6 K\theta = 253.15 \times 2^{2/7} = 308.6\,\mathrm{K}. (b) At 1000 hPa1000\,\mathrm{hPa}, θ=T=298.15 K\theta = T = 298.15\,\mathrm{K}. The potential temperature increases upwards, from 298.2 to 308.6 K308.6\,\mathrm{K}: the layer is stable for dry air. (c) It would arrive at its potential temperature, 308.6 K308.6\,\mathrm{K}, that is 35.4 ∘C35.4\,{}^{\circ}\mathrm{C}, ten degrees warmer than the air already there.

Exercise 26.6 ★★

(a) With Proposition 26.2, compute the pressure at 5 km5\,\mathrm{km} in the standard atmosphere. (b) Compute it again assuming an isothermal atmosphere at 288.15 K288.15\,\mathrm{K}. Which is lower, and why? (c) Compute the pressure at 20 km20\,\mathrm{km} from the values at 11 km11\,\mathrm{km} of Table 26.1.

Solution

Solution of Exercise 26.6.

(a) T=288.15−6.5×5=255.65 KT = 288.15 - 6.5 \times 5 = 255.65\,\mathrm{K} and p=101 325×(255.65/288.15)5.256=54 020 Pap = 101\,325 \times (255.65/288.15)^{5.256} = 54\,020\,\mathrm{Pa}, 540.2 hPa540.2\,\mathrm{hPa}. (b) H=8.43 kmH = 8.43\,\mathrm{km} and p=101 325 e−5/8.43=56 010 Pap = 101\,325\,e^{-5/8.43} = 56\,010\,\mathrm{Pa}. The real (standard) pressure is lower: the air there is colder than 288.15 K288.15\,\mathrm{K}, hence denser, and the pressure falls faster with height. (c) p=22 632 e−9/6.34=5475 Pap = 22\,632\, e^{-9/6.34} = 5475\,\mathrm{Pa}, the value of Table 26.1.

Exercise 26.7 ★★

Use the linear equation of state about (35,10 ∘C)(35, 10\,{}^{\circ}\mathrm{C}) with the coefficients of Table 26.2. (a) What change of salinity compensates, for the density, a warming of 1 ∘C1\,{}^{\circ}\mathrm{C}? (b) A water at 10 ∘C10\,{}^{\circ}\mathrm{C} is warmed by 2 ∘C2\,{}^{\circ}\mathrm{C} and made saltier by 0.30.3. Does it become denser or lighter, and by how much? (c) Answer the same question for a water at 0 ∘C0\,{}^{\circ}\mathrm{C}, with the coefficients of that temperature.

Solution

Solution of Exercise 26.7.

(a) β ΔS=α ΔT\beta\,\Delta S = \alpha\,\Delta T gives ΔS=1.67/7.61=0.22\Delta S = 1.67/7.61 = 0.22. (b) Δρ=1026.95×(−1.67×10−4×2+7.61×10−4×0.3)=−0.11 kg m−3\Delta\rho = 1026.95 \times (-1.67\times10^{-4} \times 2 + 7.61\times10^{-4} \times 0.3) = -0.11\,\mathrm{kg}\,\mathrm{m}^{-3}: the water becomes lighter. (c) Δρ=1028.11×(−0.53×10−4×2+7.85×10−4×0.3)=+0.13 kg m−3\Delta\rho = 1028.11 \times (-0.53\times10^{-4} \times 2 + 7.85\times10^{-4} \times 0.3) = +0.13\,\mathrm{kg}\,\mathrm{m}^{-3}: the same changes now make it denser, because warming barely changes the density of water near freezing.

Exercise 26.8 ★★

A layer at 260 K260\,\mathrm{K} has a lapse rate of 4 K/km4\,\mathrm{K}/\mathrm{km}. (a) Compute its buoyancy frequency and the period of the oscillation of a parcel displaced in it. (b) Compute the same quantities for an inversion at 270 K270\,\mathrm{K} where Γ=−5 K/km\Gamma = -5\,\mathrm{K}/\mathrm{km}. (c) What happens to the period as Γ\Gamma approaches Γd\Gamma_d?

Solution

Solution of Exercise 26.8.

(a) N2=(9.81/260)×(9.76−4)×10−3=2.17×10−4 s−2N^2 = (9.81/260) \times (9.76 - 4)\times10^{-3} = 2.17 \times 10^{-4}\,\mathrm{s}^{-2}, N=0.0147 s−1N = 0.0147\,\mathrm{s}^{-1}, period 2π/N=426 s2\pi/N = 426\,\mathrm{s}, about seven minutes. (b) N2=(9.81/270)×14.76×10−3=5.36×10−4 s−2N^2 = (9.81/270) \times 14.76\times10^{-3} = 5.36 \times 10^{-4}\,\mathrm{s}^{-2}, N=0.0232 s−1N = 0.0232\,\mathrm{s}^{-1}, period 271 s271\,\mathrm{s}. (c) NN tends to zero and the period grows without bound: the restoring force vanishes in a neutral layer.

Exercise 26.9 ★★

On the profile at 30.5∘S30.5{}^{\circ}\mathrm{S}, the AAIW core has θ=4.31 ∘C\theta = 4.31\,{}^{\circ}\mathrm{C}, S=34.30S = 34.30, σθ=27.20\sigma_\theta = 27.20, and the Central Water at 200 m200\,\mathrm{m} has θ=15.43 ∘C\theta = 15.43\,{}^{\circ}\mathrm{C}, S=35.54S = 35.54, σθ=26.30\sigma_\theta = 26.30. (a) Why does the AAIW lie below the Central Water, although it is fresher? (b) Using the linear equation of state about (35,10 ∘C)(35, 10\,{}^{\circ}\mathrm{C}), estimate the difference of σθ\sigma_\theta between them, and compare with the exact difference. (c) Which of the two effects, temperature or salinity, wins?

Solution

Solution of Exercise 26.9.

(a) It is denser, σθ=27.20\sigma_\theta = 27.20 against 26.30: the AAIW is fresher but much colder. (b) Δσ=1026.95×[−1.67×10−4×(4.31−15.43)+7.61×10−4×(34.30−35.54)]=1.907−0.969=0.94\Delta\sigma = 1026.95 \times [-1.67\times10^{-4} \times (4.31 - 15.43) + 7.61\times10^{-4} \times (34.30 - 35.54)] = 1.907 - 0.969 = 0.94, against 0.900.90 exactly. (c) Temperature: the cooling adds 1.91 kg m−31.91\,\mathrm{kg}\,\mathrm{m}^{-3}, the freshening removes 0.97 kg m−30.97\,\mathrm{kg}\,\mathrm{m}^{-3}.

Exercise 26.10 ★★★

A layer of dry air has a uniform potential temperature θ\theta. (a) Show that its temperature decreases with height at the rate Γd\Gamma_d. (b) Show that the pressure in it is p=ps (1−Γdz/Ts)cpM/Rp = p_s\,(1 - \Gamma_d z/T_s)^{c_pM/R}, where TsT_s and psp_s are the values at z=0z = 0. (c) Such a layer has a top: at what height would the pressure vanish, for Ts=300 KT_s = 300\,\mathrm{K}? Why is that height never reached?

Solution

Solution of Exercise 26.10.

(a) By Theorem 26.10, (g/θ) dθ/dz=(g/T)(Γd−Γ)(g/\theta)\,\mathrm d\theta/\mathrm dz = (g/T)(\Gamma_d - \Gamma); dθ/dz=0\mathrm d\theta/\mathrm dz = 0 gives Γ=Γd\Gamma = \Gamma_d. (b) T=Ts−ΓdzT = T_s - \Gamma_d z, a layer of constant gradient L=ΓdL = \Gamma_d: by Proposition 26.2, p=ps(T/Ts)gM/(RΓd)p = p_s (T/T_s)^{gM/(R\Gamma_d)}, and gM/(RΓd)=gMcp/(Rg)=cpM/RgM/(R\Gamma_d) = gMc_p/(Rg) = c_pM/R (=3.5= 3.5 for dry air). (c) pp vanishes where T=0T = 0, at z=Ts/Γd=300/9.76=30.7 kmz = T_s/\Gamma_d = 300/9.76 = 30.7\,\mathrm{km}. Real layers of uniform potential temperature are a few kilometres deep at most: they are capped by stable air, as the troposphere is by the stratosphere.

Exercise 26.11 ★★★

Two waters lie on the isopycnal σθ=27.00\sigma_\theta = 27.00: A at 2 ∘C2\,{}^{\circ}\mathrm{C}, S=33.787S = 33.787, and B at 14 ∘C14\,{}^{\circ}\mathrm{C}, S=36.050S = 36.050. (a) Compute the salinity and potential temperature of a mixture of equal parts. (b) Show that, with a linear equation of state, the mixture would have σθ=27.00\sigma_\theta = 27.00. (c) With the exact equation the mixture has σθ=27.21\sigma_\theta = 27.21. Using the coefficients α\alpha and β\beta at 2 ∘C2\,{}^{\circ}\mathrm{C} and 14 ∘C14\,{}^{\circ}\mathrm{C}, explain why it is denser. (d) Why is the effect largest for waters of very different temperatures?

Solution

Solution of Exercise 26.11.

(a) S=(33.787+36.050)/2=34.919S = (33.787 + 36.050)/2 = 34.919, θ=8.0 ∘C\theta = 8.0\,{}^{\circ}\mathrm{C}. (b) With a linear equation of state, σ\sigma is an affine function of SS and θ\theta, so the σ\sigma of the mixture is the mean of the two, 27.00. (c) α\alpha is 0.75×10−4 K−10.75 \times 10^{-4}\,\mathrm{K}^{-1} at 2 ∘C2\,{}^{\circ}\mathrm{C} and 2.07×10−4 K−12.07 \times 10^{-4}\,\mathrm{K}^{-1} at 14 ∘C14\,{}^{\circ}\mathrm{C}, while β\beta changes little (7.80×10−47.80 \times 10^{-4} and 7.53×10−47.53 \times 10^{-4}). The density falls ever faster as the temperature rises: it is a concave function of temperature. Along the mixing line, the mixture is at the mean temperature, where the density lies above the chord joining the two parents, that is above 27.00. (d) The excess comes from the curvature, a second-order term: it grows as the square of the temperature difference of the parents.

Exercise 26.12 ★★★

At the bottom of the South Atlantic at 10.5∘S10.5{}^{\circ}\mathrm{S}, the measured temperature is 0.7299 ∘C0.7299\,{}^{\circ}\mathrm{C} at 5400 m5400\,\mathrm{m} and 0.7313 ∘C0.7313\,{}^{\circ}\mathrm{C} at 5500 m5500\,\mathrm{m}. The pressures there are 5500 dbar5500\,\mathrm{dbar} and 5604 dbar5604\,\mathrm{dbar}. (a) With an adiabatic lapse rate of 0.129 ∘C0.129\,{}^{\circ}\mathrm{C} per 1000 dbar1000\,\mathrm{dbar}, how much would a parcel taken from 5500 m5500\,\mathrm{m} to 5400 m5400\,\mathrm{m} cool? (b) Is the layer statically stable? (c) Over the whole column, the adiabatic gradient averages 0.086 ∘C0.086\,{}^{\circ}\mathrm{C} per 1000 dbar1000\,\mathrm{dbar}: estimate the potential temperature of the bottom water and compare with the 0.250 ∘C0.250\,{}^{\circ}\mathrm{C} computed by the exact algorithm.

Solution

Solution of Exercise 26.12.

(a) The pressure difference is 5604−5500=104 dbar5604 - 5500 = 104\,\mathrm{dbar}: the parcel cools by 0.129×0.104=0.013 ∘C0.129 \times 0.104 = 0.013\,{}^{\circ}\mathrm{C}, to 0.718 ∘C0.718\,{}^{\circ}\mathrm{C}. (b) The water around it at 5400 m5400\,\mathrm{m} is at 0.730 ∘C0.730\,{}^{\circ}\mathrm{C}, warmer: the parcel, of nearly the same salinity, is denser than its surroundings and sinks back. The layer is stable, although the measured temperature increases downwards; its potential temperature decreases downwards. (c) θ≈0.7313−0.086×5.604=0.249 ∘C\theta \approx 0.7313 - 0.086 \times 5.604 = 0.249\,{}^{\circ}\mathrm{C}, within a thousandth of a degree of the exact value.

26.7 Problem: A radiosonde and a CTD on the same day

Problem 26.1

Weekend problem — the stratification of the ocean and of the air above it, measured on one day

A research ship in the subtropical South Atlantic lowers a CTD in the morning and launches a radiosonde in the early afternoon. Use Γd=9.76 K/km\Gamma_d = 9.76\,\mathrm{K}/\mathrm{km}, κ=2/7\kappa = 2/7, g=9.81 m s−2g = 9.81\,\mathrm{m}\,\mathrm{s}^{-2}, R/M=287 J kg−1 K−1R/M = 287\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} for dry air, and the linear equation of state of seawater about (35,10 ∘C)(35, 10\,{}^{\circ}\mathrm{C}): ρ0=1026.95 kg m−3\rho_0 = 1026.95\,\mathrm{kg}\,\mathrm{m}^{-3}, α=1.67×10−4 K−1\alpha = 1.67 \times 10^{-4}\,\mathrm{K}^{-1}, β=7.61×10−4\beta = 7.61 \times 10^{-4}. All data are exercise data.

Part I — The CTD cast. The cast gives (temperature in ∘C{}^{\circ}\mathrm{C}):

Depth (m)010255075100150200
TT22.4022.3822.3522.1019.2017.6016.0014.90
SS35.8035.8035.8035.7935.7035.6235.5535.45
Depth (m)3005008001000150025003000
TT13.009.604.803.803.002.602.50
SS35.2234.8234.3234.3634.6034.8934.90
  1. Sketch the profiles of temperature and salinity.
  2. Find the depth of the mixed layer with a threshold of 0.5 ∘C0.5\,{}^{\circ}\mathrm{C}.
  3. Compute the mean temperature gradient between 50 m50\,\mathrm{m} and 300 m300\,\mathrm{m}, in ∘C{}^{\circ}\mathrm{C} per 100 m100\,\mathrm{m}. Where is the thermocline?
  4. Compute σ\sigma at the surface with the linear equation of state. The exact equation gives σθ=24.71\sigma_\theta = 24.71: explain the difference.
  5. Identify the salinity minimum and the salinity maximum. Which water masses do they reveal, and where did these form?
  6. The exact algorithm gives a potential temperature of 2.26 ∘C2.26\,{}^{\circ}\mathrm{C} at 3000 m3000\,\mathrm{m}. What is the mean adiabatic gradient over the column, in ∘C{}^{\circ}\mathrm{C} per 1000 dbar1000\,\mathrm{dbar}, taking 3042 dbar3042\,\mathrm{dbar} at that depth?
  7. The exact σθ\sigma_\theta are 26.3526.35 at 200 m200\,\mathrm{m}, 27.1727.17 at 800 m800\,\mathrm{m} and 27.8727.87 at 3000 m3000\,\mathrm{m}. Is the column stable? Where is the pycnocline?

Part II — Two water masses meet. Take the water at 200 m200\,\mathrm{m} (call it C) and the water at 800 m800\,\mathrm{m} (call it I).

  1. Compute the salinity and potential temperature of a mixture of equal parts of C and I, taking θ≈T\theta \approx T for both.
  2. Compute σ\sigma of C, of I and of the mixture with the linear equation of state. How does the value for the mixture compare with the mean of the two parents?
  3. The exact σθ\sigma_\theta of the mixture is 26.88826.888, that of C 26.34226.342 and of I 27.15927.159. By how much is the mixture denser than the mean of its parents? Name the effect.
  4. Near which depth of the cast does a water of the mixture’s density lie? Compare its temperature and salinity with those of the mixture.
  5. What would a mixture of one part of C and three parts of I be? Is it on the cast?

Part III — The radiosonde. The sounding gives (heights above the sea, temperatures in ∘C{}^{\circ}\mathrm{C}):

LevelABCDEFGHIJ
zz (km)0.00.51.51.95.09.012.514.018.024.0
pp (hPa)10129558478075473171861477830
TT22.017.07.39.5−10.0-10.0−36.0-36.0−59.0-59.0−59.5-59.5−57.5-57.5−52.5-52.5
  1. Compute the lapse rate of each of the nine layers.
  2. Classify each layer by Method 26.11. Where is the tropopause?
  3. Taking the layer H–I as isothermal at its mean temperature, compute its scale height and check that the pressures at H and I are consistent.
  4. The sea surface is at 22.4 ∘C22.4\,{}^{\circ}\mathrm{C} (Part I). What do the temperature at A and the layer A–B say about the heat flux between the sea and the air?
  5. A parcel leaves A at 22.0 ∘C22.0\,{}^{\circ}\mathrm{C} and rises dry. Compute its temperature at B, at C and at D, and compare with the sounding.
  6. At what height does the parcel stop being warmer than the surrounding air? What stops it?
  7. Why does the layer B–C have a lapse rate so close to Γd\Gamma_d?

Part IV — Potential temperature and stability.

  1. Compute the potential temperature at every level.
  2. Which layers have θ\theta decreasing, constant and increasing with height? Compare with question 14.
  3. Compute dθ/dz\mathrm d\theta/\mathrm dz in the layers E–F and H–I.
  4. Compute the buoyancy frequency and the period of oscillation in the layer H–I, with its mean potential temperature.
  5. Compare the stability of the layer H–I with that of the thermocline of Part I, where N2=(g/ρ) dρ/dzN^2 = (g/\rho)\,\mathrm d\rho/\mathrm dz with zz the depth; use the exact σθ\sigma_\theta of 24.7924.79 at 50 m50\,\mathrm{m} and 26.3526.35 at 200 m200\,\mathrm{m}.
  6. What is the potential temperature at the tropopause?
Solution

Solution of Problem 26.1.

  1. The temperature is nearly uniform to 50 m50\,\mathrm{m}, drops fast to 300 m300\,\mathrm{m}, more slowly to 1000 m1000\,\mathrm{m}, and is nearly uniform below. The salinity decreases from 35.80 at the surface to a minimum of 34.32 at 800 m800\,\mathrm{m}, then increases to 34.90 at 3000 m3000\,\mathrm{m}.
  2. 22.40−0.5=21.90 ∘C22.40 - 0.5 = 21.90\,{}^{\circ}\mathrm{C} is reached between 50 m50\,\mathrm{m} (22.10) and 75 m75\,\mathrm{m} (19.20): 50+25×0.20/2.90=52 m50 + 25 \times 0.20/2.90 = 52\,\mathrm{m}.
  3. (22.10−13.00)/250=0.0364 ∘C/m(22.10 - 13.00)/250 = 0.0364\,{}^{\circ}\mathrm{C}/\mathrm{m}, 3.6 ∘C3.6\,{}^{\circ}\mathrm{C} per 100 m100\,\mathrm{m}. The thermocline extends from the base of the mixed layer, near 50 m50\,\mathrm{m}, to about 800 m800\,\mathrm{m}; it is sharpest between 50 m50\,\mathrm{m} and 75 m75\,\mathrm{m} (11.6 ∘C11.6\,{}^{\circ}\mathrm{C} per 100 m100\,\mathrm{m}).
  4. σ=1026.95×[1−1.67×10−4×12.4+7.61×10−4×0.80]−1000=25.45\sigma = 1026.95 \times [1 - 1.67\times10^{-4} \times 12.4 + 7.61\times10^{-4} \times 0.80] - 1000 = 25.45, 0.74 kg m−30.74\,\mathrm{kg}\,\mathrm{m}^{-3} too high. The surface water is far from the reference state, and α\alpha is much larger at 22 ∘C22\,{}^{\circ}\mathrm{C} than at 10 ∘C10\,{}^{\circ}\mathrm{C}: the linear form underestimates the lightening by the warmth.
  5. The salinity minimum at 800 m800\,\mathrm{m}, 34.32 at 4.8 ∘C4.8\,{}^{\circ}\mathrm{C}, is Antarctic Intermediate Water, formed under the rainy westerlies of the Southern Ocean. The maximum at 3000 m3000\,\mathrm{m}, 34.90 at 2.5 ∘C2.5\,{}^{\circ}\mathrm{C}, is North Atlantic Deep Water, formed in the Greenland and Labrador Seas.
  6. (2.50−2.26)/3.042=0.079 ∘C(2.50 - 2.26)/3.042 = 0.079\,{}^{\circ}\mathrm{C} per 1000 dbar1000\,\mathrm{dbar}.
  7. σθ\sigma_\theta increases with depth: the column is stable. The pycnocline is where σθ\sigma_\theta increases fastest, in the upper few hundred metres, with the thermocline: 0.820.82 between 200 m200\,\mathrm{m} and 800 m800\,\mathrm{m}, only 0.700.70 over the 2200 m2200\,\mathrm{m} below.
  8. S=(35.45+34.32)/2=34.885S = (35.45 + 34.32)/2 = 34.885; θ=(14.90+4.80)/2=9.85 ∘C\theta = (14.90 + 4.80)/2 = 9.85\,{}^{\circ}\mathrm{C}.
  9. C: 26.461; I: 27.310; mixture: 26.886, exactly the mean of the two (26.886), since the linear form is affine.
  10. The mean of the parents is (26.342+27.159)/2=26.751(26.342 + 27.159)/2 = 26.751; the mixture, at 26.888, is 0.14 kg m−30.14\,\mathrm{kg}\,\mathrm{m}^{-3} denser. This is cabbeling.
  11. With the linear form, the cast gives 26.61 at 300 m300\,\mathrm{m}, 26.88 at 500 m500\,\mathrm{m} and 27.31 at 800 m800\,\mathrm{m}: the mixture’s density lies near 500 m500\,\mathrm{m}. The water there, (34.82,9.6 ∘C)(34.82, 9.6\,{}^{\circ}\mathrm{C}), is very close to the mixture (34.885,9.85 ∘C)(34.885, 9.85\,{}^{\circ}\mathrm{C}): the cast between the Central Water and the AAIW is itself close to a mixing line between them.
  12. S=0.25×35.45+0.75×34.32=34.60S = 0.25 \times 35.45 + 0.75 \times 34.32 = 34.60, θ=0.25×14.90+0.75×4.80=7.33 ∘C\theta = 0.25 \times 14.90 + 0.75 \times 4.80 = 7.33\,{}^{\circ}\mathrm{C}. Interpolating the cast between 500 m500\,\mathrm{m} and 800 m800\,\mathrm{m} gives (34.57,7.2 ∘C)(34.57, 7.2\,{}^{\circ}\mathrm{C}) near 650 m650\,\mathrm{m}: nearly on it.
  13. In K/km\mathrm{K}/\mathrm{km}: A–B 10.0; B–C 9.7; C–D −5.5-5.5; D–E 6.29; E–F 6.50; F–G 6.57; G–H 0.33; H–I −0.50-0.50; I–J −0.83-0.83.
  14. A–B unstable; B–C neutral; C–D an inversion; D–E, E–F, F–G stable; G–H nearly isothermal, stable; H–I and I–J inversions. The tropopause is at G, 12.5 km12.5\,\mathrm{km} and −59.0 ∘C-59.0\,{}^{\circ}\mathrm{C}, where the temperature stops decreasing.
  15. The mean temperature is 214.65 K214.65\,\mathrm{K}: H=287×214.65/9.81=6.28 kmH = 287 \times 214.65/9.81 = 6.28\,\mathrm{km}, and 147 e−4/6.28=77.7 hPa147\,e^{-4/6.28} = 77.7\,\mathrm{hPa}, consistent with the 78 measured.
  16. The air at A, 22.0 ∘C22.0\,{}^{\circ}\mathrm{C}, is cooler than the sea: heat flows from the sea to the air, which is heated from below. The layer A–B, with a lapse rate above Γd\Gamma_d, is the result.
  17. At B: 22.0−9.76×0.5=17.1 ∘C22.0 - 9.76 \times 0.5 = 17.1\,{}^{\circ}\mathrm{C}, against 17.0; at C: 7.4 ∘C7.4\,{}^{\circ}\mathrm{C}, against 7.3; at D: 3.5 ∘C3.5\,{}^{\circ}\mathrm{C}, against 9.5. The parcel is slightly warmer than the air up to C, and much colder at D.
  18. At the base of the inversion, near 1.5 km1.5\,\mathrm{km}: entering the inversion, the parcel becomes colder, hence denser, than the air around it, which stops its rise.
  19. Parcels released by the unstable surface layer rise through it and mix it until its potential temperature is uniform, which makes its lapse rate Γd\Gamma_d; the inversion above stops the mixing.
  20. In K: A 294.15; B 293.99; C 294.08; D 300.51; E 312.65; F 329.29; G 346.28; H 369.50; I 446.99; J 600.91.
  21. θ\theta decreases from A to B (unstable), is nearly constant from B to C (neutral) and increases above (stable), fastest from C to D and above G: the same classification as in question 14.
  22. E–F: (329.29−312.65)/4=4.2 K/km(329.29 - 312.65)/4 = 4.2\,\mathrm{K}/\mathrm{km}; H–I: (446.99−369.50)/4=19.4 K/km(446.99 - 369.50)/4 = 19.4\,\mathrm{K}/\mathrm{km}, nearly five times more.
  23. With the mean θ=408.2 K\theta = 408.2\,\mathrm{K}: N2=(9.81/408.2)×0.01937=4.66×10−4 s−2N^2 = (9.81/408.2) \times 0.01937 = 4.66 \times 10^{-4}\,\mathrm{s}^{-2}, N=0.0216 s−1N = 0.0216\,\mathrm{s}^{-1}, a period of 291 s291\,\mathrm{s}, about five minutes.
  24. N2=(9.81/1027)×(26.35−24.79)/150=9.9×10−5 s−2N^2 = (9.81/1027) \times (26.35 - 24.79)/150 = 9.9 \times 10^{-5}\,\mathrm{s}^{-2}, N=0.0100 s−1N = 0.0100\,\mathrm{s}^{-1}, a period of 630 s630\,\mathrm{s}. The lower stratosphere is the more stable, with N2N^2 nearly five times larger than in the thermocline.
  25. At G, θ=214.15×(1000/186)2/7=346 K\theta = 214.15 \times (1000/186)^{2/7} = 346\,\mathrm{K}: the potential temperature at the tropopause is about 346 K346\,\mathrm{K}.

এই অধ্যায়ে সংজ্ঞায়িত পরিভাষা-শব্দ

পরিভাষায় সব 957 টি শব্দ দেখো