该数据集包含了2016年1月1日至2016年12月31日的黑河水文气象观测网上游垭口站涡动相关仪观测数据。站点位于青海省祁连县，下垫面是高寒草甸。观测点的经纬度是100.2421, 38.0142N，海拔4148 m。涡动相关仪的架高3.2m，采样频率是10Hz，超声朝向是正北向，超声风速温度仪（CSAT3）与CO2/H2O分析仪（Li7500A）之间的距离是15cm。 涡动相关仪的原始观测数据为10Hz，发布的数据是采用Eddypro软件处理的30分钟数据，其处理的主要步骤包括：野点值剔除，延迟时间校正，坐标旋转（二次坐标旋转），频率响应修正，超声虚温修正和密度（WPL）修正等。同时对各通量值进行质量评价，主要是大气平稳性（Δst）和湍流相似性特征（ITC）的检验。对Eddypro软件输出的30min通量值也进行了筛选：（1）剔除仪器出错时的数据；（2）剔除降水前后1h的数据；（3）剔除10Hz原始数据中每30min内缺失率大于10%的数据；（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），奥布霍夫长度L（m），感热通量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. (2018)，观测数据处理请参考Liu et al. (2011)。

该数据集包含了黑河水文气象观测网中游大满站的大孔径闪烁仪通量观测数据。中游大满站分别架设了BLS450和BLS900型号的大孔径闪烁仪，北塔为BLS900的接收端和BLS450的发射端，南塔为BLS900的发射端和BLS450接收端。观测时间为2016年1月1日至2016年12月31日。站点位于甘肃省张掖市大满灌区内，下垫面是玉米、果园和大棚，以玉米为主。北塔的经纬度是100.3785E，38.8607N，南塔的经纬度是100.3685E，38.8468N，海拔高度约1556m。大孔径闪烁仪的有效高度22.45m，光径长度是1854m，采样频率是1min。 大孔径闪烁仪原始观测数据为1min，发布的数据为经过处理与质量控制后的数据，其中感热通量主要是结合自动气象站观测数据，基于莫宁-奥布霍夫相似理论通过迭代计算得到，主要的质量控制步骤包括：（1）剔除Cn2达到饱和的数据（Cn2>1.43E-13）；（2）剔除解调信号强度较弱的数据（Average X Intensity<1000）；（3）剔除降水时刻的数据；（4）剔除稳定条件下的弱湍流的数据（u*小于0.1m/s）。在迭代计算过程中，选取Thiermann and Grassl(1992)的稳定度普适函数，详细介绍请参考Liu et al(2011, 2013)。 关于发布数据的几点说明：（1）中游LAS数据以BLS900为主，缺失时刻由BLS450观测补充，两者都缺失则以-6999标记。（2）数据表头：Date/Time ：日期/时间（格式：yyyy/m/d h:mm），Cn2 ：空气折射指数结构参数（单位：m-2/3），H_LAS ：感热通量（单位：W/m2）。数据时间的含义，如0:30代表0:00-0:30的平均；数据以*.xls格式存储。 水文气象网或站点信息请参考Li et al. (2013)，观测数据处理请参考Liu et al. (2011)。

