虞蘭蘭，張海寧，虞洋，江宇，張玉強，金海琴

(1. 國家海洋局北海信息中心 青島市 266061；2. 海軍北海艦隊海洋水文氣象中心 青島市 266003；3. 青島91049部隊 青島市 266000)

海冰觀測數(shù)據(jù)質(zhì)量控制方法研究及北海區(qū)冰情分析

虞蘭蘭，張海寧，虞洋，江宇，張玉強，金海琴

(1. 國家海洋局北海信息中心 青島市 266061；2. 海軍北海艦隊海洋水文氣象中心 青島市 266003；3. 青島91049部隊 青島市 266000)

根據(jù)海冰觀測數(shù)據(jù)的特點，分析、借鑒幾種常用的氣象、海洋數(shù)據(jù)質(zhì)量控制方法，制定海冰觀測數(shù)據(jù)的質(zhì)量控制體系，包括非法碼檢驗、全等性檢驗、極值檢驗、一致性檢驗、統(tǒng)計特性檢驗、可視化圖形繪制檢驗等方法。選取最近4個年度 (2013/2014-2016/2017年度)北海區(qū)5個海冰測站的觀測數(shù)據(jù)，分別進行人工和程序質(zhì)量控制，結(jié)果顯示，目前業(yè)務(wù)中使用的人工檢測質(zhì)量控制效率較低，海冰質(zhì)控程序則提升了海冰觀測數(shù)據(jù)質(zhì)量控制的效率和準(zhǔn)確度，在日常業(yè)務(wù)中有較高的應(yīng)用價值。同時，對質(zhì)控后的海冰數(shù)據(jù)進行北海區(qū)冰情分析發(fā)現(xiàn)，從海冰持續(xù)時間、總冰量、浮冰冰狀分布情況，均可看出2015/2016年度北海區(qū)的冰情最為嚴(yán)重；在同一年度，北海區(qū)冰情基本隨緯度增加而加重，而在相似緯度上，黃海冰情比渤海更為嚴(yán)重。

海冰；質(zhì)量控制；北海區(qū)；冰情分析

海冰是我國中、高緯度海域的重要海洋災(zāi)害之一，頻繁出現(xiàn)在渤海和黃海中北部海域 (本文簡稱北海區(qū))，主要由海水凍結(jié)及一部分江河注入的淡水冰組成，是北海區(qū)特有的水文現(xiàn)象。海冰對海面氣熱交換的影響十分顯著，直接影響地球、大氣的能量收支平衡[1]，受西伯利亞寒潮侵襲的影響，我國北海區(qū)每年冬季均出現(xiàn)不同程度的結(jié)冰[2]。據(jù)資料統(tǒng)計，北海區(qū)海冰為一年冰，發(fā)生嚴(yán)重 (或比較嚴(yán)重)海冰災(zāi)害的周期約為5年，輕冰年與重冰年的冰情差異很大[3]。輕冰年，海冰出現(xiàn)范圍小、持續(xù)時間短，厚度較薄，對航運交通和其他海上活動的影響較?。恢乇?，北海區(qū)大面積海域被厚冰覆蓋[4]，造成航道堵塞、海洋工程設(shè)施損壞，嚴(yán)重制約港口航運、海上油氣開采、漁業(yè)生產(chǎn)和水產(chǎn)養(yǎng)殖等發(fā)展，影響了人民的生產(chǎn)生活[5]。掌握有冰海域的環(huán)境條件，利用多種手段開展海冰冰情監(jiān)測對國民經(jīng)濟發(fā)展和國防建設(shè)具有十分重要的意義[6]。

對觀、監(jiān)測數(shù)據(jù)進行質(zhì)量控制，及時對發(fā)現(xiàn)的錯誤數(shù)據(jù)進行分析判斷，能夠妥善解決資料中可能存在的問題，剔除資料中的錯誤，提高資料的可靠性與準(zhǔn)確性[7]，是實測數(shù)據(jù)應(yīng)用分析的前提。而質(zhì)量控制方法的好壞則直接影響到數(shù)據(jù)質(zhì)量，是數(shù)據(jù)質(zhì)量控制的重點[8]。由此可見，對海冰資料進行質(zhì)量控制，獲取準(zhǔn)確的海冰監(jiān)測數(shù)據(jù)，從而提高海冰實測資料的可靠性，對于推進海冰觀測資料的廣泛利用、建立較為完善的海冰防災(zāi)減災(zāi)系統(tǒng)具有意義重大。因此，海冰數(shù)據(jù)的質(zhì)量控制是制約相關(guān)研究工作的技術(shù)瓶頸，是海冰冰情監(jiān)測成果應(yīng)用的基礎(chǔ)，為確保海冰數(shù)據(jù)的真實性和準(zhǔn)確性，必須對其實行準(zhǔn)確有效的數(shù)據(jù)質(zhì)量控制和檢驗評估。目前國內(nèi)外針對海冰觀測數(shù)據(jù)的質(zhì)量控制研究較少，并未產(chǎn)生一套完整可靠的海冰質(zhì)控方法體系。

