まずは静的な大気を考え、そこにはたらく力の釣り合いから静水圧平衡の式を求めます。 ある微小領域の大気にはたらく力が釣り合っていれば、大気は動かずに静止しているということになります。

　上図のような円筒状の領域と、そこにかかる力を考えます。 鉛直上向きをz軸の正の方向に取ります。 円筒下部での気圧をPとすると、円筒の下面に加わる力の大きさは圧力に面積をかけたPSとなり、向きは鉛直上向きです。 また、円筒上部での気圧をP+dPとすると、同様に円筒の上面に加わる力の大きさは(P+dP)S、向きは鉛直下向きです。 ここでは特に指定はしていませんが、大気を考えているので上空に行くほど気圧は小さくなる場合はdPは負の値となります。 次に、円筒にはたらく重力は円筒全体の質量をmとするとmgとなりますが、円筒内の気体の密度をρとするとρ・S・dz・gと書くことができます。 向きは鉛直下向きです。 以上より、円筒にはたらく力の釣り合いの式は、

　この式の両辺を密度で割ってやると

静水圧平衡の一般形

　静水圧平衡の式をより一般的にベクトルを用いて書くと

理想気体の状態方程式

　静水圧平衡の式は圧力Pについての微分方程式になっていて、これをzについて積分すれば鉛直方向の圧力分布が計算出来ることになります。 しかし、圧力、密度と温度の関係が不明であるためこのままでは微分方程式を解くことができません。 それぞれを結びつける式としては、気体の状態方程式があります。

　理想気体の状態方程式として恐らく最も知られているのは、高校物理の範囲でも取り扱う次の形だと思われます。

　気体定数Rには、次のような関係があります。

　系に含まれる全粒子数Nを体積Vで割ったものは、数密度 (number density)と呼ばれます。 文字通り、単位体積当たりに何個の粒子が含まれているかを示す値であり、[m^-3]など、体積の逆数の次元(長さのマイナス3乗の次元)を持ちます。 「個数密度」という表現もたまに見られます。 数密度を とすると、状態方程式は

モル質量を用いた表式と比気体定数

　もう一度一般的な表式に戻ります。

　物質量nは、気体の全質量をその気体のモル質量 (molar mass)で割る事で求めることができます。 モル質量とは、物質1 molあたりの質量に相当します。 例えば水分子のモル質量はおよそ18.0×10^-3 kg mol^-1あるいは18.0 g mol^-1であり、水1 molの質量は18.0 gです。 気体の全質量をM_total [kg]、その気体のモル質量をM [kg mol^-1]とすると、物質量n [mol]は で与えられます。 気体が混合物である場合は、平均したモル質量を計算に用います。 モル質量を用いて状態方程式を書き直すと以下のようになります。

比気体定数

　上の状態方程式に対して、気体定数をモル質量で割った値を比気体定数 (specific gas constant)と呼ぶ事があります。 比気体定数は次のように表されます。

　気体定数は気体の種類や混合比に関係なく一定の値を持つ定数ですが、定義より明らかなように、比気体定数は気体の種類や混合比によって異なる値を持つ量です。 気体定数Rは [J K^-1 mol^-1]の次元を持ちますが、気体定数をモル質量で割った値である比気体定数R_speの次元は[J kg^-1 K^-1]です。

　取り扱う気体の種類が変わらない場合は、気体定数よりも比気体定数を用いた方が便利な場合があります。 そのためか、特に工学の分野においては比気体定数のことを指して気体定数と呼び、単にRのみで表すことがあります。 このような場合は、(本来の)気体定数を と表記して区別します。 また、特に区別をする場合は普遍気体定数 (universal gas constant)と表現します。

　次に平均分子量を用いた表式を導出します。 先ほど導出した、統計力学でよく使われる表式は以下の通りです。

　平均分子量とm_uの積は、系に含まれる粒子の平均質量を表します。 よって、に全粒子数Nをかけた値は、系の気体の全質量となります。

水素原子質量を用いた場合

　上記の状態方程式は、統一原子質量単位m_uを用いず、水素原子の質量m_Hを用いて以下のように書かれている場合もあります。

