張玉翠，姜寒冰,3，張傳偉,3，沈彥俊*

1. 中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心/中國科學(xué)院農(nóng)業(yè)水資源重點(diǎn)實(shí)驗(yàn)室/河北省節(jié)水農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，石家莊 050022

2. 中國科學(xué)院欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗(yàn)站，石家莊 051430

3. 中國科學(xué)院大學(xué)，北京 100049

數(shù)據(jù)庫（集）基本信息簡(jiǎn)介

數(shù)據(jù)庫（集）組成 本數(shù)據(jù)集由欒城試驗(yàn)站通量觀測(cè)數(shù)據(jù)組成，含1個(gè)Excel文件，包括5個(gè)表單，分別為農(nóng)田蒸散量、降水量、凈輻射量、感熱通量、二氧化碳凈交換量數(shù)據(jù)記錄。每個(gè)表單分別包括2192條日尺度數(shù)據(jù)和6條生育季尺度數(shù)據(jù)。

引 言

陸地生態(tài)系統(tǒng)碳水循環(huán)是地球化學(xué)循環(huán)的核心過程。研究陸地生態(tài)系統(tǒng)碳水循環(huán)過程對(duì)于人類調(diào)節(jié)生物圈–水圈–大氣圈相互作用，對(duì)維持全球生態(tài)系統(tǒng)的物質(zhì)運(yùn)輸與能量循環(huán)具有重要的指導(dǎo)意義，也是人為適應(yīng)和干預(yù)全球變暖，緩解水資源供需矛盾的戰(zhàn)略需求[1-2]。經(jīng)過長(zhǎng)期的理論發(fā)展和技術(shù)改革，基于微氣象學(xué)原理的渦度相關(guān)技術(shù)已成為陸地生態(tài)系統(tǒng)與大氣間的CO2和水熱通量交換的標(biāo)準(zhǔn)觀測(cè)手段之一[2]。渦度相關(guān)系統(tǒng)的觀測(cè)原理就是將空氣流動(dòng)視為無數(shù)個(gè)渦流，通過物理量的脈動(dòng)與垂直風(fēng)速脈動(dòng)的協(xié)方差計(jì)算湍流通量[3]。目前，渦度相關(guān)系統(tǒng)已被廣泛應(yīng)用于陸地生態(tài)系統(tǒng)碳水循環(huán)過程研究，長(zhǎng)期連續(xù)的高頻率聯(lián)網(wǎng)觀測(cè)積累的通量數(shù)據(jù)可為碳水循環(huán)機(jī)理、碳源碳匯的時(shí)空分布特征以及生物地球化學(xué)循環(huán)模型的優(yōu)化等相關(guān)研究提供有力的數(shù)據(jù)支撐[4-5]。

2001 年，中國陸地生態(tài)系統(tǒng)通量觀測(cè)研究網(wǎng)絡(luò)（ChinaFLUX）依托中國生態(tài)系統(tǒng)研究網(wǎng)絡(luò)（CERN）成立，經(jīng)過十幾年的發(fā)展，ChinaFLUX 目前已建設(shè)完成100 余個(gè)通量觀測(cè)站，覆蓋森林、草地、農(nóng)田等多種類型的陸地生態(tài)系統(tǒng)，形成了全國尺度的碳水通量觀測(cè)研究體系，實(shí)現(xiàn)了我國通量觀測(cè)事業(yè)從起步到國際前沿水平的跨越式發(fā)展[6-8]。作為ChinaFLUX 重要成員之一的中國科學(xué)院欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗(yàn)站（以下簡(jiǎn)稱“欒城試驗(yàn)站”）始建于1981 年，隸屬于中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心（簡(jiǎn)稱“中科院農(nóng)業(yè)資源研究中心”）。1989 年加入CERN，成為第一批基本站。2005 年成為國家生態(tài)系統(tǒng)觀測(cè)研究網(wǎng)絡(luò)（CNERN）臺(tái)站，同年成為ChinaFLUX 農(nóng)田生態(tài)系統(tǒng)的重要觀測(cè)站點(diǎn)之一。

