CMIP6数据处理及在气候变化、水文、生态等领域中应用

气候变化对农业、生态系统、社会经济以及人类的生存与发展具有深远影响,是当前全球关注的核心议题之一。IPCC(Intergovernmental Panel on Climate Change,政府间气候变化专门委员会)的第六次评估报告明确;指出,自20世纪50年代以来,全球平均气温和海温的上升、广泛的积雪和冰川融化以及全球海平面的升高,无一不在证明气候变暖是无可争议的事实。为了对未来气候进行评估,科学家通常使用全球气候模型进行预测。

全球气候模型(Global Climate Model, GCM),亦称全球环流模型或全球大气模型,是一种数值模型,被广泛用于模拟地球的气候系统。GCM利用一系列的数学公式来描绘气候系统的各个主要组成部分,包括大气、海洋、冻土以及地表和海洋表面的生物地理过程。GCM的空间和时间精度可以根据需要进行调整。这些模型为我们提供了理解气候系统运行机制的途径,为预测气候变化趋势、评估气候变化对人类社会和生态系统的影响以及制定应对气候变化的策略提供了关键工具。

  为了进一步理解气候变化,世界气候研究计划(World Climate Research Programme, WCRP)发起了气候模型比较计划(Climate Model Intercomparison Project,CMIP)。CMIP的主要目标是收集和比较各种全球气候模型的模拟结果,以理解和预测过去、现在和未来的气候变化。

  CMIP6数据被广泛应用于全球和地区的气候变化研究、极端天气和气候事件研究、气候变化影响和风险评估、气候变化的不确定性研究、气候反馈和敏感性研究以及气候政策和决策支持等多个领域。这些数据为我们理解和预测气候变化,评估气候变化的影响和风险,以及制定有效的气候政策和决策提供了关键的信息和工具。

原文:最新CMIP6数据处理及在气候变化、水文、生态等领域中的实践技术应用

专题一  CMIP6中的模式比较计划

1.1 GCM介绍

全球气候模型(Global Climate Model, GCM),也被称为全球环流模型或全球大气模型,是一种用于模拟地球的气候系统的数值模型。这种模型使用一系列的数学公式来描述气候系统的主要组成部分,包括大气、海洋、冰冻土壤以及地表和海洋表面的生物地理过程。GCM在空间和时间上的精度可以根据需求进行调整,通常的分辨率可以从几百公里到几公里,时间步长可以从几分钟到几小时。

图片

1.2 CMIP介绍

CMIP,全称为气候模型比较计划(Climate Model Intercomparison Project),是由世界气候研究计划(World Climate Research Programme,WCRP)发起的一个国际合作项目。其目的是通过收集和比较各种全球气候模型(GCMs)的模拟结果,以理解过去的、现在的和未来的气候变化。

1.3相关比较计划介绍

图片

专题二 数据下载

2.1方法一:手动人工

图片

利用官方网站

2.2方法二:自动

利用Python的命令行工具

图片

2.3方法三:半自动购物车

图片

利用官方网站

2.4 裁剪netCDF文件

基于QGIS和CDO实现对netCDF格式裁剪

图片

QGIS中的操作

图片

裁剪效果

2.5 处理日期非365天的GCM

以BCC为例处理

专题三 基础知识

3.1 Python基础

Python 是一种高级的、解释型的编程语言,其语法简洁明了,适合快速开发。在大气科学中,Python 以其丰富的科学计算和数据分析库备受青睐。这些库如 Numpy,Scipy,Pandas 和 Xarray 等,为处理大气科学数据提供了强大的支持。

●Numpy:Numpy 是 Python 中用于科学计算的核心库,提供了高性能的多维数组对象及相关工具。对于大气科学数据的处理,例如温度、压力、风速等通常都会使用到多维数组。Numpy 提供了丰富的函数库来处理这些数组,包括数学运算、逻辑运算、形状操作、排序、选择等操作。

●Scipy:Scipy 是基于 Python 的开源软件,用于科学计算中的数值积分和微分方程数值求解,线性代数,优化,信号处理等。在大气科学中,例如对气温、气压等数据进行傅立叶分析,求解大气动力学中的偏微分方程等,都可以使用 Scipy 来实现。

●Pandas:Pandas 是基于 Numpy 构建的,使数据清洗和分析工作变得更快更简单。Pandas 是专门为处理表格和混杂数据设计的,而 Numpy 更适合处理统一的数值数组数据。在大气科学中,例如对气象站的观测数据进行时间序列分析,处理混合类型的气象数据,以及对数据进行清洗、筛选和统计等操作,Pandas 都是非常有用的工具。