　以下に、状態方程式の様々なパターンをまとめます。 物理量の単位は便宜上SI単位系(あるはMKSA単位系)で書いてあります。

大気のスケールハイト

　まずは、静水圧平衡の式を解きます。 解くべき式と、物理量を関係付ける状態方程式は次の通りです。

　理想気体の状態方程式の両辺をzで微分します。 圧力P、密度ρ、温度Tがそれぞれ高さzの関数だとすると、微分した結果は

　この式を用いて、静水圧平衡の式からPを消去すると、

　この式を見ると、密度は指数関数的に減少し、高度が z=H₀となる場所で、密度が地上の密度の 1/e倍になっていることが分かります。 このような、元の量の1/eになるまでの距離や高度 H₀の事を、スケールハイトと呼び、指数関数的に減少する際の系の典型的な長さにあたります。

　状態方程式の表現に様々な種類があるのと同様に、スケールハイトも様々な書き方が出来ます。 例えば、気体定数と平均モル質量を用いた場合のスケールハイトの表式は

圧力分布の導出

　上では圧力を消去して密度について解きましたが、密度を消去して圧力分布を導出する事も出来ます。 状態方程式を微分せず、そのまま静水圧平衡の式に代入して密度を消去すると、

　式を見ると分かるように、等温大気の場合は密度と圧力の高度依存性は同じ形で表され、密度のスケールハイトと圧力のスケールハイトは同じになります。

　大気が等温の理想気体である場合、密度分布、圧力分布は以下のように表される。

ここで、

　あるいは

は密度と圧力のスケールハイトである。

スケールハイトから分かる事

　大気のスケールハイトは、大気を扱う時のあちこちに顔を出します。 ここではその一部を紹介します。

大気の柱密度

　得られた密度を、z方向に積分します。 積分範囲は、0(地表面)から無限遠までです。

　密度をzで0から無限遠まで積分すると、地表面での密度にスケールハイトをかけた形になることが分かります。 これは、質量を面積で割った次元を持っています。 つまり、密度を[kg m^-3]、スケールハイトが[m]だとすると、これは[kg m^-2]という次元になります。 つまり、単位面積あたりの大気の質量を意味します。 これを大気の柱密度 (column density)と呼びます。 ある単位面積の上にある無限遠まで続く柱の中にある大気の質量を表しているため、このような名前で呼ばれます。

　単位面積あたりの大気質量が、地表面での密度とスケールハイトの積で書けるということは、全ての大気を表面での密度に圧縮した場合、その厚みはスケールハイトと等しくなるという事になります。

大気の重心の位置

　地表から測った時の大気の重心の位置は、次の式で計算出来ます。

　分子の部分を計算すると以下のようになります。

　よって、大気の重心の位置は

大気の全質量

　スケールハイトなどを用いて、大気の全質量を概算する事も出来ます。 大気の柱密度の部分では、単位面積あたりの大気の質量を計算しました。 そのため、柱密度に地球の表面積をかければ、大気の全ての質量になることが分かります。

　惑星の半径をR_pとすると、大気の全質量は で計算出来ることになります。 具体的な数値は後ほど紹介します。

スケールハイトを決める物理量

　スケールハイトの表式は、

　Rは定数であるため、スケールハイトに影響を及ぼすのは残りの3つの物理量になります。
　分子に温度Tが入っているため、これは大気の温度が高いほどスケールハイトは大きくなるということを意味します。 直感的には、温度が高いほど気体は膨張して大気層はより広がり、逆に温度が低いと収縮して大気層も薄くなるということに対応しています。

　分母に平均モル質量Mがありますが、平均モル質量はおおまかに言えば大気の重さを表しています。 Mが大きいほど、つまりその天体の大気が重いほどスケールハイトは小さくなり、大気層は薄くなります。 二酸化炭素主体の重い大気よりも、水素主体の軽い大気の方が大気層が大きく広がりそうですが、この感覚に対応しています。

　スケールハイトの表式 　 　を用いて、様々な天体のスケールハイトを計算してみます。

