Foundations of Physical Geography / 自然地理学基础

自然地理学基础

地球是一个复杂的多层次行星系统，其表层形态是内动力地质作用与外动力地质作用长期博弈的结果。作为自然地理学的核心研究对象，地球表层系统涵盖从地核深处涌动的岩浆活动，到大气圈中运行的全球环流，从海洋深处循环的温盐环流，到生物圈中演化的生命形态。理解这些要素之间的相互作用机制，是把握全球地理格局形成规律的关键前提。

自然地理学（Physical Geography）作为地球科学的重要分支，其研究范畴与地质学、气象学、海洋学和生态学存在显著交集，但始终保持着独特的研究视角——关注地球表层各圈层之间的界面过程与耦合关系。现代自然地理学不仅需要描述自然现象的空间分布规律，更需要揭示这些规律背后的动力机制与演化历程。本书正是基于这一学术定位，系统运用板块构造理论、年代测定技术和地形演化模型，为后续各洲的自然地理分析建立统一的技术分析框架。

本章将依次探讨地球的几何属性与坐标参照系统、地球内部的圈层结构与物质运动模式、地表各圈层的耦合作用机制、板块构造理论的技术细节、绝对与相对年代测定方法，以及地形演化的数学模型。这些内容将构成理解后续各章节区域地理特征的坚实理论基础。

---

1.1 地球的几何属性与坐标参照系统

地球形状的科学认知演进

人类对地球形状的认知经历了从直觉猜测到精确测量的漫长历程。古希腊数学家埃拉托斯特尼（Eratosthenes）早在公元前3世纪便通过塞伊尼（Syene）与亚历山大（Alexandria）两地的正午太阳高度角差异，成功估算出地球周长约为39,690公里，与现代测量值误差仅约6%。这一经典实验不仅展示了古代科学的智慧，更开创了定量研究地球几何属性的先河。

现代大地测量学对地球形状给出了更为精确的定义：地球并非完美的球体，而是一个两级略微扁平、赤道略微隆起的旋转椭球体（oblate spheroid），更精确地说是一个大地水准面（geoid）所界定的复杂曲面。大地水准面是指假设全球海洋表面完全静止、不受潮汐和海浪影响时所形成的连续等位面，它真实反映了地球内部质量分布的不均匀性。

地球的几何参数可由以下公式精确描述。对于旋转椭球体模型，其长半轴$a$（赤道半径）和短半轴$b$（极半径）的差异定义了两个关键物理量：

扁率（flattening）：$$f = \frac{a - b}{a}$$

第一偏心率（first eccentricity）：$$e^2 = \frac{a^2 - b^2}{a^2}$$

根据1984年世界大地测量系统（WGS84）的定义，地球的几何参数为：

| 参数 | 数值 | |------|------| | 长半轴 $a$ | 6,378,137.0 m | | 短半轴 $b$ | 6,356,752.314245 m | | 扁率 $f$ | 1/298.257223563 | | 椭球表面积 | 510,065,621.724 km² | | 椭球体积 | 1.08321×10¹² km³ |

扁率$f \approx 1/298$意味着地球极半径比赤道半径短约21公里，这一差异虽然相对量级微小，但在高精度定位和大地测量中不可忽视。值得注意的是，不同国家和地区采用的大地参照系可能使用略有差异的椭球参数，如苏联的克拉索夫斯基椭球体（Krasovsky ellipsoid）和中国的CGCS2000坐标系，它们在定位精度和局部拟合效果上各有侧重。

大地坐标系统的技术框架

确定地球表面任何一点的位置，需要建立科学的坐标参照系统。地理坐标系统（geographic coordinate system）以地球自转轴为基准，通过经度和纬度两个角度量来描述点在球面上的位置。

纬度（latitude） 定义为当地垂线方向与赤道平面的夹角。根据测量方法和应用场景的不同，纬度存在多种定义形式：

• **天文纬度（astronomical latitude）**：通过观测恒星高度确定的当地垂线与赤道平面的夹角，反映地球外部重力场方向
• **大地纬度（geodetic latitude）**：参考椭球体法线与赤道平面的夹角，是地图制图和导航定位的基础
• **地心纬度（geocentric latitude）**：从地心指向该点的连线与赤道平面的夹角

对于WGS84椭球体，大地纬度$\phi$与地心纬度$\phi'$之间的关系可通过以下公式转换：

$$\phi' = \arctan\left[(1 - e^2) \tan\phi\right]$$

