280 likes | 372 Views
A Numerical Approach to Model the Accretion of Icelandic Crust. Gabriele Marquart and Harro Schmeling. Bathymetry in the North Atlantic. Observations of crustal thickness. Darbyshire,2000. Thickness of the Icelandic crust from Gravity and seismic Data. Extraction. 1 cm/a. 1 cm/a.
E N D
A Numerical Approach to Model the Accretion of Icelandic Crust Gabriele Marquart and Harro Schmeling Colloquium Prague, April, 2005
Bathymetry in the North Atlantic Colloquium Prague, April, 2005
Observations of crustal thickness Colloquium Prague, April, 2005
Darbyshire,2000 Thickness of the Icelandic crust from Gravity and seismic Data Colloquium Prague, April, 2005
Extraction 1 cm/a 1 cm/a Streamlines Melting rate 1. Model Crust is simply related to extracted melt Colloquium Prague, April, 2005
Rising velocity Melting zone melting Numerical Model of a Rising Plume with Melting Anomalous temperature Melt production rate 120 - 60 km depth (T. Ruedas) Colloquium Prague, April, 2005
Texcess= 350 K (1%) Texcess= 250 K ( 0.1%) Texcess= 250 K (1%) Texcess= 250 K (3%) Texcess= 150 K (1%) Predictions for crustal thickness Colloquium Prague, April, 2005
Comparison to „observation“ Model crust Darbyshire Colloquium Prague, April, 2005
2. Model 1 cm/a 1 cm/a Streamlines Melting rate Width of emplacement zone 50 km (Gauß) Extrated material is fed back into the model Colloquium Prague, April, 2005
Kinematic model of Palmason, 1980 Colloquium Prague, April, 2005
Iceland Surface Tectonic Features Colloquium Prague, April, 2005
Structure of the Crust in Iceland Seismic findings: - Distinct upper crust 5-10 km thick - Seismically fast lower crust down to 24-50 km - Poorly constrained transition to the mantle Colloquium Prague, April, 2005
Receiver functions Low Vp-velocities (10%) beneath 40 km Schlindwein, 2001 Crustal Structure from receiver functions Colloquium Prague, April, 2005
The model concept for crustal accretion • Extrusives • fissures, magma chambers • deep dykes and sills • Underplating Colloquium Prague, April, 2005
1 cm/a 1 cm/a Streamlines Melting rate 3. Model Extracted melts are emplaced in a separate crustal model (with contstant rate...) Extraction Colloquium Prague, April, 2005
Source Functions Modeling Crustal Accretion - Equations Physical Equations Momentum conservation: Mass conservation: Energy conservation: Colloquium Prague, April, 2005
Model assumptions • 2D • Constant viscosity • Total accretion rate 2 cm/s spreading rate • T of surface lavas: 100 K • T of magma chambers: 600 K • T deep dykes: 300 K • 3 models: 1) Dominated (60%) by deep accretion 2) Dominated (60%) by magma chamber accretion 3) Dominated (60%) by shallow accretion Colloquium Prague, April, 2005
Visualization of the Accretion of Crust after 500 time steps • Accretion is traced by markers • New markers are inserted at each time step • Color indicates the source • Number of markers is according to the strength of the source • Markers are followed up for 10 Ma, • after 1Ma the color is changed • Marker positions are determined by a RK-4th order scheme Colloquium Prague, April, 2005
Accretion dominated by deep dykes (60% Mtot) Colloquium Prague, April, 2005
Accretion dominated by magma chambers (60% Mtot) Colloquium Prague, April, 2005
Accretion dominated by surface lavas (60% Mtot) Colloquium Prague, April, 2005
Comparison of Different Accretion Styles Lava flows Deep dykes Magma chamber • lateral variable crust • upper crust thinning in central region • hot in central region • vertical layering of the • middle crust • Uniformly stratified • hot crust (Gabbro, • mantle mix?) • thin seismogenic zone • cold crust • hot only in central region • downbuildung, with tilted layering Colloquium Prague, April, 2005
Crustal structure at the rift axis Krustenstruktur aus Seismik Colloquium Prague, April, 2005
Magma Chamber Deep Dykes Surface Lavas Temperature Temperature Temperature 40 C/km 30 C/km 20 C/km Location of profiles Horizontal velocity Horizontal velocity Horizontal velocity Vertical velocity Vertical velocity Vertical velocity Comparison of Different Accretion Types Colloquium Prague, April, 2005
Riftzone 0 km 50 km Comparison to the Seismogenic Crust in Iceland ~ 500°C South Iceland Seismic zone Stefanson, 1998 Deep dykes Magma Chamber Lava flows Depth: 20 km 20 km 10 km 20 km 10 km 5 km 10 km 5 km Colloquium Prague, April, 2005
Magma Chamber Deep Dykes Surface Lavas Temperature Temperature Temperature 40 C/km 30 C/km 20 C/km Location of profiles Horizontal velocity Horizontal velocity Horizontal velocity Vertical velocity Vertical velocity Vertical velocity Comparison of Different Accretion Types • Strong vertical and differential • horizontal velocities Colloquium Prague, April, 2005
Seismic Azimuthal Anisotropy from Rayleigh waves 20-40 km 50-80 km Li & Detrick, EPSL2003 Colloquium Prague, April, 2005
Preliminary Findings for the Accreton of Crust on Iceland • Thermal & geometric structure depends strongly on accretional mode • Iceland: shallow seismogenic zone, high thermal gradient suggests deep or intermediate accretion (deep dykes and magma chambers) as the dominating process (However, the seismogenic upper crust of 10-15 km is produced by shallow fissure swarm intrusions and subairial lava flows) • Then only moderate differential velocities and mixing of the different accretion zones Colloquium Prague, April, 2005