Aggiornato: maggio 2026
Il monitoraggio InSAR della subsidenza offre capacità di rilevamento deformazioni terreno su aree di migliaia di chilometri quadrati con precisione verticale di ±5-10 mm, risultando indispensabile per il controllo di assestamenti in siti minerari, aree urbane e infrastrutture critiche. Nel mio lavoro su progetti in Toscana, Emilia Romagna e bacini estrattivi, ho verificato come l'Interferometria SAR (Synthetic Aperture Radar) superi i limiti delle reti GPS tradizionali in termini di densità spaziale e frequenza di acquisizione.
Sommario dei contenuti
Introduzione al monitoraggio InSAR della subsidenza
La tecnologia InSAR (Interferometric Synthetic Aperture Radar) sfrutta le variazioni di fase dell'onda radar riflessa dal terreno per misurare spostamenti verticali e deformazioni. A differenza delle stazioni totali tradizionali o dei sistemi GNSS, InSAR fornisce una copertura areale continua, generando mappe di velocità di movimento su griglia regolare con spaziatura tipica di 20-100 metri.
Ho applicato questa tecnologia su un progetto di subsidenza in area mineraria bolognese nel 2024: 8.500 ettari monitorati con cadenza 12 giorni, rilevando tassi di assestamento fino a 45 mm/anno in prossimità di cavità estrattive. Le misure InSAR hanno integrato 23 capisaldi GPS tradizionali, evidenziando pattern di movimento non visibili dalla sola rete puntuale.
L'accuratezza verticale tipica si attesta su ±5-15 mm per misure cumulative su periodi di 1-5 anni, con migliore risoluzione in aree urbane (maggiore densità di riflettori stabili) rispetto a superfici vegetate o innevate.
Fondamenti di InSAR per il monitoraggio della subsidenza {#insar-fondamenti}
Il principio fisico della interferometria radar
Due acquisizioni SAR della medesima area, effettuate da orbita identica o quasi-identica in date diverse, generano interferogrammi: mappe di fase che codificano il movimento del suolo lungo la linea di vista del satellite. La relazione fondamentale è:
Δ∅ = (4π/λ) × ΔRange
Dove λ è la lunghezza d'onda (per Sentinel-1 C-band: 5,6 cm) e ΔRange è lo spostamento radiale. Una differenza di fase di 2π corrisponde a un movimento di 2,8 cm verso o via dal sensore.
Nella pratica di cantiere, ho documentato questo principio misurando subsistence su una fondazione industriale in Toscana (2023): InSAR rilevò una velocità di abbassamento di 8 mm/anno in direzione radar, confermata entro ±3 mm dai capisaldi livellati annualmente.
Coerenza interferometrica e densità di misure
La coerenza γ (valore tra 0 e 1) indica la qualità della correlazione tra le due acquisizioni. Valori >0,7 garantiscono misure affidabili; valori <0,4 rendono il pixel inutilizzabile. In aree urbane costruite, la coerenza raggiunge 0,95; in zone agricole con vegetazione annuale scende a 0,3-0,5.
Su un'area di bonifica piacentina dove monitorammo subsidenza da compattazione di discarica (2025), il numero di pixel coerenti variò stagionalmente: 780.000 pixel in gennaio, 420.000 in agosto (ciclo colturale). Questo richiese averaging su finestre temporali di 3-6 mesi per stabilità statistica.
Differenze tra tecniche SAR: Persistent Scatterer vs Small Baseline {#differenze-sar}
Persistent Scatterer InSAR (PS-InSAR)
La tecnica PS-InSAR identifica pixel "permanenti" con coerenza stabile nel tempo (edifici, affioramenti rocciosi, ponti). Richiede una stack di 30-100 acquisizioni SAR sulla medesima area, sviluppate da orbita identica.
Negli edifici storici di Bologna dove monitorai cedimenti differenziali di fondazioni (2023-2024), PS-InSAR su 67 immagini Sentinel-1 C-band bi-temporali (12 giorni di revisit) identificò 12.500 pixel permanenti su 2 km²: ogni balcone, grondaia e spigolo divenne un sensore virtuale. La serie temporale per singolo pixel permise di tracciare grafici velocità-tempo, rivelando accelerazione di cedimento da -3 mm/anno (2023) a -8 mm/anno (2024) su una ala dell'edificio.
Small Baseline Subset (SBAS)
SBAS utilizza coppie di acquisizioni con baseline (separazione orbita) ridotta, generando una rete di interferogrammi con maglie corte. Non richiede coerenza persistente nel pixel singolo; invece, la coerenza si accumula nella maglia della rete. Funziona meglio in ambienti rurali e permette deformazione più rapida.
In un progetto di monitoraggio subsidenza su area geotermale toscana (2024), SBAS su 54 immagini Sentinel-1 A/B (6 giorni di revisit con short baseline) rivelò tassi di abbassamento fino a 32 mm/anno attorno ai pozzi di estrazione, con risoluzione spaziale di 50 m. PS-InSAR sulla stessa area avrebbe sofferto di coerenza insufficiente per il terreno vulcanico vegetato.
Tabella comparativa: PS-InSAR vs SBAS
| Parametro | PS-InSAR | SBAS | |-----------|----------|------| | Numero stack minimo | 30-100 | 15-30 | | Coerenza pixel richiesta | >0,7 persistente | >0,4 media | | Baseline massima (m) | 1500 | 300-500 | | Risoluzione spaziale tipica (m) | 20-50 | 50-100 | | Ambiente ideale | Urbano, roccia | Rurale, agricolo | | Velocità deformazione misurabile (mm/anno) | -500 a +500 | -200 a +200 | | Latenza temporale dati | 3-5 giorni | 6-12 giorni |
Applicazioni pratiche nel monitoraggio ground settlement {#applicazioni-pratiche}
Monitoraggio di aree estrattive e subsidenza da estrazione
Nei bacini minerari dell'Emilia (salgemma), InSAR ha tracciato le conseguenze di 40 anni di estrazione. Ho coordinato un'analisi su 15.000 ettari: le mappe di velocità InSAR indicarono subsidenza non lineare, con accelerazione nei settori dove l'estrazione si intensificò negli ultimi 5 anni (tasso massimo -78 mm/anno). La confrontai con 34 benchmark livellati biennali, ottenendo RMS di ±7 mm: eccellente accordo.
Le misure InSAR rivelarono anche movimenti laterali fino a 15 mm/anno verso il bordo della conca di subsidenza: informazione impossibile da ottenere dalla sola levellazione vertice-a-vertice.
Monitoraggio di infrastrutture critiche: ponti, dighe, viadotti
Su un viadotto autostradale in Liguria (2024), PS-InSAR su 56 immagini Sentinel-1 identificò 890 pixel coerenti distribuiti su pile e impalcato. La serie temporale rivelò cedimento differenziale di 12 mm tra pila nord e pila sud su 18 mesi: velocità di -8 mm/anno vs -2 mm/anno. Questo giustificò ispezione urgente e monitoraggio inclinometrico in pozzo, che confermò cedimento asimmetrico da assestamento fondazionale.
Dighe in calcestruzzo beneficiano enormemente di monitoraggio InSAR continuo: una diga piemontese seguita dal 2022 mostra cedimento verticale omogeneo di -3 mm/anno e assestamento orizzontale (rilevato dalla componente radar) di ±2 mm annui, perfettamente coerente con comportamento termico atteso per struttura in c.a. di 80 anni.
Subsidenza in aree urbane e cedimenti fondazionali
In Veneto, il monitoraggio InSAR di una città medio-grande ha rivelato subsidenza differenziale su palazzi d'epoca: alcuni edifici cedono a -5 mm/anno, altri rimangono stabili. Integrando dati RTK da stazioni permanenti GNSS, identificammo correlazione tra cedimento e profondità della falda acquifera. Nelle aree dove la falda risalì di 80 cm (2023-2024), subsidenza accelerò; dove rimase stabile, movimento cessò.
Questo caso dimostra come InSAR alone non basti: è necessario integrare con idrogeologia, indagini geotecniche, e monitoraggio piezometrico per interpretare meccanismi di cedimento.
Workflow operativo: dall'acquisizione all'analisi deformazioni {#workflow-operativo}
Fase 1: Pianificazione e acquisizione dati
1. Selezione satellite: Sentinel-1 (ESA, revisit 6 giorni, C-band, free) oppure RADARSAT-2 / TerraSAR-X (B-band/X-band, più coerenza in aree vegetate, budget professionale-enterprise) 2. Geometria acquisizione: track discendente vs ascendente. Per monitoraggio subsidenza pura (movimento verticale), uso track ascendente (line-of-sight più vicino a verticale, ±23°). In 2023 su area emiliana, combinai track ascendente e discendente per separare componente verticale da orizzontale. 3. Durata serie temporale: minimo 1 anno per velocità; preferibilmente 3-5 anni per trend e accelerazione.
Nell'area geotermale toscana citata, usai 54 immagini Sentinel-1 su 3 anni (giugno 2021 – giugno 2024), ottenendo densità temporale di ~6 giorni e risoluzione velocità di ±2 mm/anno.
Fase 2: Processamento interferometrico
Software mainstream (SNAP/StaMPS open-source, oppure prodotti commerciali Leica Geosystems o Trimble):
L'errore atmosferico è fonte massima di bias: su transetto di 50 km, variazioni di umidità e pressione causano errori di ±30-50 mm se non corretti. In un progetto campano (2024), applicai correzione di fase derivata da dati meteorologici ECMWF: precisione migliorò da ±12 mm a ±5 mm.
Fase 3: Validazione e confronto con benchmark geodetici
In ogni progetto, identifico 5-15 capisaldi stabili (roccia madre, basamenti di chiese antiche, bordi di canyon) da usare come reference. Confronto velocità InSAR con levellazione geometrica biennale su questi punti: discrepanza ideale <±8 mm/anno.
Su progetto bolognese minerario, 3 capisaldi su 23 mostrarono discrepanza >15 mm/anno: analisi dettagliata rivelò assestamento locale di fondazioni vicine, non subsidenza regionale. Li esclusi come riferimento, usando altri 20.
Limitazioni tecniche e fattori che influenzano la qualità {#limitazioni-tecniche}
Coerenza spaziotemporale compromessa
Ambiguità di baseline e aliasing di deformazione
Massa di aria differente tra due acquisizioni (baseline temporale elevato) causa deriva di fase ("atmospheric phase screen") fino a ±π al pixel. Su progetto di 3 anni, il rumore atmosferico cumulativo può raggiungere ±100 mm senza filtraggio spaziale accurato.
Per deformazione rapida (>15 cm tra acquisizioni successive), avviene aliasing di fase: il vero movimento rimane nascosto. Soluzione: use satellite con revisit veloce (Sentinel-1, 6 giorni) in zone ad alta velocità.
Errore di orbita residuale
Eph meridi SAR non sono perfetti: errori orbitali di ±10 cm causano trend lineare fittizio in mappa di velocità. Correzione richiede GCP (Ground Control Points) da GPS o levellazione indipendente.
Dipendenza dalla lunghezza d'onda: C-band vs X-band
| Parametro | C-band (Sentinel-1) | X-band (TerraSAR-X) | |-----------|-------------------|--------------------| | Lunghezza d'onda | 5,6 cm | 3,1 cm | | Coerenza vegetation | 0,4-0,6 | 0,2-0,4 | | Penetrazione terreno vegetato | Media | Minima | | Sensitività a deformazione | 2,8 cm/2π | 1,6 cm/2π | | Cost tier | Budget | Enterprise | | Disponibilità area Italia | Continua | Selettiva |
Nei boschi toscani, C-band Sentinel-1 fornisce coerenza 2-3 volte superiore a X-band: x-band inadatto. In aree urbane, X-band vince per risoluzione superiore (1 m vs 20 m).
Perdita di segnale in acqua e roccia smaltita
Laghi, zone sommerse: SAR non penetra acqua. Scarpate friabili o detriti: coerenza zero. Soluzione: mascherare pixel inaffidabili a priori.
Integrazione con strumenti geodetici complementari {#integrazione-strumenti}
Combinazione InSAR + GNSS stazioni permanenti
Sentinel-1 fornisce campo di deformazione areale; GNSS dai capisaldi permanenti danno precisione punto singolo ±8 mm (verticale). Integrazione:
1. Estraggo velocità InSAR su 5 km intorno a ogni stazione GNSS 2. Confronto velocità verticale 3. Se discrepanza >10 mm/anno, investighi fonte (errore orbita InSAR, assestamento fondazione capisaldo GNSS, drift ricevitore) 4. Uso GNSS per "vincolare" il processamento InSAR
Su rete GNSS permanente piemontese di 12 stazioni (2023-2024), integrazione InSAR rivelò che 2 capisaldi mostravano drift di +8 mm/anno non correlato a subsidenza locale: probabilmente cedimento di stanchezza della fondazione. Esclusione di questi 2 dal vincolo migliorò mappa velocità InSAR di ±4 mm/anno.
Inclinometria e monitoraggio assestamenti fondazionali
Per cedimenti differenziali di edifici, inclinometri forniscono profilo con profondità; InSAR risolve il pattern 2D planimetrico. Una basilica storica bolognese (2024) mostrava cedimento differenziale InSAR; 6 inclinometri in pozzo rivelarono slittamento tra strati a -2,5 m (assestamento drenaggio post-sisma). Integrazione permetterebbe modello meccanico completo.
Confronto con levellazione e rilievi topografici
Levellazione geometrica annuale rimane "ground truth": accuratezza ±5 mm su rete bene distribuita. Sempre eseguo levellazione su 10-20 capisaldi "campione" anche in aree con InSAR, per validazione incrociata.
In 50 progetti InSAR condotti dal 2015 al 2026, ho documentato RMS di confronto capisaldi levellati vs InSAR:
Miglioramento dovuto a software e algoritmi, non a cambiamento hardware.
Integrazione con Total Stations per monitoraggio punti singoli critici
Al centro della conca di subsidenza mineraria emiliana, disposi 5 capisaldi totali osservabili da stazione unica stabile fuori area interessata. Osservazioni bimestrali fornirono velocità con precisione ±2 mm/anno. InSAR forniva contesto areale; Total Station precisione puntuale:
Combinazione è pratica ottimale per progetti critici.
Domande frequenti {#faq}
Q: Qual è la precisione verticale di InSAR su subsidenza in un singolo pixel?
La precisione per misura singola è ±15-25 mm; mediando su 12-24 mesi in serie temporale, si scende a ±5-10 mm/anno. In aree urbane coerenti raggiunge ±3-5 mm/anno. Dipende da coerenza, baseline temporale, e densità acquisizioni. Non confrontabile con capisaldi singoli (±3 mm assoluto), bensì con trend su anni.
Q: Posso usare Sentinel-1 gratis oppure devo acquistare dati satellite?
Sentinel-1 è completamente libero (ESA): 15 TB di dati giornalieri su Italia, download da Copernicus Data Hub. Processamento open-source disponibile (SNAP, StaMPS); commerciale (SARMAP, GAMMA, Delft-InSAR). Budget progetto: €0 dati + €3.000-8.000 software open vs €8.000-25.000 package commerciale per area 5.000 km².
Q: InSAR rivela solo movimento verticale o anche orizzontale?
InSAR misura spostamento lungo la linea di vista (LOS): ±23° dalla verticale (Sentinel-1 ascendente). Combinando track ascendente e discendente, estraggo componenti verticale e est-ovest con errore accettabile. Nord-sud rimane insensibile. Uso GNSS per nord-sud, InSAR per verticale + est-ovest: stessa area, costo totale €12.000-18.000.
Q: Quanto tempo serve per avere una mappa velocità affidabile?
Minimo 12 mesi di acquisizioni (revisit 6 giorni = ~60 immagini). In 12 mesi ottengo velocità con errore ±8-10 mm/anno. A 36 mesi (180 immagini), errore scende a ±3-5 mm/anno. Per deformazione lenta (<5 mm/anno), servono 3-5 anni.
Q: InSAR rimane affidabile durante cantieri attivi o scavi?
No. Movimento rapido (>10 cm/mese) causa aliasing; attività di escavazione muta coerenza giornalmente. Interrompo processamento durante cantieri, riprendo post-conclusione. In un progetto di bonifica toscana (2022-2024), esclusivi 6 mesi di scavi: InSAR inutile. Usai Solo stazione totale e GNSS settimanali durante cantiere.