经度（longitude） 定义为当地子午面与本初子午面之间的二面角。本初子午线（prime meridian）自1884年国际子午线会议确定通过英国格林尼治天文台以来，一直是全球经度测量的起算基准。经度值范围为0°至180°，向东为正（东经），向西为负或取正值（西经）。

由于地球自转的角速度恒定（ω ≈ 7.292115×10⁻⁵ rad/s），纬度不变的情况下，任意点的线速度$v$可表示为：

$$v = \omega \cdot R \cdot \cos\phi$$

其中$R$为该纬度处的地球曲率半径。这意味着赤道点的线速度最大（约1670 km/h），而两极的线速度为零。这一物理特性对大气环流、洋流分布以及导弹发射的轨道设计都具有重要影响。

地图投影与坐标转换

将地球球面坐标转换为平面地图坐标的过程涉及地图投影（map projection）的选择，这一选择直接影响距离、面积、角度和形状等几何属性的保真度。不存在同时保持所有几何属性不变的完美投影，任何投影都是对真实地球的某种妥协与简化。

墨卡托投影（Mercator projection） 是航海导航中最常用的投影方式，其特点是保持恒向线（rhumb line）为直线，使得航海者可沿固定方位角航行。然而，墨卡托投影的致命缺陷在于高纬度地区的面积畸变急剧增大——在60°纬度处，面积被放大至赤道的4倍；在85°纬度处，放大倍数高达33倍。因此，在使用墨卡托投影地图进行面积对比分析时，必须进行投影校正。

兰伯特等角圆锥投影（Lambert conformal conic projection） 广泛应用于中纬度地区的中等比例尺地图，其特点是角度变形小、局部形状保真度高。美国地质调查局（USGS）和中国国家测绘地理信息局均采用该投影作为1:100万地形图的基准投影。

坐标系统之间的转换是地理信息系统的核心技术之一。三维直角坐标系$(X, Y, Z)$与大地坐标$(\phi, \lambda, h)$之间的转换公式为：

$$X = (N + h)\cos\phi\cos\lambda$$ $$Y = (N + h)\cos\phi\sin\lambda$$ $$Z = \left[N(1 - e^2) + h\right]\sin\phi$$

其中$N = a / \sqrt{1 - e^2\sin^2\phi}$为卯酉圈曲率半径，$h$为大地高程。这种转换是GPS定位、卫星轨道计算和大地测量的数学基础。

地球质量与重力场特征

地球的总质量约为$5.972\times10^{24}$ kg，这一数值通过卡文迪许实验（Cavendish experiment）测定的万有引力常数$G = 6.674\times10^{-11}$ m³/(kg·s²)和地球表面的重力加速度反推得出。地球平均密度为5.51 g/cm³，明显高于地表岩石的平均密度（2.7-3.0 g/cm³），这表明地球内部存在密度随深度显著增加的物质分层。

标准重力公式（Standard Gravity Formula） 描述了地球表面理论重力加速度与纬度的关系：

$$g(\phi) = 9.780327\left(1 + 0.0053024\sin^2\phi - 0.0000058\sin^2 2\phi\right) \text{ m/s}^2$$

该公式表明，赤道处的重力加速度（约9.780 m/s²）比两极处（约9.832 m/s²）小约0.53%。这一差异源于两方面的物理机制：一是地球自转产生的惯性离心力在赤道处最大、极地处为零；二是地球两级扁平导致极地观测点到地心的距离小于赤道观测点。

重力异常（gravity anomaly）是指观测重力值与理论重力值之间的差异，可分为空间重力异常和布格重力异常。负重力异常通常指示下方存在低密度地质体，如沉积盆地、含水层或岩浆房；正重力异常则指示高密度地质体，如基性-超基性岩体或金属矿床。石油勘探中广泛应用重力测量圈定沉积凹陷和岩盐穹隆构造。

---

1.2 地球的圈层结构与内动力地质作用

地球内部圈层的物理分层

地球内部存在显著的物理和化学分层结构，这种分层是地球形成后44亿年演化历史的累积结果。根据地震波传播特征（波速不连续面和波速梯度），地球内部可划分为三大主要圈层：地壳（crust）、地幔（mantle）和地核（core）。

莫霍洛维奇不连续面（Moho discontinuity） 是分隔地壳与地幔的标志性地震波速界面，深度在全球范围内变化显著——大洋地区仅5-10 km，大陆地区则为30-70 km。这一界面的存在表明地壳与地幔具有不同的化学组成：地壳以硅铝质（sialic）岩石为主（SiO₂含量约60-70%），地幔则以镁铁质（mafic）岩石为主（SiO₂含量约45%）。