该数据集包含了2016年1月1日至2016年12月31日黑河水文气象观测网上游阿柔超级站气象要素梯度观测系统数据。站点位于青海省祁连县阿柔乡草达坂村，下垫面是高寒草地。观测点的经纬度是100.4643E,38.0473N，海拔3033m。空气温度、相对湿度、风速传感器分别架设在1m、2m、5m、10m、15m、25m处，共6层，朝向正北；风向传感器架设在10m处，朝向正北；气压计安装在2m处；翻斗式雨量计安装在阿柔超级站40m观测塔上；四分量辐射仪安装在5m处，朝向正南；两个红外温度计安装在5m处，朝向正南，探头朝向是垂直向下；光合有效辐射仪安装在5m处，朝向正南，探头朝向是垂直向上；土壤部分传感器埋设在塔体正南方向2m处，其中土壤热流板（自校正式）（3块）均埋设在地下6cm处；平均土壤温度传感器TCAV埋设在地下2cm、4cm处；土壤温度探头埋设在地表0cm和地下2cm、4cm、6cm、10cm、15cm、20cm、30cm、40cm、60cm、80cm、120cm、160cm、200cm、240cm、280cm、320cm处，其中4cm和10cm这两层有三个重复；土壤水分传感器分别埋设在地下2cm、4cm、6cm、10cm、15cm、20cm、30cm、40cm、60cm、80cm、120cm、160cm、200cm、240cm、280cm、320cm处，其中4cm和10cm这两层有三个重复。 观测项目有：风速（WS_1m、WS_2m、WS_5m、WS_10m、WS_15m、WS_25m）(单位：米/秒)、风向（WD_10m）(单位：度)、空气温湿度（Ta_1m、Ta_2m、Ta_5m、Ta_10m、Ta_15m、Ta_25m和RH_1m、RH_2m、RH_5m、RH_10m、RH_15m、RH_25m）(单位：摄氏度、百分比)、气压（Press）(单位：百帕)、降水量（Rain）(单位：毫米)、四分量辐射（DR、UR、DLR_Cor、ULR_Cor、Rn）(单位：瓦/平方米)、地表辐射温度（IRT_1、IRT_2）(单位：摄氏度)、光合有效辐射（PAR）(单位：微摩尔/平方米秒)、平均土壤温度（TCAV）(单位：摄氏度)、土壤热通量（Gs_1、Gs_2、Gs_3）(单位：瓦/平方米)、土壤水分（Ms_2cm、Ms_4cm_1、Ms_4cm_2、Ms_4cm_3、Ms_6cm、Ms_10cm_1、Ms_10cm_2、Ms_10cm_3、Ms_15cm、Ms_20cm、Ms_30cm、Ms_40cm、Ms_60cm、Ms_80cm、Ms_120cm、Ms_160cm Ms_200cm、Ms_240cm、Ms_280cm、Ms_320cm）(单位：体积含水量，百分比)、土壤温度（Ts_0cm、Ts_2cm、Ts_4cm_1、Ts_4cm_2、Ts_4cm_3、Ts_6cm、Ts_10cm_1、Ts_10cm_2、Ts_10cm_3、Ts_15cm、Ts_20cm、Ts_30cm、Ts_40cm、Ts_60cm、Ts_80cm、Ts_120cm、Ts_160cm Ts_200cm、Ts_240cm、Ts_280cm、Ts_320cm）(单位：摄氏度)。 观测数据的处理与质量控制：（1）确保每天144个数据（每10min），若出现数据的缺失，则由-6999标示；土壤热通量G1在2016.12.8-12.16之间传感器问题，数据缺失；（2）剔除有重复记录的时刻；（3）删除了明显超出物理意义或超出仪器量程的数据；（4）数据中以红字标示的部分为有疑问的数据；（5）日期和时间的格式统一，并且日期、时间在同一列。如，时间为：2016-6-10 10:30；（6）命名规则为：AWS+站点名称。 水文气象网或站点信息请参考Liu et al. (2018)，观测数据处理请参考Liu et al. (2011)。