欒城試驗(yàn)站渦度相關(guān)系統(tǒng)長(zhǎng)期連續(xù)的碳水通量聯(lián)網(wǎng)觀測(cè)開始于2007 年，本數(shù)據(jù)集包含了欒城試驗(yàn)站2007–2013 年的碳水通量觀測(cè)數(shù)據(jù)，選取了其中較為重要的組成部分，包括農(nóng)田蒸散量、降雨量、凈輻射量、感熱通量和二氧化碳凈交換量，整理形成了日尺度和生育季節(jié)尺度數(shù)據(jù)產(chǎn)品，并以此數(shù)據(jù)論文形式對(duì)數(shù)據(jù)集進(jìn)行系統(tǒng)性描述，以便于讀者溯源定位，理解數(shù)據(jù)集內(nèi)涵，推動(dòng)數(shù)據(jù)開放共享和規(guī)范化使用。

1 數(shù)據(jù)采集和處理方法

1.1 數(shù)據(jù)來源及采集方法

基于ChinaFLUX 的頂層設(shè)計(jì)，經(jīng)過觀測(cè)塔選址、觀測(cè)儀器選型和野外觀測(cè)系統(tǒng)安裝與調(diào)試等技術(shù)工作，中科院欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗(yàn)站于2007 年11 月6 日正式開始長(zhǎng)期連續(xù)的碳水通量聯(lián)網(wǎng)觀測(cè)。本數(shù)據(jù)集為2007 年10 月1 日至2013 年9 月30 日的農(nóng)田碳水通量數(shù)據(jù)，其中2007 年10 月1日至11 月5 日的數(shù)據(jù)由歷年相應(yīng)通量項(xiàng)平均值并結(jié)合欒城試驗(yàn)站氣象數(shù)據(jù)、土壤含水量和參考蒸散等插補(bǔ)得到。

欒城試驗(yàn)站全稱為中國科學(xué)院欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗(yàn)站（臺(tái)站代碼：LCA），該站地處太行山山前平原中部（37°53′N、114°41′E，海拔50.1 m），暖溫帶半濕潤(rùn)季風(fēng)氣候，年平均氣溫12.9℃，多年平均降雨量490 mm（2007–2013 年），無霜期200 天，全年平均日照時(shí)數(shù)2608 小時(shí)，≥0℃的積溫為4710℃，雨熱同期，適宜農(nóng)業(yè)發(fā)展。土壤以潮褐土為主，質(zhì)地為壤土，隨土層由淺到深依次為粉砂壤、壤土和粘壤。農(nóng)業(yè)生產(chǎn)以華北地區(qū)種植面積最為廣泛的冬小麥–夏玉米一年兩熟輪作制農(nóng)田為主。冬小麥10 月初播種，次年6 月中旬收獲；夏玉米6 月中旬播種，9 月底收獲。其中冬小麥冠層高度約為1.0 m，夏玉米冠層高度在2.0 m 左右。

結(jié)合上述植被類型、冠層高度以及通量觀測(cè)下墊面要求，欒城站架設(shè)了通量觀測(cè)塔，并將碳水分析儀安裝高度設(shè)置為3.5 m。渦度相關(guān)系統(tǒng)主要由開路式CO2/H2O 紅外分析儀和三維超聲風(fēng)速計(jì)組成，同時(shí)輔有輻射分量和常規(guī)氣象要素（降雨和空氣溫濕度等）的同步觀測(cè)，觀測(cè)設(shè)備型號(hào)及制造商等信息見表1。數(shù)據(jù)觀測(cè)和存儲(chǔ)除人工觀測(cè)降水量外由系統(tǒng)自動(dòng)完成，通常情況下每0.5–1 個(gè)月定期維護(hù)一次，每?jī)赡晷?zhǔn)一次。

表1 各觀測(cè)系統(tǒng)主要測(cè)定要素、關(guān)鍵傳感器型號(hào)、制造商信息及安裝高度