古登堡不连续面（Gutenberg discontinuity） 位于地下约2900 km深处，是地幔与地核的分界面。在该界面处，地震P波速度从约13.6 km/s骤然降低至约8.1 km/s，S波完全消失（表明外核为液态）。这一波速突变揭示了地球核心液态外核的存在。

地幔进一步可细分为上地幔（upper mantle，深度0-670 km）和下地幔（lower mantle，深度670-2890 km）。上地幔顶部存在一个重要的低速层（low velocity zone），地震波速降低5-10%，该层对应岩石圈之下的软流圈（asthenosphere）。软流圈是部分熔融状态的高温塑性层，被认为是板块构造运动的驱动源——地幔对流和火山活动的热源所在。

地核由外核（outer core，深度2890-5150 km，液态）和内核（inner core，深度5150-6371 km，固态）组成。外核的主要成分为铁镍合金，并可能含有少量轻元素（硫、氧、氢等）；内核则被认为是由纯铁或铁镍合金组成的固态晶体结构。尽管内核温度高达5000-6000°C（高于铁的熔点），但极高压力（约360 GPa）使其保持固态。

地震波速度与地球内部物质组成

地震波是研究地球内部结构最有效的探测手段，其传播速度受介质密度、弹性模量和化学组成的共同控制。纵波（P波）和横波（S波）的速度公式为：

$$V_P = \sqrt{\frac{K + 4\mu/3}{\rho}}$$ $$V_S = \sqrt{\frac{\mu}{\rho}}$$

其中$K$为体积模量，$\mu$为剪切模量，$\rho$为介质密度。S波只能在固态介质中传播，而在液态外核中完全消失；P波在液态中的传播速度显著低于在固态中的传播速度。

地震层析成像（seismic tomography）技术利用全球地震台网记录的走时残差和波形数据，重建地球内部三维速度结构。该技术揭示了若干重要发现：冷而致密的高速异常体俯冲到地幔过渡带甚至下地幔深部（如太平洋俯冲带）；热而低密度的低速异常体从核幔边界上升形成地幔柱（mantle plume）；俯冲板片在地幔过渡带（410-660 km深度）可能发生停滞或穿透。

内动力地质作用的类型与特征

内动力地质作用（endogenic geological processes）是指由地球内部能量驱动的地质作用，主要包括构造运动、岩浆作用、变质作用和地震活动。这些作用共同塑造了地表的基本构造格局，是决定区域地形骨架的控制性因素。

构造运动（tectonism） 可分为水平运动和垂直运动两种基本形式。水平运动导致地壳块体沿断裂带发生大规模位移，形成裂谷、转换断层和走滑断裂等构造形态；垂直运动表现为地壳块体的整体抬升或沉降，形成高原、山脉和沉积盆地等正向或负向地形。值得强调的是，构造运动产生的地形往往经历了长期的均衡调整——山根（mountain root）的存在是艾里均衡（Airy isostasy）理论的直接证据，该理论认为高山地区地壳厚度增加以补偿上部载荷的重量。

岩浆作用（magmatism） 包括岩浆生成、运移和固结的全过程。根据岩浆产出位置，可分为深成侵入作用（plutonism）和火山喷发作用（volcanism）。岩浆的化学组成决定了其黏度和流动性：玄武质岩浆（基性岩浆）SiO₂含量低（45-52%），温度高（1000-1200°C），黏度低，流动性好，形成广阔的熔岩台地；花岗质岩浆（酸性岩浆）SiO₂含量高（65-75%），温度相对较低（700-900°C），黏度高，流动性差，常形成陡峭的穹隆或岩钟。

变质作用（metamorphism） 指岩石在高温高压条件下发生矿物成分和结构构造改造的过程。根据变质条件的不同，可分为接触变质作用（contact metamorphism，发生在岩浆侵入体附近的热力变质）和区域变质作用（regional metamorphism，与大规模构造运动相关的温度压力联合变质）。变质岩中保存的矿物组合和结构特征是重建古构造环境的重要依据。

地震活动（seismicity）是地壳弹性应变能突然释放的表现形式。全球约90%的地震能量释放集中在环太平洋地震带（Pacific Ring of Fire）、地中海-喜马拉雅地震带（Alpine-Himalayan belt）和大洋中脊地震带（mid-ocean ridge belt）三大地震带，这些地震带与板块边界高度吻合，印证了地震活动与板块构造运动的内在联系。

地热梯度与地球内部热力学

地球内部的热状态是驱动内动力地质作用的根本能源。观测数据表明，地壳浅部的平均地热梯度（geothermal gradient）约为25-30°C/km，但在不同构造单元差异显著：年轻造山带和火山活动区的地热梯度可达50-100°C/km，而古老稳定地盾区则低至10-15°C/km。