自1969年特大冰封以來，我國海冰監(jiān)測、研究和預(yù)報工作全面展開。隨著遙感技術(shù)的發(fā)展，20世紀(jì)80年代開始大量應(yīng)用遙感技術(shù)監(jiān)測海冰。21世紀(jì)以來，北海區(qū)海洋站已形成業(yè)務(wù)化運行的海冰觀測數(shù)據(jù)，然而，在實際業(yè)務(wù)中缺乏專門的海冰質(zhì)控軟件，主要對其進行初步的人工檢驗，影響了海冰觀測資料使用的規(guī)范性，亟待開發(fā)系統(tǒng)的海冰質(zhì)控應(yīng)用程序，對海冰資料進行有效的質(zhì)量控制，以滿足業(yè)務(wù)的需求。本文借鑒氣象、海水溫鹽度、海浪等業(yè)務(wù)化觀測數(shù)據(jù)的質(zhì)量控制方法，針對海冰觀測數(shù)據(jù)的特點，初步分析制定海冰數(shù)據(jù)質(zhì)量控制的流程，通過實例進行海冰數(shù)據(jù)質(zhì)量控制的檢驗，并分析最近4個年度北海區(qū)冰情分布特征。

1 質(zhì)量控制方法簡介

根據(jù)觀測數(shù)據(jù)的實際情況進行質(zhì)量控制方法的選擇和應(yīng)用，是數(shù)據(jù)質(zhì)量控制的重點，質(zhì)量控制方法直接影響了數(shù)據(jù)的質(zhì)量。由于目前并無系統(tǒng)的海冰質(zhì)量控制方法體系，本文主要借鑒氣象及其它水文數(shù)據(jù)質(zhì)量控制的常用方法，包括非法碼檢驗、全等性檢驗、極值檢驗、一致性檢驗、海洋環(huán)境氣候特性檢驗、統(tǒng)計特性檢驗、人機交互圖形繪制檢驗等[7-14]。

1.1 非法碼檢驗

對溫鹽度、海浪、海冰、氣象數(shù)據(jù)等要素的海洋站觀測數(shù)據(jù)，必須按照規(guī)定的格式和代碼進行記錄，每一類觀測要素的標(biāo)準(zhǔn)數(shù)據(jù)格式，均具有固定的編碼方式和特定屬性。根據(jù)規(guī)定的各類要素數(shù)據(jù)編碼格式，對不同位置出現(xiàn)的符號、數(shù)字或字符進行范圍檢驗，即為非法碼檢驗。

1.2 全等性檢驗

觀測數(shù)據(jù)的部分參數(shù)相對于固定的海洋站，往往是長期不變的，如資料類型、觀測點經(jīng)緯度、觀測點水深、觀測方法、儀器名稱、觀測儀器海拔高度等參數(shù)值，質(zhì)控過程中檢驗上述參數(shù)值是否保持不變，即為全等性檢驗。

1.3 極值檢驗

極值檢驗是指結(jié)合海洋環(huán)境要素的物理征、歷史極值等統(tǒng)計經(jīng)驗，分別給定各要素的極大值和極小值，作為該要素觀測數(shù)據(jù)的上、下限，從而判定數(shù)據(jù)是否異常的方法，亦被稱為極值控制法。極值控制的關(guān)鍵是極值的確定，若觀測數(shù)據(jù)超出極值的范圍，則判定該數(shù)據(jù)可疑。選取適當(dāng)?shù)臉O值范圍可有效檢驗出極端異常的數(shù)據(jù)，然而，極值控制無法檢驗出數(shù)值在極值范圍內(nèi)的異常數(shù)據(jù)，需進一步采用更為精準(zhǔn)的控制方法。

1.4 一致性檢驗

