Hoesung Lee, voorzitter IPCC   Bron: IPCC

Het vierde en laatste deel van het AR6 rapport van het IPCC,  het syntheserapport, is uit. Het vat de belangrijkste bevindingen van de drie voorgaande hoofddelen samen voor ‘policy makers’. De eerste drie delen hadden betrekking op de fysische wetenschap van het klimaat, hoe men zich daaraan kan aanpassen (adaptatie), en manieren om de uitstoot van broeikasgassen te verminderen. Die deelrapporten werden gepubliceerd in 2021 en 2022.

Het syntheserapport bevat nog sterkere waarschuwingen dat de wereld “onomkeerbare” niveaus (1,5 °C of 2  °C) van opwarming van de aarde nadert, waarbij catastrofale gevolgen snel onvermijdelijk worden. En dat het “nu of nooit” is om drastische maatregelen te nemen om een ramp te voorkomen. Kortom: all politics.

Nu is het al héél lang 5 voor 12 in de klimaatalarmistische gemeenschap, dat maakt dus steeds minder indruk. Het fysische deel van het rapport is interessant, omdat dat een tamelijk volledig beeld geeft van wat wetenschappers wel en niet weten over de werking van het klimaat. Dat is erg dik, erg wetenschappelijk en dus moeilijk, dus dat wordt door journalisten en politici zelden  ingekeken laat staan gelezen. Het syntheserapport is een politiek stuk, vooral door ambtenaren samengesteld.

Het vorige rapport, AR5, werd in 2014 afgerond. Het technische deel van AR6 is min of meer gelijk aan dat van AR5, maar de ‘vertaling’ naar politici wordt merkwaardig genoeg wel steeds alarmistischer.  Je zou verwachten dat in de periode die AR6 beslaat, 2015 tot 2023, er wellicht klimatologische ontwikkelingen hebben plaatsgevonden die het IPCC noopten om extra hard op de alarmbel te slaan. Opvallend is dan ook dat de gemiddelde temperatuur op aarde in die periode niet is gestegen, er was zelfs sprake van een lichte temperatuurdaling.

Fig.1    Bron: NOAA

De grafiek van figuur 1 is afkomstig van NOAA en toont de temperatuuranomalie van land en oceanen op aarde ten opzichte van de referentieperiode 1901-2000 voor de periode 2015 t/m  2022. Rechts bovenin staat de trend: -0,07  °C per decennium. Dat is voor de periode 2015-2022 -0,05 °C.

Fig.2    Bron: GISS NASA

Met behulp van de website van GISS NASA kun je een kaart maken van de temperatuurtrends op aarde in de periode 2015-2022. Interessant, want de trendverschillen per regio zijn groot. Op het kaartje van figuur 2 zijn de gebieden die in de periode 2015-2022 warmer zijn geworden oranje gekleurd en de gebieden die kouder zijn geworden blauw. Wit betekent dat de temperatuur ongeveer gelijk gebleven is. De data zijn afkomstig van de Aqua satellieten, die met behulp van de AIRS (Atmospheric Infrared Sounder) sinds eind 2002 de oppervlaktetemperatuur van de planeet meten. GISS NASA komt voor de periode 2015-2022 uit op een gemiddelde trend van -0,03 °C. Ook hier dus een lichte daling van de temperatuur.

Fig.3    Data: GISS NASA

Bovenstaande figuur 3 toont de temperatuurveranderingen op het kaartje van figuur 2  per breedtegraad van de Zuidpool tot de Noordpool voor de periode 2015-2022. De sterkste opwarming vind je nabij de 75 °ZB, op de grens van het continent Antarctica  en de omringende oceanen. Veranderingen in het oppervlak drijfijs heeft sterke invloed op de albedo en daarmee op de luchttemperatuur. Opvallend is dat bij 75 °NB er nauwelijks sprake is van opwarming. Dat matcht goed met mijn observaties dat er al meer dan 10 jaar geen afname van het Arctische drijfijs waargenomen wordt.

Opvallende opwarming vindt men ook op het NH op subtropische en gematigde breedte. Over de snelle opwarming in Europa en grote delen van Azië heb ik hier al vaker geschreven. Dat is vooral het gevolg van sterk toegenomen zoninstraling (SWD). De tropen koelden in de periode 2015-2022 echter opvallend sterk af.

De eerste en laatste 4 puntjes in de grafiek van figuur 3 lopen opvallend horizontaal. Dat komt omdat de satellieten vanwege hun baan niet in staat zijn om de temperaturen in de poolgebieden van 82 tot 90 graden NB en ZB te meten. Niet op die 8 puntjes letten dus, ze geven de temperatuur op de polen niet weer.

Is dit alles een voorbode voor koelere tijden? Sommige wetenschappers zijn daar van overtuigd, anderen niet. Zolang voor velen CO2 de enige ‘driver’ is van de temperatuur op aarde past een pas op de plaats als die vanaf 2015 niet in het plaatje. Maar ik ben er van overtuigd dat die bijna ziekelijke focus op CO2 sterk overdreven is en mensen blind maakt voor natuurlijke klimaatfluctuaties, die van alle tijden zijn. De voorlopige resultaten van ons onderzoek naar de stralingsbalans op het meetstation Cabauw is daar een mooi voorbeeld van: van 2001 t/m 2020 werd daar een afname van LWD (neerwaartse langgolvige straling = broeikaseffect) gemeten, en tegelijk een toename van de SWD (zoninstraling) en een sterke toename van de gemeten luchttemperatuur.  Daar is het laatste woord nog niet over gesproken.

Fig.1    Bron CBS

Het CBS meldde op 6 maart 2023 op haar website dat Nederland  in 2022 de productie van hernieuwbare elektriciteit met 20 procent heeft zien groeien. Het aandeel in de totale stroomproductie groeide van 33 naar 40 procent.  Dat had vooral te maken met de zon die vorig jaar overvloedig scheen. Dat ziet er wel spectaculair uit, maar was dat ook zo?

Op 22 december 202 schreef ik een uitgebreid artikel getiteld “Energie in Nederland nu en later”. Daarin heb ik gepoogd de tamelijk ingewikkelde CBS tabellen terug te brengen tot eenvoudige grafieken. Een van de valkuilen in het energieverhaal is dat er een groot verschil bestaat tussen het bruto- en het netto energieverbruik, respectievelijk 3023,7 PJ en 1803,5 PJ. Het verschil, 1230 PJ, is uiteindelijk aan de verbruikerszijde niet werkelijk beschikbaar als energie. Een van de oorzaken is dat ruim 530 PJ niet als energiebron maar als grondstof gebruikt wordt in bijvoorbeeld de chemische industrie. Dat geldt voornamelijk voor olie en aardgas. Een ander deel (700 PJ) gaat verloren bij de omzetting naar andere energievormen, bij distributie, of wordt gebruikt in het productieproces van de energiesector zelf. Steenkool wordt grotendeels omgezet in elektriciteit en warmte, aardgas bijna 40%. Daarom is het beter om naar het netto energieverbruik, te kijken. Dit is de grafiek  van het netto energieverbruik 2021 uit het artikel van december 2022:

Fig.2    Data: CBS

Hernieuwbare energie (hier inclusief biomassa) is hier bijna volledig opgegaan in de sectie ‘elektriciteit’ en ‘warmte’. De belangrijkste energiebronnen in 2021 zijn aardgas (bijvoorbeeld voor verwarming) en aardolieproducten zoals benzine, diesel en dergelijke. In werkelijkheid is het aandeel van met name aardgas in het netto energieverbruik nog veel groter omdat aardgas een belangrijke rol speelt in de elektriciteitsproductie (en levering van restwarmte). Datzelfde geldt ook voor steenkool.

Maar het artikel van het CBS (figuur 1) ging alleen over elektriciteit, niet over het resterede grootste deel van het energieverbruik. Bovendien gaat het CBS artikel niet over het elektriciteitsverbruik maar over de elektriciteitsproductie. Daarover later meer.

Fig.3    Data: CBS

Figuur 3 toont de elektriciteitsproductie in 2021 naar energiebron. Van de elektriciteitsproductie in 2021 ( = 22% van het netto energieverbruik in dat jaar) was 25% afkomstig van duurzame energiebronnen (wind, zon, waterkracht). Het CBS rekent bij die duurzame energiebronnen nog steeds biomassa, waarschijnlijk omdat biomassa (= bossen kappen) door de EU bestempeld is als ‘duurzaam’. Dat is de blauwe taartpunt in figuur 3 (8%) . Zo komt het CBS aan een duurzame stroomproductie in 2021 van 33%, terwijl het in werkelijkheid 25% was.