3.2 CDO基本操作

CDO(Climate Data Operator)是大气科学领域常用的一款气候和气象数据处理工具。它是一个功能强大的命令行工具,可以处理和分析格网和无格网数据,支持多种数据格式,包括netCDF、GRIB、SERVICE, EXTRA和IEG。

CDO提供了一套丰富的函数库,可以用来进行各种常见的数据操作,包括:

●基础操作:如选择、提取和修改变量、维度、属性等。

●数值操作:如四则运算、统计运算、函数运算等。例如,可以计算数据的平均值、最大值、最小值、标准差等。

●空间操作:如重新格网、插值、汇总、选择和提取地理区域等。

●时间操作:如选择和提取时间周期、计算时间平均或累积等。

3.3 Xarray的基本操作

Xarray 是一个用于处理多维数组数据的 Python 库,它在 numpy 的基础上提供了一系列用于数据操作和分析的高级接口,并能很好地支持 netCDF 这类基于网络的自描述数据格式,因此在大气科学和气候科学中被广泛使用。

Xarray 的主要特点包括:

●基于标签的数据操作:Xarray 使用维度名称而不是轴编号进行数据选择和操作,极大地增强了代码的可读性和可维护性。

●自动对齐数据:在进行运算时,Xarray 可以自动对齐不同数据集的变量(variables)和坐标(coordinates)。

●分组运算和数据透视:Xarray 支持类似于 pandas 的分组运算(group-by)和数据透视(pivot)功能。

●l/O操作:Xarray 对多种数据格式提供了非常好的支持,尤其是对 netCDF 数据的读取和写入。

专题四  单点降尺度

4.1 Delta方法

图片

     Delta方法(Delta Change Method),也称为增量方法或差值方法,是气候模型降尺度的一种简单而常用的方法。该方法假设气候变化的幅度在未来相对于历史期间将保持恒定。因此,对于某一具体的未来时段,可以通过计算过去和现在气候的差值(即 delta),并将其应用到未来的气候预测上,来预估未来的气候状态。该方法可以应用于温度和降水等气候变量的预测。

4.2统计订正

概率分布函数(Probability Density Function, PDF)的订正。这种方法的基本思想是:通过修改大尺度模型输出的PDF,使其更符合观测数据的PDF,从而获得更准确的小尺度气候变量。

图片

4.3机器学习方法

降尺度是将粗尺度的全球气候模型(GCM)输出数据转换为地面更精细尺度的过程。机器学习方法因其在处理复杂模式识别和高维数据问题的强大能力,已经被成功应用于降尺度技术。在气候学领域,机器学习已被成功用于将粗尺度的气候模型输出(例如,温度和降水)与其他环境变量(例如,地形和土壤类型)关联,以获得更高分辨率的气候预测。

实现步骤

●建立特征

● 建立模型

●模型评估

图片

4.4多算法集成方法

多算法的集成

图片

图片

贝叶斯模型平均 (Bayesian Model Averaging, BMA)

贝叶斯模型平均是一种统计方法,用于根据观察数据确定各种模型的后验概率。与选择一个最好的模型相反,贝叶斯模型平均考虑了所有可能的模型,然后根据每个模型的后验概率进行加权平均。

Python+pymc3实现

专题五 统计方法的区域降尺度

5.1 Delta方法

图片

5.2 基于概率订正方法的

专题六 基于WRF模式的动力降尺度

      动态降尺度通常使用更高分辨率的区域气候模型(RCM),这些模型在更大尺度的全球气候模型驱动下运行。其中,WRF(Weather Research and Forecasting)模型是目前使用最广泛的区域气候模型之一。

      WRF模型是一个灵活的、大气环流模型,适合用于各种尺度的气候和气象研究。它的主要特点是具有高分辨率(可达到几公里),并且可以考虑到许多重要的地球物理过程,如云的形成、降水、陆面过程、海洋过程、边界层过程、辐射、化学过程等。

6.1制备CMIP6的WRF驱动数据

利用cdo工具对gcm的输出文件进行重新编码制备wrf的驱动数据

6.1.1针对压力坐标系的数据制备

6.1.2针对sigma坐标系GCM数据制备

6.1.3 WPS处理

图片

图片

6.2 WRF模式运行

图片

6.3 模式的后处理

● 提取变量

●变量的统计

●变量的可视化

图片

图片

专题七 典型应用案例-气候变化1

7.1针对风速进行降尺度

图片

7.2针对短波辐射降尺度

图片

专题八 典型应用案例-气候变化2

