How Uncertain Are Estimates of Marine Carbon Dioxide Removal? Insights from a Simplified Model

Technische Samenvatting: Hoe onzeker zijn de schattingen van mariene koolstofverwijdering? Inzichten uit een vereenvoudigd model

Probleemstelling
Technieken voor mariene koolstofverwijdering (mCDR), specifiek oceaanalkaliniteitverhoging (OAE) en directe oceaanverwijdering (DOR), vertrouwen op het verlagen van de partiële druk van CO2 in het oppervlaktewater om atmosferische CO2 te verminderen. Hoewel de theoretische maximale efficiëntie van deze interventies bekend is (bijv. ~0,8 mol CO2 per mol alkaliniteit voor OAE), is de tijdsverloop van koolstofverwijdering sterk gevoelig voor fysische en biogeochemische processen. Huidige schattingen van de effectiviteit van mCDR zijn sterk afhankelijk van oceanografische modellen, maar de onzekerheden die inherent zijn aan deze modellen—met name met betrekking tot verticale menging, lucht-zee gasuitwisseling en watermassa-eigenschappen—zijn slecht begrepen. Bestaande wereldwijde modellen maken vaak gebruik van vaste parametrisaties die natuurlijke variabiliteit kunnen onderschatten, en hun complexiteit maakt het moeilijk om de specifieke mechanismen te isoleren die de onzekerheid in CDR-berekeningen drijven. Dit gebrek aan gekwantificeerde onzekerheid vormt een fundamentele uitdaging voor Monitoring, Reporting, and Verification (MRV)-kaders en de timing van netto koolstofneutraliteit voor mCDR-projecten.

Methodologie
Om deze onzekerheden systematisch te onderzoeken, maakt de studie gebruik van een eendimensionaal (1D), abiotisch oceaanmodel. Deze benadering met verminderde complexiteit isoleert verticale processen (turbulente diffusie en lucht-zee gasuitwisseling) terwijl horizontale advectie, subductie en biologische interacties worden genegeerd. Het model evolueert de totale alkaliniteit (TA) en opgeloste anorganische koolstof (DIC) over de tijd, waarbij de oppervlakte-CO2flux wordt berekend op basis van de lucht-zee gasuitwisseling.

De studie maakt gebruik van twee primaire analytische kaders:

1. Gevoeligheidssimulaties: Individuele modelparameters (bijv. de diepte van de menglaag, windsnelheid, gasoverdrachtcoëfficiënt, pycnocline-diffusiviteit, achtergrond TA/DIC, injectiediepte) worden verstoord met waarden die consistent zijn met de geobserveerde natuurlijke variabiliteit. Deze worden vergeleken met een referentisimulatie om de "CDR-gevoeligheidsratio" ($r_{CDR}$) te bepalen, gedefinieerd als de ratio van de CDR in een gevoeligheidsrun ten opzichte van de referentierun.
2. Monte Carlo-ensembles: Meer dan 1.000 simulatieparen (baseline en interventie) worden gelijktijdig uitgevoerd, waarbij parameters willekeurig worden getrokken uit statistische distributies die de geobserveerde oceaanvariabiliteit vertegenwoordigen. Deze benadering beoordeelt hoe gecombineerde onzekerheden in fysische processen en watermassa-eigenschappen zich manifesteren in CDR-schattingen.

De simulaties worden uitgevoerd voor oceaan-dieptes van 200 m en 2.000 m over een periode van 50 jaar om de impact van diepte op de persistentie van onzekerheid te evalueren.

Belangrijkste Resultaten
De studie identificeert drie duidelijke temporele regimes die de CDR-onzekerheid beheersen:

1. Initiële snelle evenwichtsvorming ($t \ll \tau_{mix}$): Binnen de eerste jaren zijn CDR-schattingen zeer gevoelig voor parameters die de initiële lucht-zee gasuitwisselings-tijdsschaal ($\tau_{air-sea}$) controleren, inclusioneel de diepte van de menglaag, windsnelheid, gasoverdrachtcoëfficiënten en de achtergrond carbonaatchemie (TA, DIC, temperatuur). Tijdens deze fase kunnen CDR-schattingen variëren met een factor twee, afhankelijk van parameterbiases.
2. Verlies van alkaliniteit uit de menglaag ($t \approx \tau_{mix}$): Terwijl verticale menging alkaliniteit uit het oppervlak verwijdert, worden onzekerheden gedomineerd door de interactie tussen de lucht-zee uitwisselings-tijdsschaal en de meng-tijdsschaal ($\tau_{mix}$). Onzekerheden in de pycnocline-diffusiviteit ($\kappa_{pyc}$) worden cruciaal, waardoor de variabiliteit van CDR-schattingen piekt (tot 30%) 2–3 jaar na de injectie.
3. Diepe menging ($t \gg \tau_{mix}$): Op de lange termijn blijven onzekerheden decennia voortbestaan als de tijdsschaal voor menging in de diepe oceaan vergelijkbaar is met of sneller is dan de lucht-zee evenwichtsvorming. De studie stelt vast dat onzekerheid in de pycnocline-diffusiviteit de dominante factor is voor de langetermijn-CDR-onzekerheid, omdat deze zowel de fractie van de geëquilibreerde alkaliniteit vóór diepe subductie als de snelheid bepaalt waarmee niet-geëquilibreerde diepe wateren naar het oppervlak terugkeren.

Kwantitatieve Bevindingen:

• Omvang van de Onzekerheid: In de eerste jaren kunnen CDR-schattingen tot 100% (een factor twee) verschillen van de referentisimulatie. Na 20–50 jaar neemt dit af tot 10–20%, hoewel aanzienlijke onzekerheid blijft bestaan voor diepe oceaaninjecties.
• Injectiediepte: Het injecteren van alkaliniteit onder de menglaag (zelfs met slechts 20–40 m verschil) resulteert in aanzienlijk lagere efficiëntie en hogere onzekerheid, waarbij sommige schattingen zelfs na 50 jaar nog steeds ver onder de referentisimulatie liggen.
• Timing van Netto Neutraliteit: De tijd die nodig is voor een interventie om netto koolstofnegatief te worden (rekening houdend met de kosten van de levenscyclusanalyse) is zeer variabel. Terwijl een referentiegeval ongeveer 1,4 jaar suggereert, laat de Monte Carlo-ensemble een brede verdeling zien: ~30% van de simulaties bereikt netto neutraliteit in minder dan een jaar, 50% duurt langer dan 2 jaar, en ~4% duurt langer dan 50 jaar. Deze "long tail" is primair gerelateerd aan diepe injecties of snelle menging die vroege evenwichtsvorming verhindert.
• Diepteafhankelijkheid: Onzekerheden houden langer aan in 2.000 m diepe oceanen vergeleken met 200 m diepe oceanen, waar de gehele waterkolom binnen het simulatievenster kan evenwicht bereiken.

Betekenis en Claims
Het artikel stelt dat onzekerheden in mCDR-schattingen niet louter transient zijn, maar decennia kunnen voortduren vanwege de tijdsschalen van de menging in de diepere lagen relatief aan de lucht-zee gasuitwisseling. De auteurs beweren dat deze resultaten een fundamentele uitdaging identificeren voor regionale en wereldwijde mCDR-simulaties, waarbij de leidende fysische processen vaak slecht geconstreerd zijn.

Specifiek benadrukt de studie dat:

• MRV-implicaties: De onzekerheid in de tijdsevolutie van koolstofverwijdering heeft een directe impact op de timing van netto neutraliteit, wat cruciaal is voor koolstofkrediet-kaders. Zonder goed gekwantificeerde onzekerheden kan het toekennen van koolstofkredieten voortijdig of onnauwkeurig zijn.
• Modelbeperkingen: Huidige oceanografische modellen vertrouwen vaak op parametrisaties voor turbulentie en gasuitwisseling die onderhevidig zijn aan aanzienlijke ruimtelijke en temporele biases. De studie suggereert dat het valideren van deze specifieke fysische processen (menging en gasoverdracht) belangrijker is voor mCDR dan voor algemene klimaatprojecties.
• Injectiestrategie: De resultaten benadrukken de noodzaak om ervoor te zorgen dat alkaliniteit strikt binnen de oppervlaktemenglaag wordt geïnjecteerd om de efficiëntie te maximalen en de onzekerheid te minimalen.

De auteurs hanteren een bescheiden reikwijdte door op te merken dat hun 1D-model horizontale advectie en biologische feedback uitsluit, wat verdere onzekerheden zou kunnen introduceren. Ze concluderen echter dat zelfs in een abiotisch systeem de onzekerheden die voortvloeien uit fysische processen alleen al aanzienlijk genoeg zijn om ernstige verdere investigatie te rechtvaardigen, om te garanderen dat mCDR-interventies daadwerkelijk koolstofnegatief zijn op decenniale tijdschalen.