Dat CBS aandeel van 33% groeide volgens het CBS tot 40%, vooral door meer elektriciteit van zonnepanelen terwijl het aandeel fossiel afnam met 11%  tot 56%. Nu is ‘afname met 11%’ vervelend rekenen. Duidelijker is: het fossiel aandeel was in 2021 62% (zie figuur 3) en daalde in 2022 tot 56%. Een daling met 6% . Duurzaam steeg met 7%.

Fig.4    Data: CBS

Figuur 4 is ook afkomstig van het artikel over het energieverbruik 2021. Het aandeel van echt hernieuwbare (dus zonder hout verbranden) energiebronnen in de elektriciteitsproductie was zoals gezegd in 2021 25%, dat is dus 106 PJ. Windmolens en zonnepanelen leveren rechtstreeks elektriciteit aan het stroomnet. De productie van restwarmte is  gekoppeld aan elektriciteitsproductie met behulp van de niet hernieuwbare energiebronnen aardgas, steenkool, biomassa en kernenergie.

Kijken we naar het netto energieverbruik in 2021, 1803,5 PJ, dan was het aandeel van hernieuwbare energiebronnen dus 5,88% van het netto energieverbruik in 2021. Rekening houdend met de toename van 7% duurzaam in de elektriciteitsproductie in 2022, dan steeg het aandeel duurzaam in 2022 tot ongeveer 6,3%. Het netto energieverbruik van 2022 was bij het schrijven van dit artikel nog niet bekend.

Tot slot nog enkele woorden over de productie en het verbruik van elektriciteit. Tot voor een aantal jaren geleden ontliepen productie en verbruik elkaar niet veel. Er werd voor gezorgd dat de productie altijd wat hoger was dan het verwachte verbruik, dat goed bekend was. Op die manier was de leveringszekerheid in Nederland zeer hoog. Bijspringen en afschakelen van gascentrales was de regulator van het netwerk. De productie van elektriciteit was destijds vraaggericht  en dus goed controleerbaar.

Met de komst van met name wind- en zonne-energie werd een steeds groter deel van de elektriciteitsproductie aanbodgericht: de weersomstandigheden en het dag/nachtritme bepalen steeds meer de elektriciteitsproductie. Dat dat tot onstabiele prijsvorming leidt zag ik afgelopen 24 uur op mijn energie-app:

 
Fig. 5      Bron: ANWB app

Het uurtarief van elektriciteit daalde gisterenmiddag 4 uren lang tot € 0,00 per kWh en vannacht nogmaals 4 uren lang tot € 0,00. Fijn voor mij, ik heb (bijna) gratis een paar wasjes gedaan. Iedereen voelt wel aan dat het niet veel doller moet worden, maar dat wordt het wel vrees ik. Een toenemende onevenwichtigheid in de stroomvoorziening is te voorzien.

In het artikel van 23 februari j.l. besteedde ik voor het eerst aandacht aan de infraroodmetingen die  verricht worden op station Cabauw. Verrassend was dat in de periode 2001 t/m 2020 niet alleen de LWD (langgolvige straling van atmosfeer naar aardoppervlak) een lichte negatieve trend vertoont maar ook de LWU, de langgolvige straling van het aardoppervlak naar de atmosfeer. Overigens waren die trends statistisch niet significant.

Fig.1    Data: KNMI

Je zou verwachten dat de LWD (= broeikaseffect van broeikasgassen plus wolken) vanwege het tussen 2001 en 2020 gestegen atmosferisch CO2-gehalte gestegen zou zijn, maar blijkbaar zorgde de afname van de bewolking in die periode voor een lichte netto afname van het broeikaseffect.

De lichte daling van de LWU was verrassend, met name omdat de luchttemperatuur die op 1,5m hoogte wordt gemeten in Cabauw een forse stijging laat zien van 2001 t/m 2020:

Fig.2    Data: KNMI

De luchttemperatuur op 1,5 m hoogte is dus niet synoniem aan de LWU die uitgestraald wordt door het aardoppervlak. Die luchttemperatuur is een maat voor de gemiddelde bewegingsenergie van de atomen of moleculen op 1,5 m hoogte. Er is energie nodig om de temperatuur te laten stijgen, en om de temperatuur te laten dalen moet energie weggehaald worden. Die energie om de luchttemperatuur op 1,5m hoogte te veranderen komt in de vorm van straling (LWD, LWU, SW netto), maar ook van verticale beweging van luchtn (LH= latente warmte en SH=voelbare warmte) en horizontale beweging van lucht (advectie). De luchttemperatuur op 1,5 m hoogte is dus de uitkomst van de aan- en afvoer van diverse vormen van energie, niet alleen van LWU.

In het voorafgaande artikel liet ik al even zien hoe een grafiek op basis van de 10-minutenmetingen in Cabauw er uit  ziet. Ik was nieuwsgierig naar de verschillen in de stralingsbalans in Cabauw in de zomer en winter, en bij diverse graden van bewolking. Ik gebruikte de parameter SP (percentage van de langst mogelijke zonneschijnduur per etmaal) als proxy voor de bewolkingsgraad en filterde de etmalen met maximale relatieve zonneschijnduur (SP>90) en maximale hemelbedekking (SP=0). Dat leverde over de gehele periode van 2001 t/m 2020 94 etmalen op met SP=0 of SP>90. Vervolgens heb ik gekeken of er blokjes van 5 etmalen voorkwamen waarbij tenminste 1 etmaal PS=0 had en een ander etmaal SP>90. Dat leverde in totaal slechts 4 blokjes van 5 etmalen op over de gehele periode.

Fig.3    Data: KNMI

Figuur 3 toont de stralingsgegevens (10-minutendata) over 5 dagen in juli 2008.  Net SW is SWD minus SWU. Net SW is dus kortgolvige straling die voor 100% geabsorbeerd wordt door het aardoppervlak. Overdag is de inkomende zonnestraling evenredig met de sinus van de hoek van instraling van de zon. Te zien is dat op 4 en 5 juli de net SW nagenoeg perfect is, geen bewolking dus. Vanaf 6 juli komt er steeds meer bewolking binnendrijven, die er voor zorgt dat de grafiek van net SW niet meer perfect is. Op 8 juli is de hemel vrijwel volledig bedekt met wolken en komt de net SW nauwelijks meer boven de 100 W/m2 uit, terwijl de net SW op 4 juli nog tot 750 W/m2 reikte.

Interessant is wat dat voor consequenties heeft voor de LWU en LWD. Op 8 juli met volledig bedekte hemel liggen LWD (warmtestraling omlaag) en LWU (warmtestraling opwaarts) ongeveer op hetzelfde stralingsniveau van 400 W/m2. LWU iets hoger dan LWD. Hoe komt dat?

Laten we de Wet van Stefan Boltzmann eens toepassen op de situatie op 8 juli 2001. LWD is gemiddeld over het hele etmaal 395,3 W/m2, LWU 404,4 W/m2 Dat betekent dat de onderste straler (aardoppervlak) een gemiddelde temperatuur had van  290,61 K en de bovenste straler (wolkenbasis) een temperatuur van  288,96 K. in graden Celsius is dat respectievelijk 17,5 °C en  15,8 °C. Een klein verschil dat zo kenmerkend is voor een volledig bedekte hemel. Als we aannemen dat de atmosfeergradiënt -6,5 °C/km was, dan lag de wolkenbasis toen op ongeveer 250m hoogte. Ongeveer, want de werkelijke atmosfeergradiënt kennen we niet, en aardoppervlak en wolkenbasis zijn geen zwarte stralers.

Kijken we naar een situatie met clear sky, zoals op 4 juli 2001, dan zien we dat de LWU vrijwel onmiddellijk reageert op de toenemende net SW vanaf zonsopgang. De korte vertraging van ruim 1 uur tussen de top van net SW en LWU is het gevolg van de (geringe) warmtecapaciteit van de bodem. De LWD curve hangt af van de luchttemperatuur die zijn maximum bereikt in de late namiddag vóór zonsondergang en zijn minimum bereikt net na zonsopgang.

