【O奖解析】2022HIMCM-B题:二氧化碳与全球变暖(P1)

f29de6b2706880c7e9d6af0a768a8421.jpeg


一、问题描述

1647846e2d3bb5634640bb02e794ccc2.jpeg 1、问题背景

在工业革命之前,大气中的二氧化碳(CO2)浓度稳定在约280ppm(百万分率)。至2004年3月,大气中二氧化碳浓度达到了377.7ppm,标志着当时10年间最大的增长幅度。根据美国国家海洋和大气管理局(NOAA)以及斯克里普斯海洋研究所(SIO)的科学家们的报告,2022年5月,月平均二氧化碳浓度峰值达到了421ppm。经济合作与发展组织(OECD)的一份报告显示,预计到2050年,二氧化碳水平将达到685ppm。


《今日科学》杂志编辑部要求你们团队分析并评估当前报道及未来预测的二氧化碳浓度水平的准确性。他们提供了两组数据集(CO2数据集1和温度数据集2),以协助你们的研究。


2、问题提出

要求如下:


1、您是否同意有关二氧化碳浓度水平的说法?使用CO2数据集1分析二氧化碳的变化。

a. 您是否认同2004年3月的二氧化碳增加量比以往任何10年间的增长都要大?请说明理由。

b. 根据数据拟合多种数学模型,描述过去并预测未来大气中二氧化碳的浓度水平。

c. 利用每种模型预测2100年时大气中二氧化碳的浓度。是否有任一模型预测结果与声称的2050年达到685ppm的二氧化碳浓度水平相符?若到2050年还未达到该水平,您的模型预测何时会达到685ppm?

d. 您认为哪种模型最为精确?为什么?


2、温度与二氧化碳之间有何关系?许多科学家认为全球气温升高与大气中二氧化碳浓度之间存在关联。使用第一部分的工作成果和温度数据集2来帮助比较陆地海洋温度与二氧化碳浓度水平之间的联系。

a. 构建一个模型预测未来陆地海洋温度变化。相对于1951年至1980年的基准期,您的模型预测何时陆地海洋平均温度将分别上升1.25°C、1.50°C和2°C?

b. 构建一个模型分析自1959年以来二氧化碳浓度与陆地海洋温度之间的关系,并解释这种关系或证明二者间不存在关系的理由。

c. 将第二部分b中构建的模型扩展至未来。您的模型在多远的未来内是可靠的?对于模型预测未来二氧化碳浓度水平和/或陆地海洋温度的能力,您有哪些关注点?


3、为《今日科学》撰写一篇非技术性文章(最多1页),阐述团队的研究发现和对未来可能的建议。


3、研究思路

由以上部分可以得出,问题一主要分为以下两个小问:

1、2004年3月的二氧化碳增加量是否显著;

2、建立多个模型预测未来大气中二氧化碳的浓度水平并比较模型准确度;


O奖获得者团队的研究思路如下:

a.在2004年CO2浓度的增长情况方面,采用10年平均增长量、10年线性拟合斜率和T-test统计分析三种手段,证实2004年是否呈现出较高的增长态势。


b.在预测CO2浓度方面,团队将数据集划分成训练集和测试集,分别构建了时间序列下的线性模型、二次模型、指数模型、S型模型以及ARIMA(2,2,3)模型,并对各模型在整个数据集上的拟合优度(R-squared)进行了比较。并分别预测2100年的CO2浓度。



2004年CO2增长分析

2004年CO2增长分析旨在探究这一年是否呈现异常或显著的增加趋势。为了实现这一目标,研究团队采用了三种不同的定量方法来衡量CO2浓度增长的变化:


1、十年平均增长量计算

分析人员首先计算了从1968年至2021年间的CO2浓度(y)平均增长量。


1ce8ecd6d62a45bc423ec98d64226ab1.jpeg

结果显示,截至2004年,十年平均增量为1.86℃,相较于前一年2003年的10年平均增量1.89℃有所下降。这表明2004年与之前的十年相比,并未展现出最高的平均增长速度,但仍然处于较高水平。


f6849078193af58822e2b4cf1acddd8a.jpeg


2、十年线性拟合斜率

使用最小二乘法对每个连续十年的数据段绘制线性拟合曲线,得到相应的斜率。


8243a88046d731ebed2f8e035f31502a.jpeg