ECA极端气候指数计算

    ECA (European Climate Assessment) 是欧洲的一个气候评估项目,其在全球范围内发布了一系列的极端气候事件指数。这些指数被广泛用于气候变化研究,特别是在研究极端天气和气候事件方面。

   ECA 的极端气候指数主要包括以下几类:

温度指数:这些指数主要用于度量温度的极端情况,例如热日数(TX90p,年中最高气温超过90百分位数的天数)、冷日数(TN10p,年中最低气温低于10百分位数的天数)、热夜数(TN90p,年中最低气温超过90百分位数的天数)、冷夜数(TN10p,年中最低气温低于10百分位数的天数)等。

  降水指数:这些指数主要用于度量降水的极端情况,例如最大连续5日降水量(RX5day)、大于或等于10mm的降水日数(R10mm)、大于或等于20mm的降水日数(R20mm)、降水强度(SDII)等。

这些指数对于理解和预测极端气候事件的影响非常重要,因为极端气候事件(如热浪、干旱、洪水等)往往比平均气候变化带来更大的影响。因此,对这些指数的研究有助于我们更好地理解和适应气候变化。

lConsecutive dry days index

图片

lConsecutive frost days index per time period

图片

lConsecutive summer days index per time period

图片

lConsecutive wet days index per time period

图片

专题九 典型应用案例-生态领域

预估生长季开始和结束时间

1、建立气象数据与VIPPHEN遥感物候数据中生长季开始和结束

2、在未来气候情景下预估生长季长季开始、结束和长度

专题十 典型应用案例-水文、生态模式数据

● SWAT数据制备

●Biome-BGC数据

Biome-BGC是利用站点描述数据、气象数据和植被生理生态参数,模拟日尺度碳、水和氮通量的模型,其研究的空间尺度可以从点尺度扩展到陆地生态系统。案例中以单点模拟方式制备CMIP6的气象数据。

图片

原文:最新CMIP6数据处理及在气候变化、水文、生态等领域中的实践技术应用 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/185383.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

西班牙Wallapop是什么?原来欧洲版闲鱼也很好用!

说到国内的闲鱼大家肯定不陌生,那国外的二手闲置平台大家知道吗?在西班牙,最受欢迎的移动购物APP是Wallapop和速卖通。Wallapop是西班牙第一大二手商品网站,网站上丰富的性价比高的商品正好满足了西班牙人的需求。今天龙哥就和大家…

Yolov8实现瓶盖正反面检测

一、模型介绍 模型基于 yolov8n数据集采用SKU-110k,这数据集太大了十几个 G,所以只训练了 10 轮左右就拿来微调了 基于原木数据微调:训练 200 轮的效果 10 轮SKU-110k 20 轮原木 200 轮瓶盖正反面 微调模型下载地址https://wwxd.lanzouu.co…

网站提示不安全?

随着互联网的普及和发展,网络安全问题日益严重。黑客攻击、数据泄露、恶意软件等问题层出不穷,给企业和个人带来了巨大的损失。在这个背景下,确保网站安全显得尤为重要,而使用SSL证书是解决这些问题的有效措施。 什么是SSL证书&am…

分享超实用的软文撰写步骤!建议收藏

一想到写软文就头大,根本不知道从哪里下手,这是很多写手在创作过程中会遇到的问题。 一篇软文写得好不好,关键就要看你的创作步骤到不到位,软文创作是有一套可执行的具体方式的,跟着步骤来,你也能轻轻松松…

提升性能测试效率:JMeter中的用户自定义变量!

前言 在测试过程中,我们经常会碰到测试服务地址有改动的情况,为了方便,我们会把访问地址参数化,当访问地址变化了,我们只需要把参数对应的值改动一下就可以了。 一:添加配置元件-用户定义的变量&#xff…

LeetCode 7 整数反转

题目描述 整数反转 给你一个 32 位的有符号整数 x ,返回将 x 中的数字部分反转后的结果。 如果反转后整数超过 32 位的有符号整数的范围 [−2^31, 2^31 − 1] ,就返回 0。 假设环境不允许存储 64 位整数(有符号或无符号)。 示…

【Android】使用intent.putExtra()方法在启动Activity时传递数据

食用方法 在Android中,你可以使用Intent对象来在启动Activity时传递数据。以下是一个示例,展示了如何在startActivity时传递数据到被启动的Activity: 在启动Activity的地方,创建一个Intent对象,并使用putExtra()方法…

【PPT模板合集】关于自制内容的PPT模板合集,包括原创的PPT及改良内容的PPT,适合科研/比赛/工作

