Come si muove una calotta glaciale quando l’oceano la spinge?

Immagina una gigantesca lastra di ghiaccio che fluttua sopra l’acqua, come un surfista su un’onda. Gli scienziati hanno creato un modello matematico per capire come questa “calotta glaciale” interagisce con l’oceano. Il modello combina due mondi: da un lato, il ghiaccio si comporta come un materiale elastico, capace di deformarsi ma anche di tornare alla forma originale; dall’altro, l’acqua è un fluido viscoso, cioè scorre con una certa resistenza. Per farli “parlare” tra loro, i ricercatori hanno imposto due regole ai bordi del ghiaccio: una kinematica, che dice come la superficie del ghiaccio si muove in risposta al flusso d’acqua, e una dinamica, che descrive la pressione esercitata dall’acqua sul ghiaccio. Per semplificare i calcoli, hanno trasformato il problema in un dominio fisso usando la formulazione ALE (Arbitrary Lagrangian‑Eulerian), così da poter studiare l’evoluzione nel tempo senza dover ridisegnare continuamente la superficie. Il risultato è una serie di equazioni che garantiscono l’esistenza di soluzioni “forti” per brevi intervalli di tempo, a patto che i dati iniziali siano abbastanza regolari (appartenenti allo spazio H²). Le difficoltà principali sono nate dai termini non lineari che coinvolgono derivate verticali e dalle pressioni di ordine superiore che agiscono sull’interfaccia ghiaccio‑acqua. Questo lavoro apre la strada a simulazioni più realistiche del comportamento delle calotte glaciali in un clima in cambiamento.

Il modello proposto si basa su un accoppiamento tra le equazioni di Navier‑Stokes per il fluido viscoso e le equazioni di elasticità lineare per la calotta glaciale. La formulazione ALE consente di mappare il dominio fisico deformabile su un dominio di riferimento fissato, introducendo un campo di velocità di trasporto che separa il moto della mesh dalla velocità del fluido. Questo approccio elimina la necessità di ricostruire continuamente la geometria dell’interfaccia, riducendo la complessità computazionale. Le condizioni al contorno kinematicamente imposte garantiscono la continuità della velocità normale attraverso l’interfaccia, mentre quelle dinamiche includono il bilancio di forze: la pressione idrodinamica dell’acqua è bilanciata dalla tensione elastica del ghiaccio e dalle forze di gravità. Le stime a priori sono ottenute mediante tecniche di energia e interpolazione in spazi di Sobolev, con particolare attenzione ai termini non lineari che coinvolgono derivate verticali di ordine superiore. La presenza di pressione di ordine due sull’interfaccia richiede l’uso di operatori di proiezione per mantenere la regolarità delle soluzioni. Il risultato è una dimostrazione di esistenza locale di soluzioni forti in H², che apre la strada a studi di stabilità a lungo termine e a simulazioni numeriche più accurate di fenomeni di scioglimento e avanzamento delle calotte glaciali in risposta a variazioni di carico oceanico.