De clear sky situatie zorgt er voor dat tussen zonsopgang en zonsondergang  het verschil tussen LWU en LWD gemiddeld 100 W/m2 is, en zijn maximum van 120 W/m2 bereikt rond 12 uur ’s middags. Om 12 uur is de LWU 468,4 W/m2 en de LWD 347,5 W/m2.  Omrekenen naar °C betekent dat om 12 uur een grondtemperatuur van 28,3 °C aanwezig was en een luchttemperatuur van 6,64 °C. Dat is een verschil van 21,7 °C. Omdat een wolkenbasis ontbreekt is die luchttemperatuur hier de gemiddelde stralingshoogte van waaruit de LWD straalt. Bij een atmosfeergradiënt van -6,5 °C/km komt de LWD  gemiddeld van een hoogte van  3,3 km.

Fig.4    Data: KNMI

De grafiek van figuur 4 toont nog een zomersituatie, maar nu met 3 dagen clear sky en 2 dagen met gedeeltelijke bewolking: 1 juli lichte bewolking en 2 juli wat meer wolken zoals de net SW laat zien. Aan de snelle verhoging van de LWD rond 3 uur ‘s nachts op 2 juli is te zien dat dan meer wolken boven Cabauw verschijnen. De ruimte tussen de blauwe en bruine lijn neemt dan beduidend af.

Fig.5    Data: KNMI

Bovenstaande grafiek laat de derde zomerse situatie zien. Nu geen clear sky meer, maar 5 dagen met wat meer en wat minder bewolking. Ik denk dat de grafiek verder voor zich spreekt.

Fig.6    Data: KNMI

Fig.7    Data: KNMI

De figuren 6 en 7 geven twee wintersituaties weer.  De net SW (zon) grafieken blijven zelfs op de clear sky dagen  laag en reiken niet veel hoger dan zo’n 300 W/m2. De overgang van clear sky naar volledig bewolkt en andersom levert een prachtig niveauverschil op tussen LWD en LWU. Je kunt heel nauwkeurig aflezen hoe laat de wolken binnendrijven en wanneer ze weer verdwijnen.

Die 10-minutendata van de straling is een waar goudmijntje. Misschien een volgende keer nog wat meer gegevens uit die mooie databank van Cabauw.

Fig.1    Bron: Wikipedia

Bij de communicatiemast van Cabauw, 18 km ten ZW van Utrecht, worden naast de reguliere meteo metingen ook de langgolvige (IR) straling naar beneden (LWD) en omhoog (LWU) gemeten.  Dat gebeurt met behulp van zogenaamde pyrgeometers die op 1,5m hoogte boven het maaiveld gemonteerd zijn.

Fig.2    Bron: EPLAB

Vanaf mei 2000 wordt 6x per uur zogenaamde 10-minutenmetingen uitgevoerd, een etmaal telt dus 144 metingen. Data van gehele jaren zijn beschikbaar van 2001 t/m 2020. De data zijn beschikbaar gemaakt door het KNMI maar het bestandstype was voor mij niet leesbaar. Een vriend heeft ze om kunnen zetten zodat ik ze kon lezen en geschikt maken voor Excel. Het eerste dat ik daarna gedaan heb is de 10 minuten data omzetten naar twee extra databanken, etmaalgemiddelden en maandgemiddelden.

Fig 3    Data: KNMI

Figuur 3 toont een etmaal LWD en LWU en net SW (= SWD minus SWU) op basis van de 10 minutenmetingen op 25 mei 2001. Net SWD is wat er aan kortgolvige zonnestraling binnenvalt (global irradiation) minus het deel dat aan het aardoppervlak reflecteert.  Er is sprake van een clear sky situatie, alleen rond 8u ’s ochtends schuift er twee maal een klein wolkje voor de zon. ’s Nachts tussen 1 en 3 uur komt er blijkbaar ook bewolking voorbij, die zorgt voor een verhoging van de LWD. Wolken dragen immers bij aan het broeikaseffect. Als je goed kijkt zie je dat dat op zijn beurt weer zorgt voor een geringe stijging van de LWU. Later meer over de 10-minutenmetingen.

Fig.4    Bron: KNMI

De energiebalans van figuur 4 zijn we al enkele malen hier tegengekomen, onder andere in het recente bericht over de kortgolvige straling op station De Bilt. Bedenk dat de waarden in figuur 4 gemiddelden zijn van de gehele aarde. Vullen we de gegevens van Cabauw van figuur 3 in de stralingsbalans in, dan is de LWD ongeveer 300W/m2, LWU ongeveer 400 W/m2, en net SW (dat deel van het zonlicht dat daadwerkelijk het aardoppervlak in Cabauw verwarmt) van 0 W/m2 ’s nachts naar ruim 700 W/m2 rond 12u, gemiddeld over het etmaal 265 W/m2. Voor de laatste waarde heb ik even in de tabel van de dagdata gekeken.

Het fijne is dat nu behalve de kortgolvige stralingsdata (zonlicht SWD en SWU) ook de langgolvige warmtestralingsdata (IR LWD en LWU) voorhanden zijn. Daarmee hebben we de belangrijkste stralingscomponenten van de energiebalans van Cabauw ter beschikking. Dat kan een schat aan informatie opleveren. Vandaag kijken we naar de maanddata van kort- en langgolvige straling.

Fig.5    Data: KNMI

De blauwe lijn is de LWD, de langgolvige straling die vanuit de atmosfeer naar de aarde gericht is. Dit is het broeikaseffect, veroorzaakt door enkele broeikasgassen zoals waterdamp (de belangrijkste) en CO2. Verder spelen wolken ook een belangrijke rol in de LWD. Iedereen weet uit ervaring da een bewolkte nacht afkoeling tegengaat. De LWD volgt de seizoenen, wat logisch is: ’s zomers warmt de grond extra op, waardoor er meer LWU is die op zijn beurt weer zorgt voor een grotere LWD.

De zwarte lijn is de LWU, de infrarood straling ie het gevolg is van de opwarming van het aardoppervlak. De LWU is vrijwel altijd groter dan de LWD: het aardoppervlak wordt verwarmd door de LWD + net SW overdag + aanvoer van energie van elders door de wind.

De streepjeslijnen zijn de lineaire trendlijnen van LWD en LWU, en hun verloop wekt verbazing want er is geen toename van de trend. Sterker, de trends zijn beide licht negatief, maar dat is zo weinig dat dat de trend statistisch niet significant is. Dat er geen trend omhoog is is verbazingwekkend omdat het atmosferisch CO2 gehalte tussen 2001 en 2020 gestegen is van 371 ppm naar 414 ppm, wat goed zou moeten zijn voor een stralingsforcering van 0,5 W/m2. Mogelijke oorzaak kan een afname van de bewolking zijn, en/of veranderingen van het waterdampgehalte van de atmosfeer. In elk geval laat het wel zien dat in Cabauw CO2  van 2001 tot 2020 zeker niet een dominante factor is.

Fig.6    Data: KNMI

Figuur 6 toont de naar beneden gerichte global irradiation (SWD) en de naar boven gerichte SWU, het gedeelte van de SWD dat aan het aardoppervlak reflecteert. Opvallend is dat de SWD ’s zomers sterke verschillen kent van jaar tot jaar, vooral als gevolg van zomerse verschillen in bewolking. De onderzijde van de rode lijn ligt nagenoeg op de lineaire trendlijn: ’s winters is de instraling klein en is de bewolking in dat jaargetijde meestal prominent aanwezig. De lineaire trendlijn van de SWD toont een behoorlijk sterke stijgende trend van 8,6 W/m2, waarover in het vorige artikel van 4 februari 2023 al voldoende geschreven is. De trend van de SWU is ongeveer half zo groot als die van SWD.

Fig.7    Data: KNMI

De grafiek van figuur 7 is afgeleid van figuur 6. De oranje lijn is de net SW, dus SWD minus SWU. De paarse lijn is de albedo, de reflectiewaarde van het aardoppervlak in %. De opwaartse trend is 0,0179x, dat is 4,3 W/m2 van 201 t/m 2020. Dat is de helft van de SWD, wat hopelijk geen verbazing zal wekken. De albedo laat een viertal grote uitschieters en enkele kleine zien. Alle in de winter, dus ik ga er van uit dat dat het gevolg is van de aanwezigheid van een sneeuwdek. De albedo neemt licht toe en is gemiddeld zo’n 24%. Dat is iets hoger dan normaal voor grasland wordt aangenomen (18 a 20%).