该数据集包含了2016年1月1日至2016年12月31日的黑河水文气象观测网下游荒漠站的涡动相关仪观测数据。站点位于内蒙古额济纳旗，下垫面是荒漠。观测点的经纬度是100.9872E, 42.1135N，海拔1054m。涡动相关仪的架高4.7m，采样频率是10Hz，超声朝向是正北向，超声风速温度仪（CSAT3）与CO2/H2O分析仪（Li7500）之间的距离是15cm。 涡动相关仪的原始观测数据为10Hz，发布的数据是采用Eddypro软件处理的30分钟数据，其处理的主要步骤包括：野点值剔除，延迟时间校正，坐标旋转（二次坐标旋转），频率响应修正，超声虚温修正和密度（WPL）修正等。同时对各通量值进行质量评价，主要是大气平稳性（Δst）和湍流相似性特征（ITC）的检验。对Eddypro软件输出的30min通量值也进行了筛选：（1）剔除仪器出错时的数据；（2）剔除降水前后1h的数据；（3）剔除10Hz原始数据中每30min内缺失率大于10%的数据；（4）剔除夜间弱湍流的观测数据（u*小于0.1m/s）。观测数据的平均周期为30分钟，一天48个数据，缺失数据标记为-6999。涡动相关仪的Li7500在4月21-22日标定，数据缺失。因仪器漂移等原因引起的可疑数据用红色字体标识。 发布的观测数据包括：日期/时间Date/Time，风向Wdir（°），水平风速Wnd（m/s），侧向风速标准差Std_Uy（m/s），超声虚温Tv（℃），水汽密度H2O（g/m3），二氧化碳浓度CO2（mg/m3），摩擦速度Ustar（m/s），奥布霍夫长度L（m），感热通量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格式存储。 水文气象网或站点信息请参考Li et al. (2013)，观测数据处理请参考Liu et al. (2011)。

该数据集包含了2016年3月18日至2016年12月31日的黑河水文气象观测网中游张掖湿地站的涡动相关仪观测数据。站点位于甘肃省张掖市，下垫面是湿地。观测点的经纬度是100.44640E, 38.97514N，海拔1460.00m。涡动相关仪的架高5.2m，采样频率是10Hz，超声朝向是正北向，超声风速温度仪（Gill）与CO2/H2O分析仪(Li7500A)之间的距离是25cm。 涡动相关仪的原始观测数据为10Hz，发布的数据是采用Eddypro软件处理的30分钟数据，其处理的主要步骤包括：野点值剔除，延迟时间校正，角度订正，坐标旋转（二次坐标旋转），频率响应修正，超声虚温修正和密度（WPL）修正等。同时对各通量值进行质量评价，主要是大气平稳性（Δst）和湍流相似性特征（ITC）的检验。对Eddypro软件输出的30min通量值也进行了筛选：（1）剔除仪器出错时的数据；（2）剔除降水前后1h的数据；（3）剔除10Hz原始数据中每30min内缺失率大于10%的数据；（4）剔除夜间弱湍流的观测数据（u*小于0.1m/s）。观测数据的平均周期为30分钟，一天48个数据，缺失数据标记为-6999。因仪器漂移等原因引起的可疑数据用红色字体标识，其中2016年4月-6月涡动系统的三维超声断续出现问题，数据出现一些缺失。 发布的观测数据包括：日期/时间Date/Time，风向Wdir（°），水平风速Wnd（m/s），侧向风速标准差Std_Uy（m/s），超声虚温Tv（℃），水汽密度H2O（g/m3），二氧化碳浓度CO2（mg/m3），摩擦速度Ustar（m/s），奥布霍夫长度L（m），感热通量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格式存储。 水文气象网或站点信息请参考Li et al. (2013)，观测数据处理请参考Liu et al. (2011)。

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

