On the Mediterranean conveyor belt system

The Mediterranean Sea is a semi-enclosed sea, connected to the Atlantic Ocean through the Gibraltar Strait and subdivided in two different regions by the Sicily Strait. The geographical extension of the basin, the surface heat flux, and the water inflow from the Gibraltar Strait are some of the basic factors determining its horizontal and vertical circulation. In the Mediterranean strong salinity and temperature zonal gradients contribute to maintain the zonal-vertical circulation, while meridional-vertical cells are equally forced by winds and deep water mass formation in three regions, the Gulf of Lyon, the southern Adriatic and the Cretan Sea areas. The objective of this thesis is to study how these cells combine together to form the Mediterranean conveyor belt system. This has never been attempted before so the conclusions are necessarily preliminary. In the first part we discuss the vertical zonal and meridional circulation by reconstructing the Wust Mediterranean vertical salinity and temperature structures in an attempt to evaluate the water mass structure consistent with modern data. Our results confirm that Wust depicted vertical circulation from scarce data is reproduced by the past 27 years observations. The structure of both meridional and zonal circulations was discussed using velocity vertical streamfunctions with two different methods. The first one, eulerian, allowed us to observe vertical structures that were already reported in the literature. Recent studies in the Atlantic Ocean have shown that gyres and eddies have an important influence in the isopycnal vertical circulation. This is called the residual circulation which was computed in this study for the first time. A possible interpretation of horizontal connection between the meridional and zonal cells was discussed using horizontal streamfunction. In the last part of the thesis we have been developing an idealized numerical model to study the vertical circulation in the Mediterranean.

Assimilazione dati di precipitazione in un modello meteorologico ad alta risoluzione (MOLOCH): sviluppo ed applicazione alla previsione di eventi intensi in Liguria

L'Italia è spesso sede di eventi d'intensa precipitazione, frequentemente associati ad alluvioni con conseguente perdita di vite umane e gravi danni economici. È quindi di fondamentale importanza poter prevedere questi eventi con un adeguato anticipo. Allo stato attuale, i modelli meteorologici non permettono sempre di raggiungere tale obbiettivo e, di conseguenza, è in atto un'intensa attività di ricerca al fine di di renderne più accurata la previsione, sia attraverso il miglioramento dei modelli stessi, sia sviluppando l'assimilazione dati, la quale riduce l'incertezza della condizione iniziale da cui parte la previsione. All'interno di questo contesto, la tesi si prefigge l'obiettivo di studiare gli effetti dell'assimilazione di dati di precipitazione, effettuata mediante uno schema di nudging, nel modello non idrostatico MOLOCH. Al fine di ottimizzare lo schema e di valutarne l'impatto, sono stati simulati tutti gli eventi di maltempo di ottobre e novembre del 2014 che hanno interessato la Liguria, area frequentemente soggetta ad alluvioni. Dalla sistematica verifica dei risultati, effettuata sia qualitativamente che mediante numerosi metodi statistici (tra cui la tecnica SAL basata sull'individuazione dei nuclei di precipitazione), si riscontra un generale miglioramento della previsione della precipitazione anche se limitato alle prime ore dopo la fine del periodo di assimilazione. L'impatto dello schema di nudging varia a seconda dell'evento e, in particolare, si osserva una certa correlazione tra il miglioramento ottenuto e la tipologia di evento, come descritto dalla teoria dell'equilibrio convettivo e come riportato in alcuni studi analoghi recentemente comparsi in letteratura. La ricaduta del miglioramento della previsione meteorologica è stata valutata anche in termini di impatti al suolo accoppiando il modello meteorologico con il modello idrologico utilizzato dal Centro Funzionale della Regione Liguria, con risultati abbastanza positivi.

Analisi modellistica di un intenso ciclone di tipo tropicale nel Mediterraneo (7-8 novembre 2014): predicibilità, caratteristiche termodinamiche e confronto con osservazioni da satellite