Later wat meer details over de langgolvige stralingsdata van Cabauw.

Fig.1    Screenshot NOS journaal

Media als tv, radio en kranten waarschuwen regelmatig voor fake nieuws. Het is opmerkelijk dat dan met name misinformatie bedoeld wordt op het internet, vooral bij social media. Zelden of nooit wordt de hand in eigen boezem gestoken. Dat verbaast me, want traditionele media zijn niet vrij van smetten. Daarom ben ik zeker ook op mijn hoede bij traditionele media.

Op het gebied van klimaat is de berichtgeving vaak slecht en alarmistisch. Zo is de Volkskrant berucht vanwege de bijna activistische wijze waarop er over klimaat-items wordt bericht, maar radio en tv kunnen er ook wat van. Ik vermijd dergelijke misinformatie zo veel als mogelijk, maar soms ontkom ik er niet aan en zie bij het zappen even een flard langskomen. Dat gebeurde kort geleden ook weer, op vrijdag 17 februari j.l. bij het NOS journaal van 18 uur (vanaf 10:00) .

Het was een kort item over Antarctica. Ik heb een paar screenshots genomen en het commentaar erbij uitgeschreven, dat uit 3 zinnen bestaat: “ Sinds de eerste metingen begin jaren ’70 heeft er nog nooit zo weinig ijs op Antarctica gelegen. Het is daar nu zomer en er ligt relatief weinig ijs, maar dit is wel heel weinig. Er ontbreekt nu een gebied met zeeijs ter grootte van 21 x Nederland.”

Fig.2    Bron: NASA

Een argeloze kijker wordt hier volstrekt op het verkeerde been gezet. Antarctica is een continent, bedekt met landijs dat gemiddeld zo’n 2200 m dik is. Op de sommige plaatsen bereikt het landijspakket een dikte van bijna 4800 m. Met die hoeveelheid landijs gebeurt vrij weinig. De bewering van het NOS journaal dat er al vanaf begin jaren ’70 aan dat landijs gemeten wordt lijkt me onjuist. Pas vanaf 1992 wordt er met behulp van satellieten aan het landijs gemeten, in eerste instantie met behulp van altimetrietechnieken, vanaf 2002 met gravimetrietechnieken (satelliet Grace).

De grafiek van figuur 2 ziet er dramatisch uit, maar is het niet: met een totale ijsmassa van 27.000.000 Gt zou met het tempo van de afname sinds 2002 (169 Gt/jaar) meer dan 150.000 jaar duren voordat het ijs op Antarctica verdwenen is. Gezien de zeer korte periode dat er serieus gemeten wordt aan het Antarctische landijs is het volstrekt onzeker of de lichte afname die we zien sinds 2002 ook werkelijk een trend is die doorzet. Er lijkt in 2020 een trendbreuk waarneembaar, maar niets is zeker in dergelijke korte meetreeksen.

Fig.3    Bron: King et al 2012

Figuur 3 is afkomstig van een artikel van King et al in Nature, en toont per afstroombekken op Antarctica de veranderingen in de massabalans op basis van de Grace data tussen 2002 en 2010. Die massabalans is de uitkomst van wat er aan ijs bijkomt (door sneeuw) minus wat er door afstroming van ijs afgaat.  Een klein beetje verdwijnt door sublimatie en smelt. King et al: “We resolve 26 independent drainage basins and find that Antarctic mass loss, and its acceleration, is concentrated in basins along the Amundsen Sea coast. Outside this region, we find that West Antarctica is nearly in balance and that East Antarctica is gaining substantial mass.”

Fig.4    Bron: Singh et al 2020

Recent onderzoek van Singh et al in Nature wijst uit dat het Antarctische continent de laatste zeven decennia niet is opgewarmd, ondanks een constante toename van de atmosferische concentratie van broeikasgassen. Figuur 4 toont de gemiddelde oppervlaktetemperatuur anomalie op Antarctica over de periode 1984-2014 vergeleken met de periode 1950-1980.
De hoogtelijnen tonen de hoogte boven zeeniveau (in m).

De figuur laat zien dat slechts een klein deel van het continent, met name het schiereiland West Antarctica, opwarming vertoont. Het grootste deel van Antarctica laat afkoeling zien. De grote hoogteligging van de Antarctische ice sheet zou daar een belangrijke rol in spelen.

De eerste twee zinnen van het NOS-commentaar ( “ Sinds de eerste metingen begin jaren ’70 heeft er nog nooit zo weinig ijs op Antarctica gelegen. Het is daar nu zomer en er ligt relatief weinig ijs, maar dit is wel heel weinig.”) zijn dus volstrekte onzin. Maar in de laatste zin ( “ Er ontbreekt nu een gebied met zeeijs ter grootte van 21 x Nederland.” )  komt de aap uit de mouw: het gaat helemaal niet over het ijs op Antarctica maar over het zeeijs (drijfijs) rondom Antarctica. Dat zeeijs neemt ’s winters toe en ’s zomers af:

Fig. 4    Bron: NRLSSC

Figuur 4 is een animatie van de veranderingen van de oceaantemperatuur en het zeeijs oppervlak gedurende een jaar van eind 2017 tot eind 2018. Het continent Antarctica wordt aangegeven met de beige kleur en het zeeijs rondom het continent met een lichtgrijze kleur. Dat zeeijs kent een minimum oppervlakte rond eind februari.

Fig.5    Screenshot NOS journaal 17-2-2023

Het screenshot van figuur 5 toont het kaartje dat het NOS journaal aan het einde van het bericht liet zien. Die vreemde blauwe vlekken rond Antarctica zijn bij nadere inspectie blauwe kaartjes van Nederland , waarschijnlijk 21 😉 ( “ Er ontbreekt nu een gebied met zeeijs ter grootte van 21 x Nederland.” ).  Ik denk dat het de  bedoeling is dat die kaartjes op die plekken liggen waar nu geen zeeijs is maar ‘normaal gesproken’ wel. Maar zeker weet ik dat niet, ik heb nog nooit zo’n raar kaartje gezien. De vergelijking met 21x het oppervlak van Nederland moet waarschijnlijk de dramatiek verhogen. En wat betekent ‘normaal gesproken’?

Fig.6    Bron: NSIDC

Is het oppervlak zeeijs rond Antarctica werkelijk dramatisch laag? Figuur 6 is de meest recente grafiek en is bijgewerkt t/m 18 februari 2023. De grafiek laat zien dat er van een dramatische situatie geen sprake is. Het zeeijs-oppervlak was volgens de grafiek op 18 februari 1,846 miljoen km2 (blauwe lijntje linksonder). Het oppervlak van Nederland is 41.543 km2. Het zeeijs rond Antarctica had dus een oppervlak van ruim 44x het oppervlak van Nederland. In 2022 was dat niet veel meer.

Conclusie: het NOS journaal heeft een bericht de wereld ingestuurd dat onjuiste informatie geeft en suggestief is. Men haalt landijs en zeeijs door elkaar, zodat het lijkt alsof het ijs op Antarctica dramatisch vermindert. En als het zeeijs-oppervlak een klein beetje lager is dan afgelopen jaar spreekt men van een oppervlak ter grootte van 21x Nederland die zou ontbreken. Hoezo ‘ontbreken’, overigens?

Onlangs verscheen op deze website een stuk over de Labrijnreeks. Daarin was te zien dat vanaf het begin van de reeks in 1706 de jaartemperatuur t/m 1988 geen trend vertoonde. Figuur 1 is het laatste stuk van die Labrijnreeks en beslaat de periode vanaf 1901 t/m 2022.

Fig.1  Data: KNMI

Het eerste stuk van 1901 t/m 1988 gaat licht opwaarts maar die trend is statistisch niet significant. Vanaf 1989 t/m 2022 is de trend echter sterk (2,9 °C/eeuw) en levert van 1989 t/m 2022 in De Bilt een temperatuurstijging van 1 °C . Inclusief de temperatuursprong 1988-1989 leverde dat een temperatuurtoename op van 1,6 °C.