地球

　地球の平均気温を288 K、重力加速度を9.80665 m s^-2とします。 まずは平均モル質量を計算する必要があります。 そこで、地球大気の組成から計算します。 地球の大気組成は、
・窒素 78.08%
・酸素 20.95%
・アルゴン 0.934%
・二酸化炭素 0.04%
です。 この組成と各分子の分子量から平均モル質量を計算すると、
0.7808×28 + 0.2095×32 + 0.00934×40 + 0.0004×44 = 28.9576 g mol^-1 = 28.9576 × 10^-3 kg mol^-1 となります。

　これらの値と、気体定数 R=8.3144621 J K^-1 mol^-1を用いてスケールハイトを計算すると、8432.24867...m となるため、地球大気のスケールハイトは ~8.4 kmとなります。

　スケールハイトが 8.4 kmということは、標高 8.4 kmの地点で気圧が地上の気圧の1/e倍 (~0.358倍)になるということです。 参考として、エベレストの標高は 8848 mで、その頂上付近での気圧は地上の3分の1程度になっているため、おおむね一致していることが分かります。

　密度や圧力は

　例えば富士山の山頂(標高 3776 m)での気圧を計算すると、地上での大気圧P₀を 1013.25 hPaとして圧力の関係式に代入すると、647.50 hPa、およそ 650 hPaとなります。 これは実際の富士山山頂での平均的な気圧におおむね近い値になります。 また、エベレストの山頂での値を計算するとおよそ 355 hPaとなります。

　また、静水圧平衡の式から導出した圧力に関する微分方程式は、以下のように書けました。

火星

　次に火星大気のスケールハイトを計算します。 火星の平均気温を 210 K、重力加速度を 3.71 m s^-2とします。 火星の大気組成は、
・二酸化炭素 95.32%
・窒素 2.7%
・アルゴン 1.6%
・酸素 0.13%
となっています。 地球の場合と同様に平均モル質量を計算すると、 43.3784 × 10^-3 kg mol^-1 になります。

　スケールハイトを計算すると、11366.0458...mとなるため、火星大気のスケールハイトは ~11.4 kmとなります。

　二酸化炭素主体の火星大気は、窒素主体の地球大気よりも平均モル質量は大きくなります。 これはスケールハイトを小さくする方向にはたらきますが、火星の気温が低い事に加え、重力加速度が地球の 38%程度しかないため、火星大気のスケールハイトは地球のものよりも大きくなります。

金星

　金星の平均気温を 737 K、重力加速度は 8.87 m s^-2とします。 金星の大気組成は、
・二酸化炭素 96.5%
・窒素 3.5%
・二酸化硫黄 0.015%
となっています。 同様に平均モル質量を計算すると、 43.2296 × 10^-3 kg mol^-1 になります。 火星も金星も大気組成は二酸化炭素主体と似ているため、平均モル質量はほとんど同じ値になります。

　スケールハイトを計算すると 15899.8213...mであり、金星大気のスケールハイトは ~15.9 kmとなります。

　金星大気の平均モル質量は火星と同様に地球より大きくなります。 金星は質量・半径ともに地球よりやや小さい程度であるため、表面での重力は地球とは大差ありません。 しかし、表面温度が 737 Kと圧倒的に異なるため、スケールハイトも地球を大きく上回る値となっています。

タイタン

　次に、土星の衛星タイタンの大気での計算をしてみます。 タイタンは非常に濃い大気を持つ衛星として有名です。 タイタンの温度を 94 K、重力加速度を 1.354 m s^-2とします。 タイタンの大気組成は、
・窒素 95%
・メタン 4.9%
となっています。 よって平均モル質量は 27.384 × 10^-3 kg mol^-1となります。 窒素主体の大気であり、地球での値と近くなっています。

　スケールハイトを計算すると 21078.8281...mであり、タイタン大気のスケールハイトは ~21.1 kmとなります。

　静水圧平衡から導出した密度分布をz方向に積分する事で、柱密度が計算出来ました。 また、柱密度に天体の表面積をかけることで大気の全質量を概算することができます。 柱密度Σと大気の全質量M_atmは、それぞれ