Il lavoro di tesi riguarda lo studio dettagliato di un ciclone di tipo tropicale (tropical like cyclone, TLC) verificatosi nel Canale di Sicilia nel novembre 2014, realizzato attraverso un'analisi modellistica effettuata con i modelli BOLAM e MOLOCH (sviluppati presso il CNR-ISAC di Bologna) e il confronto con osservazioni. Nel primo capitolo è fornita una descrizione generale dei cicloni tropicali e dei TLC, indicando come la formazione di questi ultimi sia spesso il risultato dell'evoluzione di cicloni extratropicali baroclini nel Mediterraneo; sono identificate le aree geografiche con i periodi dell'anno maggiormente soggetti all'influenza di questi fenomeni, riportando un elenco dei principali TLC verificatisi nel Mediterraneo negli utlimi tre decenni e lo stato dell'arte sullo studio di questi eventi. Nel secondo capitolo sono descritte le modalità di implementazione delle simulazioni effettuate per il caso di studio e presentati i principali prodotti dell'analisi modellistica e osservazioni da satellite. Il terzo capitolo si apre con la descrizione della situazione sinottica e l'analisi osservativa con immagini Meteosat e rilevazioni radar che hanno permesso di ricostruire la traiettoria osservata del TLC. In seguito, viene dapprima fornito l'elenco completo delle simulazioni numeriche, quindi sono presentati alcuni dei più importanti risultati ottenuti, dai quali emerge che la previsione della traiettoria e intensità del TLC differisce notevolmente dalle osservazioni. Tenendo conto della bassa predicibilità che ha caratterizzato l'evento, nel quarto capitolo è descritto il metodo usato per ricostruire in maniera ottimale la traiettoria, utilizzando spezzoni da varie simulazioni, che ha permesso un confronto più realistico con i dati osservati e un'analisi dei processi fisici. Nel quinto capitolo sono riportati i principali risultati di alcuni test mirati a valutare l'impatto di aspetti legati all'implementazione delle simulazioni e ad altre forzanti fisiche.

Deterministic and probabilistic verification of multi-model meteorological forecasts on the subseasonal timescale

In questo studio, un multi-model ensemble è stato implementato e verificato, seguendo una delle priorità di ricerca del Subseasonal to Seasonal Prediction Project (S2S). Una regressione lineare è stata applicata ad un insieme di previsioni di ensemble su date passate, prodotte dai centri di previsione mensile del CNR-ISAC e ECMWF-IFS. Ognuna di queste contiene un membro di controllo e quattro elementi perturbati. Le variabili scelte per l'analisi sono l'altezza geopotenziale a 500 hPa, la temperatura a 850 hPa e la temperatura a 2 metri, la griglia spaziale ha risoluzione 1 ◦ × 1 ◦ lat-lon e sono stati utilizzati gli inverni dal 1990 al 2010. Le rianalisi di ERA-Interim sono utilizzate sia per realizzare la regressione, sia nella validazione dei risultati, mediante stimatori nonprobabilistici come lo scarto quadratico medio (RMSE) e la correlazione delle anomalie. Successivamente, tecniche di Model Output Statistics (MOS) e Direct Model Output (DMO) sono applicate al multi-model ensemble per ottenere previsioni probabilistiche per la media settimanale delle anomalie di temperatura a 2 metri. I metodi MOS utilizzati sono la regressione logistica e la regressione Gaussiana non-omogenea, mentre quelli DMO sono il democratic voting e il Tukey plotting position. Queste tecniche sono applicate anche ai singoli modelli in modo da effettuare confronti basati su stimatori probabilistici, come il ranked probability skill score, il discrete ranked probability skill score e il reliability diagram. Entrambe le tipologie di stimatori mostrano come il multi-model abbia migliori performance rispetto ai singoli modelli. Inoltre, i valori più alti di stimatori probabilistici sono ottenuti usando una regressione logistica sulla sola media di ensemble. Applicando la regressione a dataset di dimensione ridotta, abbiamo realizzato una curva di apprendimento che mostra come un aumento del numero di date nella fase di addestramento non produrrebbe ulteriori miglioramenti.

Simulazione della risposta dei radiometri satellitari nelle microonde AMSU-B e MHS per il retrieval della precipitazione nevosa e studio di sensibilità dell’algoritmo 183-WSL

