Content
This validation is based on run over 1464 time steps, the 61 days of January … February 2008. In total, 2008 consists of 8784 steps.
Table of contents
To put all variables in one plot, they are processed into relative increases c. The expectation value and variance are inferred firstly from the coarse data, yielding the reference r, and then from the downscaled data, giving f. Each 3-tuple of time step, variable and moment gets its own relative increase:
c= (f-r)/r
The reference r sometimes disappears, most commonly for solar radiation at night. The denominator then defaults to the temporal minimum of the variable excluding 0. The substitution is harmless given that the algorithm always downscales coarse fields of all 0 successfully to fine fields of all 0.
Fig. 1: The absolute value of the relative increase of the average by downscaling. Around c=1, ‘:’ symbolizes that on that time step and for that variable, there were no discernable differences in mean. R assesses its minimal discernable difference between numbers to be 2.220446e-16, given the ideal circumstances.
Fig. 2: The relative increase of the estimated variance by downscaling.
Tab. 1: The temporal range of c for the average
Variable | Minimum | Maximum |
---|---|---|
ALB_RAD | –9.96495743055829e-14 | 1.04995407175811e-13 |
ASOB_S | –1.14691715993802e-13 | 1.15130955830844e-13 |
ATHB_S | –1.6758300767169e-13 | 1.48814500692525e-13 |
PS | 1.39850203608734e-04 | 1.49032901938406e-04 |
T | –4.50791187516156e-05 | 8.44811036090759e-05 |
TD_2M | –1.07802549282895e-13 | 1.04752430230132e-13 |
T_G | –8.27152248163247e-07 | 1.21865720882416e-06 |
TOT_PREC | –7.01978457635899e-13 | 7.76570361665211e-13 |
U_10M | –1.19296404199416e-11 | 6.48444302264809e-13 |
V_10M | –3.35985657877883e-12 | 2.67334701488226e-12 |
Tab. 2: The temporal range of c for the estimated variance
Variable | Minimum | Maximum |
---|---|---|
ALB_RAD | –9.96495743055829e-14 | 1.04995407175811e-13 |
ASOB_S | –1.14691715993802e-13 | 1.15130955830844e-13 |
ATHB_S | –1.6758300767169e-13 | 1.48814500692525e-13 |
PS | 1.39850203608734e-04 | 1.49032901938406e-04 |
T | –4.50791187516156e-05 | 8.44811036090759e-05 |
TD_2M | –1.07802549282895e-13 | 1.04752430230132e-13 |
T_G | –8.27152248163247e-07 | 1.21865720882416e-06 |
TOT_PREC | –7.01978457635899e-13 | 7.76570361665211e-13 |
U_10M | –1.19296404199416e-11 | 6.48444302264809e-13 |
V_10M | –3.35985657877883e-12 | 2.67334701488226e-12 |
The validation codes source:plotting/verif.R and source:plotting/verif.F90 are independet of the algorithm.
Inference of variance
The 2nd inferred moment is the centralized 2nd moment. As the expectation value was also inferred and not given, the Bessel correction removes the bias from the estimator. Removing bias is also the reason not to infer the standard deviation, which could share a plot with the mean.
Interpretation
The increase in mean is insignificant. R’s minimal difference is just barely smaller for than the increase for most variables. The others are too small to impact any calculation by 3 orders of magnitude. By this standard, downscale succeeded.
The increase in variance is insignificant for solar and terrestial radiation. Downscaling in general leads to higher variance due to additional fine-scale variability. It usually stems from the spline interpolation, which generates much too smooth field compared to observations. The more advanced Schomburg rules increase the variance by a larger margin, as exhibited by pressure. Temperature appears to shift between these two groups of increases. The most opportune conjecture points to the switch in its Schomburg rule, which turns off the advanced downscaling in case of a high temperature lapse rate.
A subset of 6 days
This validation is based on run over 166 time steps, i.e. 166h, i.e. 6 days and 22h. The case can be made that this does not suffice as a sample for the entire year 2008 with its 8784 time steps.
Fig. 3: The relative increase of the estimator variance by downscaling.
Fig. 4: The absolute value of the relative increase of the expectation value by downscaling.