Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

A three‐dimensional numerical model of sediment transport,erosion and deposition within a network of channel belts,floodplain and hill slope: extrinsic and intrinsic controls on floodplain dynamics and alluvial architecture

A three‐dimensional numerical model of sediment transport, erosion and deposition within a network of channel belts and associated floodplain is described. Sediment and water supply are defined at the upstream entry point, and base level is defined at the downstream edge of the model. Sediment and water are transported through a network of channels according to the diffusion equation, and each channel has a channel belt with a width that increases in time. The network of channels evolves as a result of channel bifurcation and abandonment (avulsion). The timing and location of channel bifurcation is controlled stochastically as a function of the cross‐valley slope of the floodplain adjacent to the channel belt relative to the down‐valley slope, and of annual flood discharge. A bifurcation develops into an avulsion when the discharge of one of the distributaries falls below a threshold value. The floodplain aggradation rate decreases with distance from the nearest active channel belt. Channel‐belt degradation results in floodplain incision. Extrinsic (extrabasinal, allogenic) and intrinsic (intrabasinal, autogenic) controls on floodplain dynamics and alluvial architecture were modelled, and sequence stratigraphy models were assessed. Input parameters were chosen based on data from the Rhine–Meuse delta. To examine how the model responds to extrinsic controls, the model was run under conditions of changing base level and increasing sediment supply. Rises and falls in base level and increases in sediment supply occurred over 10 000 years. Rising base level caused a wave of aggradation to move up‐valley, until aggradation occurred over the entire valley. Frequency of bifurcations and avulsions increased with rate of base‐level rise and aggradation rate. Channel‐belt width varied with water discharge and the lifespan of the channel belt. Wide, connected channel belts (and high channel‐deposit proportion) occurred around the upstream inflow point because of their high discharge and longevity. Less connected, smaller channel belts occurred further down‐valley. Such alluvial behaviour and architecture is also found in the Rhine–Meuse delta. During base‐level fall, valley erosion occurred, and the incised valley contained a single wide channel belt. During subsequent base‐level rise, a wave of aggradation moved up‐valley, filling the incised valley. Bifurcation and avulsion sites progressively moved upstream. Relatively thin, narrow channel belts bordered and cut into the valley fill. These results differ substantially from existing sequence stratigraphy models. The increase in sediment supply from upstream resulted in an alluvial fan. Most bifurcations and avulsions occurred at the fan apex (nodal avulsion), and channel belts were the widest and the thickest here (giving high channel‐deposit proportion) due to their high discharge and longevity. The width and thickness of channel belts decreased down‐valley due to decreased discharge, longevity and aggradation rate. This behaviour occurs in modern alluvial fans. Intrinsic controls also affect floodplain dynamics and alluvial architecture. Variation of aggradation rate, bifurcation frequency and number of coexisting channel belts occurred over periods of 500 to 2000 years, compared with 10 000 years for extrinsic controls. This variation is partly related to local aggradation and degradation of channel belts around bifurcation points. Channel belts were preferentially clustered near floodplain margins, because of low floodplain aggradation rate and topography there.     

大陆板片多重性俯冲与折返的动力学模式——苏鲁高压-超高压变质地体的折返年龄限定

苏鲁高压超高压变质地体自南而北由高压(HP)、很高压(VHP)和超高压(UHP)变质叠覆岩片组成,前者依次叠覆在后者之上,岩片之间的界限为韧性剪切带。根据超高压变质岩片中角闪岩相岩石与高压变质岩片中绿片岩相岩石的黑云母和白云母Ar_Ar和Rb_Sr测年新结果,结合前人在该区所做的锆石SHRI MP U_Pb、全岩Sm_Nd、Rb_Sr等测年数据综合分析表明,超高压变质岩石的峰期变质年龄为240~220Ma,折返年龄为220~200Ma;而高压变质岩石的峰期变质年龄大于258Ma,起始折返年龄为258~240Ma,折返年龄比超高压变质岩石早30~40Ma。这说明扬子板片并不是整体俯冲和折返的。由于具组分和密度差异,俯冲板块的不同部位沿岩性或构造界面先后分片俯冲和折返,在北苏鲁超高压变质板片开始俯冲时,南苏鲁高压变质板片已开始折返。     

Evaluation of a numerical model of the British–Irish ice sheet using relative sea‐level data: implications for the interpretation of trimline observations

The glacial isostatic adjustment (GIA) of the British Isles is complex due to the interplay between local and non‐local signals. A number of recent studies have modelled the GIA response of the British Isles using relative sea‐level data. This study extends these previous analyses by using output from a numerical glaciological model as input to a GIA model. This is a necessary step towards more realistic GIA models, and although there have been similar studies for the major late Pleistocene ice sheets, this is the first study to do so for the British Isles. We test three reconstructions, classed as ‘minimal’, ‘median’ and ‘maximal’ in terms of their volume at maximum extent, and find it is possible to obtain good data–model fits. The minimal reconstruction is clearly preferred by the sea‐level data. The ice reconstructions tested were not constrained by geomorphological information of past ice extent (lateral and vertical). As a consequence, the reconstructions extend further than much of this information suggests, particularly in terms of ice thickness. It is notable, however, that the reconstructions produce good fits to many sea‐level data from central, mountainous regions (e.g. Scottish highlands), which lends support to the suggestion that trimlines, often used as an constraint on the palaeo ice surface location, are in fact features formed at the transition from warm‐ to cold‐based ice and so mark a minimum constraint on the ice surface altitude. Based on data–model misfits, suggestions are made for improving the ice model reconstructions. However, in many locations, the cause of the misfit could be due to non‐local, predominantly Fennoscandian ice and so interpretation is not straightforward. As a result, we suggest that future analyses of this type consider models and observations for both Fennoscandia and the British Isles. Copyright © 2012 John Wiley & Sons, Ltd.