Il lavoro di questa tesi è focalizzato sulla valutazione della sensibilità delle microonde rispetto a differenti idrometeore per le frequenze a 89 GHz e 150 GHz e nella banda di assorbimento del vapor d'acqua a 183.31 GHz. Il metodo di indagine consiste nell'utilizzo del modello di trasferimento radiativo RTTOV (Eyre, 1991) per simulare radianze dei canali dei sensori satellitari nelle microonde Advanced Microwave Sounding Unit-B (AMSU-B) e Microwave Humidity Sounder (MHS). Le simulazioni basate sul modello RTTOV si sono focalizzate su tre dataset indipendenti, forniti da ECMWF. Il primo passo tiene conto di una selezione di categorie dei profili atmosferici basato su una distinzione della fase delle idrometeore, LWP, IWP e WVP, con sottoclassi terra e oceano. La distinzione in diverse categorie permette di valutare la sensibilità di ciascuna frequenza utilizzata nelle simulazioni al variare del contenuto di acqua e ghiaccio. Un secondo approccio è usato per valutare la risposta di ciascuna frequenza nei casi di nevicate sulla terraferma. Questa indagine ha permesso lo sviluppo di un nuovo algoritmo prototipale per la stima dell'intensità di precipitazione nevosa basato su una serie di test a soglia e una equazione di combinazione lineare multipla che sfrutta una combinazione dei canali più sensibili alla snowfall: 150, 186 e 190 GHz. Una prima verifica su casi di studio pre-selezionati di snowstorm sembra fornire risultati promettenti. Infine è stato realizzato uno studio di sensibilità dell’algoritmo 183-WSL (Laviola and Levizzani, 2011) utilizzando le simulazioni di RTTOV con precipitazione/non precipitazione come predittori discreti e con le intensità di precipitazione come predittori continui. Le simulazioni RTTOV rivelano una sovrastima delle radianze in presenza di profili di pioggia e ciò potrebbe essere dovuto alle approssimazioni e parametrizzazioni adottate nel modello RTTOV-SCATT usato per la risoluzione dello scattering in presenza di precipitazione.

Particelle di cellulosa come nuclei di ghiacciamento in nube

La formazione del ghiaccio nelle nubi avviene prevalentemente per nucleazione eterogenea, grazie alla presenza di nuclei di ghiacciamento (Ice Nucleating Particles, INP), ovvero particelle di aerosol (prevalentemente polveri minerali) in grado di favorire la solidificazione di una gocciolina o il passaggio diretto dalla fase vapore alla fase ghiaccio. Recentemente, l'interesse si è esteso anche all'aerosol di tipo biologico (funghi, spore, etc.), in grado di agire come INP a temperature più elevate rispetto all'aerosol minerale. Il lavoro sperimentale di questa tesi, svolto presso il laboratorio del gruppo Nubi e Precipitazioni dell’ISAC-CNR di Bologna, ha avuto come obiettivo lo studio della capacità della cellulosa di agire come INP. Esso si inserisce in un’attività di ricerca internazionale dedicata allo studio degli INP (progetto Ice NUclei research unIT, INUIT, Germania). Il lavoro sperimentale ha riguardato la messa a punto di due diversi sistemi di generazione dell'aerosol di cellulosa, la caratterizzazione dimensionale delle particelle, la loro osservazione al microscopio elettronico e la preparazione dei filtri per le misure di INP. I risultati ottenuti hanno evidenziato che le proprietà nucleanti della cellulosa sono inferiori rispetto alle polveri minerali ma paragonabili ad altri materiali, come la cenere vulcania e l'aerosol marino. Quindi, in aree ricche di vegetazione o dove le polveri minerali non sono abbondanti, la cellulosa potrebbe costituire un importante materiale per la formazione del ghiaccio nelle nubi miste. Misure preliminari hanno evidenziato come le particelle di cellulosa con dimensioni inferiori a 0,5μm risultano meno attive nella nucleazione del ghiaccio rispetto a quelle di dimensioni maggiori. I risultati ottenuti sono stati presentati al Congresso PM2016 (Roma, 17-19 maggio 2016) e saranno presentati in un poster alla prossima European Aerosol Conference (Tours, Francia, 4-6 settembre 2016).

From the traditional verification approach to the spatial methods: a study of applicability of the SAL metric in the framework of the WMO project MesoVICT

