当前浏览: 通量观测矩阵


黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(17号点自动气象站)

该数据集包含了2012年5月12日至9月17日的通量观测矩阵中17号点的自动气象站观测数据。站点位于甘肃省张掖市大满灌区内,下垫面是果园。观测点的经纬度是100.36972E, 38.84510N,海拔1559.63m。空气温度、相对湿度传感器架设在5m处,朝向正北;气压与翻斗式雨量计安装在2m、10米处;风速与风向传感器架设在10米,朝向正北;四分量辐射仪安装在6米处,朝向正南;两个红外表面温度传感器安装在6米处,朝向正南支臂朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、60cm和100 cm处,并距离气象塔2米的正南方;土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,并距离气象塔2米的正南方;土壤热流板(3块)依次埋设在地下6 cm处,其中两块埋设在棵间,一块埋设在植株下面。红外表面温度传感器为垂直向下,在 8月6日后有一个(IRT_2)红外倾角调整为天顶角50°。观测项目有:空气温湿观测(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)(单位:瓦/平方米)、 多层土壤水分(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)确保每天1440个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2012-6-10 10:30;(5)命名规则为:AMS+站点编号 。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-10

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(戈壁站涡动相关仪)

该数据集包含了2012年5月31日至9月15日的通量观测矩阵中戈壁站涡动相关通量观测数据。站点位于甘肃省张掖市,下垫面是戈壁。观测点的经纬度是100.30420E, 38.91496N,海拔1562.00m。涡动相关仪架高4.6m,采样频率是10Hz,超声朝向是正北向,超声风速仪与CO2/H2O分析仪之间的距离是15cm。 涡动相关仪的原始观测数据为10Hz,发布的数据是采用Edire软件后处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Edire软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据每30min内缺失率大于3%的数据;(4)剔除夜间弱湍流的观测数据(u*小于0.1m/s)。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),稳定度Z/L(无量纲),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE,二氧化碳通量的质量标识QA_Fc。感热、潜热、二氧化碳通量的质量标识分为三级(质量标识0:(Δst <30, ITC<30); 1: (Δst <100, ITC<100); 其余为2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-09

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(4号点涡动相关仪)(2012年5月-9月)

该数据集包含了2012年5月31日至9月17日的通量观测矩阵中4号点的涡动相关通量观测数据。站点位于甘肃省张掖市盈科灌区内,下垫面是村庄。观测点的经纬度是100.35753E, 38.87752N,海拔1561.87m。涡动相关仪架高4.2m(8月19日后涡动相关仪架高调整为6.2m),采样频率是10Hz,超声朝向是正北向,超声风速仪与CO2/H2O分析仪之间的距离是17cm。 涡动相关仪的原始观测数据为10Hz,发布的数据是采用Edire软件后处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Edire软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据每30min内缺失率大于3%的数据;(4)剔除夜间弱湍流的观测数据(u*小于0.1m/s)。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),稳定度Z/L(无量纲),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE,二氧化碳通量的质量标识QA_Fc。感热、潜热、二氧化碳通量的质量标识分为三级(质量标识0:(Δst <30, ITC<30); 1: (Δst <100, ITC<100); 其余为2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-09

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(13号点自动气象站)

该数据集包含了2012年5月6日至9月20日的通量观测矩阵中13号点的自动气象站观测数据。站点位于甘肃省张掖市大满灌区农田内,下垫面是玉米。观测点的经纬度是100.37852E, 38.86074N,海拔1550.73m。空气温度、相对湿度传感器架设在5m处,朝向正北;气压与翻斗式雨量计安装在2米、10米处;风速与风向传感器架设在10米,朝向正北;四分量辐射仪安装在6米处,朝向正南;两个红外表面温度传感器安装在6米处,朝向正南支臂朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、60cm和100 cm处,并距离气象塔2米的正南方; 土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,并距离气象塔2米的正南方;土壤热流板(3块)依次埋设在地下6 cm处,其中两块埋设在棵间,一块埋设在植株下面。观测项目有:空气温湿观测(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)(单位:瓦/平方米)、 多层土壤水分(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)确保每天1440个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2012-6-10 10:30;(5)命名规则为:AMS+站点编号 。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-09

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(5号点自动气象站)

该数据集包含了2012年6月4日至9月18日的通量观测矩阵中5号点的自动气象站观测数据。站点位于甘肃省张掖市盈科灌区农田内,下垫面是玉米。观测点的经纬度是100.35068E, 38.87574N,海拔1567.65m。空气温度、相对湿度传感器架设在5m处,朝向正北;气压与翻斗式雨量计安装在塔顶10米处;风速(10米)与风向传感器(10米)架设在10米,朝向正北;四分量辐射仪安装在6米处,朝向正南;两个红外表面温度传感器安装在4米处,支臂朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、60cm和100 cm处,并距离气象塔2米的正南方; 土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,并距离气象塔2米的正南方;土壤热流板(3块)依次埋设在地下6 cm处,其中两块埋设在棵间,一块埋设在植株下面。观测项目有:空气温湿观测(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)(单位:瓦/平方米)、 多层土壤水分(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)确保每天1440个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2012-6-10 10:30;(5)命名规则为:AMS+站点编号 。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-09

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(5号点涡动相关仪)(2012年6月-9月)

该数据集包含了2012年6月3日至9月18日的通量观测矩阵中5号点的涡动相关通量观测数据。站点位于甘肃省张掖市盈科灌区农田内,下垫面是玉米。观测点的经纬度是100.35068E, 38.87574N,海拔1567.65m。涡动相关仪架高3m,采样频率是10Hz,超声朝向是正北向,超声风速仪与CO2/H2O分析仪之间的距离是17cm。 涡动相关仪的原始观测数据为10Hz,发布的数据是采用Edire软件后处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Edire软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据每30min内缺失率大于3%的数据;(4)剔除夜间弱湍流的观测数据(u*小于0.1m/s)。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),稳定度Z/L(无量纲),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE),二氧化碳通量的质量标识(QA_Fc)。感热、潜热、二氧化碳通量的质量标识分为三级(质量标识0:(Δst <30, ITC<30); 1: (Δst <100, ITC<100); 其余为2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-09

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集 (16号点涡动相关仪) (2012)

该数据集包含了2012年6月6日至9月17日的通量观测矩阵中16号点的涡动相关通量观测数据。站点位于甘肃省张掖市大满灌区农田内,下垫面是玉米。观测点的经纬度是100.36411E, 38.84931N,海拔1564.31m。涡动相关仪架高4.9m,采样频率是10Hz,超声朝向是正北向,超声风速仪与CO2/H2O分析仪之间的距离是20cm。 涡动相关仪的原始观测数据为10Hz,发布的数据是采用Edire软件后处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,角度订正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Edire软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据每30min内缺失率大于3%的数据;(4)剔除夜间弱湍流的观测数据(u*小于0.1m/s)。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),稳定度Z/L(无量纲),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE,二氧化碳通量的质量标识QA_Fc。感热、潜热、二氧化碳通量的质量标识分为三级(质量标识0:(Δst <30, ITC<30); 1: (Δst <100, ITC<100); 其余为2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-09

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(9号点涡动相关仪)

该数据集包含了2012年6月4日至9月17日的通量观测矩阵中9号点的涡动相关通量观测数据。站点位于甘肃省张掖市盈科灌区农田内,下垫面是玉米。观测点的经纬度是100.38546E, 38.87239N,海拔1543.34m。涡动相关仪架高3.9m,采样频率是10Hz,超声朝向是正北向,超声风速仪与CO2/H2O分析仪之间的距离是20cm。 涡动相关仪的原始观测数据为10Hz,发布的数据是采用Edire软件后处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,角度订正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Edire软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据每30min内缺失率大于3%的数据;(4)剔除夜间弱湍流的观测数据(u*小于0.1m/s)。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),稳定度Z/L(无量纲),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE,二氧化碳通量的质量标识QA_Fc。感热、潜热、二氧化碳通量的质量标识分为三级(质量标识0:(Δst <30, ITC<30); 1: (Δst <100, ITC<100); 其余为2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-08

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(宇宙射线土壤水分)

该数据集包含了2012年6月1日至9月20日的通量观测矩阵中两台宇宙射线仪器(crs_a和crs_b)的观测数据。站点位于甘肃省张掖市盈科灌区农田内,下垫面是玉米地。crs_a观测点的经纬度是100.36975E, 38.85385N,海拔1557.16m,crs_b观测点的经纬度是100.37225E, 38.85557N,海拔1557.16m,仪器探头底部距地面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。具体见数据引用文献. 在此,根据两台仪器源区内的Soilnet土壤水分数据对仪器进行率定,并根据实际情况之间建立土壤体积含水量θv和快中子之间的关系,即将公式(2)中的θm换作θv处理。分别选取干湿状况差异比较明显的6月26日-27日和7月16日-17日四天的数据,其中6月26日-27日率定数据显示土壤水分较小,因此选取4厘米、10厘米和20厘米的三个值平均值作为率定数据,其变化范围为22%-30%,而7月16日-17日率定数据显示土壤水分较大,因此选取4厘米、10厘米的两个值平均值作为率定数据,其变化范围为28%-39%,最后平均得到crs_a和crs_b的值N0分别为3252、3597。 4) 土壤水分计算 计算得到每小时的土壤含水量数据。具体计算公式见数据引用文献. 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Zhu et al.(2015)。

2016-07-08

黑河生态水文遥感试验:非均匀下垫面地表蒸散发的多尺度观测试验-通量观测矩阵数据集(3号点涡动相关仪)

该数据集包含了2012年6月3日至9月18日的通量观测矩阵中3号点的涡动相关通量观测数据。站点位于甘肃省张掖市盈科灌区农田内,下垫面是玉米。观测点的经纬度是100.37634E, 38.89053N,海拔1543.05m。涡动相关仪架高3.8m,采样频率是10Hz,超声朝向是正北向,超声风速仪与CO2/H2O分析仪之间的距离是20cm。 涡动相关仪的原始观测数据为10HZ,发布的数据是采用Edire软件后处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,角度订正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Edire软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据每30min内缺失率大于3%的数据;(4)剔除夜间弱湍流的观测数据(u*小于0.1m/s)。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),稳定度Z/L(无量纲),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE,二氧化碳通量的质量标识QA_Fc。感热、潜热、二氧化碳通量的质量标识分为三级(质量标识0:(Δst <30, ITC<30); 1: (Δst <100, ITC<100); 其余为2)。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 多尺度观测试验或站点信息请参考Liu et al. (2016),观测数据处理请参考Xu et al.(2013)。

2016-07-08