A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Growth response of an early successional assemblage of coralline algae and benthic diatoms to ocean acidification

The sustained absorption of anthropogenically released atmospheric CO2 by the oceans is modifying seawater carbonate chemistry, a process termed ocean acidification (OA). By the year 2100, the worst case scenario is a decline in the average oceanic surface seawater pH by 0.3 units to 7.75. The changing seawater carbonate chemistry is predicted to negatively affect many marine species, particularly calcifying organisms such as coralline algae, while species such as diatoms and fleshy seaweed are predicted to be little affected or may even benefit from OA. It has been hypothesized in previous work that the direct negative effects imposed on coralline algae, and the direct positive effects on fleshy seaweeds and diatoms under a future high CO2 ocean could result in a reduced ability of corallines to compete with diatoms and fleshy seaweed for space in the future. In a 6-week laboratory experiment, we examined the effect of pH 7.60 (pH predicted to occur due to ocean acidification just beyond the year 2100) compared to pH 8.05 (present day) on the lateral growth rates of an early successional, cold-temperate species assemblage dominated by crustose coralline algae and benthic diatoms. Crustose coralline algae and benthic diatoms maintained positive growth rates in both pH treatments. The growth rates of coralline algae were three times lower at pH 7.60, and a non-significant decline in diatom growth meant that proportions of the two functional groups remained similar over the course of the experiment. Our results do not support our hypothesis that benthic diatoms will outcompete crustose coralline algae under future pH conditions. However, while crustose coralline algae were able to maintain their presence in this benthic rocky reef species assemblage, the reduced growth rates suggest that they will be less capable of recolonizing after disturbance events, which could result in reduced coralline cover under OA conditions.