Quality check

Quality check procedure

    Quality Check - Temporal aggregation - FAQ - Bibliography

  Prior to the comparison of a dataset to reference data (validation protocol), you need to be sure that these reference data are of sufficiently good quality to serve as reference.
This page aims at giving a large insight of the procedure adopted by MINES ParisTech and Transvalor in term of quality check of ground station measurements. In particular, the first section is dedicated to the quality control steps applied onto the 1 minute Baseline Surface Radiation Network station data. Please refer to the selected publications for more information.

Quality Check (full description for sub-hourly values)

The quality check procedure is derived from the one published in the deliverable 3.2 on data harmonization provided within the framework of the European FP7 ENDORSE project. The table available page 60, summarizes the quality check procedure for all the meteorological components. Please note that this summarized quality check procedure is also well described in the poster associated with the paper of Espinar et al. 2012. All the differences between our implementation and this document are highlighted in orange.

NB: This quality check procedure described below is for sub-hourly measurements, which is typically the procedure adopted onto the 1 min BSRN station data. If you need to handle ground measurements with other time steps (such as hourly, or daily values), you need to switch to the adequate column in the above p60 table.

Let's define GHIGRD, DHIGRD and BNIGRD respectively the sub-hourly GHI, DHI and BNI data for the ground station. SZA is the Solar Zénithal angle and SEA the Solar Elevation angle (please note that SEA=90°-SZA). Let's name GHITOA and BNITOA respectively the Global and the Direct component at the Top of Atmosphere.

Step Quality Check Equation
Step 1 - Night
  • BNI night values are set to 0
  • GHI and DHI dark night values are set to 0
  • Sunrise and sunset: if GHI orDHI value is negative, set it to 0
  • for SEA , BNIGRD=0
  • for SEA -3°, GHIGRD=0 and DHIGRD=0
  • for -3° SEA ≤ 0° & GHIGRD < 0, GHIGRD=0 (idem for DHI)
Step 2 - Missing values
  • Set the missing values to NaN
if GHIGRD is missing, GHIGRD=NaN (idem for DHI and BNI)
Step 3 - Outliers
  • Discard the values which exceed "Extremely Rare Limits (ERL)" and "Physical Possible Limits (PPL)" criteria.

 

NB: the maximum step for two following GHIGRD measures of 1000 Wm-2 hasn't been implemented.

 

NB2: All the GHIGRD below 10 W/m² and all the BHIGRD below 4 W/m² are discarded.
NB3: Pour les validation reports, et pour PT15M, Lucien recommandait 30 W/m2

 

 

GHI:

  • PPL: 0.03 GHITOA < GHIGRD < min(1.2 BNITOA,1.5 BNITOA cos(SZA)1.2+100)
  • ERL: 0.03 GHITOA < GHIGRD < 1.2 BNITOA cos(SZA)1.2+50

BNI:

  • PPL: 0 < BNIGRD < BNITOA
  • ERL: 0 < BNIGRD < 0.95BNITOAcos(SZA)0.2+10

DHI:

  • PPL: 0.03 GHITOA < DHIGRD < min(0.8 BNITOA,0.95 BNITOA cos(SZA)1.2+50)
  • ERL: 0.03 GHITOA < DHIGRD < 0.75 BNITOA cos(SZA)1.2+30
Step 4 - Consistency
  • If the ground station provides several components of radiation, apply the consistency checks (cross-comparisons)

For 2 components (GHIGRD and DHIGRD)

For GHIGRD > 50 Wm-2

  • for SZA < 75°, DHIGRD/GHIGRD < 1.05
  • for 75° < SZA < 93°, DHIGRD/GHIGRD < 1.10

For 3 components (GHIGRD, DHIGRD and BNIGRD): compute first BHIGRD* from BNIGRD, and then compute GHIGRD*=BHIGRD*+DHIGRD and compare it to GHIGRD.

 For GHIGRD* > 50 Wm-2

  • for SZA < 75°, 0.92 GHIGRD/GHIGRD* ≤ 1.08
  • for 75° < SZA < 93°, 0.85 GHIGRD/GHIGRD* ≤ 1.15
 

Top of page

Temporal aggregation procedure

The below temporal aggregation procedure corresponds to the comparison of a HelioClim-3 (HC3) dataset against 1 min ground station measurements (typically BSRN).

Time step Steps
1 MIN
  1. Apply the Quality Check procedure (see above section)
from 1 to 15 min
  1. Generate the 15 min data if at least 85% of the 1 min slots are available, by applying an "intelligente interpolation" which takes into account the sun position at each instant.
15 MIN
  1. Set 15 min HC3 to NaN where 15 min GRD is NaN, and reverse.
  2. Compute the stats at 15 min