通过对2004年前后的数据进行比对,可以了解增长趋势的变化情况。2004年对应的10年线性拟合斜率为1.853,比以往任何10年都要大,此后斜率仍呈上升趋势。


03cd93881ca019c7ce67c7640fd24f31.jpeg


3、十年增长序列的T检验与秩和检验

T-Test(独立样本双尾检验),是一种用于检验两组数据平均值是否有显著差异的统计方法,适用于数据满足正态分布的情况。研究人员首先对年度CO2浓度增量序列进行了正态性检验和拟合,结果表明年度CO2增长数据符合正态分布的特点。

eb8991c42fcce549f3025a88717baa1d.jpeg

与此同时,考虑到原始数据分布未知,团队也运用了非参数检验——秩和检验(Rank-sum test),它不受数据总体分布的具体形式影响,适用于各种分布类型的数据。


然后,将2004年的十年平均增量作为参考序列,从1968年至2003年共36个十年平均增量序列作为比较序列。

经过两种统计检验后,得到以下结果:


a5bab5d7f7659f47c1acb69f032350c3.jpeg

T-Test和秩和检验结果显示在H矩阵中,大部分情况下H值为1,意味着2004年的十年平均增量与对比样本存在显著差异,即2004年序列相对于该年序列的增长速度显著·;

虽然也有部分H值为0,但总体来看,2004年的增长仍表现出了显著性增强的特征。



CO2浓度预测模型的建立与比较

研究团队探究和预测未来二氧化碳(CO2)浓度的变化。为了检验和对比模型,他们将CO2浓度数据集划分为训练集和测试集,其中80%的数据用于训练模型,剩余20%的数据用于检验模型。

接下来研究团队,建立了多种模型来模拟和预测CO2浓度随时间的变化情况。

1、线性模型

简单直观地描述CO2浓度随时间的线性增长趋势。


769abbb718acce2e3cc670afe9be051b.jpeg


2、多项式模型

考虑了非线性的增长模式,可能涉及二次或更高次项,以便更好地捕捉潜在的加速增长现象。


93e773562d599b9db02cf4d5efe7b398.jpeg


3、指数模型

假设CO2浓度随时间呈指数增长,反映出可能出现的快速累积效应。


54938d953a4b75faed01c248a82d43a7.jpeg


4、S型模型(sigmoid函数)

这类模型通常用于描述生长过程先快后慢的饱和状态,用来推测CO2浓度是否接近某种极限值。


ff9e8ed1dfc5b7ee7e13b77f75c6e8d4.jpeg


5、ARIMA模型

一种广泛应用于时间序列分析中的自回归积分滑动平均模型,它可以捕获时间序列数据中的季节性和趋势信息。


df1f61639c303c023943153992a681e3.jpeg


6、评价指标

研究小组针对每种模型分别在测试集上评估它们的性能。评价标准包括模型拟合优度R-squared、均方根误差(RMSE)、平均绝对误差(MAE)和平均相对误差(MRE)等多个统计指标。


e4a8f08c4ecaab434506504dd02e31fd.jpeg

通过在测试集上比较不同模型的预测效果,经过模型性能评估和误差分析,发现指数模型在训练集和测试集中均表现出最佳准确度。


a2c375dd9ccb543931bcf92a59ebeadc.jpeg


7、未来预测

最终,研究团队使用各个模型对未来CO2浓度进行预测。预测到2100年时,线性模型给出的预测值为534.8825 ppm,而指数模型则预测出688.3268 ppm的浓度水平。而对于达到685ppm的时间点,不同模型预测结果也不尽相同,其中指数模型、二次模型和ARIMA模型预测的时间点都在21世纪中叶附近,而S型模型则表明在当前趋势下,CO2浓度可能永远无法达到685ppm这一水平。


d7ef50552ad68650dae3b74d942a50f4.jpeg



资料获取

今天分享就到这里,后续将持续更新优秀论文解析。近20年优秀论文与HIMCM备考软件、资料已备好,需要的小伙伴可关注微信公众“思谋EDU”,添加助理微信免费领取。


78G超全数学建模资料:

db3ff6014c3811b2c6cde4b7c03f3df0.jpeg

历年HIMCM优秀论文——O奖:

890afea181bb82cd8af32c21ce55e6ee.jpeg



  • 7
    点赞
  • 29
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值