　スケールハイトは既に計算しましたが、柱密度を計算するためには表面での大気の密度が必要です。 そこで、理想気体の状態方程式より、

地球

　地球の場合、地表での大気密度と、大気の柱密度は
ρ₀ = 1.23 kg m^-3
Σ = 1.03 × 10⁴ kg m^-2
となります。 ここで計算した柱密度の値から、1平方メートルの地面の上に乗っている大気の質量はおおよそ 10000 kg、すなわち 1平方メートルあたり10トンの大気が乗っているということが分かります。

　これに地球の表面積をかけると、地球大気の全質量は

M_atm = 5.27 × 10¹⁸ kg

と求まります。

　NASAのEarth Fact Sheetに掲載されている地球大気の全質量は5.1 × 10¹⁸ kgとなっており、おおむね一致していると言えます。 地球の質量は 5.9726 × 10²⁴ kgなので、地球大気の質量は地球質量のわずか 0.000088%ということになります。 大気層がいかに薄くわずかなものであるかが分かります。

火星

　火星の表面での大気圧は 6.36 hPaです。 地球の大気圧が 1013.25 hPaであることを考えると、非常に希薄な大気です。 ここから、火星の地表での大気密度と、大気の柱密度は
ρ₀ = 1.58 × 10^-2 kg m^-3
Σ = 1.80 × 10² kg m^-2
と求められます。 火星の場合は、1平方メートルあたりの大気質量は 180 kgであり、10トンあった地球大気よりも遥かに小さい値になっています。

　これに火星の表面積をかけると、火星大気の全質量は

M_atm = 2.59 × 10¹⁶ kg

と求まります。 火星大気の質量は地球大気の質量の 0.49%、おおよそ200分の1程度の値になることが分かります。

　NASAのMars Fact Sheetに掲載されている火星大気の全質量は2.5 × 10¹⁶ kgとなっており、こちらもおおむね一致していると言えます。 また、火星の質量は 0.64174 × 10²⁴ kgであるため、火星大気の質量は火星質量の 0.000004%ということになります。

金星

　金星の表面での大気圧は 92000 hPaです。 地球の大気圧のおよそ 90倍の値であり、非常に高圧の環境になっています。 ここから、金星の地表での大気密度と、大気の柱密度は
ρ₀ = 6.52 × 10 kg m^-3
Σ = 1.04 × 10⁶ kg m^-2
と求められます。 金星の場合は、1平方メートルあたりの大気質量はおよそ 1000000 kg、つまり1000トンになります。 地球の場合は10トンだったため、およそ 100倍となっています。

　柱密度に金星の表面積をかけると、金星大気の全質量は

M_atm = 4.77 × 10²⁰ kg

と求まります。 金星大気の質量は地球大気の質量の90.5倍となり、金星は地球の90倍の大気を持っていると言うことができます。

　NASAのVenus Fact Sheetに掲載されている金星大気の全質量は 4.8 × 10²⁰ kgとなっており、こちらもおおむね一致していると言えます。 また、金星の質量は 4.8676 × 10²⁴ kgであるため、金星大気の質量は金星質量の 0.0098%ということになります。 地球や火星に比べて大気の量が遥かに多いため、大気質量の割合も他の天体に比べてずっと大きな値になっています。

タイタン

　タイタンの表面での大気圧はおよそ 1600 hPaです。 地球の大気圧よりも大きい値になっています。 ここから、タイタンの地表での大気密度と、大気の柱密度は
ρ₀ = 5.26 kg m^-3
Σ = 1.11 × 10⁵ kg m^-2
と求められます。 タイタンの場合は、1平方メートルあたりの大気質量は 111000 kg、おおむね100トンになるため、単位面積あたりの大気の量は地球の10倍だと分かります。

　タイタンの表面積をかけると、タイタン大気の全質量は

M_atm = 9.23 × 10¹⁸ kg

と求まります。 タイタンは地球よりずっと小さい衛星ですが、大気の全質量は地球の1.75倍となり、地球より多くの大気を持っていることが分かります。

