Geophysikalisches Geländepraktikum SS 2004 Gravimetrie-Messungen in Vulkaneifel.

1. Geophysikalisches Geländepraktikum SS 2004 Gravimetrie-Messungen in Vulkaneifel. Leitung: PD Dr. N. Bagdassarov Inhalt: 1. Schwerefeld: Messungen und Inversion 2. Messungen von Richtungen und Strechen mittels ein Lasertachymeter

2. Gravimetrie Zum Versuch der Höhenbestimmung mit Altimetern bleibt zu bemerken, daß die Schwankungen (bei Wiederholung der Messung) bis zu einigen Metern betrugen. Dabei ist ein Teil der Ungenauigkeit mit Sicherheit auf die Luftdruckschwankungen durch eine aufziehende Regenfront zurückzuführen. Da für die Gravimetrie jeder Meter einen Fehler von 0.3 mGal bewirkt und die Meßgenauigkeit bei ca. 0.01 mGal liegt, läßt sich als Ergebnis der Bemühungen sagen, daß das Verfahren für unsere Anwendung wohl untauglich ist. Die Grundlage der Angewandten Gravimetrie ist die Massenanziehung nach dem Newtonschen Gravitationsgesetz. Danach ist die Kraft zwischen zwei punktförmigen Massen M und m im Abstand r

3. Meßgröße ist die Schwerebeschleunigung ( im Mittel ca. 980 Gal). Einheit von ist . Wir werden mit LaCoste-Romberg Gravimeter messen. b b b Interessant für uns sind nur die Änderungen der Schwere durch Dichteinhomogenitäten im Untergrund. Effekte durch verschiedene geographische Breiten, unterschiedliche Meereshöhen oder Topographie müssen reduziert werden, um vergleichbare Werte zu erlangen.

4. Datenkorrektur • Die Breiten korrigiert man mit der Normalschwereformel (für das Referenzellipsoid). • Für die Höhe wählt man eine Bezugsfläche, etwa das Meeresniveau, legt die Stationen rechnerisch auf dieses Niveau und korrigiert außerdem für die Massen dazwischen. • Topographische Reduktion

