Slab Detachment - Desprinderea slabului in urma incetarii subductiei

S-a folosit pachetul de elemente finite CitcomS.py Version 3.0.1 (Tan et al., 2006) disponibil de la Computational Infrastructure for Geodynamics (CIG http://geodynamics.org). Calculele au fost realizate in cadrul unui domeniu sferic, 3D (θ, φ, r), unde θ este latitudinea, φ este longitudinea si r este raza. Raza interioară corespunde unei adancimi de 1900 km iar raza exterioara suprafetei Pamantului. Dimensiunes pe latitudine si longitudine este de 23° si 28°. Acest domeniu este impartit in mod egal in 129 elemente pe dimensiunea radiala si 150 x 257 in planul latitudine x longitudine. Aceasta configuratie corespunde unei rezolutii de 15 x 17 x 12 km. Conditiile la frontiera sunt urmatoarele: la limitele de sus exista viteze impuse ce corespund miscarilor tectonice inregistrate pentru regiunea Panonica acum 30 Ma. Placa oceanica se subduce cu 1.5cm/yr, iar placa de deasupra se misca cu aceeasi viteza. Limitele superioara si inferioara sunt presupuse isotermice. La partea inferioara se considera "free-slip" iar lateral se considera pereti reflectori. Structura termica initiala este descrisa de un strat termic cu varsta variabila, in care litosfera veche are 200 Myr iar litosfera tanara are 10Myr. Mantaua este separata in 4 straturi: litosfera (cu grosime variabila), mantaua superioara (cu grosime variabila), zona de tranzitie (410 - 670 km) si mantaua inferioara (670 - 1300 km). Litosfera prezinta grosime variabila, intre 100km pentru litosfera tanara si intre 100 - 200 km pentru litosfera cea mai veche. Suprafata slabului este modelata folosind particule trasatoare (particulele magenta). Vitezele particulelor trasatoare sunt calculate prin interpolarea a opt viteze nodale folosind functii bilineare. TRacerii sunt advectati folosind metoda de punc mediu (cu o precizie de ordin 2).