地热梯度$\frac{dT}{dz}$与热流密度$q$的关系由傅里叶定律描述：

$$q = -k \frac{dT}{dz}$$

其中$k$为岩石热导率（typical value ~2.5 W/m·K for continental crust）。全球平均热流密度约为65 mW/m²，但不同构造单元差异悬殊：年轻海洋地壳热流高达100-300 mW/m²，而前寒武纪地盾区热流仅约30-40 mW/m²。

地球内部热源主要来自三个方面：

1. 放射性衰变热（radiogenic heat）：铀（²³⁸U、²³⁵U）、钍（²³²Th）和钾（⁴⁰K）等长寿命放射性同位素在衰变过程中释放热量。据估计，当前地球内部约50%的热流来自放射性衰变热。

2. 原始残留热（primordial heat）：地球形成初期星子吸积和小天体撞击产生的热量仍有部分残留在深部地幔和地核中。

3. 核幔分异热（core-mantle differentiation heat）：地核与地幔边界处可能存在的轻元素化学分异释放潜在热量。

---

1.3 岩石圈-大气圈-水圈-生物圈耦合系统

地球表层系统的圈层构成

地球表层系统（Earth's surface system）由岩石圈（lithosphere）、大气圈（atmosphere）、水圈（hydrosphere）和生物圈（biosphere）四大圈层组成，各圈层在空间上相互叠置、在物质和能量上高度交换，形成了复杂而精密的耦合关系。

岩石圈指地球固体外壳，包括地壳和上地幔顶部刚性层，厚度约100-200 km。岩石圈分为大洋岩石圈（oceanic lithosphere，平均厚度5-10 km）和大陆岩石圈（continental lithosphere，厚度50-200 km）。大洋岩石圈密度较大（~3.3 g/cm³），年龄相对年轻（0-200 Ma），参与板块循环；大陆岩石圈密度较小（~2.8 g/cm³），年龄老（可达4.0 Ga），具有浮力优势，理论上可长期保存。

大气圈是包围地球的气态外壳，质量约5.15×10¹⁸ kg，平均海平面气压101.325 kPa。按化学组成，大气分为干空气（氮气78.08%、氧气20.95%、氩气0.93%、二氧化碳0.04%）和可变组分（水汽、臭氧、悬浮颗粒物）。按垂直结构，大气分为对流层（troposphere，0-12 km，涵盖大气质量80%）、平流层（stratosphere，12-50 km）、中间层（mesosphere，50-85 km）和热层（thermosphere，85-600 km）。

水圈包括海洋、湖泊、河流、冰川、地下水、大气水汽和生物水，总质量约1.4×10²¹ kg。其中97.5%为海洋咸水，2.5%为淡水，而淡水中约68.7%封存于冰川和冰盖中，仅约0.3%为地表可利用的淡水资源。这种水资源分布格局对人类生存和社会经济活动构成根本性约束。

生物圈是地球系统中生命存在和活动的空间范围，涵盖从深海热泉口（~110 MPa，400°C）到高山冻原带（~6000 m）的广阔环境。生物圈通过光合作用、呼吸作用和分解作用参与全球碳循环，通过固氮作用和养分循环参与全球氮、磷循环，对大气组成和岩石风化过程产生深刻影响。

圈层界面的物质能量交换

各圈层之间的界面是物质迁移和能量转换的关键地带。这些界面过程（interface processes）决定了地球表层系统的演化方向和速率。

岩石圈-大气圈界面的相互作用主要表现为物理风化（physical weathering）和化学风化（chemical weathering）两大过程。物理风化包括温度变化引起的热胀冷缩、岩石孔隙中水分冻结膨胀、冰劈作用以及生物根劈作用等，导致岩石崩解为碎屑物质；化学风化则涉及矿物与水、氧气、二氧化碳和有机酸的化学反应，如硅酸盐矿物的水解、长石的高岭石化等。化学风化过程中消耗大气CO₂，具有全球气候调节功能。

据估计，全球河流每年向海洋输送约15-20亿吨溶解物质和100-150亿吨悬浮物质，这些物质主要源自岩石风化过程。化学剥蚀速率（chemical erosion rate）与区域气候、构造抬升速率和岩性密切相关：在季风气候区的碳酸盐岩分布区，化学剥蚀速率可达100-200 t/km²/yr；而在干旱寒冷的高山地区，化学剥蚀速率可能低于10 t/km²/yr。

岩石圈-水圈界面涉及沉积物的