该数据集包含了2016年1月1日至2016年12月31日黑河水文气象观测网中游黑河遥感站气象要素观测数据。站点位于甘肃省张掖市党寨镇东侧，下垫面是人工草地。观测点的经纬度是100.4756E, 38.8270N，海拔1560m。空气温度湿度传感器架设在1.5m处，朝向正北；气压计在防水箱内；翻斗式雨量计安装在0.7 m处；风速风向传感器架设在10m处，朝向正北；四分量辐射仪安装高度为1.5m，朝向正南；两个红外温度计安装高度为1.5m，朝向正南，探头朝向是垂直向下；土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、80cm、120cm、160cm处；土壤水分探头埋设在2cm、4cm、10cm、20cm、40cm、80cm、120cm、160cm处；平均土壤温度探头埋设在2cm和4cm；土壤热流板（3块）依次埋设在地下6cm处；两个光合有效辐射仪分别架设在冠层上方1.5m（探头垂直向上和向下方向各一个），朝向正南。 观测项目有：空气温湿度（Ta_1.5m、RH_1.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_2cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_80cm、Ts_120cm、Ts_160cm）(单位：摄氏度)、土壤水分（Ms_0cm、Ms_2cm、Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_80cm、Ms_120cm、Ms_160cm）(单位：%)、向上与向下光合有效辐射（PAR_U_up、PAR_U_down）(单位：微摩尔/平方米秒)、平均土壤温度（TCAV）（单位：摄氏度）。 观测数据的处理与质量控制：（1）确保每天144个数据（每10min），若出现数据的缺失，则由-6999标示；2016.1.01-1.29由于采集器问题，很多观测要素出现较多的错误值；（2）剔除有重复记录的时刻；（3）删除了明显超出物理意义或超出仪器量程的数据；（4）数据中以红字标示的部分为有疑问的数据；（5）日期和时间的格式统一，并且日期、时间在同一列。如，时间为：2016-6-10 10:30；（6）命名规则为：AWS+站点名称。 水文气象网或站点信息请参考Liu et al. (2018)，观测数据处理请参考Liu et al. (2011)。

样地调查数据为，于2013年8月份，在天涝池流域设置森林样地30块，样地规格为10 m×20 m，样地长边与山坡走向平行，其中青海云杉林26块，祁连圆柏林2块，云杉圆柏混交林2块，在样地内，采用围尺测量每株树木的胸径（树干1.3 m高度处的直径），采用手持超声波测高器测量每株树木的树高、枝下高（树冠下端第一活枝的高度），采用皮尺测量南北方向和东西方向冠幅，利用差分GPS对样地进行定位.以样地碳储量数据为优化控制条件，以Kriging插值得到的生物量空间分布图驱动场，采用HASM算法模拟天涝池森林生物量空间分布图，模拟结果符合研究区的植被分布规律，获得较好的效果。

EC-EARTH-Heihe是采用EC-EARTH全球模式输出结果作为驱动场模拟1980-2005年和RCP4.5情景下2006-2080黑河流域6小时数据。空间范围：模拟区域的网格中心位于（40.30N，99.50E）, 水平分辨率为3 km，模式的模拟网格点数为161（经向）X 201（纬向）。 投影方式：LAMBERT正形投影，两个标准纬度为30N 和60N。 时间范围：1980年1月1日－2010年12月31日，时间间隔为6小时。 文件内容说明：采用grads无格式月存储。除最高、最低温度为日尺度以外，其他变量都是6小时数据。 可以采用MATLAB进行读取,可见tmax_erain_xiong_heihe.m文件说明。 黑河流域数据说明： 1) Anemometer west wind(m/s)
 简称usurf 2) Anemometer south wind(m/s) 简称vsurf 
3) Anemometer temperature (deg K) 简称tsurf 
4)maximal temperature (deg K) 简称tmax 5) minimal temperature (deg K) 简称tmin 
6) Anemom specific humidity (g/kg)
 简称qsurf 7) Accumulated precipitation (mm/hr) 简称precip 
8) Accumulated evaporation (mm/hr) 简称evap 9) Accumulated sensible heat (watts/m**2/hr) 简称sensible 
10) Accumulated net infrared radiation (watts/m**2/hr) 简称netrad 
文件名定义: 简称-ec-earth-6hour.年月 例如：precip-ec-earth-6hour.198001 为1980年1月6小时降水资料 （1） 采用EC-EARTH全球气候模式1980-2005年驱动模拟的历史6小时一次数据。 （2） 2006-2080年为EC-EARTH全球气候模式在RCP 4.5情景下生产黑河流域6小时一次数据资料。