一致性檢驗法通過分析要素之間的內(nèi)在關(guān)系和時空分布聯(lián)系，對要素的正確性進行判定，包括內(nèi)部一致性檢驗、空間一致性檢驗及時間持續(xù)性檢驗。

內(nèi)部一致性檢驗主要分析各要素自身屬性及要素之間的聯(lián)系是否正確，如對于風(fēng)向這一要素，出現(xiàn)大于360°或小于0°的風(fēng)向、風(fēng)速為0而風(fēng)向未記為C的情況，則風(fēng)向未通過內(nèi)部一致性檢驗；要素之間的聯(lián)系，如總云量≥低云量、氣溫≥露點溫度、風(fēng)速與風(fēng)浪的浪高成正比等，若出現(xiàn)有悖于上述關(guān)系的要素時，則該要素未通過內(nèi)部一致性檢驗。

由于不同空間的數(shù)據(jù)分布特征具有較大的差異和不確定性，故空間一致性檢驗使用四分位距檢查法，又稱四分位法[15-17]。四分位法主要用于描述觀測數(shù)據(jù)的偏態(tài)分布、分布不明確及末端無確定值等特征，可消除異常大值的影響，相對于均值、標(biāo)準(zhǔn)差等更加穩(wěn)定可靠。本文根據(jù)四分位法的實質(zhì)內(nèi)涵，對其現(xiàn)有的計算公式進行了歸納、精簡。將觀測數(shù)據(jù)整理為升序排列的樣本X=[x1， x2，…，xn]， n≥3， 選取 3 個數(shù)據(jù)點將 X 劃分為四等分，數(shù)據(jù)點對應(yīng)數(shù)值即為四分位數(shù)，分別用Q1、Q2、Q3表示，計算步驟如下所示:

另外，由于月經(jīng)血內(nèi)會混有一些具有抗凝作用的活性酶，所以月經(jīng)為非凝固性的，這樣有利于經(jīng)血的排出。但是當(dāng)月經(jīng)量過多的時候，抗凝血酶原的量會出現(xiàn)相對不足，就會使經(jīng)血內(nèi)混有血塊。此時，由于經(jīng)血里一些細(xì)碎甚至大塊的血塊，所以顏色也會偏黑。

1.4.1 計算第1個四分位數(shù)Q1

Q1在X樣本中的下標(biāo)為進行mQ向下取整 (即截斷取整)得到整數(shù) ，則1可由用xm和xm+1插值得到:

1.4.2 計算第3 個四分位數(shù)Q3

Q3在X樣本中的下標(biāo)為進行pQ向下取整 (即截斷取整)得到整數(shù)，則1可用xp和xp+1插值得到:

1.4.3 計算四分位距IQR

1.4.4 確定樣本X中異常值的界限為:

公式(4)中的Fl、Fu分別為異常值的下限和上限，在[Fl， Fu]以外的數(shù)據(jù)認(rèn)定為異常值。然而四分位法也有一定的局限性，該方法必須保證異常數(shù)據(jù)占總數(shù)據(jù)量的比例較小，否則將無法有效識別異常數(shù)據(jù)，故使用時應(yīng)結(jié)合其它檢驗方法對數(shù)據(jù)進行質(zhì)量控制。

時間持續(xù)性檢驗主要檢測數(shù)據(jù)是否不隨時間發(fā)生改變，應(yīng)根據(jù)要素的特征進行檢查時間的設(shè)置，當(dāng)某要素在檢查時間內(nèi) (如24 h)樣本的標(biāo)準(zhǔn)差小于給定的最低值，則該要素未通過時間持續(xù)性檢驗。

1.5 海洋環(huán)境氣候態(tài)特性檢驗

根據(jù)各要素的海洋環(huán)境氣候態(tài)分布特性進行資料檢驗，將某海洋站某要素的月統(tǒng)計值與該站該要素的多年歷史統(tǒng)計值進行比較，分析其超出或低于給定值的原因，排除異常天氣現(xiàn)象的影響因素后，即認(rèn)為該要素未通過海洋環(huán)境氣候態(tài)特性檢驗。氣候態(tài)檢驗必須保證已有的觀測數(shù)據(jù)序列足夠長 (至少20年)，否則不具備進行該項檢驗的條件。

1.6 統(tǒng)計特性檢驗

