当前浏览: 水文


祁连山综合观测网:黑河流域地表过程综合观测网(景阳岭站自动气象站-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网上游景阳岭站气象要素观测数据。站点位于青海省祁连县景阳岭垭口,下垫面是高寒草甸。观测点的经纬度是101.1160E, 37.8384N,海拔3750m。空气温度、相对湿度传感器架设在5m处,朝向正北;气压计安装在地面上的防撬箱内;翻斗式雨量计安装在10m处;风速与风向传感器架设在10m,朝向正北;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下4cm、10cm、20cm、40cm、80cm、120cm、160cm处,在距离气象塔2m的正南方;土壤水分探头埋设在地下4cm、10cm、20cm、40cm、80cm、120cm、160cm处,在距离气象塔2m的正南方;土壤热流板(3块)依次埋设在地下6cm处,在距离气象塔2m的正南方。 观测项目有:空气温湿度(Ta_5m、RH_5m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、风速(WS_10m)(单位:米/秒)、风向(WD_10m)(单位:度)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤温度(Ts_0cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_80cm、Ts_120cm、Ts_160cm)(单位:摄氏度)、土壤水分(Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_80cm、Ms_120cm、Ms_160cm)(单位:百分比)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;由于雪覆盖太阳能板引起供电不足,导致2.3-3.18间数据缺失;由于传感器问题,6.6-6.20、6.22-7.30风速风向出现较多NAN无效值,9-10月四分量向上短波辐射数据出错。(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2019-9-10 10:30;(6)命名规则为:AWS+站点名称。 黑河综合观测网或站点信息请参考Liu et al. (2018)和Che et al. (2019),观测数据处理请参考Liu et al. (2011)。

2020-06-24

祁连山综合观测网:黑河流域地表过程综合观测网(阿柔超级站大孔径闪烁仪-2019)

该数据集包含了黑河流域地表过程综合观测网上游阿柔站的大孔径闪烁仪通量观测数据。上游阿柔站分别架设了BLS900和RR-RSS460型号的大孔径闪烁仪,北塔为RR-RSS460的接收端和BLS900的发射端,南塔为RR-RSS460的发射端和BLS900的接收端。观测时间为2019年1月1日至2019年12月31日。站点位于青海省祁连县阿柔乡草达坂村,下垫面是高寒草地。北塔的经纬度是100.4712E,38.0568N,南塔的经纬度是100.4572E,38.0384N,海拔高度约3033m。大孔径闪烁仪的有效高度13.0m,光径长度是2390m,采样频率是1min。 大孔径闪烁仪原始观测数据为1min,发布的数据为经过处理与质量控制后的数据,其中感热通量主要是结合自动气象站观测数据,基于莫宁-奥布霍夫相似理论通过迭代计算得到,主要的质量控制步骤包括:(1)剔除Cn2达到饱和的数据(BLS900:Cn2>7.25E-14,RR-RSS460:Cn2>7.84 E-14);(2)剔除解调信号强度较弱的数据(BLS900:Average X Intensity<1000;RR-RSS460:Demod>-20mv);(3)剔除降水时刻的数据;(4)剔除稳定条件下的弱湍流的数据(u*小于0.1m/s)。在迭代计算过程中,对于BLS900,选取Thiermann and Grassl(1992)的稳定度普适函数;对于RR-RSS460,选取Andreas(1988)的稳定度普适函数,详细介绍请参考Liu et al. (2011, 2013)。由于仪器调整和供电不足,大孔径闪烁仪数据由于信号和供电等问题缺失的日期为:2019.04.20-2019.04.30;2019.05.07-2019.05.13;2019.10.21-2019.10.25;2019.12.30-2019.12.31。 关于发布数据的几点说明:(1)上游LAS数据以BLS900为主,缺失时刻由RR-RSS460观测补充,两者都缺失则以-6999标记。(2)数据表头:Date/Time :日期/时间(格式:yyyy/m/d h:mm),Cn2 :空气折射指数结构参数(单位:m-2/3),H :感热通量(单位:W/m2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xlsx格式存储。 黑河流域地表过程综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。

2020-06-23

祁连山综合观测网:黑河流域地表过程综合观测网(四道桥超级站大孔径闪烁仪-2019)

该数据集包含了黑河流域地表过程综合观测网下游四道桥站的大孔径闪烁仪通量观测数据。下游四道桥站分别架设了BLS450、BLS900和RR-RSS460型号的大孔径闪烁仪,北塔为接收端,南塔为发射端。观测时间为2019年1月1日至2019年12月31日。站点位于内蒙古额济纳旗,下垫面是柽柳、胡杨、裸地和耕地。北塔的经纬度是101.137E,42.008N,南塔的经纬度是101.131E,41.987N,海拔高度约873m。大孔径闪烁仪的有效高度25.5m,光径长度是2350m,采样频率是1min。 大孔径闪烁仪原始观测数据为1min,发布的数据为经过处理与质量控制后的数据,其中感热通量主要是结合自动气象站观测数据,基于莫宁-奥布霍夫相似理论通过迭代计算得到,主要的质量控制步骤包括:(1)剔除Cn2达到饱和的数据(BLS900&BLS450:Cn2>7.25E-14,RR-RSS460:Cn2>7.84 E-14);(2)剔除解调信号强度较弱的数据(BLS900&BLS450:Average X Intensity<1000;RR-RSS460:Demod>-20mv);(3)剔除降水时刻的数据;(4)剔除稳定条件下的弱湍流的数据(u*小于0.1m/s)。在迭代计算过程中,对于BLS900&BLS450,选取Thiermann and Grassl(1992)的稳定度普适函数;对于RR-RSS460,选取Andreas(1988)的稳定度普适函数,详细介绍请参考Liu et al. (2011, 2013)。由于仪器故障,大孔径闪烁仪数据缺失的日期为:2019.04.14-2019.04.25;2019.11.02-2019.11.12。 关于发布数据的几点说明:(1)下游LAS数据以BLS900为主,其次为BLS450,最后为RR-RSS460,最终缺失时刻以-6999标记。(2)数据表头:Date/Time :日期/时间(格式:yyyy/m/d h:mm),Cn2 :空气折射指数结构参数(单位:m-2/3),H:感热通量(单位:W/m2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xlsx格式存储。 黑河流域地表过程综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。

2020-06-23

祁连山综合观测网:黑河流域地表过程综合观测网(四道桥超级站气象要素梯度观测系统-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网下游四道桥超级站气象要素梯度观测系统数据。站点位于内蒙古额济纳旗达来呼布镇四道桥,下垫面是柽柳。观测点的经纬度是101.1374E, 42.0012N,海拔873m。空气温度、相对湿度、风速传感器分别架设在5m、7m、10m、15m、20m、28m处,共6层,朝向正北;风向传感器架设在15m处,朝向正北;气压计安装在防水箱内;翻斗式雨量计安装在28m处;四分量辐射仪安装在10m处,朝向正南;两个红外温度计安装在10m处,朝向正南,探头朝向是垂直向下;两个光合有效辐射仪安装在10m处,朝向正南,探头垂直向上和向下方向各一个;土壤部分传感器安装在塔体南侧2m处,其中土壤热流板(自校正式)(3块)依次埋设在地下6cm处;平均土壤温度传感器TCAV埋设在地下2cm、4cm处;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、80cm、120cm、160cm和200cm处;土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、80cm、120cm、160cm和200cm处。 观测项目有:风速(WS_5m、WS_7m、WS_10m、WS_15m、WS_20m、WS_28m)(单位:米/秒)、风向(WD_15m)(单位:度)、空气温湿度(Ta_5m、Ta_7m、Ta_10m、Ta_15m、Ta_20m、Ta_28m和RH_5m、RH_7m、RH_10m、RH_15m、RH_20m、RH_28m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、向上和向下光合有效辐射(PAR_U_up、PAR_U_down)(单位:微摩尔/平方米秒)、平均土壤温度(TCAV)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤水分(Ms_2cm、Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_80cm、Ms_120cm、Ms_160cm、Ms_200cm)(单位:体积含水量,百分比)、土壤温度(Ts_0cm、Ts_2cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_80cm、Ts_120cm、Ts_160cm、Ts_200cm)(单位:摄氏度)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;4-10月由于传感器的问题,2m深度土壤温度数据出错;11.13日后红外温度2数据出错;9.15日后土壤温湿度断续出现错误;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2019-9-10 10:30;(6)命名规则为:AWS+站点名称。 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。

2020-06-22

祁连山综合观测网:黑河流域地表过程综合观测网(四道桥超级站宇宙射线观测系统-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网下游四道桥超级站宇宙射线观测系统数据。站点位于内蒙古额济纳旗四道桥,下垫面是柽柳。观测点的经纬度是101.1374E, 42.0012N,海拔873 m。仪器探头底部距地面0.5m,采样频率是1小时。 宇宙射线仪器的原始观测项目包括:电压Batt(V)、温度T(℃)、相对湿度RH(%)、气压P(hPa)、快中子数N1C(个/小时)、热中子数N2C(个/小时)、快中子采样时间N1ET(s)及热中子采样时间N2ET(s)。发布的数据为经过处理计算后的数据,数据表头包括:Date Time(日期 时间)、P(气压 hPa)、N1C(快中子数 个/小时)、N1C_cor(气压订正的快中子数 个/小时)和VWC(土壤体积含水量 %),其处理的主要步骤包括: 1) 数据筛选 数据筛选共四条标准:(1)剔除电压小于和等于11.8伏特的数据;(2)剔除空气相对湿度大于和等于80%的数据;(3)剔除采样时间间隔不在60±1分钟内的数据;(4) 剔除快中子数较前后一小时变化大于200的数据。剔除及缺失数据用-6999补充。 2) 气压订正 根据仪器说明手册中提到的快中子气压订正公式,对原始数据进行气压订正,得到订正后的快中子数N1C_cor。 3) 仪器率定 在计算土壤水分的过程中需要对计算公式中的N0进行率定。N0为土壤干燥条件下的快中子数,通常使用测量源区内的土样得到实测土壤水分(或者通过比较密集的土壤水分无线传感器获取)θm(Zreda et al. 2012)和对应时间段内的快中子校正数据N,再通过公式反求得到N0。 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Zhu et al. (2015)。

2020-06-22

祁连山综合观测网:黑河流域地表过程综合观测网(大满超级站宇宙射线观测系统-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网中游大满超级站宇宙射线观测系统数据。站点位于甘肃省张掖市大满灌区农田内,下垫面是玉米田。观测点的经纬度是100.3722E, 38.8555N,海拔1556m。仪器探头底部距地面0.5m,采样频率是1小时。 宇宙射线仪器的原始观测项目包括:电压Batt(V)、温度T(℃)、相对湿度RH(%)、气压P(hPa)、快中子数N1C(个/小时)、热中子数N2C(个/小时)、快中子采样时间N1ET(s)及热中子采样时间N2ET(s)。发布的数据为经过处理计算后的数据,数据表头包括:Date Time(日期 时间)、P(气压 hPa)、N1C(快中子数 个/小时)、N1C_cor(气压订正的快中子数 个/小时)和VWC(土壤体积含水量 %),其处理的主要步骤包括: 1) 数据筛选 数据筛选共四条标准:(1)剔除电压小于和等于11.8伏特的数据;(2)剔除空气相对湿度大于和等于80%的数据;(3)剔除采样时间间隔不在60±1分钟内的数据;(4) 剔除快中子数较前后一小时变化大于200的数据。剔除及缺失数据用-6999补充。 2) 气压订正 根据仪器说明手册中提到的快中子气压订正公式,对原始数据进行气压订正,得到订正后的快中子数N1C_cor。 3) 仪器率定 在计算土壤水分的过程中需要对计算公式中的N0进行率定。N0为土壤干燥条件下的快中子数,通常使用测量源区内的土样得到实测土壤水分(或者通过比较密集的土壤水分无线传感器获取)θm(Zreda et al. 2012)和对应时间段内的快中子校正数据N,再通过公式反求得到N0。 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Zhu et al. (2015)。

2020-06-21

祁连山综合观测网:黑河流域地表过程综合观测网(花寨子荒漠站自动气象站-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网中游花寨子荒漠站气象要素观测数据。站点位于甘肃省张掖市花寨子,下垫面是盐爪爪山前荒漠。观测点的经纬度是100.3201E, 38.7659N,海拔1731m。空气温度、相对湿度传感器架设在5m、10m处,朝向正北;气压计安装防水箱内;翻斗式雨量计安装在10m处;风速风向传感器架设在5m、10m处,朝向正北;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,在距离气象塔2m的正南方;土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,在距离气象塔2m的正南方;土壤热流板(3块)依次埋设在地下6cm处。 观测项目有:空气温湿度(Ta_5m、RH_5m、Ta_10m、RH_10m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、风速(WS_5m、WS_10m)(单位:米/秒)、风向(WD_5m、WD_10m)(单位:度)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤水分(Ms_2cm、Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_60cm、Ms_100cm)(单位:体积含水量,百分比)和土壤温度(Ts_0cm、Ts_2cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_60cm、Ts_100cm)(单位:摄氏度)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;风速风向偶尔出现一些错误值,由于传感器的问题,2.28-3.10间20cm土壤温度,5-7月土壤热通量2观测数据出现错误值;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2019-6-10 10:30;(6)命名规则为:AWS+站点名称 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。

2020-06-21

祁连山综合观测网:黑河流域地表过程综合观测网(大沙龙站自动气象站-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网上游大沙龙站气象要素观测数据。站点位于青海省祁连县西侧沙龙滩地区,下垫面是沼泽化高寒草甸。观测点的经纬度是98.9406E, 38.8399N,海拔3739m。空气温度、相对湿度传感器架设在5m处,朝向正北;气压计安装在地面上的防撬箱内;翻斗式雨量计安装在10m处;风速与风向传感器架设在10m,朝向正北;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下4cm、10cm、20cm、40cm、80cm、120cm、160cm处,在距离气象塔2m的正南方;土壤水分探头埋设在地下4cm、10cm、20cm、40cm、80cm、120cm、160cm处,在距离气象塔2m的正南方;土壤热流板(3块)依次埋设在地下6cm处,并距离气象塔2m的正南方。 观测项目有:空气温湿度(Ta_5m、RH_5m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、风速(WS_10m)(单位:米/秒)、风向(WD_10m)(单位:度)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤温度(Ts_0cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_80cm、Ts_120cm、Ts_160cm)(单位:摄氏度)、土壤水分(Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_80cm、Ms_120cm、Ms_160cm)(单位:体积含水量,百分比)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2019-9-10 10:30;(6)命名规则为:AWS+站点名称。 黑河综合观测网或站点信息请参考Liu et al. (2018)和Che et al. (2019),观测数据处理请参考Liu et al. (2011)。

2020-06-21

祁连山综合观测网:黑河流域地表过程综合观测网(荒漠站自动气象站-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网下游荒漠站气象要素观测数据。站点位于内蒙古额济纳旗荒漠滩,下垫面是红砂荒漠。观测点的经纬度是100.9872E, 42.1135N,海拔1054m。空气温度、相对湿度传感器架设在5m、10m处,朝向正北;气压计安装防水箱内;翻斗式雨量计安装在10m处;风速传感器架设在5m、10m处,风向传感器架设在10m处,朝向正北;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,在距离气象塔2m的正南方;土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,在距离气象塔2m的正南方;土壤热流板(3块)依次埋设在地下6cm处。 观测项目有:空气温湿度(Ta_5m、RH_5m、Ta_10m、RH_10m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、风速(WS_5m、WS_10m)(单位:米/秒)、风向(WD_10m)(单位:度)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤水分(Ms_2cm、Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_60cm、Ms_100cm)(单位:体积含水量,百分比)和土壤温度(Ts_0cm、Ts_2cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_60cm、Ts_100cm)(单位:摄氏度)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2019-6-10 10:30;(6)命名规则为:AWS+站点名称 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。

2020-06-21

祁连山综合观测网:黑河流域地表过程综合观测网(阿柔超级站宇宙射线观测系统-2019)

该数据集包含了2019年1月1日至2019年12月31日黑河流域地表过程综合观测网上游阿柔超级站宇宙射线观测系统数据。站点位于青海省祁连县阿柔乡草达坂村,下垫面是亚高山山地草甸。观测点的经纬度是100.4643E, 38.0473N,海拔3033m。仪器探头底部距地面0.5m,采样频率是1小时。 宇宙射线仪器的原始观测项目包括:电压Batt(V)、温度T(℃)、相对湿度RH(%)、气压P(hPa)、快中子数N1C(个/小时)、热中子数N2C(个/小时)、快中子采样时间N1ET(s)及热中子采样时间N2ET(s)。发布的数据为经过处理计算后的数据,数据表头包括:Date Time(日期 时间)、P(气压 hPa)、N1C(快中子数 个/小时)、N1C_cor(气压订正的快中子数 个/小时)和VWC(土壤体积含水量 %),其处理的主要步骤包括: 1) 数据筛选 数据筛选共四条标准:(1)剔除电压小于和等于11.8伏特的数据;(2)剔除空气相对湿度大于和等于80%的数据;(3)剔除采样时间间隔不在60±1分钟内的数据;(4) 剔除快中子数较前后一小时变化大于200的数据。剔除及缺失数据用-6999补充。 2) 气压订正 根据仪器说明手册中提到的快中子气压订正公式,对原始数据进行气压订正,得到订正后的快中子数N1C_cor。 3) 仪器率定 在计算土壤水分的过程中需要对计算公式中的N0进行率定。N0为土壤干燥条件下的快中子数,通常使用测量源区内的土样得到实测土壤水分(或者通过比较密集的土壤水分无线传感器获取)θm(Zreda et al. 2012)和对应时间段内的快中子校正数据N,再通过公式反求得到N0。 黑河综合观测网或站点信息请参考Liu et al. (2018)和Che et al. (2019),观测数据处理请参考Zhu et al. (2015)。

2020-06-19