Increasing in resolution of numerical weather prediction models has allowed more and more realistic forecasts of atmospheric parameters. Due to the growing variability into predicted fields the traditional verification methods are not always able to describe the model ability because they are based on a grid-point-by-grid-point matching between observation and prediction. Recently, new spatial verification methods have been developed with the aim of show the benefit associated to the high resolution forecast. Nested in among of the MesoVICT international project, the initially aim of this work is to compare the newly tecniques remarking advantages and disadvantages. First of all, the MesoVICT basic examples, represented by synthetic precipitation fields, have been examined. Giving an error evaluation in terms of structure, amplitude and localization of the precipitation fields, the SAL method has been studied more thoroughly respect to the others approaches with its implementation in the core cases of the project. The verification procedure has concerned precipitation fields over central Europe: comparisons between the forecasts performed by the 00z COSMO-2 model and the VERA (Vienna Enhanced Resolution Analysis) have been done. The study of these cases has shown some weaknesses of the methodology examined; in particular has been highlighted the presence of a correlation between the optimal domain size and the extention of the precipitation systems. In order to increase ability of SAL, a subdivision of the original domain in three subdomains has been done and the method has been applied again. Some limits have been found in cases in which at least one of the two domains does not show precipitation. The overall results for the subdomains have been summarized on scatter plots. With the aim to identify systematic errors of the model the variability of the three parameters has been studied for each subdomain.

Bottom pressure time series study: the ANTARES deep sea observatory case

Questo elaborato riguarda l'analisi delle variazioni di pressione sul fondale del Mar Mediterraneo. Il primo obiettivo era quello di identificare la variabilità a basse frequenze della pressione sul fondo. I dati osservati sono stati raccolti presso l'osservatorio ANTARES che si trova 40 km al largo di Tolone e ricopre un'area di 180x180 m. Sono state analizzate tre serie temporali, L1, L3, L8, e da esse è stata ricostruita una serie temporale dall'1 gennaio 2009 al 31 dicembre 2013. Sono state applicate delle procedure di despiking, media su 30 minuti e gap-filling. Sono stati proposti due metodi di detiding, una media su tre giorni e il filtro di Doodson. Successivamente è stato eseguito un detrending. Un confronto tra l'analisi armonica e il filtro di Doodson mostra una similarità nei risultati. Per questa ragione, per continuare il lavoro, è stata utilizzata la serie a cui era stato tolto il trend e le maree tramite il filtro di Doodson. Togliere le maree ha reso possibile la rimozione dei segnali ad alta frequenza per potersi focalizzare sui segnali rimasti a bassa frequenza. Uno spettro sulla serie ha mostrato i segnali stagionali dominanti insieme alla cosiddetta Chandler wobble con una frequenza di 14 mesi.

Fluid dynamic analysis of trees influence in dispersion of pollutant in urban street canyon

degli elementi vegetali nella dinamica e nella dispersione degli inquinanti nello street canyon urbano. In particolare, è stato analizzata la risposta fluidodinamica di cespugli con altezze diverse e di alberi con porosità e altezza del tronco varianti. Il modello analizzato consiste in due edifici di altezza e larghezza pari ad H e lunghezza di 10H, tra i quali corre una strada in cui sono stati modellizati una sorgente rappresentativa del traffico veicolare e, ai lati, due linee di componenti vegetali. Le simulazioni sono state fatte con ANSYS Fluent, un software di "Computational Fluid Dynamics"(CFD) che ha permesso di modellizare la dinamica dei flussi e di simulare le concentrazioni emesse dalla sorgente di CO posta lungo la strada. Per la simulazione è stato impiegato un modello RANS a chiusura k-epsilon, che permette di parametrizzare i momenti secondi nell'equazione di Navier Stokes per permettere una loro più facile risoluzione. I risultati sono stati espressi in termini di profili di velocità e concentrazione molare di CO, unitamente al calcolo della exchange velocity per quantificare gli scambi tra lo street canyon e l'esterno. Per quanto riguarda l'influenza dell'altezza dei tronchi è stata riscontrata una tendenza non lineare tra di essi e la exchange velocity. Analizzando invece la altezza dei cespugli è stato visto che all'aumentare della loro altezza esiste una relazione univoca con l'abbassamento della exchange velocity. Infine, andando a variare la permeabilità delle chiome degli alberi è stata trovatta una variazione non monotonica che correla la exchange velocity con il parametro C_2, che è stata interpretata attraverso i diversi andamenti dei profili sopravento e sottovento. In conclusione, allo stadio attuale della ricerca presentata in questa tesi, non è ancora possibile correlare direttamente la exchange velocity con alcun parametro analizzato.

Analisi meteorologica di eventi con precipitazione intensa nell'Italia Centrale