　上のスケールハイト導出のときは、大気は等温であるという仮定をおきました。 そのため、状態方程式を微分した際に出てくる温度のz微分の項は0になりました。 しかし実際には、大気の温度は高度に依存して変化します。 また、そのほかの物理量も常に一定である保証はありません。 そのため、一般的なスケールハイトは複雑な形になります。

　ここまでは、高度によって変化するのは密度と圧力だけだとして計算して来ました。 しかしより一般的には、理想気体の状態方程式は

　静水圧平衡の式は、より一般的には

　平均モル質量(または平均分子量)と重力加速度の変化は無視出来るとした場合、スケールハイトは温度依存性を持つ値である事が分かります。 等温大気の仮定では、この温度依存性を無視して温度を定数として扱っているため、スケールハイトは高度依存性の無い定数として書くことができたということになります。

　大気の温度、大気の組成などが高度依存性を持つ一般的な場合、大気のスケールハイトは

　あるいは

で与えられ、高度によって変化する関数となる。

　静水圧平衡の式と、理想気体の状態方程式は同じものを用います。 ただし、温度の高度依存性を残したものを使います。

　ここで温度の高度依存性を、次のような形で与えます。

　理想気体の状態方程式と温度の高度依存性の式を、静水圧平衡の式に代入することで以下の式が得られます。

　ここで、

　ここで境界条件から積分定数を決定する事を考えます。 上で行ったスケールハイトの導出の時と同様に、地表面、すなわちz=0で圧力はP₀であるとすると、

　一方、理想気体の状態方程式に圧力構造の式を代入する事によって、密度構造を以下のように求めることができます。

　以上より、次の結果が得られます。

　静水圧平衡にある大気において、温度が

で与えられる高度依存性を持つとき、圧力と密度構造は以下の式で与えられる。

ただし、 とする。

等温大気との比較

高度依存性

　等温大気の場合は、圧力構造は

　しかし等温でない場合は、形は似ているものの指数が1だけ違っています。 そのため、非等温大気の場合は圧力と密度の高度依存性は異なることが分かります。

　また、等温大気の場合は指数関数的に減衰していますが、この非等温大気の場合は指数関数ではなくべき関数的になっています。

大気上端の存在

　等温大気の場合は指数関数的な減衰になっているため、圧力や大気は無限遠の高度で0になる構造になっています。 しかし非等温大気の場合は、式の形より の時に圧力・密度ともに0になります。

　従って、非等温大気では有限の高度で圧力・密度が0になるということが分かります。 今考えている問題設定の下では、地表での温度と気温減率から大気の上端が決まります。

　地表面の温度が高ければ大気の上端は高くなり、気温減率が大きく温度の減少が大きい場合は大気の上端は低くなります。

スケールハイト

　一般的な条件でのスケールハイトは、

　このスケールハイトの式に、z=0を代入すると

上空での圧力の具体値

　非等温大気での圧力構造の式

　地表での大気圧をP₀ = 1013.25 hPa、T₀ = 288 Kとし、平均モル質量や重力加速度もこれまで用いたものと同じ値を代入します。 このとき、富士山山頂のz = 3776 mでの気圧は 640 hPa、エベレスト山頂のz = 8848 mでの気圧は 324 hPaとなります(等温大気の場合はそれぞれ 648 hPa、355 hPa)。 大きな差ではありませんが、温度の依存性を考えることによって値が変わっていることが分かります。 特にエベレスト頂上での値は、一般的に言われている「地上の3分の1」という値により近くなっています。

　なお、この仮定のもとで、圧力が0になる高度 を計算すると、 288 [K] / 0.006 [K m^-1] = 48 [km] となります。 つまり高度 48 kmの地点で大気が途切れるということになります。 もちろん実際はそのような事にはなっていません。 地球大気はいくつかの層に分かれており、高度が上がるに従って気温が下がっていく対流圏は、高度 11 km程度まで続いています。 そこより上は成層圏ですが、そこでは気温は高度につれて低下せず、むしろ上昇する傾向にあります。 そのため、温度が高度に従って一定の割合で低下して行くという仮定は、高度 11 km程度より上空では間違っているということになります。