通量觀測(cè)塔下墊面為充分灌溉管理下的典型冬小麥–夏玉米一年兩熟輪作制農(nóng)田生態(tài)系統(tǒng)，下墊面均勻且地勢(shì)平坦，風(fēng)浪區(qū)（＞350 m）足夠大。觀測(cè)系統(tǒng)利用CR5000 型數(shù)據(jù)采集器以10 Hz 的頻率進(jìn)行原始數(shù)據(jù)采樣和存儲(chǔ)，連續(xù)自動(dòng)監(jiān)測(cè)生態(tài)系統(tǒng)的碳水通量，每30 min 輸出一組平均值。常規(guī)氣象要素的觀測(cè)記錄輸出頻度同樣為30 min，由相應(yīng)的數(shù)據(jù)采集系統(tǒng)自動(dòng)完成數(shù)據(jù)獲取和存儲(chǔ)。通過U 盤下載原始數(shù)據(jù)并存儲(chǔ)備份，然后根據(jù)ChinaFLUX 數(shù)據(jù)管理要求進(jìn)行后續(xù)的質(zhì)量控制和標(biāo)準(zhǔn)化處理。

1.2 數(shù)據(jù)處理方法

通量數(shù)據(jù)的原始觀測(cè)數(shù)據(jù)統(tǒng)一遵循ChinaFLUX 技術(shù)體系進(jìn)行標(biāo)準(zhǔn)化的質(zhì)量控制和處理，處理流程如圖1。為規(guī)范和便于讀者對(duì)數(shù)據(jù)的使用，本數(shù)據(jù)集將10 Hz 的原始數(shù)據(jù)處理計(jì)算形成了日尺度和生育季節(jié)尺度數(shù)據(jù)產(chǎn)品[9]。

圖1 通量數(shù)據(jù)質(zhì)量控制與處理流程

數(shù)據(jù)質(zhì)量控制過程主要包括原始數(shù)據(jù)分析、超聲虛溫校正、坐標(biāo)軸旋轉(zhuǎn)、WPL 校正、穩(wěn)態(tài)測(cè)試、湍流積分特性、野點(diǎn)剔除、能量閉合度評(píng)價(jià)以及夜間摩擦風(fēng)速閾值篩選[10-15]。野點(diǎn)剔除指原始數(shù)據(jù)中的異常值剔除，剔除方法包括絕對(duì)值法和相對(duì)值法[16]。絕對(duì)值法即依據(jù)觀測(cè)指標(biāo)理論上可能出現(xiàn)的數(shù)值范圍設(shè)置上下限，將明顯超出正常范圍的數(shù)值直接剔除；相對(duì)值法即計(jì)算數(shù)據(jù)的平均值和方差，以4 倍方差作為檢驗(yàn)標(biāo)準(zhǔn)，剔除偏離過大的數(shù)據(jù)。坐標(biāo)軸旋轉(zhuǎn)也稱為傾斜校正，其目的是保證超聲風(fēng)速儀與地表平行，減少由于儀器傾斜導(dǎo)致觀測(cè)的垂直風(fēng)速中包含水平風(fēng)速分量而造成的通量數(shù)據(jù)誤差。本研究采用中傾斜校正，主要包括三次坐標(biāo)軸旋轉(zhuǎn)（Tiple Coordinate Rotation，TR）、WPL（Webb-Pearman-Leuning）校正、頻率響應(yīng)校正。其中三次坐標(biāo)軸旋轉(zhuǎn)可使x 軸平行于平均風(fēng)速方向，側(cè)向平均風(fēng)速與垂直平均風(fēng)速及兩者的協(xié)方差為零；WPL 校正針對(duì)由于水熱通量輸送引起氣體密度變化而導(dǎo)致的潛熱通量和感熱通量觀測(cè)誤差；頻率響應(yīng)校正針對(duì)儀器的系統(tǒng)誤差，即由于儀器對(duì)高頻和低頻信號(hào)響應(yīng)存在的缺陷。根據(jù)平均值檢驗(yàn)法（Average Values Test，AVT）對(duì)夜間CO2通量數(shù)據(jù)進(jìn)行校正[17]。本數(shù)據(jù)集前期主要利用編寫的程序代碼（2007–2011 年）以及軟件EddyPro（LICOR 公司提供，2011–2013 年）對(duì)通量數(shù)據(jù)進(jìn)行了以上質(zhì)量控制和處理。2017 年后期利用ChinaFLUX統(tǒng)一程序代碼對(duì)原始數(shù)據(jù)又進(jìn)行了一次處理，由于所選計(jì)算方法基本一致，數(shù)據(jù)一致性達(dá)99.5%。針對(duì)缺失數(shù)據(jù)的插補(bǔ)，短時(shí)段（小于2 h）的缺失數(shù)據(jù)采用線性內(nèi)插法，較長(zhǎng)時(shí)段（大于2 h）的缺失數(shù)據(jù)利用平均晝夜變化法（Mean Diurnal Variation，MDV）進(jìn)行插補(bǔ)[18]。對(duì)于蒸散量（ET）連續(xù)多日（7 日以上）的數(shù)據(jù)缺失時(shí)段或MDV 插值后仍存在的異常值，我們根據(jù)蒸散發(fā)生的機(jī)制對(duì)數(shù)據(jù)處理方法進(jìn)行了改進(jìn)，利用土壤含水量與ET/ET0的定量關(guān)系進(jìn)行重新插補(bǔ)，這也是我們對(duì)數(shù)據(jù)插補(bǔ)方法的一點(diǎn)創(chuàng)新。具體方法如下：