In this work we analyze five events of intense precipitation that affected Central Italy between 2011 and 2015. These events caused important flood phenomena, with a significant impact on both the economy and people. This study mainly uses data from ERA5, the latest reanalysis produced by the ECMWF (European Center for Medium-range Weather Forecast), in addition to data from the meteorological network managed by Centro Funzionale of Protezione Civile d’Abruzzo. In the first part is reported an overview of the general issues that characterize the study of intense precipitation events, with a deepening of some recently analyzed cases. Subsequently the data and tools used are reported, with a description of the events studied and the method of investigation. We then go to analyze the synoptic framework in which the event developed and to investigate the causes that generated it. Finally, the precipitation pattern was studied and the degree of reproducibility by the ERA5 reanalysis was tested through a comparison with the meteorological network.

Caratterizzazione della struttura e delle misure di precipitazione da una rete di ponti commerciali a microonde in Emilia Romagna

Lo scopo di questa tesi consiste nell’analisi della rete di ponti a microonde, o CML commercial microwave link, gestita da Lepida Spa; la rete è costituita da un gruppo di connessioni telefoniche, collocate sul territorio dell’Emilia Romagna, ed utilizzata all’interno di questo lavoro come strumento di misura della precipitazione. In particolare è stata elaborata una descrizione delle caratteristiche di misura, trasmissione di segnale e geometria delle antenne coinvolte ed eseguita un’analisi statistica di verifica della misura, sviluppata a partire dal confronto tra le osservazioni estrapolate dalla rete, ottenute per mezzo dell’applicazione dell’algoritmo RAINLINK, e le misure effettuate da radar e pluviometri presenti sullo stesso territorio monitorato. Le analisi iniziali hanno evidenziato la peculiarità delle caratteristiche dei ponti all’interno della rete: processi di controllo hanno dimostrato la presenza di vuoti d’informazione relativamente ad alcuni ponti ed uno studio legato alla trasmissione del segnale ha mostrato la presenza di effetti di attenuazione di cui non è possibile offrire una parametrizzazione, anche in assenza di precipitazione. La misura di precipitazione da rete CML è stata confrontata separatamente con i dati estratti da radar e pluviometro, cercando in entrambi i casi di presentare la migliore strategia possibile di comparazione tra misure con diverse caratteristiche spaziali; dal confronto emerge come l’accordo tra le diverse osservazioni risulti maggiore durante i mesi estivi. I risultati hanno evidenziato la presenza di un sottogruppo di strumenti, collocati generalmente nelle zone in cui la rete risulta meno densa, le cui misure presentano un grado di accordo inferiore con il dato osservato da pluviometro e radar, relativamente al resto della rete. É stato inoltre osservato che, la maggior parte dei ponti della rete, tende a sottostimare la presenza di precipitazione.

Long-period oscillations in GPS Up time series. Study over the European/Mediterranean area

The surface of the Earth is subjected to vertical deformations caused by geophysical and geological processes which can be monitored by Global Positioning System (GPS) observations. The purpose of this work is to investigate GPS height time series to identify interannual signals affecting the Earth’s surface over the European and Mediterranean area, during the period 2001-2019. Thirty-six homogeneously distributed GPS stations were selected from the online dataset made available by the Nevada Geodetic Laboratory (NGL) on the basis of the length and quality of the data series. The Principal Component Analysis (PCA) is the technique applied to extract the main patterns of the space and time variability of the GPS Up coordinate. The time series were studied by means of a frequency analysis using a periodogram and the real-valued Morlet wavelet. The periodogram is used to identify the dominant frequencies and the spectral density of the investigated signals; the second one is applied to identify the signals in the time domain and the relevant periodicities. This study has identified, over European and Mediterranean area, the presence of interannual non-linear signals with a period of 2-to-4 years, possibly related to atmospheric and hydrological loading displacements and to climate phenomena, such as El Niño Southern Oscillation (ENSO). A clear signal with a period of about six years is present in the vertical component of the GPS time series, likely explainable by the gravitational coupling between the Earth’s mantle and the inner core. Moreover, signals with a period in the order of 8-9 years, might be explained by mantle-inner core gravity coupling and the cycle of the lunar perigee, and a signal of 18.6 years, likely associated to lunar nodal cycle, were identified through the wavelet spectrum. However, these last two signals need further confirmation because the present length of the GPS time series is still too short when compared to the periods involved.

Sviluppo e validazione di mappe di copertura nevosa per lo studio della variabilità del manto nevoso in Emilia Romagna nel periodo 2000-2020