De temperaturen in figuur 1 zijn gebaseerd op de dagelijkse Tg van station De Bilt. De gemeten temperaturen van 1901 tot 1 september 1951 zijn door het KNMI in 2016 gecorrigeerd (gehomogeniseerd). Door Dijkstra et al (2022) is aangetoond dat er daarbij sprake was van overcorrectie. Het KNMI heeft als reactie hierop aangegeven dat die homogenisatie in 2023 opnieuw bekeken gaat worden. Ik gebruik hier derhalve de gemeten (niet gehomogeniseerde) temperaturen.

Als het over opwarming of klimaatverandering gaat dan denken veel mensen automatisch dat die sterke toename wel het gevolg zal zijn van de uitstoot van CO2 door de mens. Gezien de eenzijdige berichtgeving over dat onderwerp is dat overigens niet verbazingwekkend. Nu is het zeker juist dat CO2 een broeikasgas is, het is na waterdamp zelfs het sterkste broeikasgas en het staat vast dat  broeikasgassen de temperatuur aan het aardoppervlak verhogen. De cruciale vraag is echter hoe groot dat CO2-effect is en welke andere factoren mogelijk een rol spelen bij die opmerkelijke opwarming in De Bilt. Anders dan radio, tv en kranten vaak melden is daar in de wetenschappelijke wereld nog veel onzekerheid over.

Ik ga proberen meer zicht te krijgen op die opvallende stijging van de temperatuur in De Bilt vanaf 1989 en kijk eerst naar de straling.

Straling

Fig.2  Naar: KNMI

De zon levert voor bijna alle processenaan het aardoppervlak de benodigde energie, direct of indirect. De energiestromen die daarmee gemoeid zijn worden vaak zichtbaar gemaakt in de zogenaamde stralingsbalans (figuur 2). ‘Balans’ slaat op het evenwicht tussen wat het aarde-atmosfeersysteem binnenkomt en wat het systeem verlaat. Binnen komt kortgolvige straling van de zon (341 W/m2), er uit gaan een deel van het zonlicht dat gereflecteerd wordt plus langgolvige (warmte-)straling. Die warmtestraling ontstaat door absorptie van zonlicht door het aardoppervlak en deeltjes in de atmosfeer.

De straling van de zon (kortgolvige straling) is gemiddeld 341 W/m2 aan de top van de atmosfeer (TOA) en wordt op zijn neerwaartse weg door de dampkring gedeeltelijk teruggekaatst naar de ruimte door wolken en aërosolen in de atmosfeer (79 W/m2) en een deel ervan wordt in de atmosfeer geabsorbeerd (78 W/m2). Wat overblijft (184 W/m2) valt in op het aardoppervlak en wordt globale straling genoemd. Een deel (23 W/m2) wordt door het aardoppervlak gereflecteerd en verlaat ongebruikt de dampkring, samen met de door wolken en aerosolen weerkaatste kortgolvige straling.

Fig.3    Bron: ECMWF

De globale straling minus de reflectie aan het aardoppervlak is het deel van de kortgolvige straling die het aardoppervlak verwarmt. Dat wordt de netto kortgolvige straling genoemd en is in de stralingsbalans 161 W/m2 groot. Dat is de waarde van de gemiddelde netto kortgolvige straling op aarde. De netto kortgolvige straling op een bepaalde plaats is onder andere afhankelijk van breedteligging, het seizoen en de albedo (reflectiewaarde) van het aardoppervlak ter plaatse. Figuur 3 laat de ruimtelijke verschillen in waarde zien van de netto kortgolvige straling. Het zal niet verwonderlijk zijn dat de waarden gemiddeld het hoogst zijn in de tropen en het laagst op de polen. De smalle lichtere strook boven de oceanen nabij de evenaar toont de ITCZ (intertropische convergentiezone) waar door botsende passaten een sterke opwaartse luchtstroom ontstaat met enorme onweerswolken die instraling belemmeren.

Aan de rechterzijde van de stralingsbalans (fig. 2) is de langgolvige straling zichtbaar die het gevolg is van de opwarming van het aardoppervlak en de atmosfeer door de kortgolvige straling van de zon. Door wolken en broeikasgassen wordt warmtestraling (IR) die vanaf het aardoppervlak omhoog gaat grotendeels geabsorbeerd en direct naar alle kanten uitgestraald, deels gericht op het aardoppervlak. Daardoor neemt de temperatuur aan het aardoppervlak toe, het zogenaamde broeikaseffect. Zonder dat broeikaseffect zou de gemiddelde temperatuur op aarde -18 °C zijn en dus onleefbaar koud. Een klein deel van de langgolvige straling vanaf het aardoppervlak, 40 W/m2, vindt overigens ongehinderd zijn weg naar de ruimte door het zogenaamde atmosferisch venster. Wie nieuwsgierig is hoe de langgolvige straling buiten het atmosferisch venster toch zijn weg naar de ruimte weet te vinden leest boven in het menu onder “Feiten” de pagina “Broeikastheorie”.

Fig.4    Bron: ECMWF

In de stralingsbalans van figuur 2 is te zien dat de neerwaartse langgolvige straling gemiddeld 333 W/m2 is en de opwaartse langgolvige straling 396 W/m2. Het verschil, de netto langgolvige straling vanaf het aardoppervlak is dus 63 W/m2. Samen met de energiestromen van voelbare en latente warmte verlaat dus 160 W/m2 energie in de vorm van langgolvige straling (warmtestraling) het aardoppervlak, (bijna) gelijk aan wat er aan netto kortgolvige straling (zonlicht) binnenkwam. Figuur 4 toont de netto langgolvige straling aan het aardoppervlak. Het  – teken betekent dat er vanaf het aardoppervlak gezien meer uitgaat dan inkomt. Woestijnen hebben de hoogste langgolvige straling ‘up’ door het nagenoeg ontbreken van wolken en de hele lage waterdampconcentraties in woestijnlucht.

De Wet van Stefan Boltzmann geeft aan hoe de stralingskracht van een voorwerp afhangt van de temperatuur en de grootte van het oppervlak van een hemellichaam:

P =σAT⁴

Waarbij P = vermogen in Watt, σ = constante 5,670373·10^-8  , A = oppervlak in m² en T = temperatuur in K.

Vullen we deze formule in met P van de uitgaande langgolvige straling ( 239 W/m2), bekend door satellietmetingen), dan is de temperatuur van de aarde 255 K, dus -18 °C. Dat is de stralingstemperatuur van het systeem aarde, dus inclusief de atmosfeer. We schatten de gemiddelde temperatuur aan het aardoppervlak op 15 °C, dus de temperatuur van de dampkring op grotere hoogte moet dan veel lager zijn dan het aardoppervlak om gemiddeld op -18 °C uit te komen. Dat klopt ook: op 10 km hoogte zijn temperaturen van ongeveer -50 °C heel gewoon. Dat verschil tussen de berekende temperatuur van -18 °C en de gemeten temperatuur van +15 °C aan het aardoppervlak is dus het gevolg van het broeikaseffect. In werkelijkheid is het verschil nog groter, maar dat wordt deels gedempt door de stromen van voelbare (convectie, 17 W/m2) en latente warmte (verdamping, 80 W/m2) vanaf het oppervlak. Leven op aarde is mogelijk dankzij het broeikaseffect.

Fig.5    Bron: ECMWF

Figuur 5 toont de ruimtelijke verdeling van de temperatuur aan het aardoppervlak, gemiddeld over een jaar. De temperatuur is uiteraard afhankelijk van de waarden van de kort- en langgolvige straling die hiervoor besproken zijn. Maar de temperatuurverdeling op aarde wordt ook sterk beïnvloed door transport van warmte door zeestromen en winden. Zonder dat transport zou het in de tropen aanzienlijk warmer zijn en op de polen aanzienlijk kouder. De figuur laat ook het effect zien van de Golfstroom/Noord-Atlantische Drift, de warme zeestroom die vanaf de Caraïbische Zee  langs de West-Europese kust stroomt en warmte vanuit de tropen tot in de Arctische Zee brengt. Die zeestroom zorgt voor een extra opwarming van Europa van ongeveer 10 °C, waardoor voedsel verbouwen hier mogelijk is.

De Bilt

Fig.6    Data: KNMI

Terug naar De Bilt, waar ik bekijk welke invloed de zon vanaf 1950 heeft gehad op de temperatuurontwikkeling aldaar. Alle energiestromen van de stralingsbalans uit figuur 2 zijn natuurlijk ook in De Bilt werkzaam, alleen de waarden verschillen. Omdat de cijfers van de globale straling in De Bilt pas vanaf 1958 beschikbaar zijn analyseer ik de periode van 1950 t/m 2022. De lineaire trend van 1950 t/m 1988 is klein en statistisch niet significant. Dat geldt niet voor de periode 1989-2022. De temperatuursprong van 1988-1989 was 0,66 °C en de temperatuurstijging van 1989 t/m 2022 was 1 °C.