大気の柱密度と大気の全質量

　密度構造の式をz方向に積分すると大気の柱密度が計算でき、そこから大気の全質量が計算出来ることは既に紹介しました。 先程で計算した非等温大気の場合も大気の密度構造が導出出来たため、これをz方向に積分する事で同じ議論が出来るはずです。 ここまでの議論から、大気は で途切れることが分かっています。 よって積分の上端は無限遠ではなく までという事になり、柱密度は以下の積分から計算出来ることになります。

　 を代入すれば、柱密度は

　大気が等温である事を仮定すると、密度と圧力は指数関数的な減少をし、大気が非等温であり、一定の割合で温度が下がって行く大気を仮定すると、密度と圧力はべき関数的に下がって行くことが分かりました。 また、適切と思われる気温減率を仮定して実際の圧力を計算した結果、小さなずれはあるもののおおむね近い値になる事も分かりました。

　両者の違いは温度構造のみであり、等温大気は温度は地上の温度で常に一定、非等温大気の場合はここでは という構造を仮定しました。 この場合、気温減率が0のときが等温大気であるということになります。 そのため気温減率が小さい極限(Γ→0)では、非等温大気の構造は等温大気の構造に漸近するはずです。 ここではそれを確かめてみます。

　大気の温度が で与えられる場合の圧力構造は

　ここで、

　この極限の形は、高校数学で扱った通り、以下のように計算出来ます。

　この結果を用いると、

　よって、気温減率が0となる極限では、非等温大気の構造は等温大気の構造と一致することが証明出来ました。

　等温の大気構造と、非等温の大気構造をグラフにして比較をしてみます。

　横軸が地表面からの高度 [km]、縦軸が気圧 [hPa]です。 地表面での圧力はP₀ = 1013.25 hPaとし、地上での温度T₀ = 288 K、気温減率Γ = 0.006 K m^-1としています。 平均モル質量や重力加速度は地球大気の計算の際に上で用いたものと同じ値を使っています。 赤で示したのが等温大気で、指数関数的な減衰をするものです。 緑で示したものが非等温大気で、べき関数的な減衰をするものです。

　地表付近では両者の違いはあまり大きくありません。 富士山山頂とエベレスト山頂での違いを比較しましたが、どちらも高度 10 km未満であるためこのグラフ上では違いがほとんど見えない領域です。 しかし高度が上がるにつれて違いが顕著になって行きます。

　また、等温大気の場合は圧力は0に漸近しますが、有限の高度では0になりません。 それに対して、非等温大気の場合は途中で0になっています。 先程の議論から、表面での温度が 288 K、気温減率が 0.006 K m^-1の場合は高度 48 kmで圧力が0になるのでした。 そのあたりでちゃんと圧力が0になっている事が確認出来ます。

　上のグラフでは減少が急激で見づらいため、縦軸の対数を取ったグラフを表示します。

　上のものと同じデータですが、縦軸が対数目盛になっています。 そのため、指数関数である等温大気(赤)は直線になっています。 また、非等温大気(緑)が 48 km の所で0になっている(縦軸が対数のため、グラフ上では無限小に発散する)ことが分かります。
　地上付近では等温大気と非等温大気のずれは非常に小さいですが、高度が高くなると一気にずれが大きくなります。 逆に言うと、地上付近では等温大気の仮定はおおむね正しいと言えるということでもあります。

　気温減率が0の極限では、非等温大気と等温大気は一致することを示しました。 そこで、気温減率の値を変えた場合の大気の圧力構造を比較してみます。

　赤が等温大気です。 緑で示したのが先程と同じで気温減率がΓ = 0.006 K m^-1の場合、青がΓ = 0.003 K m^-1、マゼンタがΓ = 0.001 K m^-1の場合です。 気温減率が小さい場合は温度が0になる高度が高いため、大気が途切れる位置も高くなっていることが分かります。 また、気温減率を小さくするほど、赤の等温大気の線に漸近して行く事も分かります。

参考文献