基于Python常见地球科学数据（ERA5、雪深、积雪覆盖、海温、植被指数、土地利用）处理实践技术应用

基于Python常见地球科学数据ERA5、雪深、积雪覆盖、海温、植被指数、土地利用处理一、全球大气再分析数据-ERA5数据ECMWF中心推出的ERA5全球大气再分析数据提供了大量大气、陆地和海洋气候变量的逐小时数据。这些数据在30km网格上覆盖了全球在时间跨度上从1979至今。该数据能够提供全球范围的格点气象数据。课程中将针对该数据介绍以下内容a、多年数据的读取一般在下载过程中会将按照最小的时间单位获取文件。因此需要将数据进行后处理进行合并之后开展后续的操作。在本案例中将对气温变量进行读取。b、趋势分析在多年气温数据的基础上使用mann kendall趋势检验方法判断趋势使用线性回归方法计算变化倾向率。c、多时间尺度统计利用小时数据计算日、月、年尺度数据。二、北半球雪深数据-加拿大雪深数据加拿大气象中心CMC制作了全球每日雪深业务分析数据。数据时间分辨率为1天空间分辨率为0.25°。数据包括雪深和估计雪水当量SWE2个变量。该数据可以为全球和区域尺度冰冻圈积雪相关水文和生态研究提供支持。课程中将针对该数据介绍以下内容A、将ASCII数据转化为netCDF数据将日数据进行整理。使用multiprocess工具并行转化netCDF格式数据。B、分析北半球雪深的气候态空间分布使用Xarray工具计算多年、年和月的均值。C、分析时间序列的趋势变化及与气温的关系使用线性回归分析月尺度上趋势计算雪深变化序列与气温之间的相关系数。三、积雪覆盖数据-IMS积雪覆盖数据美国国家冰雪中心NSIDC从1997年2月至今的北半球雪盖和海冰的地图。这些数据以ASCII文本和GeoTIFF格式提供并提供有三种不同的分辨率1公里、4公里和24公里。该数据可以为全球和区域尺度冰冻圈积雪相关水文和生态研究提供支持。课程中将针对该数据介绍以下内容A、将ASCII数据转化为netCDF数据将日数据进行整理。使用multiprocess工具并行转化netCDF格式数据。B、分析新疆北疆、青藏高原和东北地区气候态积雪分布使用Xarray工具计算多年、年和月的均值。使用basemap进行可视化。C、分析新疆北疆、青藏高原和东北地区积雪面积变化规律使用Xarray工具进行汇总绘制不同子区的面积图时间序列。四、海温数据-NOAA OISSTNOAA 1/4°每日最佳内插海面温度OISST是一个长期的气候数据记录它将来自不同平台卫星、船舶、浮标和Argo浮标的观测数据纳入一个定期的全球网格。该数据集经过插值处理以填补网格上的空白并创建一个空间上完整的海面温度图。卫星和船舶的观测数据以浮标为参照以补偿平台的差异和传感器的偏差。OISST数据空间分辨率0.25°时间跨度从1981至今。该数据可以为气候变化、海~气关系等研究提供支持。课程中将针对该数据介绍以下内容A、读取多年月海温数据使用Xarray工具读取海温数据B、分析海温变化趋势C、全球海温与长江流域夏季降水关系分析CN05中长江流域夏季降水与全球海温之间的相关关系五、植被指数数据-GIMMS NDVI植被指数是对植被冠层绿叶中叶绿素吸收的光合有效辐射的辐射测量因此是衡量一个地区生理功能的植被覆盖水平的指数。最新版本的GIMMS NDVI数据集横跨1981年7月至2015年12月。课程中将针对该数据介绍以下内容A、读取多年生长季NDVI数据B、计算青藏高原生长季NDVI与降水关系六、土地利用数据-MODIS土地利用数据地覆盖动态产品包括标志着季节性周期的植被生长、成熟和衰老的时间层。MODIS TerraAqua组合土地覆盖产品包含了五个不同的土地覆盖分类方案通过监督决策树分类方法得出。主要的土地覆盖方案确定了IGBP定义的17个等级包括11个自然植被等级3个人为改变的等级以及3个非植被等级。课程中将针对该数据介绍以下内容A、读取多年土地利用数据B、裁剪长江流域土地利用C、分析长时间序列长江流域土地利用变化