①利用聯(lián)合國糧農(nóng)組織（Food and Agricultural Organization，F(xiàn)AO）推薦的Penman-Montieth 公式計(jì)算參考作物蒸散量ET0[19]：

式中：ET0為潛在蒸散發(fā)（mm），Rn為表層凈輻射（MJ/(m2·day)），G 為土壤熱通量（MJ/(m2·day)），Ta為日平均氣溫（oC），U 為2 m 高度處風(fēng)速（m/s），VPD 為飽和水汽壓差（kPa），Δ 為飽和水汽壓曲線斜率（kPa/oC），γ 為干濕表常數(shù)（kPa/oC），ET0計(jì)算用到的所有氣象數(shù)據(jù)均來自欒城試驗(yàn)站標(biāo)準(zhǔn)氣象場(chǎng)。

②篩選出連續(xù)7 日以上的數(shù)據(jù)缺失時(shí)段或異常值時(shí)段。

③篩選出存在缺值的日期及其前后2–3 天的數(shù)據(jù)。

④從篩選結(jié)果中選擇與缺值處土壤含水量相近的日期（日期相同，年份不同），假設(shè)每個(gè)生育期的作物系數(shù)在同時(shí)期不變，將各ET/ET0取平均值作為缺值日的ET/ET0。其中土壤含水量數(shù)據(jù)由IHII 型中子儀（英國Didcot 公司生產(chǎn)）和自動(dòng)傳感器（Stevens 公司生產(chǎn)的Hydra Probe II）共同測(cè)定，觀測(cè)深度為0–2 m。

⑤將平均所得的ET/ET0乘以當(dāng)天對(duì)應(yīng)的ET0，得到插補(bǔ)后的ET。

2 數(shù)據(jù)樣本描述

2.1 數(shù)據(jù)集命名規(guī)則及數(shù)據(jù)量

本數(shù)據(jù)集為欒城站2007 年10 月1 日至2013 年9 月30 日的農(nóng)田生態(tài)系統(tǒng)日尺度碳水通量觀測(cè)數(shù)據(jù)，共計(jì)1 個(gè)EXCEL 數(shù)據(jù)文件，大小為295 KB。文件命名為“欒城典型灌溉農(nóng)田生態(tài)系統(tǒng)2007–2013 年日通量數(shù)據(jù)集”，其中包括5 個(gè)表單，分別命名為“蒸散量”“降水量”“凈輻射量”“感熱通量”和“二氧化碳凈交換量”。每個(gè)表單分別包括2192 條日尺度數(shù)據(jù)和6 條生育季尺度數(shù)據(jù)，其中降雨量數(shù)據(jù)包括了同時(shí)期的人工觀測(cè)和自動(dòng)觀測(cè)數(shù)據(jù)。