Il telerilevamento satellitare costituisce una delle tecniche di osservazione maggiormente efficace nel monitoraggio e nello studio dell'estensione del manto nevoso. Il manto ricopre un ruolo di estrema importanza quale risorsa idrica per gli ecosistemi e per le applicazioni agricole e industriali. Inoltre, è un indicatore climatologico in un contesto di cambiamenti climatici regionali e globali. In questo senso, la copertura nevosa è da considerarsi come una importante risorsa economica e sociale. Lo scopo della presente tesi è di produrre mappe di copertura nevosa giornaliere per la stima dell'estensione del manto nevoso della regione Emilia-Romagna e di analizzare la sua variabilità spazio-temporale nel periodo 2000-2020. Le mappe sono state sviluppate sulla base dei prodotti di neve, M*D10A1, del sensore MODIS, che consistono in mappe di classificazione della copertura in funzione dell'indice NDSI. Inizialmente, è stato costruito un albero decisionale con criteri a soglia multipla in grado di rielaborare la classificazione della superficie del sensore e produrre mappe di copertura nevosa sulla base di tre classi: neve, no-neve e nube. L'accuratezza delle mappe è stata validata tramite test statistici effettuati confrontandole con i dati di altezza del manto nevoso in situ di 24 stazioni meteorologiche, per il periodo di compresenza dei dati 2000-2012. I risultati della procedura di validazione hanno mostrato come, in generale, vi sia buon accordo tra il dato satellitare e la rispettiva osservazione al suolo, soprattutto nei pressi di stazioni lontane da vegetazione sempreverde e/o di ambiente urbano. Infine, è stata valutata la variabilità climatologica dell'estensione del manto nevoso tramite l'elaborazione degli indici di neve SCF, SCD e diversi indici SCA. L'attenzione è stata particolarmente focalizzata sugli indici di massima estensione invernale del manto, del valore mediano e del numero di giorni con estensione superiore al 39.5% della regione.

New assessments of global sea level rise from tide gauges

Estimates of global sea-level change rates based on observations from Tide Gauges (TGs) show a long-term global mean sea-level rise (GMSLR) of 1÷2 mm/yr for the 20th century. The considerable scatter in these estimates is mainly attributable to the uneven distribution of the TG sites and to several physical phenomena that cause local sea level to deviate from the global mean, or to affect the TG record through land subsidence or uplift. The main cause of vertical ground motion on a regional space scale is the response of the Earth to past ice loads, called Glacial Isostatic Adjustment (GIA), which is often modelled and corrected for. In this work, a simple average approach was used to revisit two past estimates based on small sets of long, high-quality TG records in view of the longer record available, employing a newer GIA model (ICE-6G) from Peltier et al. [2015]. The value of GMSLR obtained from both sets is (1.5±0.4) mm/yr. In addition, a much larger set of TGs was used to estimate the contemporary (post 1993) GMSLR using satellite estimates from Cazenave et al. [2018] as a benchmark, in an attempt to understand how a simple average approach could perform for larger sets. The resulting estimate of (3.4÷3.5)±0.2 mm/yr (depending on the GIA correction applied) is comparable to the satellite result of (3.1±0.3) mm/yr.

Anisotropic seismic tomography with the reversible jump Markov chain Monte Carlo

The established isotropic tomographic models show the features of subduction zones in terms of seismic velocity anomalies, but they are generally subjected to the generation of artifacts due to the lack of anisotropy in forward modelling. There is evidence for the significant influence of seismic anisotropy in the mid-upper mantle, especially for boundary layers like subducting slabs. As consequence, in isotropic models artifacts may be misinterpreted as compositional or thermal heterogeneities. In this thesis project the application of a trans-dimensional Metropolis-Hastings method is investigated in the context of anisotropic seismic tomography. This choice arises as a response to the important limitations introduced by traditional inversion methods which use iterative procedures of optimization of a function object of the inversion. On the basis of a first implementation of the Bayesian sampling algorithm, the code is tested with some cartesian two-dimensional models, and then extended to polar coordinates and dimensions typical of subduction zones, the main focus proposed for this method. Synthetic experiments with increasing complexity are realized to test the performance of the method and the precautions for multiple contexts, taking into account also the possibility to apply seismic ray-tracing iteratively. The code developed is tested mainly for 2D inversions, future extensions will allow the anisotropic inversion of seismological data to provide more realistic imaging of real subduction zones, less subjected to generation of artifacts.