Fig.7    Data: KNMI

Zoals ik al enkele malen eerder op deze website heb laten zien  heeft de invallende kortgolvige straling zeker veel -zo niet alles- te maken met de temperatuursprong en de stijging na de sprong. Globale straling is het totaal van direct en diffuus zonlicht dat in De Bilt binnenvalt. Diffuus licht ontstaat door verstrooiing van zonlicht door gassen, deeltjes en wolken. Ruwweg is zo’n 60% van het invallende licht in De Bilt diffuus en de rest direct.

De grafiek van figuur 7 toont het verloop van de globale straling Q in De Bilt van 1958 t/m 2022. Vergelijking tussen figuur 6 (temperatuur) en figuur 7 (globale straling) laat zien dat de ligging van de breuk in beide grafieken identiek is. De sterke opwaartse trend van 1989 t/m 2022 levert een toename van de globale straling op van maar liefst 16 W/m2.

Fig.8    Data: KNMI

Figuur 8 toont het verloop van het aantal zonuren per jaar vanaf 1950. Ook hier de opvallende breuk eind jaren ’80. De sterke trend in het tweede deel van de reeks laat zien dat vanaf 1989 het aantal zonuren is toegenomen van 1540 uren tot 1893 zonuren eind 2022. Een formidabele toename met 353 zonuren, dat is maar liefst 23% meer zonuren in de periode na de breuk.

Fig.9    Data: KNMI

Bovenstaande grafiek toont zowel Q (globale straling) als SQ (zonuren) op station De Bilt. R2=0,9, dus de correlatie tussen SQ en Q is zeer sterk. Dat is niet verwonderlijk, directe zonnestraling levert (veel) meer instraling op waardoor de globale straling Q automatisch groter wordt. De globale straling ‘volgt’  dus de zonuren.

Die zonuren worden berekend uit de data van de globale straling. De WMO (World Meteorological Organization) definieert zonneschijnduur als de periode waarin de directe zonne-instraling een drempelwaarde van 120 W/m2 overschrijdt.

Q (globale straling) is een belangrijke parameter als het gaat om de temperatuur op een bepaalde plaats. Het aantal zonuren per jaar is een mooie proxy voor de bewolking. Meer bewolking betekent immers minder direct zonlicht. Het is daarom niet moeilijk om in te zien dat de spectaculaire toename sinds 1989 van het aantal zonuren met maar liefst 33% het gevolg is van een sterke afname van de bewolking.

Fig10    Bron: KNMI

De lucht boven Nederland wordt met name in het zomerhalfjaar steeds zonniger en droger, schrijft het KNMI. Bij lagere relatieve vochtigheid raakt de lucht minder snel verzadigd met waterdamp, ontstaan er minder snel wolken en is er meer direct zonlicht. De grafiek van figuur 10 laat zien dat er gemiddeld 42 zonuren méér zijn als de relatieve vochtigheid 1 procent daalt.

Fig.11    Data: KNMI

Figuur 11 laat zien dat het verloop van de relatieve luchtvochtigheid in De Bilt aardig in de pas loopt met de toename van het aantal zonuren. De relatieve vochtigheid kent gemiddeld over de hele periode een aflopende trend, maar vooral vanaf eind jaren ’90 is de daling spectaculair.

Veranderende luchtcirculatie

De sterke stijging van de temperatuur in De Bilt vanaf 1989 is dus vooral het gevolg van toegenomen globale straling. En die toegenomen globale straling kan verklaard worden door afnemende bewolking. Hierboven zagen we dat de afgenomen relatieve vochtigheid daarbij een rol kan spelen. Maar wat is de oorzaak achter die afname van de bewolking en relatieve vochtigheid? Het KNMI zoekt het in en een tweetal zaken. Deels komt dit doordat er deze eeuw  in het zomerhalfjaar vaker hogedrukgebieden waren, met weinig wolken, veel zon en drogere lucht. De andere factor is de opwarming van het land, waardoor de bodem ‘s zomers droog wordt en door de afwezigheid van verdamping de temperatuur extra stijgt. Vochtige oceaanlucht warmt dan boven het land op en de relatieve vochtigheid daalt (KNMI).

Maar die extra opwarming ’s zomers boven het land door het ontstaan van een droge bodem lijkt  vooral het gevolg van meer hogedrukgebieden. Immers, in hogedrukgebieden valt minder neerslag en is de verdamping met name door de sterke toename van direct zonlicht groter. Natuurlijk spelen bij verdamping ook de luchttemperatuur en de windsnelheid ook een rol, maar direct zonlicht heeft een sterke invloed op de verdamping. De verdamping is vooral afhankelijk van de beschikbare hoeveelheid zonlicht en is daarmee in de zomermaanden belangrijk hoger dan gedurende de winter. De achterliggende oorzaak lijkt dus een verandering te zijn in de luchtcirculatie boven ons  deel van Europa.

Fig. 12    Bron: Pik Report No 119

Eind 19^e eeuw begon de Duitse weerdienst met het registreren van luchtcirculatiepatronen boven Europa. Bovenstaande tabel is afkomstig van een rapport van het Potsdam Institute Climate Impact Research (PIK) dat de classificatie van luchtcirculatietypen beschrijft. Die circulatiepatronen vormt de basis van de classificatie  die de Duitse meteorologische Dienst (DWD) vanaf het eind van de 19^e eeuw hanteert. De luchtcirculatie boven NW Europa wordt dan ingedeeld op basis van 29 zogenaamde Grosswetterlagen (GWL). Grosswetterlagen zijn een bron van nuttige informatie over luchtcirculatieveranderingen in de loop van de tijd. Ik heb daar in het verleden al vaker aandacht aan besteed, zie onder andere hier.

DWD vat de 29 GWL samen in een achttal Grosswettertypen en drie Zirkulationsformen. Die laatste zijn zonaal (west-oost gericht), meridionaal (noord-zuid of zuid-noord gericht) en mengvormen van die twee. Er is sprake van een zonaal circulatietype wanneer er een min of meer gelijkmatige west-oostelijke stroming bestaat tussen een subtropisch hogedrukgebied boven de Noord-Atlantische Oceaan en een subpolair lagedrukgebied. Stationaire, blokkerende hogedrukgebieden tussen 50 en 65 graden NB zijn kenmerkend voor het meridionale circulatietype, alsmede alle troggen met een noord-zuid noord-zuidas. De lucht beweegt dan van N naar Z en andersom.  Bij het gemengde circulatietype zijn de zonale en meridionale stromingscomponenten ongeveer even groot.

Om een en ander overzichtelijk te houden gebruik ik alleen indeling in de drie circulatietypen en bekijk de ontwikkeling daarvan in de loop der jaren voor elk van de seizoenen. Omdat de signalen nogal volatiel zijn heb ik ze ‘gesmoothed’ met een loess filter om langjarige trends goed te kunnen zien.

Fig.13    Data: DWD

Figuur 13 laat voor de lente zien dat de laatste decennia de meridionale luchtcirculatie sterk gestegen is ten koste van de gemengde luchtcirculatie. In de praktijk betekent dat meer dagen met warme lucht uit het zuiden. De meridionale aanvoer van lucht uit het N blijft gelijk. De zonale circulatie blijft min of meer stabiel.

Fig.14    Data: DWD

De zomer laat in de afgelopen halve eeuw een toename zien van het meridionale circulatietype ten koste van het gemende type. Dat betekent meer dagen met aanvoer van warme lucht vanuit het Z. De zonale (westelijke) circulatie neemt ’s zomers vanaf 1950 geleidelijk af.

Fig.15    Bron: DWD

De herfst vertoont vanaf ongeveer 1980 een sterk wisselend beeld,  vergeleken met de stabiele situatie in de 100 jaren daarvoor. De meridionale circulatie neemt vanaf de jaren ’80 sterk toe, vooral ten koste van de zonale circulatie.

Fig.16    Bron: DWD

