6 jaar aan data van GFS DMO's en hun correcties

Er is weer een extra jaar aan data het rekenmodel in gegaan en daarmee komt het totaal nu op 6 jaar.

Om niet elk jaar het zelfde te hoeven schrijven, verwijs ik naar het overzicht van vorig jaar, waar ik heb uitgelegd hoe de input van het rekenmodel tot stand komt.

Zie: 5 jaar aan data van GFS DMO's en hun correcties

Al die data gaat in een groot totaal overzicht, maar wordt tevens gerangschikt per maand.

In onderstaande overzichten (afbeelding 1 en 2) staan alle afwijkingen van de aflopen 6 jaar.

Afbeelding 1:

Afbeelding 2: 

Onderaan bovenstaand overzicht zijn de gemiddelde afwijkingen te zien over de gehele dataset.

Voor de Tn is dat:

- een gemiddelde positieve afwijking van 2,2°C

- een gemiddelde negatieve afwijking van -1,1°C

- een gemiddelde afwijking van 1,1°C

Voor de Tx is dat:

- een gemiddelde positieve afwijking van 1,1°C

- een gemiddelde negatieve afwijking van -1,6°C

- een gemiddelde afwijking van -0,9°C

Samengevat geeft dat aan dat GFS gemiddeld de Tn 1,1°C te hoog inschat en de Tx 0,9°C te laat inschat.

Omdat GFS bij stralingscondities (weinig/geen wind (of wind die weg valt) en geen bewolking of opklaringen die ontstaan) de plank qua Tn behoorlijk mis slaat, heb ik daar een stralingscorrectie voor bedacht.

Alle afwijkingen van de daadwerkelijke Tn ten opzichten van de gemiddelde DMO die meer dan 4°C zijn, tellen mee voor deze stralingscorrectie.

Op een totale dataset van 2191 dagen, is dit 214x gebeurd, wat neerkomt op een percentage van 9,77

De gemiddelde afwijking gedurende die 214x bedraagt 4,9°C zoals in onderstaand overzicht te zien is:

Omdat de ene maand de andere niet is, heb ik dit jaar voor de stralingscondities een opsplitsing per maand gemaakt.

In dat overzicht is per maand de zien de maximaal opgetreden stralingscorrectie, de gemiddelde stralingscorrectie, aantal keer dat het in een maand is opgetreden en dit vergeleken met het aantal keer dat het had kunnen optreden.

Dat geeft dan een percentage dat de statistische kans uitdrukt dat er in een maand een stralingscorrectie nodig is.

Met 6 jaar aan data is die kans in juni het kleinst (1,67%) en in september het grootst (20,56%).

De gemiddeld grootste correctie is nodig in december (5,3°C) terwijl deze in juni het kleinst is (4,2°C).

In de laatste kolom is zichtbaar wanneer een stralingscorrectie voor het laatst is opgetreden.

Daaruit blijkt dat er in 2023 geen enkele correctie was in januari, mei en oktober en dat de laatste correctie tot nu toe (data is bijgewerkt tot en met 30-11-2023) was op 20 november van dit jaar.

Tot op heden corrigeer ik de gemiddelde DMO met de gevonden gemiddelde afwijkingen, maar uiteindelijk is het mijn doel om correcties per seizoen of misschien zelfs per maand te gaan gebruiken.

Om inzichtelijk te krijgen of dat een meerwaarde zou hebben, heb ik de afwijkingen per seizoen én per maand gerangschikt.

Voor de Tn ziet dat er als volgt uit (zie afbeelding 3 en 4):

Afbeelding 3:

Afbeelding 4:

Dit levert een interessant beeld op.

Als we nog even in gedachten nemen dat de gemiddelde afwijking voor de Tn 1,1°C bedraagt en je dit dan vergelijkt met allereerst de seizoenen, dan zie je dat de lente (1,1°C) en zomer (1,1°C) exact gelijk zijn aan die gemiddelde afwijking.

Voor die seizoenen zou het vooralsnog geen toegevoegde waarde hebben om een andere correctie te gaan gebruiken dan de gemiddelde afwijking over de hele dataset.

Kijk je echter naar de winter (0,8°C) en herfst (1,4°C) dan zou dit weldegelijk een toegevoegde waarden hebben.

Bekijk je vervolgens binnen een seizoen naar de afwijkingen per maand dan valt  in de lente op dat april een afwijking heeft van 1,5°C , waardoor voor die maand een maandcorrectie een toegevoegde waarde zou zijn.

Waar de zomer als geheel niet gecorrigeerd hoeft te worden, zie je bij de individuele maanden wel grote verschillen (jun: 0,5°C , jul: 1,2°C , aug: 1,5°C)

In de winter zijn de verschillen tussen de verschillende maanden klein, maar in de herfst zie je dat september er met een afwijking van 1,9°C duidelijk uit springt.

Ook voor de Tx is een opsplitsing per seizoen en per maand gemaakt, zoals te zien is in afbeelding 5 en 6:

Afbeelding 5:

Afbeelding 6:

Met een gemiddelde afwijking voor de Tx van -0,9°C valt op dat de winter en herfst hierop geen afwijkingen laten zien.

In de lente (-1,4°C ) en zomer (-0,7°C) zou een correctie per seizoen dan weer wel een meerwaarde kunnen hebben.

Kijk je echter naar de afzonderlijke maanden binnen de seizoenen, dan zijn daar wel maanden te vinden waarbij een maandcorrectie een meerwaarde zou kunnen hebben.

Omdat ik me goed kan voorstellen dat voor sommigen al deze cijfers voor de ogen duizelen, heb ik dit jaar twee grafieken gemaakt, waarin al die cijfers zijn samengevat en worden weergegeven per maand.

Dit is te zien in afbeelding 7 en 8:

Afbeelding 7:

Afbeelding 8:

Al met al ontstaat er ieder jaar een beter beeld van waar en hoeveel GFS-uitvoer gecorrigeerd moet worden om tot een betrouwbare verwachting voor Tn en Tx te komen (op één cijfer achter de komma).

Het behoeft waarschijnlijk geen toelichting dat alles in beginsel valt of staat met wat GFS als uitvoer geeft.

Als de Opers voor een bepaalde dag er run na run helemaal naast zitten dan kan ook een rekenmodel zoals hierboven toegelicht, de meubels niet meer redden.

Op naar volgend jaar december, wanneer ik een toelichting kan geven op 7 jaar aan data.

[edit 30-12-2023] Door de hack van Weerwoord, waren de afbeeldingen verloren gegaan. Deze zijn weer toegevoegd.