2.2 數(shù)據(jù)文件示例

表2 為數(shù)據(jù)表頭示例，詳細(xì)介紹了數(shù)據(jù)項(xiàng)含義、數(shù)據(jù)類型、計(jì)量單位等信息。

表2 通量觀測(cè)數(shù)據(jù)表頭說明

數(shù)據(jù)項(xiàng) 數(shù)據(jù)類型 計(jì)量單位 數(shù)據(jù)項(xiàng)說明 示例 P-Year 數(shù)字 mm 生育年尺度降水量 597.5 Rn 數(shù)字 MJ/m2/day 日尺度凈輻射量 5.9 Rn-wheat 數(shù)字 MJ/m2 冬小麥季凈輻射量 1232.2 Rn-maize 數(shù)字 MJ/m2 夏玉米季凈輻射量 1003.0 Rn-Year 數(shù)字 MJ/m2/year 生育年尺度凈輻射量 2235.2 Hs 數(shù)字 MJ/m2/day 日尺度感熱通量 2.4 Hs-wheat 數(shù)字 MJ/m2 冬小麥季感熱通量 270.7 Hs-maize 數(shù)字 MJ/m2 夏玉米季感熱通量 202.6 Hs-Year 數(shù)字 MJ/m2/year 生育年尺度感熱通量 473.3 NEE 數(shù)字 gC/m2/day 日尺度CO2 凈交換量 0.9 NEE-wheat 數(shù)字 gC/m2 冬小麥季CO2 凈交換量 ?309.2 NEE-maize 數(shù)字 gC/m2 夏玉米季CO2 凈交換量 ?393.1 NEE-Year 數(shù)字 gC/m2/year 生育年尺度CO2 凈交換量 ?702.3

數(shù)據(jù)表頭說明：

（1）Year 表示生育年，10 月1 日至次年9 月30 日，即從小麥播種到玉米收獲的一個(gè)輪作周期。

（2）P（Precipitation）表示降水量（mm），表示從天空降落到地面上的液態(tài)或固態(tài)（經(jīng)融化后）水，未經(jīng)蒸發(fā)、滲透、流失，而在水平面上積聚的深度。

（3）ET（Evapotranspiration）表示農(nóng)田蒸散量（mm），指單位時(shí)間單位面積農(nóng)田土壤蒸發(fā)量和植物蒸騰量的總和。

（3）Rn（Net Radiation）表示凈輻射量（MJ/m2），是射入到地表的輻射能與地表射出的輻射能的差。

（4）Hs（Sensible Heat Flux）：表示感熱通量（MJ/m2），指由于溫度變化引起的大氣與下墊面之間發(fā)生的湍流形式的熱交換。

（5）NEE（Net ecosystem exchange of CO2）表示生態(tài)系統(tǒng)CO2凈交換量（gC/m2）。NEE 為正值表示農(nóng)田生態(tài)系統(tǒng)向大氣釋放CO2；NEE 為負(fù)值表示農(nóng)田生態(tài)系統(tǒng)從大氣吸收CO2。

3 數(shù)據(jù)質(zhì)量控制和評(píng)估