In tegenstelling tot de drie andere seizoenen is de winter vanaf 1900 een toonbeeld van gelijkmatigheid. De drie circulatietypen tonen de laatste decennia geen opvallend verloop. Wel gaan de ups en downs van de meridionale circulatie – net zoals in de herfst- ten koste van de zonale circulatie.

Al met al is het beeld op basis van de GWL zodanig dat er de afgelopen decennia sprake is van sterke verschuivingen in de luchtcirculatie boven West Europa. Dat was overigens ook – maar in mindere mate- het geval aan het begin van de 20^e eeuw.

Fig.17    Bron: Pik Report No 119

Figuur 17 zet twee trends tegenover elkaar: de afnemende invloed van GWL met een noordcomponent en de toenemende invloed van GWL met een zuidcomponent. Veel duidelijker kan het niet: de luchtcirculatie in Europa wordt vanaf de jaren ’50 zuidelijker en dus warmer.

Wiskundestudent Jippe Hoogeveen heeft niet zo lang geleden een klimaatmodel gemaakt dat net zoals de GWL van de Duitse meteorologische dienst rekening houdt met de brongebieden en  aanvoerroutes van binnenstromende lucht, maar dat heeft hij weten te vatten in een wiskundig model. Samen met zijn vader (wiskundige aan de universiteit van Utrecht) heeft hij zijn klimaatmodel kunnen publiceren.

Zijn model is sterk afwijkend van de methodiek van de GWL, daarom is vergelijken lastig. Maar het model van Jippe laat zien dat er  de afgelopen decennia sprake is van grote verschuivingen in de brongebieden van met name de maritieme luchtsoorten:  meer W en SW en minder NW en N. De continentale luchtsoorten laten vooral in de zomer een toename uit de warmere hoek zien (E, SE, S). Dat alles betekent meer aanvoer van warmte dan voorheen.

Fig.18    Bron: Hoogeveen et al 2022

Dat het model van Jippe en zijn vader deugt laat figuur 18 zien. In rood de voorspelling van de jaarlijkse temperatuur in De Bilt door het model, in zwart de gemeten jaartemperatuur in De Bilt. De grafiek toont een zeer goede fit met een R kwadraat van 0,85. Volgens Jippe verklaart luchtcirculatie het grootste deel van de recente opwarming in Nederland. Opvallend was dat in het onderzoek van Hoogeveen et al CO2 in de recente opwarming van Nederland geen rol van betekenis lijkt te spelen.

Zoals we al zagen leverde in De Bilt de sterke opwaartse trend van de globale straling  van 1989 t/m 2022 een toename van de straling op van 16 W/m2. Hoe zit dat met CO2? Verdubbeling van CO2 veroorzaakt volgens het IPCC een stralingsforcering van 3,7 W/m2 . Recente gedetailleerde radiative transfer calculations kwamen tot een vergelijkbare schatting. (Myhre et al 2017). Dus 16 W/m2 erbij door extra zonlicht is heel veel.

Fig.19    Naar: Myhre et al 2017

Sterker nog: in de periode 1989-2022 steeg het atmosferisch CO2-gehalte op Mauna Loa van 353,2 ppm naar 418,6 ppm (bron: NOAA). De grafiek van figuur 19 laat zien dat die stijging van het CO2-gehalte van 1989 t/m 2022 met 65,4 ppm een extra radiative forcing teweeg bracht van 0,9 W/m2. Vergelijk dat eens met de extra radiative forcing van de toegenomen globale straling in De Bilt in de periode 1989-2022 met 16 W/m2!

Conclusie: de opvallende grote stijging van de temperatuur op station De Bilt tussen 1989 en 2022 is grotendeels het gevolg van de sterk toegenomen globale straling, die op zijn beurt vooral het gevolg is van wijzigingen in de grootschalige luchtcirculatie boven Europa en het noordelijk deel van de Atlantische Oceaan.

Een volgende keer meer aandacht voor de langgolvige kant van de stralingsbalans en de situatie in Europa.

Psycholoog dr. Jordan B. Peterson bevraagt klimaatonderzoeker dr. Richard Lindzen. Ze duiken in de feiten en onzin over klimaatverandering, de modellen die worden gebruikt om die te voorspellen, de sombere toestand van de academische wereld en de gepolitiseerde wereld van de “professionele” wetenschap.

Jordan Peterson kan heel goed luisteren en toont een verbazingwekkend vermogen om kennis en begrip te krijgen over een onderwerp dat niet zijn vakgebied is. Richard Lindzen is een erudiete wetenschapper die meer dan 30 jaar professor was aan het MIT en zijn sporen meer dan verdiend heeft met honderden publicaties. Hij neemt geen blad voor de mond en prikt met gemak door onzinnige beweringen over klimaatverandering. Het is een genoegen om naar iemand te luisteren die weet waar hij het over heeft als het over klimaat gaat.

Een juweeltje van een gesprek tussen twee hoog opgeleide en intelligente mannen. Niet te missen!

Hoe werkt dat eigenlijk, een overheid die af wil van heel veel boerenbedrijven en een bevolking die daar in meerderheid weinig voor voelt? Bekijk de documentaire van Marijn Poels en verbaas je over de macht van bepaalde ngo’s die samen met de overheid de dienst lijken uit te maken. En verbaas je ook over de Postcodeloterij die daar een opvallende rol in speelt.

 Meten is weten; een ouderwets gezegde, maar in deze tijd van klimaatmodellen meer dan ooit waar.  Er is een grote behoefte aan meer duidelijkheid over de invloed van de mens op klimaatveranderingen. Kennis van de natuurlijke variabiliteit van klimaat is daarbij onontbeerlijk, en dus zijn lange homogene meetreeksen nodig. Maar die zijn schaars.

De Labrijnreeks, ook wel de Zwanenburg/De Bilt reeks genoemd,  is een van de langste temperatuurreeksen ter wereld en begint in 1706. Het is een samengestelde reeks, een samenvoeging van reeksen van een zestal stations: Delft, Rijnsburg, Zwanenburg, Haarlem, Utrecht en De Bilt:

Fig.1    Bron: KNMI

Het laatste deel van de reeks zijn de KNMI data vanaf 1898 in De Bilt, vrijwel volledig de ‘moderne’  data vanaf 1901 bevattend. Er is al het nodige werk verzet om van deze heterogene datareeksen een min of meer consistente reeks te maken. Dat heeft niet zozeer te maken met de kwaliteit van de thermometers (die over het algemeen goed was) maar vooral met de wisselende omstandigheden op de diverse locaties. Ook de tijdstippen van meten waren niet voortdurend gelijk, er moest dus ook gestandaardiseerd worden. Al met al kan gesteld worden dat je de oudere Labrijndata met enige voorzichtigheid moet hanteren. Voor gemiddelde jaartemperaturen geldt dat voorbehoud natuurlijk minder dan voor maand- en vooral etmaaldata.

Fig.2    Bron: KNMI

Bovenstaande grafiek geeft de Labrijn jaartemperatuur (verticaal) weer t/m 2022 met loessfilter van 10 jaar. Opvallend is dat de reeks eind jaren ’80 van de vorige eeuw een trendbreuk laat zien.  Van 1706 tot eind jaren ‘80 is er nauwelijks sprake van een langjarige trend van 0,1 °C/100 jaar, daarna schiet de trend omhoog naar bijna 3 °C/100 jaar.

Jammer is dat de Labrijnreeks niet nog wat langer is, want vooral in het begin van de reeks is er in het temperatuursignaal sprake van behoorlijk grote schommelingen. Gelukkig beschikt Metoffice in GB over een reeks die nog langer is dan de Labrijnreeks, de zogenaamde CET (Central England Temperature). De CET  dataset is de langste instrumentele temperatuurreeks  in de wereld. De reeks van gemiddelde dagelijkse gegevens begint in 1772 en die van gemiddelde maandelijkse gegevens al in 1659. Deze dag- en maandtemperaturen zijn representatief voor het gebied dat in figuur 3 omgeven is door de rode lijn.

Fig.3    Bron: Parker et al 1991

Natuurlijk zijn ook in de CET-reeks de oudere data onbetrouwbaarder dan recente data. Maar de CET reeks is- anders dan de Labrijnreeks waar de afzonderlijke stationreeksen aan elkaar geplakt zijn- een ensemble van stations. Dat verkleint de afwijkingen doordat er gemiddeld wordt. Bovendien maken we gebruik van jaargemiddelde temperaturen.

Fig.4    Data: Metoffice

