查看原文
其他

GSOD全球逐日气象站点数据介绍与下载

走天涯徐小洋 走天涯徐小洋地理数据科学 2023-10-24

前面给大家分享过用R语言获取NOAA ISD气象数据,那个是逐小时的,想要做逐日或者逐月的研究需要的合成比较复杂。在这里给大家推荐一套逐日的气象数据。

Global Surface Summary of the Day - GSOD数据,本套数据衍生自The Integrated Surface Hourly (ISH)数据集,从1929年至今,包括全球超9000个站点的气象观测数据。

  • GSOD官网:https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C00516#Documentation

GSOD包含的主要内容

从这套数据的内容来看,和ISD数据差异还是比较大的,这个字段比较少,字段单位比较奇葩,是英制单位,不是常见的国际公制单位,使用时需要进行单位换算

字段说明

  • STATION 气象站代码,WMO/DATSAV3合并w/WBAN码
  • DATE 日期,mm/dd/yyyy格式
  • LATITUDE 小数点格式的纬度,南半球负值
  • LONGITUDE 小数点格式的经度,西半球负值
  • ELEVATION 高程,单位米
  • NAME, 气象站名称
  • TEMP, Mean temperature平均温度,单位Fahrenheit华氏度,精确到一位小数,缺失值9999.9
  • TEMP_ATTRIBUTES 参与平均温度计算的观测值个数
  • DEWP,Mean dew point平均露点,单位Fahrenheit华氏度,精确到一位小数,缺失值9999.9
  • DEWP_ATTRIBUTES 参与平均露点计算的观测值个数
  • SLP, Mean sea level pressure平均海平面压强,单位mb毫巴mbar,1mb=1hPa百帕斯卡,缺失值9999.9
  • SLP_ATTRIBUTES,参与平均海平面压强计算的观测值个数
  • STP,Mean station pressure平均测站压强,单位mb 毫巴mbar,缺失值9999.9
  • STP_ATTRIBUTES,参与平均测站压强计算的观测值个数
  • VISIB,Mean visibility平均能见度,单位miles英里,缺失值999.9
  • VISIB_ATTRIBUTES,参与平均能见度计算的观测值个数
  • WDSP,Mean wind speed平均风速,单位knots节(海里/小时),缺失值999.9
  • WDSP_ATTRIBUTES,参与平均风速计算的观测值个数
  • MXSPD,Maximum sustained wind speed最大持续风速,单位knots节,缺失值999
  • GUST,Maximum wind gust最大阵风,单位knots节,缺失值999.9
  • MAX,Maximum temperature最高温度,单位Fahrenheit华氏度,缺失值9999.9
  • MAX_ATTRIBUTES
    • 空值,自动观测到的最高温度
    • *星号,从逐小时数据中提取得到的最高温度
  • MIN, Minimum temperature最低温度,单位Fahrenheit华氏度
  • MIN_ATTRIBUTES
    • 空值,自动观测到的最低温度
    • *星号,从逐小时数据中提取得到的最低温度
  • PRCP,Precipitation amount降水量(雨或融化的雪),单位inches英寸,缺失值99.99,“0”值表示当日没有测量到降水
  • PRCP_ATTRIBUTES
    • A,1 report of 6-hour precipitation amount.
    • B,Summation of 2 reports of 6-hour precipitation amount.
    • C,3次报告6小时降水量的和
    • D,4次报告6小时降水量的和
    • E,1次报告12小时降水量
    • F,2次报告12小时降水量的和
    • G,1次报告24小时降水量
    • H,气象站报告当日降水量为0,但是在逐小时观测中报告了至少一次降水,可能是当日数据不完整
    • I,气象站当日没有报告任何降水,逐日观测中也没有报告任何降水,但是仍然可能存在降水只是没有被报告
  • SNDP,Snow depth积雪深度,单位inches英寸,取当日最后一次观测值(如果观测了多次),缺失值999.9
  • FRSHTT,指示器,当日报告的天气情况,是否有雾Fog ('F' - 1st digit)、雨Rain or Drizzle ('R' - 2nd digit)、雪Snow or Ice Pellets ('S' - 3rd digit)、冰雹Hail ('H' - 4th digit)、雷Thunder ('T' - 5th digit)、龙卷风或漏斗云Tornado or Funnel Cloud ('T' - 6th digit)。
    • 1,有
    • 0,没有

单位换算

  • 压强单位:1hPa = 1mb
  • 温度单位: ,常见温度对照: 
  • 距离单位:
    • 英里转公里
    • 英寸转毫米
  • 速度单位: 1节=1海里/小时=1.852km/h

数据下载

  • 全球数据下载地址:https://www.ncei.noaa.gov/data/global-summary-of-the-day/access/
这里是所有的数据
  • 按条件筛选数据下载:https://www.ncei.noaa.gov/access/search/data-search/global-summary-of-the-day
按条件筛选数据并下载

R语言下载GSOD数据

既然有处理ISD数据的R包,我猜肯定也有处理GSOD数据的R包,还真让我找到了,目前CRAN上面公开的就是这个GSODR。

和其他CRAN发布的R包一样,直接安装即可

加载GSODR包,包中内置isd_history.rda数据,内含气象站基本情况。

library(GSODR)
library(tidyverse)
load(system.file("extdata""isd_history.rda", package = "GSODR"))
as.tibble(isd_history)
isd_history

关于GSODR更多内容可以查看GSODR网站:

  • https://ropensci.github.io/GSODR/articles/GSODR.html

或者等下期内容我再介绍如何使用GSODR和tidyverse下载和整理气象数据。

多多点赞再看,或者加入培训班会员据说可以让下期推文更快出来。

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存