目前全球通量觀測(cè)研究領(lǐng)域普遍使用的數(shù)據(jù)質(zhì)量評(píng)價(jià)體系主要包括一貫性檢驗(yàn)和完全湍流假設(shè)的檢驗(yàn)。本數(shù)據(jù)集基于此對(duì)數(shù)據(jù)質(zhì)量進(jìn)行了系統(tǒng)評(píng)價(jià)。功率譜和協(xié)方差譜檢驗(yàn)結(jié)果表明，三維風(fēng)速、CO2濃度和H2O 濃度的功率譜變化模態(tài)符合慣性副區(qū)?2/3 斜率理論值，表明渦度相關(guān)系統(tǒng)設(shè)備響應(yīng)正常，不存在系統(tǒng)性的相移或失真[20]。結(jié)合湍流穩(wěn)定性檢驗(yàn)和積分統(tǒng)計(jì)特性檢驗(yàn)對(duì)數(shù)據(jù)進(jìn)行總體質(zhì)量評(píng)價(jià)和等級(jí)劃分，結(jié)果顯示觀測(cè)樣本符合湍流的方差相似性規(guī)律，數(shù)據(jù)質(zhì)量較高[21]。能量閉合度是通量數(shù)據(jù)質(zhì)量評(píng)價(jià)的重要參考標(biāo)準(zhǔn)之一，Wilson 等通過統(tǒng)計(jì)分析了全球通量觀測(cè)網(wǎng)絡(luò)（FLUXNET）22 個(gè)站點(diǎn)50 年通量數(shù)據(jù)，發(fā)現(xiàn)所有植被類型和氣候條件下的觀測(cè)均存在不同程度的能量不閉合現(xiàn)象，平均不閉合度達(dá)20%[22]。為反映生態(tài)系統(tǒng)通量實(shí)際變化情況，本數(shù)據(jù)集中通量觀測(cè)數(shù)據(jù)未進(jìn)行強(qiáng)制能量閉合處理。但是借助于改進(jìn)的ET（LE）插補(bǔ)方法，能量閉合度可達(dá)91%。其中小麥季觀測(cè)結(jié)果的能量閉合度可達(dá)95%，玉米季可達(dá)88%，遠(yuǎn)高于全球通量觀測(cè)網(wǎng)絡(luò)臺(tái)站匯交數(shù)據(jù)能量閉合程度的平均水平，數(shù)據(jù)質(zhì)量較好。此外，為保證降水?dāng)?shù)據(jù)質(zhì)量，本數(shù)據(jù)集不僅包含了人工觀測(cè)降水?dāng)?shù)據(jù)，也收錄了同時(shí)期的自動(dòng)觀測(cè)降水?dāng)?shù)據(jù)。

本數(shù)據(jù)集中，日尺度蒸散量和二氧化碳凈交換量有效觀測(cè)數(shù)據(jù)為90%和93%；人工降水量有效觀測(cè)數(shù)據(jù)為100%，自動(dòng)觀測(cè)為96%；凈輻射量和感熱通量有效觀測(cè)數(shù)據(jù)約為95%。數(shù)據(jù)缺失的原因主要分為兩類：一是觀測(cè)儀器運(yùn)行故障導(dǎo)致的數(shù)據(jù)缺失，包括供電故障、設(shè)備維護(hù)等；二是數(shù)據(jù)處理過程導(dǎo)致的數(shù)據(jù)缺失，比如夜間和降水天氣數(shù)據(jù)的篩選等。

4 數(shù)據(jù)使用方法和建議

本數(shù)據(jù)集為欒城試驗(yàn)站自通量觀測(cè)系統(tǒng)布設(shè)完成至2013 年積累的通量數(shù)據(jù)，其觀測(cè)、處理和質(zhì)量控制與評(píng)估均采用國際通用方法，并根據(jù)自身?xiàng)l件和觀測(cè)情況進(jìn)行了改進(jìn)，可靠性高。數(shù)據(jù)跨度從2007 年10 月至2013 年9 月，覆蓋6 個(gè)連續(xù)作物年，觀測(cè)樣地屬于華北地區(qū)典型的潮褐土高產(chǎn)農(nóng)業(yè)生態(tài)類型，具有較強(qiáng)的代表性。本數(shù)據(jù)集可為農(nóng)業(yè)節(jié)水理論與水文水資源研究等相關(guān)領(lǐng)域的發(fā)展提供堅(jiān)實(shí)的數(shù)據(jù)支撐，包括農(nóng)田耗水特征與蒸散結(jié)構(gòu)分離[23-24]，土壤水利用層次與地下水補(bǔ)給途徑解析[23,25-26]，農(nóng)田生態(tài)系統(tǒng)碳氮水循環(huán)過程[27-28]，農(nóng)田水平衡和適水種植制度調(diào)整[16-17,29-30]，農(nóng)業(yè)生產(chǎn)力與水資源可持續(xù)利用程度評(píng)價(jià)[31]，作物模型的驗(yàn)證與改進(jìn)[29-30,32]等。此外，本數(shù)據(jù)集可與其他觀測(cè)臺(tái)站的數(shù)據(jù)進(jìn)行綜合集成，服務(wù)于陸地生態(tài)系統(tǒng)對(duì)全球變化的響應(yīng)、區(qū)域尺度陸地生態(tài)系統(tǒng)物質(zhì)循環(huán)與能量流動(dòng)過程以及生態(tài)系統(tǒng)管理政策制定等相關(guān)領(lǐng)域。

