Top Menu

Jump to content
  • downscale
    • View all projects
Home
    • Work packages
    • News
    • Getting started
    • Introduction video
    • Welcome to OpenProject

      Get an overview

      Get a quick overview of project management and team collaboration with OpenProject.
      You can restart this video from the help menu

    • Help and support
    • Upgrade to Enterprise Edition
    • User guides
    • Videos
    • Shortcuts
    • Community forum
    • Professional support

    • Additional resources
    • Data privacy and security policy
    • OpenProject website
    • Security alerts / Newsletter
    • OpenProject blog
    • Release notes
    • Report a bug
    • Development roadmap
    • Add and edit translations
    • API documentation
  • Sign in
      Forgot your password?

Side Menu

  • Overview
  • Activity
  • Work packages
  • News
  • Repository
  • You are here:Wiki main page
    • Table of contents
      • Hierarchy leafGravimetric humidity bias on 2008-04-07, 05 UTC
      • Expanded. Click to collapseCollapsed. Click to showWiki main page
        • Expanded. Click to collapseCollapsed. Click to showCoarse data
          • Hierarchy leafCoarse external data
          • Hierarchy leafIntroduction of new data
          • Hierarchy leafTemperature
        • Hierarchy leafDependencies
        • Expanded. Click to collapseCollapsed. Click to showExchange of variable fields
          • Hierarchy leafnamcouple
          • Expanded. Click to collapseCollapsed. Click to showSchomburg scheme
            • Expanded. Click to collapseCollapsed. Click to showSchomburg rules
              • Hierarchy leafSchomburg rule
        • Expanded. Click to collapseCollapsed. Click to showFine data
          • Hierarchy leafEstimated resolutions
          • Expanded. Click to collapseCollapsed. Click to showFine external data from the Climate Limited-area Modeling community
            • Hierarchy leafDivided by 7
            • Hierarchy leafEmulation of the test data resolution
          • Hierarchy leafSRTM data
        • Hierarchy leafIntroduction
        • Hierarchy leafLiterature
        • Expanded. Click to collapseCollapsed. Click to showModel1
          • Hierarchy leafDeprecated dummy data
          • Hierarchy leafDummy coordinates
          • Hierarchy leafProcessing
          • Expanded. Click to collapseCollapsed. Click to showRead_grib1
            • Hierarchy leafEccodes_handler
        • Hierarchy leafNaming schemes
        • Hierarchy leafOASIS3
        • Hierarchy leafPrivacy policy
        • Hierarchy leafTitle page
        • Expanded. Click to collapseCollapsed. Click to showValidation
          • Expanded. Click to collapseCollapsed. Click to showPreservation of the sub-scale average
            • Expanded. Click to collapseCollapsed. Click to showBias anecdote for gravimetric humidity
              • Hierarchy leafGravimetric humidity bias on 2008-04-04, 18 UTC
            • Expanded. Click to collapseCollapsed. Click to showBias anecdote for surface pressure
              • Hierarchy leafSurface pressure bias on 2008-04-04, 18UTC
              • Hierarchy leafSurface pressure bias on 2008-04-07, 05UTC
            • Expanded. Click to collapseCollapsed. Click to showBias anecdote for temperature
              • Hierarchy leafTemperature bias on 2008-04-07, 05UTC
              • Hierarchy leafTemperature bias on 2008-04-04, 18UTC
            • Hierarchy leafgrib2bin.pl
          • Hierarchy leafSpatial average and variance before and after downscaling
        • Expanded. Click to collapseCollapsed. Click to showWhat the Community Land Model requires
          • Hierarchy leafsva_DWD_forcing.ncl
You are here:
  • Wiki main page
  • Validation
  • Spatial average and variance before and after downscaling

Content

Spatial average and variance before and after downscaling

  • More
    • Table of contents

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

  • Inference of variance
  • Interpretation
  • A subset of 6 days

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.

Relative differences in mean
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.

relative differences in estimated variance
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.

The relative increase in variance by downscaling
Fig. 3: The relative increase of the estimator variance by downscaling.

The relative change in mean by downscaling
Fig. 4: The absolute value of the relative increase of the expectation value by downscaling.

Loading...