扫一扫
关注中图网
官方微博
本类五星书更多>
-
>
宇宙、量子和人类心灵
-
>
(精)BBC地球故事系列-星际旅行
-
>
从一到无穷大
-
>
图说相对论(32开平装)
-
>
一本有趣又有料的化学书
-
>
刘薰宇的数学三书:原来数学可以这样学全3册
-
>
光学零件制造工艺学
气象卫星资料:变分同化和质量保证 版权信息
- ISBN：9787030749123
- 条形码：9787030749123 ; 978-7-03-074912-3
- 装帧：一般胶版纸
- 册数：暂无
- 重量：暂无
- 所属分类：>
气象卫星资料:变分同化和质量保证 本书特色
必要的数学基础知识配合详尽的实际案例和推导，给读者描绘了大气资料同化的全貌
气象卫星资料:变分同化和质量保证 内容简介
本书简述如何从物理观测资料中提取大气变量有用信息的大气资料同化方法。内容包括资料同化在数值天气预报中的作用、大气控制方程组、函数拟合插值方法、逐步订正方法、**插值方法、滤波器、极小化算法、伴随技术、三维和四维变分资料同化方法、极轨和静止环境气象卫星资料、全球定位系统无线电掩星资料和多种卫星观测资料中的热带气旋结构。书中还包含了微积分、线性代数、概率论和泛函分析中的一些基础数学理论，论述深入浅出，数学推导步骤完整，图示优美简洁，方便读者理解相关内容。
气象卫星资料:变分同化和质量保证 目录
目录
前言
第1章 资料同化概述 1
1.1 资料同化是干什么的? 1
1.2 热力学变量和大气状态方程 6
1.3 数值天气预报中的大气控制方程 12
第2章 线性代数和统计方法 15
2.1 引言 15
2.2 内积和线性映射的伴随映射 15
2.3 *小二乘拟合 16
2.4 高斯分布 17
2.5 *大似然估计 19
2.6 拟合误差 20
2.7 线性回归 21
2.8 拉格朗日乘子 23
2.9 *小方差估计 24
2.10 结束语 26
第3章 插值 27
3.1 引言 27
3.2 多项式函数拟合 27
3.3 局部多项式函数拟合 30
3.4 区域多项式函数拟合 37
3.5 引入背景场 44
3.6 逐步订正 46
3.7 *优插值 48
3.8 位势场的背景误差协方差 51
3.9 结束语 54
第4章 滤波 55
4.1 引言 55
4.2 傅里叶变换 55
4.3 不同尺度的方差贡献 56
4.4 混淆 57
4.5 非递归滤波器 58
4.6 滤波器设计 61
4.7 递归滤波器 65
4.8 经验集合模态分解 66
4.9 结束语 69
第5章 极小化 70
5.1 引言 70
5.2 泛函极值 70
5.3 标量函数、梯度和Hesse矩阵 73
5.4 极小化迭代、线搜索与收敛速度 75
5.5 *速下降法和牛顿法 77
5.6 共轭梯度法 77
5.7 搜索方向的秩一更新公式 79
5.8 搜索方向的秩二更新公式 80
5.9 L-BFGS方法 81
5.10 线搜索方法 82
5.11 结束语 84
第6章 伴随模式 85
6.1 引言 85
6.2 伴随常微分方程 85
6.3 非线性模式 86
6.4 切线模式 87
6.5 伴随模式 89
6.6 伴随变量与拉格朗日乘子之间的等价关系 91
6.7 解析伴随模式方程 92
6.8 数值伴随模式的计算程序编写 99
6.9 伴随模式在敏感性和相对敏感性研究中的应用 104
6.10 结束语 107
第7章 微波温度计探测资料 109
7.1 引言 109
7.2 携带微波温度计的极轨环境卫星计划 110
7.3 极轨卫星轨道特点 112
7.4 辐射吸收和发射 117
7.5 两点定标方程和观测资料误差 119
7.6 辐射吸收和发射权重函数 122
7.7 快速辐射传输模式 124
7.8 通道特点、扫描模式、视场、权重函数 129
7.9 ATMS条带噪声分析与减噪 138
7.10 风云三号卫星微波温度计 147
7.11 AMSU-A亮温的相对敏感性 149
7.12 结束语 154
第8章 三维变分资料同化 155
8.1 引言 155
8.2 确定性数学公式 155
8.3 统计数学公式 158
8.4 构造B矩阵的NMC方法 160
8.5 构造B矩阵的递归滤波器 165
8.6 3D-Var和卡尔曼滤波比较 168
8.7 AMSU-A液态水路径反演和云检测 169
8.8 偏差估计和偏差订正 172
8.9 AMSU-A资料同化对定量降水预报的影响 183
8.10 附加说明 188
第9章 四维变分同化 189
9.1 引言 189
9.2 4D-Var公式和利用伴随模式的梯度计算 189
9.3 控制重力波振荡的惩罚方法 193
9.4 含有“开-关”过程的物理参数化方案的伴随算子 194
9.5 完整物理过程4D-Var全球同化系统研发 195
9.6 完整物理过程区域伴随模式系统研发 196
9.7 参数估计 198
9.8 增量4D-Var资料同化方法 200
9.9 4D-Var与扩展卡尔曼滤波和集合卡尔曼滤波方法的比较 202
9.10 附加说明 203
第10章 GPS无线电掩星技术 204
10.1 引言 204
10.2 GPS和LEO卫星轨道特征以及GPS掩星计划 205
10.3 超相位延迟和超多普勒频移 208
10.4 弯角和影响参数 211
10.5 折射率反演 213
10.6 折射率和弯角的局地观测算子 214
10.7 弯角的非局地射线跟踪观测算子 218
10.8 超相位延迟的非局地切线连接观测算子 221
10.9 多路径现象的发生与检测 222
10.10 观测误差源 226
10.11 液态云和冰云的影响 228
10.12 温度、气压和水汽反演 232
10.13 卫星微波温度计探测资料的发射后校准 235
10.14 结束语 242
第11章 GOES成像仪 244
11.1 引言 244
11.2 GOES卫星高度 244
11.3 GOES计划和瞬时几何视场特性 245
11.4 AHI和ABI仪器特征 250
11.5 云检测 253
11.6 ABI偏差特征 261
11.7 GOES成像仪资料同化 267
11.8 GOES和POES遥感资料的同时同化 268
11.9 结束语 277
第12章 卫星资料中的热带气旋及其环境 278
12.1 引言 278
12.2 热带气旋概述 279
12.3 AHI观测到的台风玛莉亚（2018） 282
12.4 AMSU-A和ATMS台风暖核反演 296
12.5 锥形扫描微波辐射计AMSR2和MWRI 305
12.6 微波湿度计MHS、ATMS、MWHS和MWHS2 319
12.7 TOMS和OMPS总柱臭氧卫星资料 329
12.8 红外高光谱大气垂直探测仪AIRS、IASI和CrIS 338
12.9 涡旋初始化 345
12.10 附加说明 352
参考文献 353
索引 381
前言
第1章 资料同化概述 1
1.1 资料同化是干什么的? 1
1.2 热力学变量和大气状态方程 6
1.3 数值天气预报中的大气控制方程 12
第2章 线性代数和统计方法 15
2.1 引言 15
2.2 内积和线性映射的伴随映射 15
2.3 *小二乘拟合 16
2.4 高斯分布 17
2.5 *大似然估计 19
2.6 拟合误差 20
2.7 线性回归 21
2.8 拉格朗日乘子 23
2.9 *小方差估计 24
2.10 结束语 26
第3章 插值 27
3.1 引言 27
3.2 多项式函数拟合 27
3.3 局部多项式函数拟合 30
3.4 区域多项式函数拟合 37
3.5 引入背景场 44
3.6 逐步订正 46
3.7 *优插值 48
3.8 位势场的背景误差协方差 51
3.9 结束语 54
第4章 滤波 55
4.1 引言 55
4.2 傅里叶变换 55
4.3 不同尺度的方差贡献 56
4.4 混淆 57
4.5 非递归滤波器 58
4.6 滤波器设计 61
4.7 递归滤波器 65
4.8 经验集合模态分解 66
4.9 结束语 69
第5章 极小化 70
5.1 引言 70
5.2 泛函极值 70
5.3 标量函数、梯度和Hesse矩阵 73
5.4 极小化迭代、线搜索与收敛速度 75
5.5 *速下降法和牛顿法 77
5.6 共轭梯度法 77
5.7 搜索方向的秩一更新公式 79
5.8 搜索方向的秩二更新公式 80
5.9 L-BFGS方法 81
5.10 线搜索方法 82
5.11 结束语 84
第6章 伴随模式 85
6.1 引言 85
6.2 伴随常微分方程 85
6.3 非线性模式 86
6.4 切线模式 87
6.5 伴随模式 89
6.6 伴随变量与拉格朗日乘子之间的等价关系 91
6.7 解析伴随模式方程 92
6.8 数值伴随模式的计算程序编写 99
6.9 伴随模式在敏感性和相对敏感性研究中的应用 104
6.10 结束语 107
第7章 微波温度计探测资料 109
7.1 引言 109
7.2 携带微波温度计的极轨环境卫星计划 110
7.3 极轨卫星轨道特点 112
7.4 辐射吸收和发射 117
7.5 两点定标方程和观测资料误差 119
7.6 辐射吸收和发射权重函数 122
7.7 快速辐射传输模式 124
7.8 通道特点、扫描模式、视场、权重函数 129
7.9 ATMS条带噪声分析与减噪 138
7.10 风云三号卫星微波温度计 147
7.11 AMSU-A亮温的相对敏感性 149
7.12 结束语 154
第8章 三维变分资料同化 155
8.1 引言 155
8.2 确定性数学公式 155
8.3 统计数学公式 158
8.4 构造B矩阵的NMC方法 160
8.5 构造B矩阵的递归滤波器 165
8.6 3D-Var和卡尔曼滤波比较 168
8.7 AMSU-A液态水路径反演和云检测 169
8.8 偏差估计和偏差订正 172
8.9 AMSU-A资料同化对定量降水预报的影响 183
8.10 附加说明 188
第9章 四维变分同化 189
9.1 引言 189
9.2 4D-Var公式和利用伴随模式的梯度计算 189
9.3 控制重力波振荡的惩罚方法 193
9.4 含有“开-关”过程的物理参数化方案的伴随算子 194
9.5 完整物理过程4D-Var全球同化系统研发 195
9.6 完整物理过程区域伴随模式系统研发 196
9.7 参数估计 198
9.8 增量4D-Var资料同化方法 200
9.9 4D-Var与扩展卡尔曼滤波和集合卡尔曼滤波方法的比较 202
9.10 附加说明 203
第10章 GPS无线电掩星技术 204
10.1 引言 204
10.2 GPS和LEO卫星轨道特征以及GPS掩星计划 205
10.3 超相位延迟和超多普勒频移 208
10.4 弯角和影响参数 211
10.5 折射率反演 213
10.6 折射率和弯角的局地观测算子 214
10.7 弯角的非局地射线跟踪观测算子 218
10.8 超相位延迟的非局地切线连接观测算子 221
10.9 多路径现象的发生与检测 222
10.10 观测误差源 226
10.11 液态云和冰云的影响 228
10.12 温度、气压和水汽反演 232
10.13 卫星微波温度计探测资料的发射后校准 235
10.14 结束语 242
第11章 GOES成像仪 244
11.1 引言 244
11.2 GOES卫星高度 244
11.3 GOES计划和瞬时几何视场特性 245
11.4 AHI和ABI仪器特征 250
11.5 云检测 253
11.6 ABI偏差特征 261
11.7 GOES成像仪资料同化 267
11.8 GOES和POES遥感资料的同时同化 268
11.9 结束语 277
第12章 卫星资料中的热带气旋及其环境 278
12.1 引言 278
12.2 热带气旋概述 279
12.3 AHI观测到的台风玛莉亚（2018） 282
12.4 AMSU-A和ATMS台风暖核反演 296
12.5 锥形扫描微波辐射计AMSR2和MWRI 305
12.6 微波湿度计MHS、ATMS、MWHS和MWHS2 319
12.7 TOMS和OMPS总柱臭氧卫星资料 329
12.8 红外高光谱大气垂直探测仪AIRS、IASI和CrIS 338
12.9 涡旋初始化 345
12.10 附加说明 352
参考文献 353
索引 381
展开全部
气象卫星资料:变分同化和质量保证 节选
第1章资料同化概述 1.1资料同化是干什么的？ 本书读者可能对以下问题感兴趣：什么是资料同化？为什么需要做资料同化？资料同化主要关心什么问题？资料同化的理论基础是什么？怎么做资料同化？资料同化的关键技术是什么？需要储备多少资料特点和同化方法方面的知识才能解决不同资料同化情景下的关键技术？对解决某些问题为什么某个同化方法比其他方法更具优势？资料同化的主要挑战是什么？ 大气资料同化的目的之一是为数值天气预报（NWP）提供初始条件。数值天气预报中的控制方程可以形式上表达为 （1.1） 其中，x0是模式变量向量的初始条件，譬如地球大气在初始时刻（t0）和三维空间的风矢量（M，V，W）、密度（o）、位温（60和比湿（g）；x（t）是描述未来时刻（t，t>to）大气状态的预报向量，它的分量所表达的大气变量与xa相同；F（x（t））包含所有影响时间倾向项的大气动力过程、外源强迫（地形、太阳辐射）、显式和隐式参数化物理过程（积云对流参数化、行星边界层过程、微物理过程、辐射参数化等）。 用欧拉和蛙跳格式来近似方程（1.1）左边的时间导数（Haltiner和Williams，1980），我们得到不同时间积分步长％，上的控制方程表达式： （1.2） 其中，At是数值积分步长。一旦能确定变量初始条件X0，通过式（1.2）中的时间积分步骤，就能得到未来时刻的大气状态。这样的一个求解过程在数学上叫作正问题。因此，数值天气预报（NWP）得到的是正问题答案。 数值天气预报中需要解决的一个逆问题例子如下：给定一组观测资料和一个根据模式变量得到观测量（Fn）的观测算子（Hn） （1.3） 其中，满足式（1.2）。我们需要根据这些观测资料估计初始条件向量X0的值，资料同化要解决的是类似这样的逆问题。资料同化的目的是要根据式（1.2）、式（1.3）和已知资料，得到真实大气在某个时间和给定空间分辨率 上的一个“*优”估计值。这里的“*优”是由*大似然或*小方差这两个统计估计方法来定义的（见第2章）。 大部分逆问题比相应的正问题难解得多。数学中的逆问题理论是关于求解逆问题的理论（Tarantola，1987）。它可以为模式中未知参数提供信息，验证模式的正确性，在几个模式中筛选出一个*佳模式，为决策层设计外场观测试验，还可以从大量数据中获取关键信息。值得强调的是，逆问题理论和逆问题方法是对一个已有模式进行改进并提供更深入的分析，但不能提供一个崭新的模式。 我们用以下三个例子来说明为什么想做和需要做资料同化。图1.1展示了欧洲中期天气预报中心（ECMWF）的全球大尺度再分析资料ERA5（Hoffmann等，2018）给出的2018年9月6日0000UTC时飓风“佛罗伦斯”（Florence）的海平面气压分布。飓风中心的海平面气压比实际观测值高了46hPa。根据飓风Florence在2018年9月6日0000UTC的中心气压、*大风速、*大风速半径和34kt风速半径（用4bkt=255km），利用Fujita经验公式（Fujita，1952），我们可以得到对应飓风Florence的人造涡旋海平面气压分布，如图1.1b所示。利用包含物理过程参数化的四维变分同化系统同化人造热带气旋，可以改进飓风路径和强度预报水平（Zou和Xiao，2000;Xiao等，2000;Park和Zou，2004;Tlan和Zou，2019b）。 第二个例子比较了2000年9月17日0000UTC飓风“戈登”（Gordon）附近的QulkSCAT二维海面风分布（图1.2a）（Rlcciardulli和Wentz，2015）和美国国家环境预报中心（NCEP）FNL再分析资料（图1.2b）（NCEP，2000），ECMWF Interim再分析资料（图1.2c）（Simmons等，2007）、ECMWFERA5再分析资料（图1.2d）（Hoffmann等，2018）。飓风Gordon在这个时间达到一级飓风强度。NCEPFNL再分析资料中的海面风速比QuikSCAT观测的海面风速弱很多。这是因为NCEPFNL再分析资料的水平分辨率（1.0°x1.0°）比卫星QuikSCAT海面风、Interim再分析资料和ERA5再分析资料的水平分辨率（0.25°x0.25°）低。QuikSCAT海面风卫星观测资料在飓风Gordon内呈非对称分布，*强风速分布在飓风中心东北方向，*大风速半径略小于100km，*小风速在*佳观测路径飓风中心（图1.2a）。然而，NCEPFNL再分析资料和ERA5再分析资料中的海面风*小风速中心在*佳观测路径飓风中心的西面（图1.2b，d）。由于分辨率较低，NCEPFNL再分析资料*强风速与*小海平面气压中心的距离约为200km（图1.2b），比观测的*大风速半径大了一倍。虽然与QuikSCAT海面风观测资料的水平分辨率相同，Interim再分析资料的海面风几乎完全没有反映飓风Gordon的反气旋环流特征（图1.2c），而ERA5再分析资料中飓风Gordon的海面风反气旋环流特征明显，接近QuikSCAT海面风观测特征。此外，以上三个全球再分析资料（图1.2b～d）都没有显示QuikSCAT观测到的位于飓风中心东北方向的一个海面风速次大值（图1.2a）。因此，我们有理由尝试把QuikSCAT海面风卫星观测资料应用到飓风涡旋初始化、资料同化和数值天气预报中去。 第三个例子显示了飓风“邦尼”（Bonnie）附近由特殊传感器微波成像仪（SSM/I）给出的中心频率分别为19GHz和89GHz的垂直极化通道的二维亮温观测分布图，如图1.3a、b所示（Raytheon，2000），观测时间是1998年8月24日2230UTC。19GHz低频通道观测亮温显示出飓风Bonnie的一个非对称分布特征，飓风Bonnie中的*高观测亮温区位于飓风中心的东北侧，在275K左右，比飓风环境观测亮温高了约55K。然而，在19GHz低频通道*高观测亮温区，89GHz高频通道观测亮温值低得多，约225K，远比飓风环境区域的观测亮温低。除温度和水汽外，SSM/I观测亮温的这些分布特征与飓风Bonnie中云的分布特征有关。飓风Bonnie云中液态水路径和冰水路径的分布如图1.3c～d所示，液态水路径是1998年8月24日2230UTC由SSM/I反演得到，冰水路径是1998年8月24日1700UTC由热带降水测量任务微波成像仪TMI反演得到。液态云辐射增加了卫星微波成像仪接收到的19GHz频率附近的红外辐射量；由于云冰散射作用，卫星微波成像仪接收到89GHz频率附近的红外辐射量减小。卫星微波成像仪观测到的这些特征促成了飓风卫星资料同化的一系列相关研究（Amerault和Zou，2003，2006;Amerault等，2008，2009）。同化受云影响的卫星高频通道亮温观测资料需要考虑辐射传输模式中的冰粒子分布与中尺度数值预报模式中微物理过程参数化方案的一致性（Amerault和Zou，2003），以及云水和冰水变量的背景场误差协方差矩阵（Amerault和Zou，2006），还要发展中尺度数值预报模式和微物理过程参数化方案的伴随算子（Amerault等，2008）。Amerault等（2009）利用美国海军实验室（ONR）海洋和大气耦合中尺度预报系统COAMPS（Coupled Ocean/Atmosphere Mesoscale Prediction System）及其伴随模式（包含COAMPS微物理过程伴随算子），同化受云冰影响的SSM/I亮温观测资料后，改进了飓风预报水平。 观测资料是进行资料同化的原因。大气资料种类繁多，包括GPS无线电掩星弯角和折射率，搭载在极轨环境卫星上的微波和红外辐射仪辐射亮温，静止卫星成像仪红外通道辐射亮温，地基GPS总降水量资料，卫星总柱臭氧，卫星海面风，地面雷达径向风和反射率，卫星水汽反演风矢量，无线电探空仪观测、地面站观测、飞机观测、机载雷达飓风观测、外场观测、下投式探空仪飓风观测数据，等等。即使是在已有的同化系统中加入这些资料，科学家们仍需要完成以下工作：开发合适的观测算子，定量估计资料和模式的误差偏差和方差，开发物理上合理的偏差订正、云检测、资料稀疏化和质量控制方法，保证极小化程序收敛，进行资料同化影响评估，分析资料同化改进或不改进数值天气预报水平的原因。 在过去的80年里，随着计算机、通信、卫星和雷达遥感观测仪器科学和技术的迅猛发展，资料同化从简单的函数拟合主观方法逐步发展成为越来越依赖于计算机能力且更复杂的客观方法。这些能在计算机上实现的方法包括逐步订正法、*优插值、三维变分、四维变分、增量四维变分、卡尔曼滤波、扩展卡尔曼滤波、集合卡尔曼滤波、混合集合卡尔曼滤波和三维变分、混合集合卡尔曼滤波和四维变分。重要的是，我们不仅要知道每个方法的数学公式，也要清楚每个方法中使用的近似假设的实施细节以及各自的含义。 在实际应用中选择某一个同化方法需要考虑该方法的优缺点，以及与其他同化方法的相似性，确定所选同化方法能服务于特定的应用目的。事实上，无论用什么方法，资料同化所产生的大气状态变量的分析场（/）可以近似写成下面的一般形式： （1.4） 其中，xb是背景场向量，矩阵W被称为后验权重。所以，不同资料同化方法之间的*大差别是怎样获得后验权重，以及式（1.4）隔多长时间计算一次。当然，在实际应用中，只有在向量xb的维数较小的情况下，才直接利用式（1.4）得到分析场/。对大部分数值天气预报问题，向量xb或/的维数远大于106。在这种情况下，大气状态变量的分析场/是通过别的数学计算方法间接得到的，式（1.4）近似成立。 资料同化需注意以下几点：观测值被拟合到观测误差之内；②保证观测算子得到的模拟值与真实资料具有物理一致性；③充分利用已知背景场信息；④若模式和观测资料有偏差，要合理估计这些偏差，并在模拟值和观测资料值中减去偏差；⑤合理估计模式和观测资料的统计误差协方差特征；⑥在同化过程中加入在所关注大气尺度上成立的大气动力和物理限制条件；⑦确保计算误差和数据噪声被抑制；⑧计算速度和计算内存可承受；⑨定量估计并给出资料同化得到的大气状态变量分析场/的误差统计特征。总之，通过把各类观测资料与数值预报模式结合，大气资料同化得到在指定分辨率上的、在统计意义上“*优”的大气状态变量估计值。大气资料同化不只是求解一个单纯的数学逆问题。掌握观测物理量和天气系统及其误差结构等知识，对大气资料同化至关重要。当然，计算条件的约束也很重要。 1.2热力学变量和大气状态方程 大气状态是由三维空间的风速矢量v和数个热力学变量决定的。这一节，我们简要介绍以下一些热力学变量的定义：温度（D、气压&）、密度（p）、比容（a）、比湿（g）、水汽混合比（w）、位温（的、内能（M）、热焓（幻、等容热容（Ca）、等压热容（Cp）。不同形式的大气状态方程和热力学定律可以反映大气中不同热力学变量之间的一些约束关系。 地球大气可以近似为理想大气，大气内分子间的相互作用可以忽略。即使一个气团作为整体在某个时间是静止的，气团内的分子却不停地向四面八方运动着，分子运动的动能之和决定了该气团的温度。理想气体温度定义为（Bohren和Albrecht，1998） 其中，mm是单个分子质量，是第i个分子运动速度，k是玻尔兹曼常数，表示对气团内所有分子求和的平均算子。
书友推荐
- >
企鹅口袋书系列·伟大的思想20:论自然选择(英汉双语)
企鹅口袋书系列·伟大的思想20:论自然选择(英汉双语)
¥6.3¥14.0 - >
【精装绘本】画给孩子的中国神话
【精装绘本】画给孩子的中国神话
¥20.9¥55.0 - >
随园食单
随园食单
¥20.6¥48.0 - >
苦雨斋序跋文-周作人自编集
苦雨斋序跋文-周作人自编集
¥6.9¥16.0 - >
巴金－再思录
巴金－再思录
¥15.6¥46.0 - >
推拿
推拿
¥12.2¥32.0 - >
山海经
山海经
¥20.4¥68.0 - >
新文学天穹两巨星--鲁迅与胡适/红烛学术丛书(红烛学术丛书)
新文学天穹两巨星--鲁迅与胡适/红烛学术丛书(红烛学术丛书)
¥9.9¥23.0
本类畅销
-
-
时间简史:插图版
¥34.3¥70 -
勒维特之星-大发现系列丛书
¥6.1¥16 -
宇宙与原子
¥14.8¥39 -
天空之光如何照亮文明:中国早期天文学与文明若干专题研究
¥15.4¥48 -
气象学与气候学教程
¥29.4¥42