5. Z= w²·d=w²·R(j)·cos(j) g=b+Z Konstante Schwerepotetialwert = -6,264·107 m²/s² Referenzellipsoid: a=6378,16 km c=6356,775 km Abplattung: f=(a-c)/a=1/298,25 983,3 Normalschwerewerte g =978,049 ·(1+ +0,0052884·sin2(j) -0,0000059·sin²(2j) j 978 Eifel: j=50°14‘51,2“ N l=6°44‘24.4“E

6. Die Reduktion gemessener Schwerewerte auf eine Bezugsfläche Gemessene Schwerewert setzt aus folgenden Beiträgen zusammen: g=g(j)+dgF+dgB+dgT+dgG+dgt g(j) ist die Normaleschwere, Breiabhängigkeit von g dgj(j) ~0,81·sin(j)·y(km) = 0,81 mGal/km im Richtung Norden. dgF ist Freiluftreduktion = -2·g(j)·h/RE = -0,3086·h mGal, wobei h ist positive nach oben in m. dgB Gesteinsplatten- oder Bougersche Reduktion = 2·p·G·r·h = +0,0419 r·h mGal, wobei r ist Dichte in g/cm³, h ist in m.

7. dgT Geländereduktion: Wegen einer Erhebung des Volumelements DV mit einer Grundfläche DF und der Hohe h im Schwerpunktabstand s: Masse des Elements = r·DV DdgT =-G· r·DV/s²·sin(f)=- =-1/2·G· r·DF/s³ ·h² DdgT=-1/2·G· r·(DFi/si³ ·hi²) dgG Massenunregelmessigkeiten f Sin(f) =h/(2·s)

9. Gezeiten Dgt zeitlichen Änderungen Das Verhältnis zwischen zwei Amlituden von Sonnengezeiten und Mondgezeiten ergibt sich durch: Amlitude der Gezeiten vom Mond Amlitude der Gezeiten von der Sonne

10. 1. Ebenfalls korrigiert werden muß der Gerätegang und der Einfluß der Gezeiten. Beides erreicht man durch wiederholte Messung eines Basispunktes oder, wie in unserem Fall, durch überlappende Messung von Punkten. Die dadurch entstehende zeitliche Verlaufskurve wird von den Werten abgezogen. Das Verfahren heißt Gangreduktion. 2. Die danach verbleibende Anomalie nennt man Bouguer-Anomalie und ist die Grundlage für den Rückschluß auf die Dichteverteilung im Untergrund unter der Bezugsfläche. 3. Als besonders wichtig für eine genaue Messung ist uns dabei die genaue Horizontierung und ein stabiler Stand des Gravimeters aufgefallen.

11. Inversion x1/2=z(22/3-1) N=Dm · G/g · 1/(x²+z²)1/2 Dgmax=Dm · G/z² Die Eingrenzung auf eine Lösungen ist möglich aufgrund der geologischen Situation: tektonische Stufe, magmatische Spalteintrusion, sowie die Annahme realistischer Dichtewerte. Dm=4/3 · p ·R³ ·(r-ro) s |b|=Dm·G/s² Dg=bz=b ·sin(j)=Dm · G·z/(z²+x²)3/2

12. Inversion Betrachten wir eine dünne horizontale Platte mit der Dicke a und der Dichte r. Schwerewirkung des Plattenelements mit der Oberfläche df dDg=dbz=G·r·a/s²·df·sin(j)=G·r·a·df“= G·r·a·dW wobei W ist der Raumwinkel. Für die gesamte Platte Dg= G·r·a·W Für eine allseitig unbegrenzte horizontale PlatteDgB= 2p·G·r·a·W

13. Inversion Beiderseits unbergrenzte horizontale Massenstreifen erscheinen auf der Einheitskugel unter dem Flachenwinkel 2y W=4py/2p=2y DgB= G·r·a·2y Dg=2G·(r-ro) ·a·y= 2G·a· (r-ro) ·[arctg(x/z)-p/2]

14. Laser-Tachymeter ELTA R55 Das Laser-Tachymeter ist ein Peilgerät mit dem man vertikale und horizontale Winkel sowie Abstände bestimmen kann. An jeden Punkt wird eine Meßlatte mit einem Retroreflektor (Prisma) gestellt. Das Meßgerät wird an einem geeigneten Punkt dazwischen aufgestellt. Die Entfernung wird über die Laufzeit eines im Retroreflektor reflektierten Pulses gemessen. Der Vertikalwinkel durch Anpeilen des Reflektormittelpunktes. Aus diesen Werten für beide Punkte kann die Höhendifferenz ermittelt werden. Genauigkeiten von einigen cm sind erreichbar. Sind die Koordinaten zweier anpeilbarer Punkte bekannt, können über die Messung des Horizontalwinkels die Koordinaten eines dritten Punktes bestimmt werden. Größtes praktisches Problem ist es, einen geeigneten Standpunkt für das Gerät zu finden, von dem aus möglichst viele Punkte gut einsehbar sind.

15. Lotrichtungokular

16. Bestimmung des Koordinatssystems

17. Änderung des Bezugpunkts

18. B D C A

19. S Gemessene Profile W O N

20. Topographie, m Dg, mgal

22. Maare: Bei phreatomagmatischen Eruptionen ensteht als typischer Vulkanbau in der Regel ein Maar. Ein Maar ist eine in den prävulkanischen Untergrund eingesenkte Hohlform, die von einem Ringwall, der aus Auswurfmaterial besteht, umgeben ist. Das Volumen des Ringwalls entspricht nicht selten in etwa dem Volumen der Hohlform und in der Tat wird der Ringwall häufig von Bruchstücken der ehemaligen Maar"füllung" aufgebaut.

23. Wasser w=vertikale Geschwindigkeit der Perkolation des Wasserdampfes a= thermische Duffusionskoffezient h = Viskosität vom Dampf rg=Dichte vom Dampf k=Permeabilität der Gesteine Dampf H DP=(rf-rDamf)•g•(H-z) z w=k•rf•g/h x k= 10-12 m² z= 2·10² m h= 3 10-5 kg/m/s rf= 10³ kg/m³ a=10-3 m²/s T~800°C 2d~ 50 m

24. S Geländeplan Basispunkt W O N

25. Messdaten und Modell des Störkörpers