https://www.metoffice.gov.uk/hadobs/hadcet/data/meantemp_monthly_totals.txt

Dat de Labrijnreeks van goede kwaliteit is wordt duidelijk als we de Labrijnreeks (loess met 10-jarig filter) vergelijken met CET  loess smoothing.

Fig.5    Data: KNMI en Metoffice

CET en Labrijn volgen elkaar heel mooi, met hier en daar wat opvallende afwijkingen zoals de laatste 10 jaren en begin 19^e eeuw, waarschijnlijk als gevolg van het feit dat De Nederlandse stations wat meer landinwaarts liggen. De correlatiecoëfficiënt R2 tussen beide lijnen is hoog: 0.81 .

Bij vergelijking van beide reeksen in figuur 5 blijkt dat het eerste deel van de CET reeks tot 1706 (dat ontbreekt in de Labrijnreeks) laat zien dat er vanaf 1695 gedurende een periode van ongeveer 40 jaar er sprake is van een opvallende stijging, vergelijkbaar met die aan het einde van de grafiek:

Fig.6    Data: Metoffice

De eerste sprong (t1) laat een temperatuurstijging zien van 1,5 °C, t2 van ongeveer 1,2  °C.  Beide periodes beslaan ongeveer 40 jaar. Gezien het synchrone verloop van Labrijn en CET lijkt het aannemelijk dat de temperatuursprong van 1695-1734 ook in Nederland aanwezig geweest moet zijn.

Fig.7    Data: Metoffice

Minstens zo interessant is de grafiek van figuur 7, die de 30-jarige moving trend weergeeft. Van elk jaar is de lineaire trend (in °C/100y) weergegeven van het voorafgaande tijdvenster van 30 jaren. Zo is in het laatste jaar 2022 de trend (1,81 °C/100y) weergegeven van de periode 1993-2022, et cetera. Daarbij vallen 2 toppen op: 1720 met een moving trend van 4,92 °C/100y en 2007 met een moving trend van 5,23 °C/100y. De dalen zijn 1756 met een negatieve moving trend van -3,90 °C/100y en 1892 met een moving trend van -3,98 °C/100y.

Fig.8    Bron: Luterbacher et al 2016

Wil je nog verder terug in de geschiedenis dan moet je je toevlucht zoeken in zogenaamde proxies, zoals Luterbacher et al deden. De grafiek in figuur 8 laat  hun reconstructie zien van o.a. de zomertemperatuur (blauwe lijn) in Europa vanaf 800 n.Chr.  De gemiddelde Europese zomertemperaturen in de relatief warme MCA (‘Medieval Climate Anomaly’, 900–1200 CE) waren bijna net zo hoog als in de laatste helft van de 20^e eeuw.

De relatief koele periode Little Ice Age bestrijkt de periode 1250 tot 1800. Na het Last Maunder Minimum (rond 1700) veert de zomertemperatuur weer op. Zowel  de Labrijnreeks als de CET volgen dus grotendeels de opwarming die plaatsvindt na de LIA. Interessante vraag blijft natuurlijk in welke mate het temperatuurverloop van de Labrijnreeks (sinds 1706) en de CET (sinds 1659) het gevolg zijn van natuurlijke factoren en in welke mate van antropogene factoren.

Fig.9    Bron: Wikipedia

Tot de 20^e eeuw is de invloed van de mens op het klimaat beperkt geweest. De schommelingen van de zomertemperatuur in Europa zoals gereconstrueerd door Luterbacher et al (fig.8) volgen opvallend goed de schommelingen in de zonneactiviteit. Paleoklimatologen als Bas van Geel van UvA waren al vroeg op de hoogte van de correlatie tussen de zonneactiviteit en het klimaat, zie onder andere hier.

Die correlatie loopt via schommelingen in de hoeveelheid 14C zoals te zien figuur 9. Horizontaal is de tijd uitgedrukt in jaren Before Present. BP = 1950 in de christelijke jaartelling. Verticaal zijn de schommelingen in het 14C isotoop. Minder zonneactiviteit (zonnevlekken) betekent een zwakker zonnemagnetisch veld. Daardoor kan kosmische straling beter de aarde bereiken en wordt in de hogere atmosfeer meer van het onstabiele 14C gevormd. Dat verschijnsel is niet alleen bruikbaar om organisch materiaal te dateren (aan de hand van de halfwaardetijd van 14C) maar lijkt ook het klimaat te sturen. Onder andere de Deense natuurkundige Svensmark heeft daar in het recente verleden onderzoek naar gedaan.

Vanaf het begin van de 20^e eeuw wordt de emissie van CO2 in de atmosfeer een factor, alhoewel het nog vele jaren zou duren voordat CO2 in klimaatverandering de status kreeg die het momenteel heeft.. Het IPCC is van mening dat het zeker is dat vanaf 1950 de dominante oorzaak voor klimaatverandering de CO2-emissie door de mens is. Merkwaardig is dat daarmee automatisch de natuurlijke oorzaken van klimatologische variabiliteit terzijde worden geschoven en genegeerd.

Fig.10    Data: European Environment Agency  en NOAA (vanaf 1959)

De grafiek van figuur 10 laat de groei zien van CO2 in de atmosfeer sinds 1750 (278 ppm) t/m 2022 (419 ppm). Een groot deel van die groei wordt toegeschreven aan de uitstoot van CO2 door menselijke activiteiten. Een kleiner deel is afkomstig van de oceanen die bij stijgende watertemperaturen CO2 afstaan aan de atmosfeer.

Fig.11    Bron: ourworldindata

De grafiek van figuur 11 toont de berekende CO2 emissies (voornamelijk afkomstig van fossiele brandstoffen en cementindustrie). De CO2-emissies komen pas goed op gang vanaf 1950. Het dipje rechts bovenaan is het gevolg van verminderde emissies als gevolg van de coronacrisis. Het dunne verticale zwarte lijntje toont de hoogte van CO2-emissies in 1988.

Lutenbacher et al toonden in figuur 8 de reconstructie van de zomertemperatuur in Europa vanaf 800 tot 2000 n. Chr. Daar is te zien dat vanaf het Late Maunder Minimum begin 18^e eeuw er vanaf 1900 er sprake is van een sterke stijging van de zomertemperatuur. Ik was benieuwd of dat in de Labrijnreeks ook te zien is.

Fig.11    Data: KNMI

In de grafiek van figuur 11 valt een paar dingen op. Zo is de zomer van 2018 de warmste van de hele reeks, maar ik tel sinds 1706 maar liefst 19 jaren met een zomertemperatuur in Nederland hoger dan 18 °C. De nummers 2 en 3 zijn de zomers van 1826 en 1947, beide 18,7 °C. 1947 uiteraard op basis van de niet gehomogeniseerde data van het KNMI.

Verder is opvallend dat de zomertemperaturen van 1706 t/m 1988 geen enkele trend vertonen.  Vanaf 1989 is er sprake van een duidelijke trendbreuk en wordt de trend 3,55 °C/eeuw. Dat is  opvallend, zeker als je dat vergelijkt met de ontwikkeling in de tijd van het atmosferisch CO2 (figuur 10) en van de CO2-emissies (figuur 11). Er is veel meer aan de hand dan toename van het CO2-gehalte. Met CO2 als dominante factor verklaar je het temperatuurverloop in Engeland en Nederland niet zoals we dat hiervoor zagen bij de CET-reeks en Labrijnreeks. De lange vlakke trend en de sterke trendbreuk rond 1988-1989  kunnen niet verklaard worden met behulp van het broeikasgas CO2. Zie o.a. hier.

Paleoklimatoloog Bas van Geel schreef al in 2004:  “In feite is het klimaat nooit werkelijk stabiel geweest en het is dan ook maar de vraag of de nadruk die de laatste jaren gelegd wordt op het versterkt broeikaseffect wel terecht is. Het is namelijk niet onwaarschijnlijk dat de temperatuurstijging van de afgelopen decennia voornamelijk veroorzaakt werd door de toegenomen zonneactiviteit. De aanwijzingen die we – via paleo-onderzoek – hebben voor de zon als sturende factor bij natuurlijke klimaatveranderingen zijn in elk geval zo sterk dat het hedendaagse en het toekomstige klimaat waarschijnlijk ook voor een belangrijk deel door fluctuaties in zonneactiviteit worden, respectievelijk zullen worden bepaald.”