08美国数学建模比赛A题.docVIP

不可忽视的事实：一个海平面上升预测模型 杜克大学：Jason Chen，Joonhahn Cho，Brian Choi 目 录 目录1 问题介绍2 II．模型建立4 海平面上升模型4 温度数据分析5 冰原模型5 物质平衡积累模型6 物质平衡消融模型7 物质平衡与海平面上升模型9 热膨胀模型9 局域化说明9 III. 结果11 海平面上升模型输出数据11 淹没模型仿真结果12 IV. 讨论与结论17 V．建议18 参考文献20 附录A 海平面上升仿真代码21 附录B 地质图点矩阵脚本TOPEX/Poseidon卫星测高仪观测出，在1993-1998年期间，海平面每年上涨3.2 ± 0.2毫米。泰特斯等人估计海平面上升可能会导致270475亿美元的赔偿因温度变化热膨胀海平面上升的重大家普遍相反浮冰没有发挥重大的作用ρwater）的浮冰融化成密度为ρocean 的水体，其体积的增量为ΔV，即： （1） 海水的密度大约是1024.8kg/m3 ；北极海冰的重量大约是2x103 kg。因此，如果北极所有的海冰都融化，体积变化将由下式给出： （2） 近360吨水海平面1毫米（3） 海平面如此微小的变化对我们的模型来说是无关紧要的，因为精度已经低于千分之一毫米。 我们也忽略了南极冰原的影响，因为其总体影响海平面上升是极少，难以量化。卡瓦列利等人预测最低限度模型概述积累/消融纽曼条件高分辨率空间温度数据 图1：模拟流程图 温度数据分析 温度数据在我们的模型中是唯一使用的数据，所以必须仔细分析考虑。因为需要建立不同情景下的模型，我们的温度数据就必须包含这些受控制的且只有一个不同变量的情景。另外，温度数据必须是高质量的，并且能保证为我们的模拟提供正确的临时分析。为了到达这个目的，我们决定使用一个全球气候模型（GCM），通过输入便于我们控制的变量来处理我们自己的温度数据。由于计算能力和时间的限制，我们选择EdGCM模型。EdGCM模型是一个为教学设计的快速模型。该模型是基于美国航天局GISS气候变化模型设计出来的。此程序可以满足我们所有的需要，特别是在模拟的速度方面（大约10小时模拟一个50年的气候阶段），使我们可以分析若干个不同的温度情景。 温度情景分析中纳入了IPCC第三次评估报告得到的关于碳排放量的三个因素在IS92系列中的结果，低、高和中等。IS92e(高)，IS92a(中)以及IS92c(低)情景都是在EdGCM中非常近似的结果。这些近似的碳作用在下图2中表示出来。其他所有的参量根据美国航天局GISS模型，被默认保存。用这种方式，在三个时间序列中，获得了全球地表空气温度。 图2：二氧化碳在EdGCM模型中的影响 对于EdGCM模型，一个不利的影响是它只能输出全球气温变化。虽然能计算区域气温变化，但是却很难到达并且所得结果空间精度较低。然而，根据Chylek等人的研究，格陵兰岛的气温和全球气温之间的关系可以近似的表示为： （4） Chylek等人用未受NAO影响的区域证明，并用气候模型的输出预测了这一结果。 冰原模型 将冰原简化假设为一个矩形，且冰原上层的每一处温度假设为恒定的Ta。这是因为我们的气候模型没有达到格陵兰岛的空间精度，所以微小的温度变化被忽略，不考虑。冰原下层的冻土层有恒定的温度Tl。冰原模型的描述如图3所示。 图3：冰原模型的侧面图 为了计算冰原融化和升华的热通量，我们假设其为无限多的微分单元，如下图4。 图4：冰原微分单元 最初，高度h的计算是利用Williams等人提供的数据。 我们的模型中主要的计算方法是依据物质平衡，物质平衡即计算减去消融量后所积累的量。 积累，是指冰原上增加冰块，其主要形式是降雪。冰体的消融主要是两个过程的结果，即升华和融化。 物质平衡——积累 首先，我们建立积累的模型。Huybrechts等人证明格陵兰岛的温度还没有高达可以融化大面积的冰雪。此外，以Knight的经验推断，冰原的积累率和时间存在非常近似的线性关系，覆盖格陵兰大陆的冰雪积累量是0.3m/year。因此，积累率是0.025m/month。根据物质平衡得出： （5） LD的乘积指冰原的表面积。 物质平衡——消融 其次，我们对消融的升华和融化这两个过程进行建模。升华率（质量变化量）由下式给出：