【PPT模板合集】关于自制内容的PPT模板合集,包括原创的PPT及改良内容的PPT,适合科研/比赛/工作 零、前言一、校园层面的PPT模板1.1 各种毕业答辩1.2 夏令营答辩1.3 奖学金答辩1.4 比赛/项目答辩 二、学术层面的PPT模板2.1 学术汇报2.2 会议海报类型 三、…

滚珠丝杆在各种自动化设备中的作用

滚珠丝杆因其具有高精度、高刚度和长寿命等特性,成为许多设备中的重要组成部分,在许多行业中都有广泛的应用,接下来我们看看滚珠丝杆的具体应用有哪些? 1、打孔机:提供精确的导向,使打孔机的滑块能够沿固定…

【算法每日一练]-图论(保姆级教程篇9 最小生成树 ,并查集篇)#道路修建 #兽径管理

目录 题目:道路修建 思路: 题目:兽径管理 思路: 题目:道路修建 思路: “让这些点全部连在一起的最小代价”很明显是最小生成树。绝对不能kruskal,存边一定会超内存。所以只能prim。 但是…

解读免费化潮流:为何数据可视化软件向免费迈进?

近年来,我们见证了数据可视化软件呈现出明显的免费化趋势。这个趋势的背后隐藏着许多关键原因,影响着整个数据行业的发展和走向。为何数据可视化软件开始朝着免费方向发展?让我们一同深入探讨。 普及数字化需求: 数字化已经深入到…

力扣283题 移动零 双指针解法

移动零 给定一个数组 nums,编写一个函数将所有 0 移动到数组的末尾,同时保持非零元素的相对顺序。 请注意 ,必须在不复制数组的情况下原地对数组进行操作。 示例 1: 输入: nums [0,1,0,3,12] 输出: [1,3,12,0,0]示例 2: 输入: nums [0] 输出…

windows下使用iperf3

提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档 文章目录 前言一、使用步骤1.下载2.解压3.测试1.作为服务器使用2.作为客户端使用 总结 前言 windows使用iperf3和linux是一样的,而且可以相互操作。 一、使用步骤…

Linux系统iptables扩展

目录 一. iptables规则保存 1. 导出规则保存 2. 自动重载规则 ①. 当前用户生效 ②. 全局生效 二. 自定义链 1. 新建自定义链 2. 重命名自定义链 3. 添加自定义链规则 4. 调用自定义链规则 5. 删除自定义链 三. NAT 1. SNAT 2. DNAT 3. 实验 ①. 实验要求 ②. …

mybatis关于namespace以及id以及Mapper接口命名的说明(了解)

1、建库建表 CREATE DATABASE mybatis-example;USE mybatis-example;CREATE TABLE t_emp(emp_id INT AUTO_INCREMENT,emp_name CHAR(100),emp_salary DOUBLE(10,5),PRIMARY KEY(emp_id) );INSERT INTO t_emp(emp_name,emp_salary) VALUES("tom",200.33); INSERT INTO…

C/C++转义符:\x

文章目录 什么是转义符使用"\x"定义char数组宏定义中的\ 什么是转义符 在C语言中,转义符用于将一些特殊字符表示为单个字符,常用的转义符有: \\:反斜杠符号\:单引号\":双引号\a&#xff1…

如何获取1688的订单详情

获取1688订单详情需要申请 第一步:先去1688-开放平台申请,申请不一定能通过,审批很严。 第二步:首次登录请先注册,注册成功后即可登录。 第三步:选择:我是第三方开发者 第三方是封装好的ap…

解析编程中的技术迷题:常见挑战与应对策略

前言 在数字化时代的浪潮中,编程已经成为了一项至关重要的技能。无论是在软件开发、数据分析、人工智能还是互联网领域,编程都扮演着不可或缺的角色。作为一种创造性的活动,编程不仅仅是代码的书写,更是一种解决问题和创新的思维方…

.Net 8 Blazor下 Auto交互渲染模式试用

一、环境 C:\Users\zhuji>dotnet --version 8.0.100C:\Users\zhuji>dotnet --list-sdks 5.0.403 [C:\Program Files\dotnet\sdk] 6.0.404 [C:\Program Files\dotnet\sdk] 8.0.100 [C:\Program Files\dotnet\sdk] Microsoft Visual Studio Enterprise 2022 (64 位) - Cu…

Jmeter进阶使用:BeanShell实现接口前置和后置操作!

一、背景 我们使用Jmeter做压力测试或者接口测试时,除了最简单的直接对接口发起请求,很多时候需要对接口进行一些前置操作:比如提前生成测试数据,以及一些后置操作:比如提取接口响应内容中的某个字段的值。举个最常用…