一定時間序列的觀測資料在理論上往往符合某些概率統(tǒng)計的特性，數(shù)據(jù)對應(yīng)的隨機變量和隨機過程相互獨立并服從一定的分布。對觀測資料進行數(shù)據(jù)的統(tǒng)計特性檢驗，可選用的檢驗方法包括3δ方法 (萊茵達檢驗法)、卡方擬合優(yōu)度檢驗、輪次檢驗、格拉布斯檢驗以及狄克遜檢驗等，應(yīng)根據(jù)數(shù)據(jù)的觀測特性和實際情況選取統(tǒng)計檢驗方法。對時間序列為一個月的觀測數(shù)據(jù)一般采用萊茵達檢驗法進行質(zhì)量控制。

1.7 人機交互圖形繪制檢驗

國內(nèi)外成功的質(zhì)量控制體系表明，數(shù)據(jù)質(zhì)量控制的完全自動化無法解決所有的資料質(zhì)量問題。質(zhì)量控制程序計算出的異常數(shù)據(jù)未必是錯誤數(shù)據(jù)，若不加分析直接算作缺測或內(nèi)插得出替代值，則會造成數(shù)據(jù)誤判。因而，使用人工檢查和計算機判定相結(jié)合的方法進行數(shù)據(jù)的統(tǒng)計和分析，是觀測資料質(zhì)量控制的主要方法之一，而便捷友好的人機交互界面是實現(xiàn)該方法的關(guān)鍵步驟。在人機交互過程中，工作人員根據(jù)自己的專業(yè)知識和經(jīng)驗，核查相關(guān)要素的變化情況，對質(zhì)控軟件判定出的異常觀測數(shù)據(jù)進行最終的分析判斷。

2 海冰觀測數(shù)據(jù)質(zhì)量控制體系的建立與檢驗

觀測數(shù)據(jù)質(zhì)量控制的核心是通過分析數(shù)據(jù)的特征，選擇適當(dāng)?shù)馁|(zhì)量控制方法，以提高可疑數(shù)據(jù)的檢出率。海冰觀測要素包括觀測冰量、冰型、冰表面特征、冰狀、最大浮冰塊水平尺度、浮冰密集度、浮冰漂流方向和速度、沿岸冰堆積量、沿岸冰堆積高度等。本文根據(jù)海冰各要素的特性，全面考慮海冰資料的各種誤差，選用非法碼檢驗、全等性檢驗、極值檢驗、一致性檢驗、統(tǒng)計特性檢驗、人機交互圖形繪制檢驗等方法(表1)，初步建立海冰觀測數(shù)據(jù)質(zhì)量控制的方法體系。由于觀測數(shù)據(jù)時間序列的限制，海洋環(huán)境氣候態(tài)特性檢驗暫未考慮，應(yīng)在時間序列達到20年以上時將該檢驗加入海冰質(zhì)控方法體系。

表1 海冰觀測數(shù)據(jù)質(zhì)量控制方法體系

現(xiàn)有的海冰質(zhì)量控制方法主要為人工檢驗，極易造成異常數(shù)據(jù)的漏檢。選取最近4年的北海區(qū)海冰觀測數(shù)據(jù) (2013/2014-2016/2017)，對本文建立的海冰質(zhì)控體系的準(zhǔn)確性進行驗證與分析。首先，對海冰觀測數(shù)據(jù)分別進行人工檢測與程序質(zhì)量控制的結(jié)果對比，將發(fā)現(xiàn)的問題條數(shù)列出，如表2所示:人工檢測主要針對標(biāo)題記錄和說明記錄，最核心的數(shù)據(jù)記錄的質(zhì)量無法得到保證；而質(zhì)控程序輸出的質(zhì)控報告顯示，其數(shù)據(jù)記錄亦存在較多的問題，需要資料產(chǎn)出單位進行檢查和修改。經(jīng)人工復(fù)核，質(zhì)控程序檢測出的異常數(shù)據(jù)無錯檢現(xiàn)象，充分證明了海冰質(zhì)控程序的可行性，提升了海冰觀測數(shù)據(jù)質(zhì)控過程的效率和準(zhǔn)確度，可在業(yè)務(wù)運行中推廣使用。

表2 近4年(2013/2014-2016/2017年度)海冰觀測數(shù)據(jù)質(zhì)量控制結(jié)果對比

3 北海區(qū)海冰冰情分析

選取2013/2014-2016/2017年度 (下文全部以13/14-16/17年度代替)北海區(qū)最北部5個海冰測站的觀測數(shù)據(jù)，進行質(zhì)量控制后作相應(yīng)的修正并剔除可疑值，以得出最近4年北海區(qū)海冰(圖1)的冰情分布特征。

