Modelling of Continental Lithosphere Breakup and Rifted Margin Formation in Response to an Upwelling Divergent Flow Field Incorporating a Temperature Dependent Rheology
Abstract
We numerically model continental lithosphere deformation leading to breakup and sea floor spreading initiation in response to an imposed upwelling and divergent flow field applied to continental lithosphere and asthenosphere. The model is used to predict rifted continental margin lithosphere thinning and temperature structure. Model predictions are compared with observed rifted margin structure for four diverse case studies. Prior to application of the upwelling divergent flow field the continental lithosphere is undeformed with a uniform temperature gradient. The upwelling divergent flow field is defined kinematically using boundary conditions consisting of the upwelling velocity Vz at the divergence axis and the half divergence rate Vx . The resultant velocity field throughout the continuum is computed using finite element (FE) code incorporating a Newtonian temperature dependent rheology. The flow field is used to advect the continental lithosphere material and lithospheric and asthenospheric temperatures. Viscosity structure is hence modified and the velocities change correspondingly in a feedback loop. We find the kinematic boundary conditions Vz and Vx to be of first order importance. A high Vz/Vx (greater than10), corresponding to buoyancy assisted flow, leads to minimal mantle exhumation and a well defined continent ocean transition consistent with observations at volcanic margins. For Vz/Vx near unity, corresponding to plate boundary driven divergence, mantle exhumation over widths of up to 100 km is predicted which is consistent with observations at non-volcanic margins. The FE method allows the upwelling velocity Vz to be propagated upwards from the top of the asthenosphere to the Earth's surface without the requirement of imposing Vx. When continental breakup is achieved the half divergence velocity Vx can be applied at the lithosphere surface and the upwelling velocity Vz left free. We find this time and space dependent set of boundary conditions is more plausible than a constant corner flow type solution and predicts levels of depth dependent stretching and continent ocean transitions consistent with observation. Depth dependent lithosphere stretching, which is observed at rifted continental margins, is predicted to occur before continental breakup and sea-floor spreading initiation. The model may be used to predict surface heat flow and bathymetry, and to provide estimates of melt production rates and cumulative thickness. We compare model predictions with observed margin structure for four diverse rifted margins: the Lofoten Margin (a mature volcanic margin), Goban Spur (a mature non-volcanic margin), the Woodlark Basin (a neotectonic young ocean basin) and the Faroe-Shetland Basin (a failed attempt at continental breakup). This work forms part of the NERC Margins iSIMM project. iSIMM investigators are from Liverpool and Cambridge Universities, Badley Geoscience & Schlumberger Cambridge Research supported by the NERC, the DTI, Agip UK, BP, Amerada Hess Ltd, Anadarko, Conoco¬Phillips, Shell, Statoil and WesternGeco. The iSIMM team comprises NJ Kusznir, RS White, AM Roberts, PAF Christie, A Chappell, J Eccles, R Fletcher, D Healy, N Hurst, ZC Lunnon, CJ Parkin, AW Roberts, LK Smith, V Tymms & R Spitzer.
- Publication:
-
AGU Spring Meeting Abstracts
- Pub Date:
- May 2005
- Bibcode:
- 2005AGUSM.T53A..10T
- Keywords:
-
- 8105 Continental margins and sedimentary basins;
- 8120 Dynamics of lithosphere and mantle: general