from 15 min to hourly data
  1. From now on, always perform partial sums. Generate hourly values by summing up the 15 min irradiation slots if at least 75% of the slots are available.
HOUR
  1. Compute the stats at the hourly time step
from hourly to daily values
  1. Generate the partial daily values by summing up the 15 min values if at least 65% of the 15 min values are available.
DAY
  1. Compute the stats at the daily time step
from daily to monthly data
  1. Generate the partial monthly values by summing up the 15 min values if at least 50% of the of the 15 min values are available.
MONTH
  1. Compute the stats at the monthly time step
 

Top of page

Discussion: If the time step of your ground station measurements or of your estimation is different, this procedure needs to be adapted.
For instance, if we collected hourly in-situ measurements to carry out a comparison with HC3, we must generate hourly HC3 values from the 15 min HC3 values using an intelligent interpolation like in step 2.

FAQ

I_ "What are 'partial sums'? Why did you opt for partial sums instead of averages?"

If you apply any kind of averaging on measurements, you introduce an information which doesn't come from a real measurement. That is why we decided to avoid the distorsion of the values by summing up only the available values. This leads to partial sums for the day or the month. This means also that the resulting value is not representive of an actual daily or monthly value.

II_ "How did you choose the thresholds in your temporal aggregation procedure?"

Originally, we didn't set any threshold, and created aggregated values even if only one slot was available. Even if partial sums are not representative of, for instance, a daily or a monthly value, we decided that at least 75% of the 15 min slots were necessary to generate the hourly value, and 75% of the 15 min data to generate the daily or monthly value.

NB: discussion for the hour from 15 min slots: when the sun is rising, you could have less than 4 valid 15 min slots to generate the hourly value. As a consequence, the 75% threshold accepts 1 available slot over 1 valid slot, 2 over 2, 3 over 3 and 4 over 4, plus 3 over 4 valid slots. The first discarded is 2 available over 3 valid slots.

III_ "Why do you sum up 15 min slots to generate the daily and monthly values?"

After several tests trying to sum up the hourly values to generate the daily value, we realized that two many days were discarded. Summing up the 15 min slots returns more valid days for the comparisons.

IV_ "Why did you discard all the GHIGRD values below 50 Wm-2?"

This decision was motivated by the fact that most instruments present a high uncertainty in low radiation conditions. Moreover, as the consistency checks are not applied for values below this threshold, we finally decided to discard the values.

Top of page

Bibliography Quality

Bibliography Quality

(In alphabetic and reverse chronological order)

Espinar et al. 2012 Espinar B., P. Blanc, L. Wald, C. Hoyer-Klick, M. Schroedter Homscheidt, and T. Wanderer 2012. "Controlling the quality of measurements of meteorological variables and solar radiation. From sub-hourly to monthly average time periods" EGU General Assembly 2012, Apr 2012, European Geosciences Union 2012, Vienna (Austria) 22–27 April 2012. Poster . Quality check

Content: Poster recapitulative of QC procedures

ISO Guide 1995 ISO Guide to the Expression of Uncertainty in Measurement: first edition, International Organization for Standardization, Geneva, Switzerland, 1995 Validation protocol

 

Ohmura et al. 1998 Ohmura, A., Gilgen, H., Hegner, H., Mueller, G., Wild, M., Dutton, E. G., Forgan, B., Froelich, C., Philipona, R., Heimo, A., Koenig-Langlo, G., McArthur, B., Pinker, R., Whitlock, C. H., and Dehne, K.: Baseline Surface Radiation Network (BSRN/WCRP): New precision radiometry for climate research. B. Am. Meteorol. Soc. 1998, 79, 2115–2136 BSRN

Content: Table 3 p 2119 gives the uncertainty on the GHI, DHI and BHI components measured by the BSRN stations: +/-5 W/m² for the GHI and the DHI components, and +/-2 W/m² for the BHI. This is why a threshold of respectively 10 W/m² and 4 W/m² is applied on the measurements in the adopted quality check procedure.

Roesch et al. 2011 Roesch, A., Wild, M., Ohmura, A., Dutton, E.G., Long, C.N., and Zhang, T., 2011. 'Assessment of BSRN radiation records for the computation of monthly means', Atmos. Meas. Tech., 4, 339–354, doi:10.5194/amt-4-339-2011, 2011. Corrigendum, Atmos. Meas. Tech., 4, 973-973, doi:10.5194/amt-4-973-2011. BSRN

 

WMO 2008 WMO: Guide to meteorological instruments and methods of observation, World Meteorological Organization, WMO - No 8, 7th Edn., Geneva, Switzerland, 2008. Quality check recommen-dations

 

WMO 1981 WMO, 1981. Technical Note N° 172, WMO-No. 554, World Meteorological Organization, Geneva, Switzerland, pp. 121-123. Validation protocol

Top of page