圖1 北海區(qū)海冰測站的位置示意圖

統(tǒng)計各海冰測站能夠觀測到冰量的有冰日數(shù)，即初冰日期與終冰日期之間的天數(shù)，可表征海冰持續(xù)時間的長短，如圖2所示:從時間上分析，15/16年度北海區(qū)海冰持續(xù)的時間最長，5個測站總計309天，其余3個年度分別為249、234、260天。從地域上分析，黃海北部近岸 (東港)海冰持續(xù)的時間最長，遼東灣北部近岸 (葫蘆島、鲅魚圈)次之，遼東灣南部 (芷錨灣、溫坨子)海冰持續(xù)時間最短。

圖2 北海區(qū)13/14-16/17年度各海冰測站有冰日數(shù)統(tǒng)計(單位: 天)

統(tǒng)計各海冰測站觀測到的總冰量，即各觀測時刻海冰覆蓋面積占整個能見海面的成數(shù)之和?？偙吭酱?，表明海冰面積越大，如圖3所示:從時間上分析，15/16年度北海區(qū)海冰總冰量最大，5個測站總計3 074，其余3個年度分別為1 566、1 782、1 746。從地域上分析，基本上緯度較高的葫蘆島、鲅魚圈總冰量最大，東港次之，芷錨灣、溫坨子最小。

圖3 北海區(qū)13/14-16/17年度各海冰測站總冰量統(tǒng)計

由于13/14-16/17年度觀測得出的固定冰冰量均為0，為研究海冰水平尺度的特征，圖4繪制出浮冰的冰狀統(tǒng)計情況。假設(shè)浮冰最大水平尺度為L(單位:m)，L與浮冰冰狀的關(guān)系為:巨冰盤 (L≥2000)、大冰盤 (500≤L＜2000)、中冰盤 (100≤L＜500)、小冰盤 (20≤L＜100)、冰塊 (2≤L＜20)、碎冰 (L＜2)。

浮冰冰狀分布特征為:北海區(qū)海冰主要以冰塊、碎冰這兩類較小水平尺度的形式存在，小冰盤出現(xiàn)次數(shù)約為冰塊出現(xiàn)次數(shù)的一半；除15/16年度出現(xiàn)40次大冰盤的統(tǒng)計外，巨冰盤和大冰盤的統(tǒng)計次數(shù)均為0。從冰狀的時間分布來看，15/16年度浮冰各冰狀出現(xiàn)次數(shù)均遠(yuǎn)高于其他年度。

圖4 北海區(qū)13/14-16/17年度浮冰冰狀統(tǒng)計次數(shù)(單位: 次)

綜上所述，從海冰持續(xù)的時間、總冰量大小、浮冰冰狀分布情況，均可看出15/16年度北海區(qū)的冰情在最近4個年度中最為嚴(yán)重。在同一年度，北海區(qū)冰情隨緯度增加而加重；在相似緯度上，黃海冰情比渤海更為嚴(yán)重。

4 結(jié)論與展望

(1)北海區(qū)海冰觀測數(shù)據(jù)在實際質(zhì)控業(yè)務(wù)中主要進行人工檢驗，影響了海冰觀測資料使用的規(guī)范性。本文借鑒氣象、海水溫鹽度、海浪等業(yè)務(wù)化觀測數(shù)據(jù)的質(zhì)量控制方法，針對海冰觀測要素的特點，分析其質(zhì)量控制過程，初步建立了北海區(qū)海冰質(zhì)控方法體系，并將其應(yīng)用于海冰的實測數(shù)據(jù)。與目前的人工檢驗結(jié)果相比，本文建立的海冰質(zhì)控程序更加方便快捷，質(zhì)檢正確率亦得到了提升。

(2)從海冰持續(xù)的時間、總冰量大小、浮冰冰狀分布情況可看出，15/16年度北海區(qū)的冰情在13/14-16/17 4個年度中最為嚴(yán)重?；旧显谕荒甓龋焙^(qū)冰情隨緯度增加而加重；在相似緯度上，黃海冰情比渤海更為嚴(yán)重。

(3)在海冰質(zhì)量控制體系的應(yīng)用過程中，由于實測海冰數(shù)據(jù)時、空序列的限制，無法進行海洋環(huán)境氣候態(tài)特性和空間一致性的檢驗過程，在時間序列達到20年以上、空間觀測點分布達到相應(yīng)密度后，應(yīng)在海冰質(zhì)控體系中分別加入長期觀測數(shù)據(jù)的氣候態(tài)統(tǒng)計信息和四分位距檢查法。