目前通量觀測(cè)數(shù)據(jù)仍普遍存在不同程度的能量不閉合現(xiàn)象，數(shù)據(jù)處理和質(zhì)量控制也存在一定的不確定性。因此，本數(shù)據(jù)集在使用過程中需注意以下幾個(gè)方面：

（1）隨著ChinaFLUX 數(shù)據(jù)管理技術(shù)體系的更新和完善，不同年份匯交的通量數(shù)據(jù)處理方法不完全相同，因此本數(shù)據(jù)集可能與早期發(fā)表的數(shù)據(jù)存在一定差異（誤差允許范圍內(nèi)）。

（2）本數(shù)據(jù)集的數(shù)據(jù)管理和質(zhì)量控制方法詳細(xì)信息可參考Zhang 等[16,26]、Shen 等[17]、于貴瑞等[33]、Wen 等[20]和Yu 等[6]已發(fā)表的文獻(xiàn)。

（3）本數(shù)據(jù)集各項(xiàng)觀測(cè)數(shù)據(jù)受儀器運(yùn)行狀態(tài)或數(shù)據(jù)處理過程影響，存在不同程度的數(shù)據(jù)缺失。數(shù)據(jù)文件中對(duì)插補(bǔ)數(shù)據(jù)進(jìn)行了標(biāo)注，建議優(yōu)先采用未插補(bǔ)數(shù)據(jù)，以減少不確定性。

（4）針對(duì)降水?dāng)?shù)據(jù)，人工觀測(cè)可能會(huì)在降水量較大的連續(xù)降水日無法及時(shí)觀測(cè)導(dǎo)致數(shù)據(jù)偏小，而自動(dòng)觀測(cè)可能對(duì)小雨0.2 mm 的降水難以捕捉，有時(shí)也會(huì)由于儀器問題出現(xiàn)缺測(cè)，這是目前降水觀測(cè)普遍存在的問題，所以建議根據(jù)研究實(shí)際情況選擇所需降雨數(shù)據(jù)。

致 謝

衷心感謝CERN 綜合中心和ChinaFLUX 數(shù)據(jù)資源管理人員以及欒城站觀測(cè)人員的支持！感謝中國科學(xué)院地理所張雷明、陳智老師在數(shù)據(jù)質(zhì)量控制和處理程序方面給予的指導(dǎo)！感謝中科院農(nóng)業(yè)資源研究中心王玉英、程一松、閔雷雷、齊永青和郭英博士在數(shù)據(jù)收集、處理與監(jiān)管方面給予的幫助！

數(shù)據(jù)作者分工職責(zé)

張玉翠（1984—），女，山東省濱州市人，博士，副研究員，研究方向?yàn)樯鷳B(tài)水文學(xué)與同位素水文學(xué)。主要承擔(dān)工作：觀測(cè)系統(tǒng)的長(zhǎng)期連續(xù)運(yùn)行、數(shù)據(jù)質(zhì)量控制與整合。

姜寒冰（1993—），女，河北省泊頭市人，博士研究生，研究方向?yàn)檗r(nóng)田耗水機(jī)理。主要承擔(dān)工作：數(shù)據(jù)處理與分析。

張傳偉（1991—），男，山東省棗莊市人，碩士研究生，研究方向?yàn)槎喑叨绒r(nóng)田耗水特征及其模擬。主要承擔(dān)工作：數(shù)據(jù)處理與質(zhì)量控制。

沈彥俊（1971—），男，河北省康?？h人，博士，研究員，研究方向?yàn)樗膶W(xué)與水資源。主要承擔(dān)工作：觀測(cè)系統(tǒng)的總體運(yùn)行和數(shù)據(jù)集的整體架構(gòu)。