• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Data Quality Control Method of a New Drifting Observation Technology Named Drifting Air-Sea Interface Buoy

    2024-03-12 11:14:52LIShuoWANGBinDENGZenganWUBaoqinZHUXiandeandCHENZhaohui
    Journal of Ocean University of China 2024年1期

    LI Shuo , WANG Bin , , DENG Zeng’an , WU Baoqin, ZHU Xiande,and CHEN Zhaohui

    1) National Ocean Technology Center, Tianjin 300112, China

    2) School of Marine Science and Technology, Tianjin University, Tianjin 300072, China

    3) Key Laboratory of Ocean Observation Technology, Ministry of Natural Resources, Tianjin 300112, China

    4) People’s Liberation Army 92859, Tianjin 300061, China

    5) Physical Oceanography Laboratory, Ocean University of China, Qingdao 266100, China

    Abstract An integral quality control (QC) procedure that integrates various QC methods and considers the design indexes and operational status of the instruments for the observations of drifting air-sea interface buoy was developed in the order of basic inspection followed by targeted QC. The innovative method of combining a moving Hampel filter and local anomaly detection complies with statistical laws and physical processes, which guarantees the QC performance of meteorological variables. Two sets of observation data were used to verify the applicability and effectiveness of the QC procedure, and the effect was evaluated using the observations of the Kuroshio Extension Observatory buoy as the reference. The results showed that the outliers in the time series can be correctly identified and processed, and the quality of data improved significantly. The linear correlation between the quality- controlled observations and the reference increased, and the difference decreased. The correlation coefficient of wind speed before and after QC increased from 0.77 to 0.82, and the maximum absolute error decreased by approximately 2.8 m s?1. In addition, air pressure and relative humidity were optimized by 10?3 – 10?2 orders of magnitude. For the sea surface temperature, the weight of coefficients of the continuity test algorithm was optimized based on the sea area of data acquisition, which effectively expanded the applicability of the algorithm.

    Key words drifting air-sea interface buoy; quality control; oceanic variables; meteorological variables; continuity test

    1 Introduction

    The drifting air-sea interface buoy (DrIB), a small-sized,lightweight, and cost-effective mini-buoy, is mainly used to measure basic variables at the air-sea interface by expendable and drifting methods in real time. The large-scale application of the DrIB can realize the high-frequency acquisition of key meteorological (wind speed, wind direction, and air pressure at approximately 3 m and air temperature and relative humidity at approximately 2.8 m above the sea surface) and oceanic (sea surface temperature (SST)and wave parameters) variables in the global ocean (Centurioniet al., 2019), which is beneficial to understanding ocean dynamics, air-sea interaction, and global climate change and forecasting of modal accuracy.

    The external structure of the DrIB consists of two parts:the above-water and underwater units (Fig.1). The above water unit includes a hemispherical float and a mast, and the underwater unit consists of two mutually perpendicular water sails and an underwater counterweight. The DrIB is composed of a meteorological observation module, a hydrographic observation module, a measurement and control system, a communication module, and a power supply module. After the DrIB is deployed, it carries out observations of key meteorological and hydrological variables at the air-sea interface by drifting with the currents.After preliminary online calibration, the observations are transmitted to the shore/ship-based data management system, which can receive, process, store, and distribute observations in real time through the satellite communication system. The DrIB has a sampling frequency of 1 Hz,and the output data are the average of 120 samples with optional sampling intervals of 10 min / 30 min / 1 h. Table 1 shows the observational performance of the DrIB and the instruments installed on it.

    Table 1 Observation indexes and associated instruments of the DrIB

    Fig.1 (a) Dimensional view of the drifting air-sea interface buoy and (b) the site of drifting air-sea interface buoy observation in the Western Pacific. The abbreviations AT and RH in (a) represent air temperature and relative humidity,respectively.

    Marine science is an observation-based discipline (Wuet al., 2020). Observations are the critical foundation for scientific research, which is directly affected by data quality. However, operation errors, marine environment, unstable communication transmission,etc., can interfere with the observation results and cause data anomalies, and most of the observations cannot be applied directly. Therefore,rigorous and effective quality control (QC) of the received data and removal of the error and abnormal values hidden in the time series are necessary to improve the accuracy and reliability of observations.

    For the QC research of marine environmental observations, many institutions and experts have shared their valuable experiences. Foreign research has been carried out earlier, and after continuous improvement and development,QC technology is more mature nowadays. For example,the drifting buoy of the Surface Velocity Program (SVP)can transmit back its position and SST data, and the SVP data collection center conducts QC based on whether the value violates the set continuity criterion (Hansen and Herman, 1989; Hansen and Poulain, 1996). The National Data Buoy Center of America manages 50 ocean stations and more than 100 moored buoys. It has developed an automated QC program to perform transmission error checks,range inspections, time continuity checks, correlation checks,etc.(Jarraud, 2008; National Data Buoy Center, 2009). The Kuroshio Extension Observatory (KEO) buoy was first deployed in 2004. Currently, the time series and statistical properties of hydrometeorological variables, including wind,air temperature, air pressure, relative humidity, precipitation, long- and short-wave radiation, seawater temperature,salinity,etc., can be independently checked and identified for data quality by a complete QC program (Andersonet al., 2021). In addition, the glider toolbox based on Matlab is specially used to correct the thermal lag bias of the salinity observed by underwater gliders (Troupinet al.,2015).

    Comparatively, China’s research is relatively late. Verified by the observations of weather stations and ocean observatories, real-time meteorological and hydrological data can be quality-controlled by determining the range threshold and physical relationships between elements (Li and Li, 1997). The QC software developed by the National Marine Data and Information Service realizes the automatic QC of the data collected by ocean observatories (Yanget al.,2017); Wuet al.(2019)used a laboratory high-precision salinometer to compare and evaluate the observations of the ship-borne Conductivity, Temperature, Depth (CTD)instrument for complete QC. Liuet al.(2020)proposed a real-time QC procedure including thermal lag correction and 12 tests based on the temperature-salt profile obtained by the domestic Sea Wing underwater glider. Renet al.(2021) innovatively integrated the spike inspection and PauTa criterion, effectively eliminating the abnormal values in sea temperature and salinity observations.

    In summary, extensive research has been carried out on QC technologies of weather and ocean observatories, fixed buoy observations, and data from underwater sensors on mobile apparatus. However, studies on simultaneous observations of key meteorological and hydrological variables at the air-sea interface through drifting ways are relatively scarce. Moreover, with the rapid development of the DrIB observation network, observations have accumulated,and QC research has become urgent. In addition, drifting or mobile apparatuses, such as wave gliders and unmanned ships, are currently emerging, representing a significant trend in the development of marine observation technologies (Grareet al., 2021). However, no systematic data QC method has been formed yet, and the reliability of observations is uncertain. The research on the QC of DrIB observations and related theories and methods also applies to other mobile observation equipment, which can effectively compensate for shortcomings and solve the problem of poor quality of original data.

    With a focus on the air-sea interface scenarios, an integral QC procedure for DrIB observations was established in conjunction with the instrument’s design indexes and operational records. Basic inspection was performed first,and the observed variables were categorized for targeted QC. Two sets of observations from the Northwest Pacific were used as examples to analyze the effect of QC and verify the applicability, effectiveness, and advancement of QC methods and the related procedure.

    2 QC Procedure

    QC is mainly used for two types of errors in long-term continuous hydrometeorological observations. The gross error exceeds the observation capability of the equipment or normal variation range of the observed variables and is intuitive and easy to identify and eliminate; meanwhile,the outlier is within a reasonable range but deviates from the overall or local trend and thus needs to be screened by a variety of means (Tanet al., 2022).

    This paper applied various QC methods, such as land inspection, repetition check, incremental inspection, range inspection, and continuity inspection. Applicable continuity tests were used to carry out targeted QC on location data and meteorological and oceanic variables. The abnormal data were then eliminated or linearly interpolated based on the nearest samples. Each variable was considered as an independent one-dimensional time series. Fig.2 shows the specific steps of the QC procedure.

    Fig.2 Schematic chart of the QC procedure.

    2.1 Basic QC

    Previously, QC of observations tended to detect anomalies only from the data themselves, which might have affected efficiency. For DrIB, the observed variables are relatively abundant, and the operation of the buoy and changes in the air-sea environment are recorded well. QC should consider the instrument’s design indexes and operating attitude and maximize the advantages of operation records to identify abnormal data quickly and accurately. The basic QC mainly includes the following steps.

    1) Land inspection. Unlike mooring buoy, the DrIB may start operating before entering the sea due to land or shipbased functional test and other reasons, sending land information and directly causing abnormal SST observations.Such errors usually only exist at the initial stage of a continuous-time series. Therefore, land inspection needs to be performed first to eliminate invalid data when the DrIB is not in the sea and ensure that the observations are of the real ocean environment. By assessing whether the swing angle conforms to its in-water motion and whether a fault difference exists between the initial and stable periods of SST to make an accurate judgment of the working environment, the observations on land are removed, and a separate file of land observations is generated to facilitate subsequent inspection and verification.

    2) Repetition check. In cases of unstable satellite transmission and failed data storage modules, more than one piece of data may be stored in the land-based data receiving station simultaneously, which results in duplicate errors. The repetition check is applied to eliminate duplications and ensure a ‘one-to-one’ relationship between time and observation. A file of duplications is generated at the same time.

    3) Time increment test. Observations should be arranged in chronological order. ‘Time reversal’ is prevented by performing an incremental time test on the time records to ensure correct timing. The data detected with chaotic time are detected. If necessary, the order is adjusted, and deletion and adjustment actions are recorded in real time in the timing-related log.

    4) Range inspection. Range inspection is a practical judgment of whether values are reasonable based on the basic cognition of geography, general laws of hydrometeorological variables, and the observation capability of the equipment itself. For example, the SST generally varies within ?4℃ to 44℃ (Fenget al., 1999), and the maximum relative humidity does not exceed 100%. The design observation indexes of DrIB are shown in Table 1, and the data beyond the range are considered as set W.

    2.2 Targeted QC

    After basic QC, evident errors in the time series are eliminated. The working characteristic of ‘drifting with the tide’ determines the continuity and regularity of the observed variables. Next, the observations are divided into three categories: location data, oceanic variables, and meteorological variables. Then, the corresponding continuity criterion is set to identify the values that violate the continuing trend of the series for further targeted QC.

    1) Location data. For the sudden abnormal change (spike)in the drifting trajectory of DrIB, the spike test is used to determine whether the latitude and longitude at the spike are correct (Gouretski, 2018), which is determined by the following:

    whereln?1,ln, andln+1represent three consecutive longitude or latitude values, andαis a critical parameter, which is considered as 0.1 when the data sampling interval is 1 h.If the formula holds, thenlnis abnormal, and the abnormal values screened are set to G.

    2) Oceanic variables. The SVP data collection center integrates the spike and continuity tests and sets a weighted averagemiof iterative changes as the SST ‘spike’ detector to check the continuity of the SST time series (Hansen and Poulain, 1996). The determination ofmiis based on the recent changes in SST experienced by the buoy, and the nearest samples of the detected data are given the highest weight. The test is as follows: the first valid SST ism1, and the observations are tested sequentially by the following:

    The SST is anomalous if the formula cannot be formed andmiis not updated. The SVP data collection letsCbe 1/4 and ? be 0.5 K and applies the method to the detection of SST data in tropical Pacific Ocean.

    Compared with the relatively balanced and stable environment in tropical oceans, mid-latitude oceans exhibit more complex dynamics and a stronger air-sea nonlinear instability interaction. Thus, the SST observations are more dependent on the actual change processes of seawater at the time of acquisition. Moreover, for the improved applicability of the method, the contribution weight ofin situobservations should be appropriately increased in the calculation of the SST ‘spike’ detectormifor DrIB in mid-latitude oceans (Yuet al., 2010). In other words, the value of parameterCshould be increased. After experimental verification, the detection performance is the best when it is raised to 1/3. Correspondingly, ? is set to 0.5℃.The QC results may be incorrect when the DrIB crosses a strong front or when the sampling interval is large. Thus,the data are tested in both forward and reverse directions.The values identified as anomalous in both directions are taken as set S.

    3) Meteorological variables. Compared with the ocean,the atmosphere is more active and fast-changing, and the changing patterns of key meteorological variables, such as temperature, air pressure, wind speed, and relative humidity, are more significant. When the QC method mentioned above is used for oceanic variables to test meteorological variables, a growing gap exists between the observations and the detectormi, which results in most normal data being judged as abnormal and prevented from being directly applied.

    According to statistical theory, meteorological variables observed continuously at fixed positions tend to obey a probability distribution; therefore, statistical tests, such as the PauTa criterion (3δ criterion), which considers data with a residual error of more than three times the standard deviation from the arithmetic mean as outliers, can be used to check data (Luet al., 2019). However, if the arithmetic mean and standard deviation are used as judgment criteria,the extreme outliers are prone to biased results. Hampel filter is a filtering algorithm for identifying anomalies in a time series, and its operating principle is similar to that of the PauTa criterion, except that the criteria are the median and converted median absolute deviation, which can effectively avoid the interference of extreme values (Gholamiet al., 2022). In addition, as the spatial location of DrIB is dynamically changing, and the atmospheric environment varies significantly in different seas, the study of the statistical characteristics of its long-term time series has little significance. Based on the continuity and regularity of temporal and spatial changes of drifting observation, a moving Hampel filter has been proposed to assess the degree of data departure from daily variability using days as the time window unit to find outliers in the time series of meteorological variables observed by DrIB (Liu, 2017; Renet al., 2021). Suppose a time series of observations,q1,q2, ···,qn, with a data acquisition interval of 1 h. The time window for sampling is one day; that is,q1,q2, ···,q25are taken as a group, and the median isqme. Then, the following are obtained:

    wheredmstands for the local conversion ofMAD, and theqithat is more than three times ofdmfromqmeis abnormal.After the detection of the first group, the time window moves one piece of information in sequence and conducts detection again until the last group of the time series has been identified. The detected anomalies are set as Q1.

    The data in set Q1 are anomalous in terms of mathematical statistics. If the atmospheric environment is relatively stable, the fluctuation of variables is balanced, the local conversionMADis small, and correct data can be easily misclassified as abnormal. Local anomaly detection is further applied to avoid the risk of misjudgment by the moving Hampel filter. Theqsof any observations in the set Q1 is correct if the difference between it and the observation of any adjacent moment does not exceed the set maximum change threshold. The data misjudged as abnormal are considered set Q2. The change threshold should be larger than the design observation error of the DrIB to prevent the instrument’s error from interfering with the judgment.

    Given the rapid and instantaneous change, poor data continuity, and complex influencing factors of wind direction,only basic QC is performed. For the abnormal location data set Wg + G, abnormal oceanic observations set Ws + S,and abnormal meteorological observations (except wind direction) set Wq + Q1 ? Q2, for the retention of data continuity based on not breaking the statistical characteristics and changing trends of the time series, the results of linear interpolation of the nearest observations are used to replace one or two separate abnormal data. Wind direction anomalies and other variables with more than two consecutive anomalies are deleted.

    In summary, the QC procedure of DrIB observations developed in this study focuses on data classification for QC based on the full use of measurements. For oceanic variables, by increasing the weight ofin situobservations in the SST continuity test algorithm, the applicability of the algorithm is extended from tropical waters to higher latitudes, and the accuracy of outlier detection is improved.For meteorological variables, a QC method combining the moving Hampel filter and local anomaly detection is proposed. The median andMADare chosen as determination parameters to avoid interference by extreme values when applying the PauTa criterion. Moving detection performed in one-day steps is more consistent with the working characteristics of DrIB rather than using long-term statistical characteristics of fixed observations. The local anomaly detection method considers the actual physical process of meteorological variables, which effectively reduces the probability of misjudgment in identifying anomalies based on mathematical features only.

    3 Application and Validation of QC Procedure

    In this study, two sets of DrIB observations, namely,5840 and 8990, were used as examples to systematically analyze the QC performance and verify the methods’ effectiveness and advancement. The results showed that the QC scheme is applied well in the processing of DrIB observations.

    3.1 Data Description

    DrIB 5840 continued its drift and observation in the Kuroshio Extension from August 30, 2018, to January 1,2019. The sampling interval was 1 h. A total of 2909 records, including longitude, latitude, SST, air temperature,air pressure, wind speed, wind direction, relative humidity,and buoy operation records, were obtained. DrIB 8990 continued to drift and observe in eastern Japan from November 4, 2019, to November 6, 2019. The sampling interval was 10 min, and 251 pieces of data were collected with the same categories as DrIB 5840.

    3.2 Analysis of QC Results

    Table 2 shows the QC results of DrIB 5840 observations. The roll and pitch angles of the first six observation records were all zero, which was inconsistent with the swaying attitude of the buoy in seas. Fig.3(b) shows a steep descent at the beginning of the SST curve, and the significant difference in SST was more than 3℃. Therefore,these data were judged to be sent by the DrIB on land by mistake. Usually, DrIB functions well, and the observations are reliable and do not exceed the normal range. Thus, the range check focuses on the missing data stored in ?999.999 format. Figs.3 and 4(a) show the time series of the ob-served variables and drifting trajectory plotted using the original observations, respectively. The abrupt straight lines on the figures are precisely due to the missing data,that is, the anomalies detected by the range inspection, in which all observed variables were missed twice on December 18, 2018.

    Table 2 QC results of observations from DrIB 5840

    Fig.3 Time series of DrIB 5840 observations (before QC).

    After the basic check, the DrIB 5840 observations were classified for targeted QC. After the spike test, two abnormal data were screened out from the longitude and latitude time series, and the two spikes near 32?N, 144?E on the drifting trajectory were removed. The drifting trajectory after QC (Fig.4(b)) was clear compared with the previous observation, and the curve was naturally connected without abnormal spikes and easily visualized the dynamic changes of the observed sea areas. In addition, the QC performance was evident and intuitive. The SST observations were of high quality with a few anomalies, and seven records, which included both true anomalies and misjudgment results, were detected by the method proposed by the SVP data collection center (Fig.5(a)). According to the local diagram of SST QC results (Fig.5(b)), three anomalies from 2:00 to 4:00 on September 7, 2018 were consistent with the overall and local developing trends of the time series, which were misjudgment results. After adjusting the coefficients weights in the algorithm, three anomalies were detected and manually checked to confirm that the results were correct. The findings proved that the improved method effectively improved the accuracy of discrimination while expanding the applicability. For meteorological variables, with air pressure used as an example,the air pressure QC results were obtained by checking the SST (Fig.6(a)). The outliers were numerous and primarily distributed in patches, and all the observations from November 22 to 24, 2018 were abnormal, which was inconsistent with the facts. A total of 73 anomalies violating the daily variability of air pressure were detected using the moving Hampel filter, and the percentage of anomalies was 2.5%. By comparison, seven anomalies exceeding 10 hPa of air pressure at adjacent moments were detected from these data after the addition of local anomaly detection (Fig.6(b)). The percentage of anomalies was reduced to 0.24%, which indicates that as the physical processes significantly reduced the probability of misjudgment, the effect of QC was more guaranteed.

    Fig.4 Drifting trajectory of DrIB 5840.

    Fig.5 QC results of SST observed by DrIB 5840.

    Fig.6 QC results of air pressure observed by DrIB 5840.

    Fig.7 shows the DrIB 5840 observations after QC. Abrupt values and outliers out of the developing trends were identified and deleted or replaced by interpolation. The continuity and regularity of the time series were intact, the characteristics of the developing trends were evident, the data quality was significantly improved, and the data can be used for subsequent applications. The DrIB 8990, which had a short observation time, fewer observation data, and better quality, showed no anomaly in the detection by the basic QC. All anomalies in Table 3 came from the targeted QC. Given that the observation interval of DrIB 8990 was 10 min, the criteria of the continuity test for variables were adjusted adaptively. Fig.8 shows the comparison of the time series before and after QC, in which abnormal values were detected from air pressure, wind speed, SST, and relative humidity observations, and the sharp points on the curve were eliminated after QC.

    Table 3 QC results of observations from DrIB 8990

    Fig.7 Time series of DrIB 5840 observations (after QC).

    Fig.8 QC results of DrIB 8990 observations.

    3.3 Validation of QC Effect

    To better assess the effect of QC, we compared and evaluated the observations of DrIB 8990 before and after QC using the internationally recognized KEO buoy observations as the evaluation standard. The observation of DrIB 8990 covered the area of Kuroshio Extension and was adjacent to the KEO buoy (Fig. 9(a)). The KEO buoy was deployed in the southern area of the Kuroshio Extension, and the nominal location was 32.3?N, 144.6?E. Its observations were widely used in the research related to the Kuroshio Extension.

    Given the drifting characteristics of DrIB, spatiotemporal matching of the evaluated data was necessary. Fig.9(b)shows the distance variation between the DrIB 8990 and the KEO buoy with time. The paired data were the observations obtained at the same time within 10 km of the two buoys, that is, 164 records from 18:50 on November 4,2019, to 22:00 on November 5, 2019. A 10 km distance was considered because the KEO buoy observes in a circular scope with a diameter of more than 10 km, and the air-sea environment within its 10 km can be considered relatively consistent (Zhanget al., 2016; Andersonet al.,2020). As the KEO buoy did not provide 10 min observations of SST for the selected periods, the comparison of variables only included air pressure, wind speed, and relative humidity. In addition, DrIB and KEO buoys observe meteorological variables at different heights. The former observes relative humidity at approximately 2.8 m and air pressure and wind at approximately 3 m, and the latter pro-vides wind at a nominal height of 4 m and air pressure and relative humidity at a nominal height of 2.5 m. Therefore,the data must be unified to the same height first. By applying the COARE 3.0 formula and inputting variables, including SST, air pressure, air temperature, relative humidity, wind speed, and the corresponding observation heights,the meteorological field at the height of 10 m above the sea surface can be outputted, thus unifying the heights of the two sets of meteorological variables (Fairallet al.,2003). Given the lack of 10 min SST of the KEO buoy, the quality-controlled SST of DrIB 8990 was used to perform the height transformation of the KEO buoy observations instead.

    Fig.9 Position (a) and distance (b) of DrIB 8990 relative to the KEO buoy.

    Evaluation of the effect of QC mainly relied on two indicators, namely, data correlation and data difference, and the statistics used included correlation coefficientr, mean absolute error (AE), maximum absolute error (MAE), and root mean square error (RMSE), which are defined as follows.

    wherenis the number of paired data,biis the DrIB observation,kiis the KEO buoy observation,represents the arithmetic mean of the DrIB observations, anddenotes the arithmetic mean of the KEO buoy observations.The comparison results are shown in Fig.10 and Table 4.The solid curves in Fig.10 correspond to the left ordinate,and dashed and dotted curves refer to the right ordinate.

    Table 4 Comparison between the observations of DrIB 8990 and KEO buoy

    Fig.10 QC effect of DrIB 8990 observations.

    The results showed that the paired air pressure was the most consistent among the three variables. The correlation coefficient can intuitively reflect the degree of linear correlation between the two air pressures. The correlation coefficient between the original air pressure of DrIB 8990 and the KEO buoy was over 0.98, which was highly correlated at the 95% confidence level, and the changes in the two curves were almost the same. All three difference statistics were less than 0.5 hPa, which shows that the quality of DrIB observations reached the advanced international level. The quality of air pressure can only be slightly improved due to the high correlation and presence of a few anomalies. Ther, AE, andRMSEafter QC were optimized by about 10?3orders of magnitude, andMAEwas reduced by about 0.1 hPa. Contrary to the air pressure, the correlation coefficient between the relative humidity of DrIB 8990 and the KEO buoy was the smallest (approximately 0.4), and the difference statistics also showed weak data conformity. The relative humidity observed by the DrIB 8990 was approximately 3% higher than that of the KEO buoy on average, and the maximum exceeded 9%. Fig. 10(c)shows a few anomalies of relative humidity data detected by the QC procedure, and the data quality improved overall.The correlation coefficient increased, the deviation decreased, and the optimizations were mostly between 10?3– 10?2orders of magnitude. The time series of wind speed fluctuated strongly, but the developing trends of the two were generally consistent. Before QC, ther, AE, andRMSEwere approximately 0.77, 1 m s?1, and 1.2 m s?1, respectively. After QC, the extreme outliers out of the time series were eliminated. As a result, the correlation coefficient increased to greater than 0.8, andMAEwas reduced from approximately 5.6 m s?1to 2.8 m s?1. QC of the wind speed was the most effective, and the reliability of wind speed observations significantly improved.

    4 Conclusions

    Focusing on the air-sea interface, this study integrated various QC methods and established a scientific and practical QC procedure for the new drifting buoy observations in the order of an initial basic inspection and then targeted QC. The main conclusions and innovations are as follows.

    1) The QC procedure considers the swaying attitude of the buoy in the water, design observation indexes, and the error of the instrument itself, and it is more reasonable and more efficient than most QC methods that only rely on data.

    2) For different observational variables, respective methods were adopted to reduce the probability of misjudgment and significantly improve the effect of QC.

    3) The coefficient weight of the SST continuity check algorithm was optimized based on the differences in the sea area where the data came from. Thus, the algorithm was consistent with the characteristics of the actual air-sea environment, which expanded the application of the algorithm and effectively improved the accuracy of the QC.

    4) The innovative method combining moving Hampel filter and local anomaly detection ensured the reliability of QC results of meteorological variables in terms of mathematical statistics and actual physical processes.

    Although the evaluation of local and short-term observations was less reliable in demonstrating the effect of QC,overall, the outliers in the time series were effectively handled. After QC, the linear correlation between the observations from DrIB 8990 and KEO buoy as the control subjects increased, and the deviation decreased. The wind speed exhibited the best QC performance, the correlation coefficient before and after QC improved from around 0.77 to 0.82,andMAEwas reduced by approximately 2.8 m s?1. In addition, the quality of air pressure and relative humidity was optimized by 10?3– 10?2orders of magnitude.The results showed that the QC procedure for DrIB observations established in this study is scientific, effective,and practical, and it can provide references for the processing of observations from mobile platforms, promote the operational application of the new air-sea interface observation technologies, and offer data basis for marine environmental protection, ocean forecasting, air-sea interaction, and climate research. Marine observations are affected by many factors, such as the instruments themselves and changes in the ocean environment, and the data collected by DrIB are numerous and heterogeneous. Anomalous data also exhibit complex morphological characteristics. Therefore, finding a general QC method that can accurately identify and deal with all the anomalies is difficult. The QC procedure established in this study considers the circumstances that may lead to abnormalities in DrIB observations and has been experimentally verified to be feasible and practical. The QC procedure effectively detects independent and individual outliers in the time series.However, for continuous anomalous values, the QC performance needs further improvement, which will be the focus of our subsequent research.

    Acknowledgements

    This study is supported by the Natural Resources Development Special Fund Project of Jiangsu Province (No.JSZRHYKJ202009), the Taishan Scholar Funds (No. tsqn 201812022), the Fundamental Research Funds for the Central Universities (No. 202072001), the Guangxi Key Laboratory of Marine Disaster in the Beibu Gulf, Beibu Gulf University (No. 2021KF03), and the National Natural Science Foundation of China (No. 42176020). The authors are grateful for the financial support.

    在线免费观看不下载黄p国产| 看免费成人av毛片| 国产探花极品一区二区| 亚洲乱码一区二区免费版| 日本猛色少妇xxxxx猛交久久| 亚洲五月天丁香| 婷婷色av中文字幕| 亚洲国产精品国产精品| 亚洲色图av天堂| 午夜福利在线观看免费完整高清在| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美另类精品又又久久亚洲欧美| 一边摸一边抽搐一进一小说| 亚洲精品乱码久久久v下载方式| 亚洲人成网站高清观看| 国产精品一区二区性色av| 成人毛片a级毛片在线播放| 日韩人妻高清精品专区| 色综合站精品国产| 非洲黑人性xxxx精品又粗又长| 久久精品熟女亚洲av麻豆精品 | 熟妇人妻久久中文字幕3abv| 天堂中文最新版在线下载 | 亚洲成人av在线免费| 国产一区二区在线观看日韩| 亚洲国产精品成人综合色| 一级av片app| 久久亚洲国产成人精品v| 在线观看av片永久免费下载| 亚洲伊人久久精品综合 | 伦精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 亚洲自偷自拍三级| 免费av不卡在线播放| 成人毛片60女人毛片免费| 精品一区二区三区人妻视频| 成人毛片60女人毛片免费| 热99在线观看视频| 又粗又硬又长又爽又黄的视频| 久久综合国产亚洲精品| 国产在视频线精品| 男人狂女人下面高潮的视频| 国产v大片淫在线免费观看| 老师上课跳d突然被开到最大视频| 亚洲欧美精品综合久久99| 婷婷色av中文字幕| 内地一区二区视频在线| 午夜爱爱视频在线播放| 色视频www国产| 久久精品国产鲁丝片午夜精品| 亚洲精品一区蜜桃| 国产色婷婷99| 亚洲欧美成人精品一区二区| 日本熟妇午夜| 日本黄色视频三级网站网址| 久久久午夜欧美精品| 精品久久久久久电影网 | 午夜免费男女啪啪视频观看| 国产精品无大码| 亚洲精品国产成人久久av| 性插视频无遮挡在线免费观看| 欧美bdsm另类| 成人毛片60女人毛片免费| 欧美日韩在线观看h| 国产欧美另类精品又又久久亚洲欧美| 一级二级三级毛片免费看| av在线播放精品| 特大巨黑吊av在线直播| 九草在线视频观看| 亚洲自偷自拍三级| 久久人人爽人人片av| kizo精华| 在线观看一区二区三区| 日韩高清综合在线| 一级毛片aaaaaa免费看小| 成人毛片60女人毛片免费| 精品久久国产蜜桃| 超碰97精品在线观看| 欧美日本视频| 日本熟妇午夜| 国产91av在线免费观看| 国产免费又黄又爽又色| av专区在线播放| 精品酒店卫生间| 国产av码专区亚洲av| 男人舔奶头视频| 久久精品久久久久久噜噜老黄 | 国产中年淑女户外野战色| 欧美成人免费av一区二区三区| 国产一区有黄有色的免费视频 | 精品人妻一区二区三区麻豆| 欧美人与善性xxx| 波多野结衣高清无吗| 国产视频首页在线观看| 国产单亲对白刺激| 最近中文字幕2019免费版| 日日撸夜夜添| 日韩 亚洲 欧美在线| 少妇的逼水好多| 只有这里有精品99| 简卡轻食公司| videos熟女内射| 日本免费在线观看一区| 91久久精品国产一区二区三区| 国产精品国产三级国产专区5o | 一本久久精品| a级一级毛片免费在线观看| 亚洲精品日韩在线中文字幕| 69av精品久久久久久| 国产精品永久免费网站| 中文天堂在线官网| 激情 狠狠 欧美| 国产精品伦人一区二区| 精品熟女少妇av免费看| 国产精品国产三级国产av玫瑰| 一本一本综合久久| 国产伦一二天堂av在线观看| 成人国产麻豆网| 最后的刺客免费高清国语| 欧美97在线视频| 国产免费福利视频在线观看| 青春草亚洲视频在线观看| 极品教师在线视频| 身体一侧抽搐| 免费观看人在逋| 亚洲美女搞黄在线观看| 男的添女的下面高潮视频| 激情 狠狠 欧美| 亚洲国产成人一精品久久久| 91久久精品电影网| 69av精品久久久久久| 午夜久久久久精精品| 国产免费一级a男人的天堂| 国产视频首页在线观看| 又爽又黄a免费视频| 国产成人精品久久久久久| 亚洲国产欧美人成| 成人鲁丝片一二三区免费| 大香蕉久久网| 亚洲内射少妇av| 色视频www国产| 波野结衣二区三区在线| 91久久精品国产一区二区三区| 免费看日本二区| 久久久久久久久久久免费av| 日韩三级伦理在线观看| 蜜桃亚洲精品一区二区三区| 亚洲av日韩在线播放| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合久久99| 精品国产露脸久久av麻豆 | 国产高清有码在线观看视频| 国产亚洲av片在线观看秒播厂 | 亚洲欧美成人综合另类久久久 | 91aial.com中文字幕在线观看| 一二三四中文在线观看免费高清| 国产 一区 欧美 日韩| 日本免费一区二区三区高清不卡| 男女视频在线观看网站免费| 中文字幕av在线有码专区| www日本黄色视频网| 午夜老司机福利剧场| 国产精品电影一区二区三区| 免费播放大片免费观看视频在线观看 | 久久久国产成人精品二区| 国产伦一二天堂av在线观看| www.色视频.com| 亚洲欧美中文字幕日韩二区| 白带黄色成豆腐渣| 九色成人免费人妻av| 中文乱码字字幕精品一区二区三区 | 美女被艹到高潮喷水动态| 深夜a级毛片| 亚洲国产最新在线播放| 久久人妻av系列| 日韩av不卡免费在线播放| 69人妻影院| 尤物成人国产欧美一区二区三区| 亚洲真实伦在线观看| 99久久无色码亚洲精品果冻| 久久久欧美国产精品| 免费播放大片免费观看视频在线观看 | 色尼玛亚洲综合影院| 午夜福利网站1000一区二区三区| 男女视频在线观看网站免费| 国产午夜精品论理片| 亚洲伊人久久精品综合 | 亚洲最大成人中文| 色哟哟·www| 国产v大片淫在线免费观看| 日本-黄色视频高清免费观看| 一级毛片电影观看 | 国产精品综合久久久久久久免费| 色播亚洲综合网| 青春草视频在线免费观看| 亚洲国产精品sss在线观看| 插阴视频在线观看视频| 成年版毛片免费区| 麻豆国产97在线/欧美| 建设人人有责人人尽责人人享有的 | 日韩制服骚丝袜av| 成年免费大片在线观看| 国产老妇女一区| 嫩草影院精品99| 亚洲图色成人| 亚洲乱码一区二区免费版| 精品免费久久久久久久清纯| 国产v大片淫在线免费观看| 国产午夜精品一二区理论片| 亚洲精品aⅴ在线观看| 99久久中文字幕三级久久日本| 一二三四中文在线观看免费高清| 国产精品久久久久久久久免| 日韩国内少妇激情av| 亚洲欧美精品专区久久| 乱人视频在线观看| a级毛色黄片| 久久久欧美国产精品| 男人舔奶头视频| 噜噜噜噜噜久久久久久91| 午夜激情欧美在线| 18+在线观看网站| 又粗又爽又猛毛片免费看| 精品久久久久久久久av| 亚洲国产色片| 国产色婷婷99| av在线蜜桃| 最近最新中文字幕免费大全7| 一区二区三区乱码不卡18| 人人妻人人看人人澡| 久久99热这里只有精品18| 日韩 亚洲 欧美在线| 亚洲成色77777| a级毛片免费高清观看在线播放| 高清毛片免费看| 国产在线一区二区三区精 | 久久久成人免费电影| 亚洲欧美日韩高清专用| a级一级毛片免费在线观看| 2021天堂中文幕一二区在线观| 一个人免费在线观看电影| 蜜桃久久精品国产亚洲av| 你懂的网址亚洲精品在线观看 | 国产午夜精品一二区理论片| 免费看a级黄色片| 国产黄片视频在线免费观看| 亚洲不卡免费看| 亚洲第一区二区三区不卡| 国产成人免费观看mmmm| 亚洲欧美精品专区久久| 婷婷色麻豆天堂久久 | a级毛色黄片| 中国美白少妇内射xxxbb| 国产一区二区三区av在线| 18禁在线无遮挡免费观看视频| 中文字幕免费在线视频6| 日韩欧美国产在线观看| 22中文网久久字幕| 18禁在线无遮挡免费观看视频| 午夜精品一区二区三区免费看| 边亲边吃奶的免费视频| 久久韩国三级中文字幕| 国产成人一区二区在线| av又黄又爽大尺度在线免费看 | 日韩欧美 国产精品| 久久99热这里只频精品6学生 | 欧美日韩在线观看h| 美女cb高潮喷水在线观看| 亚洲av中文字字幕乱码综合| 久久久久久久久大av| 国产视频内射| ponron亚洲| 欧美成人a在线观看| 亚洲人成网站高清观看| 亚洲精品aⅴ在线观看| 成人漫画全彩无遮挡| 女的被弄到高潮叫床怎么办| 夫妻性生交免费视频一级片| 一级二级三级毛片免费看| 国产高清三级在线| 搡女人真爽免费视频火全软件| 三级毛片av免费| 最近中文字幕2019免费版| 观看美女的网站| 欧美97在线视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲人与动物交配视频| 日本黄色视频三级网站网址| 一个人观看的视频www高清免费观看| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 69人妻影院| 亚洲丝袜综合中文字幕| 国产亚洲91精品色在线| 乱码一卡2卡4卡精品| 日本与韩国留学比较| 精品欧美国产一区二区三| 色综合色国产| 少妇裸体淫交视频免费看高清| 91久久精品国产一区二区成人| 夜夜爽夜夜爽视频| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 日本五十路高清| 国产淫片久久久久久久久| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 日韩一区二区视频免费看| 久久久成人免费电影| 女人被狂操c到高潮| 91久久精品电影网| 最近手机中文字幕大全| 免费av不卡在线播放| 久久综合国产亚洲精品| kizo精华| 嫩草影院精品99| av专区在线播放| 国产精品,欧美在线| 中文字幕亚洲精品专区| 亚洲欧美日韩无卡精品| 久久久久性生活片| 天天躁夜夜躁狠狠久久av| 欧美xxxx黑人xx丫x性爽| 亚洲在久久综合| 爱豆传媒免费全集在线观看| 级片在线观看| 久久99热这里只有精品18| 成人鲁丝片一二三区免费| 欧美一区二区精品小视频在线| 最后的刺客免费高清国语| 天天一区二区日本电影三级| 欧美性感艳星| 久久亚洲国产成人精品v| 国产av不卡久久| 日韩欧美国产在线观看| 高清日韩中文字幕在线| 国产在线一区二区三区精 | 亚洲成人久久爱视频| 国产精品久久久久久久电影| 能在线免费看毛片的网站| 色综合色国产| 一个人观看的视频www高清免费观看| 搞女人的毛片| 视频中文字幕在线观看| 国产成人一区二区在线| 一个人免费在线观看电影| 看黄色毛片网站| 长腿黑丝高跟| 成人综合一区亚洲| 精品午夜福利在线看| 国产乱人偷精品视频| 精品人妻视频免费看| 在线免费观看不下载黄p国产| 国产免费福利视频在线观看| 亚洲四区av| 国语自产精品视频在线第100页| 久久99热这里只频精品6学生 | 色综合站精品国产| 天堂av国产一区二区熟女人妻| 国产成年人精品一区二区| 国产在线一区二区三区精 | 久久精品国产99精品国产亚洲性色| 欧美又色又爽又黄视频| a级毛色黄片| 午夜福利在线在线| 久久久久性生活片| 国产一区二区三区av在线| 日本黄色片子视频| 亚洲成人精品中文字幕电影| 亚洲一区高清亚洲精品| 村上凉子中文字幕在线| 亚洲一区高清亚洲精品| 五月伊人婷婷丁香| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 免费黄网站久久成人精品| 精华霜和精华液先用哪个| 蜜桃亚洲精品一区二区三区| 亚洲成人中文字幕在线播放| 晚上一个人看的免费电影| 91av网一区二区| 成人鲁丝片一二三区免费| 亚洲成人精品中文字幕电影| 老司机福利观看| 久久国产乱子免费精品| 欧美成人免费av一区二区三区| www日本黄色视频网| 男人狂女人下面高潮的视频| 五月伊人婷婷丁香| av线在线观看网站| 免费av毛片视频| 99久久人妻综合| 热99在线观看视频| 又爽又黄a免费视频| 99久久精品热视频| 日韩欧美国产在线观看| 亚洲av福利一区| 99热精品在线国产| 久久久久久久午夜电影| 18+在线观看网站| 亚洲精品国产成人久久av| 99久久精品一区二区三区| 高清av免费在线| 国产伦精品一区二区三区四那| 精品久久久久久久久av| 小蜜桃在线观看免费完整版高清| 免费看美女性在线毛片视频| 两性午夜刺激爽爽歪歪视频在线观看| av在线播放精品| 亚洲国产精品合色在线| 天天躁日日操中文字幕| 人体艺术视频欧美日本| 国产老妇伦熟女老妇高清| 午夜福利网站1000一区二区三区| 国产精品美女特级片免费视频播放器| 久久精品人妻少妇| 国产精品一二三区在线看| 亚洲精品,欧美精品| 日本黄色片子视频| 午夜a级毛片| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄 | eeuss影院久久| 老司机福利观看| 婷婷六月久久综合丁香| 亚洲av成人精品一区久久| av又黄又爽大尺度在线免费看 | 日韩欧美三级三区| 大又大粗又爽又黄少妇毛片口| 欧美精品国产亚洲| 波多野结衣巨乳人妻| 一级二级三级毛片免费看| 亚洲最大成人中文| 性色avwww在线观看| 亚洲国产高清在线一区二区三| 99久国产av精品国产电影| 三级国产精品欧美在线观看| 51国产日韩欧美| 国产精品国产三级专区第一集| 青春草国产在线视频| 亚洲成人精品中文字幕电影| 高清毛片免费看| 97超视频在线观看视频| av专区在线播放| av在线天堂中文字幕| 一个人看视频在线观看www免费| 99久久无色码亚洲精品果冻| 国产午夜福利久久久久久| 亚洲,欧美,日韩| 日本免费a在线| www.av在线官网国产| 插阴视频在线观看视频| 男女啪啪激烈高潮av片| 特大巨黑吊av在线直播| 日日撸夜夜添| 国产午夜精品一二区理论片| 嫩草影院新地址| 亚洲色图av天堂| 男插女下体视频免费在线播放| 国内少妇人妻偷人精品xxx网站| 少妇熟女欧美另类| 黄片wwwwww| 99久久成人亚洲精品观看| 五月玫瑰六月丁香| 色5月婷婷丁香| 国产精品久久久久久av不卡| 纵有疾风起免费观看全集完整版 | 草草在线视频免费看| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区久久| 日韩欧美在线乱码| 成人午夜高清在线视频| 全区人妻精品视频| 最近2019中文字幕mv第一页| 亚洲精品日韩av片在线观看| 国产在视频线在精品| 又粗又硬又长又爽又黄的视频| 赤兔流量卡办理| 国产黄色小视频在线观看| .国产精品久久| 在现免费观看毛片| 国产av在哪里看| 欧美极品一区二区三区四区| 精品午夜福利在线看| 亚洲国产日韩欧美精品在线观看| 亚洲国产精品成人综合色| 永久网站在线| 人妻系列 视频| 看非洲黑人一级黄片| 国产又黄又爽又无遮挡在线| 久久久精品欧美日韩精品| 亚洲国产精品合色在线| av专区在线播放| 69人妻影院| 久久精品91蜜桃| 久久精品夜色国产| 欧美bdsm另类| 尤物成人国产欧美一区二区三区| 国国产精品蜜臀av免费| 少妇熟女aⅴ在线视频| 一级毛片我不卡| 国产爱豆传媒在线观看| 极品教师在线视频| 免费不卡的大黄色大毛片视频在线观看 | 不卡视频在线观看欧美| 2021少妇久久久久久久久久久| 亚洲真实伦在线观看| 九九爱精品视频在线观看| 国产av不卡久久| 日韩国内少妇激情av| 18禁裸乳无遮挡免费网站照片| 午夜福利高清视频| 99久国产av精品| 精品午夜福利在线看| 久久久久久久久久成人| 午夜免费激情av| 亚洲av中文av极速乱| 中文字幕免费在线视频6| 国产国拍精品亚洲av在线观看| 亚洲性久久影院| av国产久精品久网站免费入址| 久久久久久九九精品二区国产| 一级av片app| 免费电影在线观看免费观看| 国产精品人妻久久久久久| 午夜福利在线观看吧| 国内少妇人妻偷人精品xxx网站| 老司机福利观看| 亚洲av成人av| 黄色配什么色好看| 简卡轻食公司| 亚洲美女视频黄频| 看片在线看免费视频| 久久6这里有精品| 国内精品美女久久久久久| 日本黄色片子视频| 亚洲自偷自拍三级| 天天一区二区日本电影三级| 国产淫片久久久久久久久| 欧美日韩精品成人综合77777| 日韩视频在线欧美| 精品免费久久久久久久清纯| 精品熟女少妇av免费看| 久久精品久久久久久久性| 久久久久久久久大av| 91精品一卡2卡3卡4卡| 性色avwww在线观看| 人人妻人人澡欧美一区二区| 男女那种视频在线观看| 免费看美女性在线毛片视频| 永久网站在线| 中文字幕制服av| 久久精品影院6| 国产探花极品一区二区| 亚洲av免费高清在线观看| 精华霜和精华液先用哪个| 成人一区二区视频在线观看| av视频在线观看入口| 永久免费av网站大全| 18禁动态无遮挡网站| 男女边吃奶边做爰视频| 看十八女毛片水多多多| 少妇人妻精品综合一区二区| 精品一区二区免费观看| 国产精品人妻久久久久久| 午夜免费激情av| 久久国产乱子免费精品| 可以在线观看毛片的网站| 日韩av在线大香蕉| 在线观看av片永久免费下载| 亚洲不卡免费看| 欧美性猛交黑人性爽| .国产精品久久| 色视频www国产| 免费不卡的大黄色大毛片视频在线观看 | 在线a可以看的网站| 久久人人爽人人片av| 亚洲综合色惰| 91狼人影院| 大香蕉97超碰在线| 日本午夜av视频| 亚洲丝袜综合中文字幕| 欧美xxxx黑人xx丫x性爽| 村上凉子中文字幕在线| 青春草视频在线免费观看| 国产精品久久久久久av不卡| 人人妻人人澡欧美一区二区| 夜夜爽夜夜爽视频| 天天躁日日操中文字幕| 国产三级在线视频| 一级爰片在线观看| 两个人视频免费观看高清| 又粗又爽又猛毛片免费看| av黄色大香蕉| 超碰av人人做人人爽久久| 91久久精品电影网| 精品久久久久久久久久久久久| 成人漫画全彩无遮挡| 成人毛片60女人毛片免费| 别揉我奶头 嗯啊视频| 免费大片18禁| 国产免费一级a男人的天堂| 精品99又大又爽又粗少妇毛片| 男女国产视频网站| 国产av码专区亚洲av| 亚洲精品456在线播放app| 国产成人福利小说| 人体艺术视频欧美日本| 欧美成人免费av一区二区三区| 欧美性猛交黑人性爽| 国产极品精品免费视频能看的| 九九爱精品视频在线观看| 超碰av人人做人人爽久久| 久久热精品热| 男女那种视频在线观看|