(4)本文僅以萊茵達檢驗法進行了海冰觀測數(shù)據(jù)的統(tǒng)計特性檢驗，為更全面地檢驗數(shù)據(jù)的隨機過程，應(yīng)在進一步的質(zhì)控方法研究工作中考慮卡方擬合優(yōu)度檢驗、輪次檢驗、格拉布斯檢驗、狄克遜檢驗等方法。此外，應(yīng)考慮更加便捷友好的人機交互界面，以便于工作人員根據(jù)程序質(zhì)控結(jié)果進行人工復(fù)核校驗的分析判斷。

[1] 侯海倩， 李家國， 米曉飛等. 基于資源 3號衛(wèi)星的海冰監(jiān)測研究[J]. 遙感信息， 2014， 29(6): 36-40.

[2] 呂曉東. 海冰基礎(chǔ)知識及船舶冰區(qū)航行的注意事項[J].學(xué)術(shù)， 2015， (2): 77-79.

[3] 周群， 魏立新， 黃煥卿. 秋季巴倫支海海溫異常對冬季我國渤海冰情的可能影響[J]. 海洋學(xué)報， 2016， 38(3):40-48.

[4] 鄭冬梅， 張書穎， 周志強等. 逐步回歸分析在渤海海冰等級預(yù)報中的應(yīng)用[J]. 海洋預(yù)報， 2015， 32(2): 57-61.

[5] 吳鑫， 王曉東， 王召進等. 東營市海冰發(fā)展規(guī)律及防范措施淺析[J]. 海洋開發(fā)與管理， 2014， (10): 54-57.

[6] 朱亞平， 程周杰. 基于FY-2C的渤海海冰監(jiān)測分析研究[J]. 海洋技術(shù)學(xué)報， 2015， 34(6): 6-12.

[7] 于婷， 劉玉龍， 楊錦坤等. 實時和延時海洋觀測數(shù)據(jù)質(zhì)量評估方法研究[J]. 海洋通報， 2013， 32(6): 610-625.

[8] 史靜濤， 周智海. 海洋站數(shù)據(jù)質(zhì)量控制技術(shù)探討[J]. 海洋技術(shù)， 2011， 30(1): 114-117.

[9] 耿愛菊， 孟楠. 地面氣象觀測軟件質(zhì)控異常數(shù)據(jù)的處理[J]. 山東氣象， 2013， 33(4): 64-65.

[10] 郭春迓， 李天然， 胡東明等. 南海北部測站風(fēng)速質(zhì)量的控制方法[J]. 廣東氣象， 2016， 38(1): 44-52.

[11] 全國海洋標(biāo)準(zhǔn)化技術(shù)委員會中國標(biāo)準(zhǔn)出版社第五編輯室.海洋觀測工作執(zhí)行標(biāo)準(zhǔn)匯編[S]. 北京: 中國標(biāo)準(zhǔn)出版社， 2010: 15-81.

[12] 任芝花， 趙平， 張強等. 適用于全國自動站小時降水資料的質(zhì)量控制方法[J]. 氣象， 2010， 36(7): 123-132.

[13] 汪寧， 張曉慧. 海洋水文氣象實時數(shù)據(jù)質(zhì)量控制[J]. 中國科學(xué)技術(shù)綜述， 2014， (19): 223-223.

[14] 許銘， 楊國輝， 董海英. 區(qū)域自動站雨量資料質(zhì)量控制方法及應(yīng)用[J]. 北京農(nóng)業(yè)·下旬刊， 2013， (6): 129-130.

[15] GonzMez-Rouco J F， JimSnez J L， Quesada V， et a1. Quality control and homogeneity of precipitation data in the southwest of europe. J Climate， 2000， 14: 964-978.

[16] 何宇飛， 楊冬梅， 朱榮等. 四分位法在地磁日變化觀測數(shù)據(jù)分析中的應(yīng)用[J]. 地震地磁觀測與研究， 2011， 32(1): 65-70.

[17] 趙永寧， 葉林， 朱倩雯. 風(fēng)電場棄風(fēng)異常數(shù)據(jù)簇的特征及處理方法[J]. 電力系統(tǒng)自動化， 2014， 38(21): 39-46.

2017-08-23

2015年北海分局海洋科技項目--海洋站海冰觀測資料的質(zhì)控方法研究與初步應(yīng)用 (2015B14)。