黑河流域1km/5天合成叶面积指数（LAI）数据集提供了2010-2014年的5天LAI合成结果，该数据利用Terra/MODIS、Aqua/MODIS、以及国产卫星FY3A/MERSI和FY3B/MERSI传感器数据构建空间分辨率1km、时间分辨率5天的多源遥感数据集。多源遥感数据集可在有限时间内提供比单一传感器更多的角度和更多次的观测，但是，由于传感器的在轨运行时间及性能差异，多源数据集的观测质量参差不齐。因此，为更有效的利用多源数据集，算法首先对多源数据集进行了质量分级，根据观测合理性分为一级数据、二级数据、三级数据。三级数据为受薄云污染的观测，不用于计算。质量评估及分级的目的是为LAI反演时最优数据集的选择及反演算法流程设计提供依据。叶面积指数产品反演算法设计为区分山地平地、区分植被类型使用不同模型的神经网络法反演。基于全球DEM图和地表分类图，针对草地和农作物等连续植被采用PROSAIL模型，针对森林和山地植被采用坡面GOST模型。利用黑河上游森林和中游绿洲的地面实测数据生成的参考图，并将对应的高分辨率参考图升尺度到1km分辨率，与LAI产品进行比较，产品在农田和森林区域与参考值间均具有良好的相关性，总体精度基本满足GCOS规定的误差不超过 (0.5, 20%)的精度阈值。将本产品与MODIS、GEOV1和GLASS等LAI产品进行交叉对比，相比较参考值而言，本LAI产品精度优于同类产品。总之，黑河流域1km/5天合成LAI数据集综合利用多源遥感数据以提高LAI参数产品的估算精度、时间分辨率等，更好的服务于遥感数据产品的应用。

该套数据是新发展的陆面生态水文模式CLM_LTF的模拟结果。该模式在NCAR发展的陆面过程模式CLM4.5之上，耦合了地下水侧向流动模块，并且考虑了人类灌溉作用。 模式运行时间是1981-2013年，空间分辨率为30弧秒(0.0083度)，时间步长为1800秒，模拟范围为黑河流域。1981-2012年大气强迫使用的是由中国科学院青藏高原研究所青藏高原多圈层数据同化与模拟中心开发的中国区域高时空分辨率地面气象要素驱动数据集，2013年大气强迫使用的是国家气象信息中心制作的高分辨率的风压湿温降水辐射数据集。地表覆盖数据为MICLCover黑河流域1公里土地覆盖格网数据集，灌溉数据见寒区旱区科学数据中心“黑河流域1981-2013年30弧秒分辨率月尺度地表水及地下水灌溉量数据集”。模式输出为月平均数据。 文件说明如下： 地下水埋深数据：Heihe_ZWT.nc 2cm土壤湿度数据：Heihe_H2OSOI_2CM.nc 100cm 土壤湿度数据：Heihe_H2OSOI_100CM.nc 蒸散发数据：Heihe_evaptanspiration.nc 数据为netcdf格式。有3个维度，依次为month, lat, lon. 其中month为月份，数值为0-395，代表1981-2013年逐个月份，lat为网格纬度信息，lon为网格经度信息。 数据储存在data变量中，地下水埋深数据单位为m, 土壤湿度数据单位为m^3/m^3, 蒸散发数据单位为mm/month

2008年全国遥感年平均地表温度和冻结指数是冉有华等（2015）基于MODIS Aqua/Terra逐日四次的5公里瞬时地表温度数据产品，发展了新的年平均地表温度和冻结指数估计方法，该方法利用上下午LST观测的平均获取日平均地表温度，方法的核心是如何恢复LST产品的缺失数据，该方法有两个特点：（1）将遥感观测到的日地表温度变幅进行了空间插值，利用插值获取的空间连续的日地表温度变幅，使一天只有一次的卫星观测数据得到应用；（2）利用了一个新的缺失数据时间序列滤波方法，即基于离散余弦变换的惩罚最小二乘回归方法。 验证表明，年平均地表温度与冻结指数的精度只与原始MODIS LST的精度有关，即保持了MODIS LST产品的精度。可用于冻土制图及相关资源环境应用。