Nina Kukowski

The impact of analogue material properties on the geometry, kinematics, and dynamics of convergent sand wedges

Abstract Simulation of geodynamic processes in sandbox experiments requires analogue materials with a deformation behaviour that reproduces the deformation mechanisms of typical crustal rocks. We present data on the frictional strength of different sand types employing static and dynamic shear tests. The sand types analysed are characterised by an elastic/frictional plastic mechanical behaviour with a transient strain-hardening and strain-softening phase prior to transition to stable sliding. This is in conflict with the standard assumption of an ideal cohesionless Coulomb-material with constant frictional properties. The influence of the identified transient material properties on the kinematics, growth mechanisms, and internal deformation patterns of convergent sand wedges results in characteristic wedge segments which vary—depending on material compaction—between wedges with well defined segments (i.e. frontal-deformation zone, frontal-imbrication zone and internal-accumulation zone) with straight slopes and wedges with a continuous convex topographic profile. For most materials, only the frontal part of the wedge is critical during experimental runs. Taper and strength of the wedge segments can be shown to be controlled by the frictional properties of active faults. Wedge segmentation is controlled by a bulk-strength increase toward the rear of the wedge due to fault rotation in mechanically less-favourable orientations and plastic material hardening. The limit between the frontal critical parts of a wedge and internal stable parts is largely controlled by a critical state of stress upon which either renewed failure or fault inactivation occurs. On this basis, we suggest that critical-taper analysis of wedges must be restricted to specific kinematic segments. Comparison of the experimental results with the Nankai accretionary wedge suggests that our interpretation also applies to natural convergent wedges. Moreover, we provide constraints for the selection of adequate granular analogue materials to simulate typical crustal rocks in natural convergent wedges.

Backstop geometry and accretionary mechanics of the Sunda margin

The convergent Sunda margin off Indonesia displays all geological features characteristic of an accretion-dominated subduction zone. A combined interpretation of prestack depth-migrated seismic reflection data and velocity information gained from refraction studies is supplemented by high-resolution bathymetric data and for the first time allows the exact mapping of backstop regimes. Initially, the outer high evolved as material was pushed against a static rigid arc framework backstop underlying a forearc basin. Increasing material strength of the outer high due to lithification formed a dynamic backstop, which controls accretion today. An out-of-sequence thrust marks the transition from the recent active frontal accretionary prism to the outer high and may be traced in the seismic and bathymetric data over the whole extent of the study area. The existence of a static as well as a dynamic backstop controls the forearc geometry and is associated with the segmentation of the forearc, which is observed in regimes of frontal as well as of oblique subduction. Mass balance calculations, which account for porosity changes and metamorphism, indicate a subduction history dominated by accretionary processes since the late Eocene. Accretion is associated with the low values of basal friction inferred for the Sunda margin. Structural investigations of conjugate fault planes indicate a very weak basal detachment. Effective stress analyses reveal that intrinsically weak material causes the high strength ratio of the detachment to the overlying sediments, whereas overpressuring within the frontal accretionary prism is negligible.

New seismic images of the Cascadia subduction zone from cruise SO108 — ORWELL

Abstract In April and May 1996, a geophysical study of the Cascadia continental margin off Oregon and Washington was conducted aboard the German R/V Sonne. This cooperative experiment by GEOMAR and the USGS acquired wide-angle reflection and refraction seismic data, using ocean-bottom seismometers (OBS) and hydrophones (OBH), and multichannel seismic reflection (MCS) data. The main goal of this experiment was to investigate the internal structure and associated earthquake hazard of the Cascadia subduction zone and to image the downgoing plate. Coincident MCS and wide-angle profiles along two tracks are presented here. The plate boundary has been imaged precisely beneath the wide accretionary wedge close to shore at ca. 13 km depth. Thus, the downgoing plate dips more shallowly than previously assumed. The dip of the plate changes from 2° to 4° at the eastern boundary of the wedge on the northern profile, where approximately 3 km of sediment is entering the subduction zone. On the southern profile, where the incoming sedimentary section is about 2.2 km thick, the plate dips about 0.5° to 1.5° near the deformation front and increases to 3.5° further landwards. On both profiles, the deformation of the accretionary wedge has produced six ridges on the seafloor, three of which represent active faulting, as indicated by growth folding. The ridges are bordered by landward verging faults which reach as deep as the top of the oceanic basement. Thus the entire incoming sediment package is being accreted. At least two phases of accretion are evident, and the rocks of the older accretionary phase(s) forms the backstop for the younger phase, which started around 1.5 Ma ago. This documents that the 30 to 50 km wide frontal part of the accretionary wedge, which is characterized by landward vergent thrusts, is a Pleistocene feature which was formed in response to the high input of sediment building the fans during glacial periods. Velocities increase quite rapidly within the wedge, both landward and downward. At the toe of the deformation front, velocities are higher than 4.0 km/s, indicating extensive dewatering of deep, oceanic sediment. Further landward, considerable velocity variation is found, which indicates major breaks throughout the accretionary history.

Analogue benchmarks of shortening and extension experiments

Abstract We report a direct comparison of scaled analogue experiments to test the reproducibility of model results among ten different experimental modelling laboratories. We present results for two experiments: a brittle thrust wedge experiment and a brittleviscous extension experiment. The experimental set-up, the model construction technique, the viscous material and the base and wall properties were prescribed. However, each laboratory used its own frictional analogue material and experimental apparatus. Comparison of results for the shortening experiment highlights large differences in model evolution that may have resulted from (1) differences in boundary conditions (indenter or basal-pull models), (2) differences in model widths, (3) location of observation (for example, sidewall versus centre of model), (4) material properties, (5) base and sidewall frictional properties, and (6) differences in set-up technique of individual experimenters. Six laboratories carried out the shortening experiment with a mobile wall. The overall evolution of their models is broadly similar, with the development of a thrust wedge characterized by forward thrust propagation and by back thrusting. However, significant variations are observed in spacing between thrusts, their dip angles, number of forward thrusts and back thrusts, and surface slopes. The structural evolution of the brittle-viscious extension experiments is similar to a high degree. Faulting initiates in the brittle layers above the viscous layer in close vicinity to the basal velocity discontinuity. Measurements of fault dip angles and fault spacing vary among laboratories. Comparison of experimental results indicates an encouraging overall agreement in model evolution, but also highlights important variations in the geometry and evolution of the resulting structures that may be induced by differences in modelling materials, model dimensions, experimental set-ups and observation location.

Ridge subduction at an erosive margin: The collision zone of the Nazca Ridge in southern Peru

The 1.5-km-high, obliquely subducting Nazca Ridge and its collision zone with the Peruvian margin have been imaged by wide-angle and reflection seismic profiles, swath bathymetry, and gravity surveying. These data reveal that the crust of the ridge at its northeastern tip is 17 km thick and exhibits seismic velocities and densities similar to layers 2 and 3 of typical oceanic crust. The lowermost layer contributes 10–12 km to the total crustal thickness of the ridge. The sedimentary cover is 300–400 m thick on most parts of the ridge but less than 100 m thick on seamounts and small volcanic ridges. At the collision zone of ridge and margin, the following observations indicate intense tectonic erosion related to the passage of the ridge. The thin sediment layer on the ridge is completely subducted. The lower continental slope is steep, dipping at ∼9°, and the continental wedge has a high taper of 18°. Tentative correlation of model layers with stratigraphy derived from Ocean Drilling Program Leg 112 cores suggests the presence of Eocene shelf deposits near the trench. Continental basement is located <15 km landward of the trench. Normal faults on the upper slope and shelf indicate extension. A comparison with the Peruvian and northern Chilean forearc systems, currently not affected by ridge subduction, suggests that the passage of the Nazca Ridge along the continental margin induces a temporarily limited phase of enhanced tectonic erosion superposed on a long-term erosive regime.

Upward delamination of Cascadia Basin sediment infill with landward frontal accretion thrusting caused by rapid glacial age material flux

The Cascadia convergent margin is a first-order research target to study the impact of rapid sedimentation processes on the mechanics of frontal subduction zone accretion. The near-trench part of the accretionary prism offshore Washington is affected by strongly increased glacial age sedimentation and fan formation that led to an outstanding Quaternary growth rate with landward vergent thrust faulting that is rarely observed elsewhere in accretionary wedges. Multichannel seismic reflection data acquired on the ORWELL project allows us to study the structure and dynamics of the atypical frontal accretion processes. We performed a kinematical and mechanical analysis of the frontal accretion structures, and developed a dynamic Coulomb-wedge model for the landward-verging backthrust formation. Backthrusting results from heterogeneous diffuse strain accumulation in the mechanically heterogeneous Cascadia basin sediment succession entering the subduction zone, and strain partitioning along a midlevel detachment that is activated by gravitational loading caused by rapid glacial age sedimentation. These complex deformation processes cause the passive “upward” delamination of the upper turbidite beds from the basal pelagic carbonate section similar to triangle-zone formation and passive backthrust wedging in foreland thrust belts caused by rapid burial beneath syntectonic sediment deposits. The deformation mechanism at the tectonic front of the Cascadia margin is an immediate response to the strongly increased late Pleistocene sediment flux rather than to atypical physical boundary conditions as generally thought.

Sediment accretion against a buttress beneath the Peruvian continental margin at 12 ° S as simulated with sandbox modeling

Reflection seismic data from the Peruvian continental margin at 12° S clearly reveal an accretionary wedge and buttress. Sandbox experiments applying the physical concept of the Coulomb theory allow the systematic investigation of the growth and deformation of such an accretionary structure. The style of deformation of the buttress and the internal structure of the wedge is observed in the sandbox models. The possibility of underplating material beneath the buttress and the amount of tectonic erosion depend on the physical properties of the materials, mainly internal friction, cohesion and basal friction. Boundary conditions such as the height of the subduction gate and the thickness of incoming sand also constrain the style of growth of the model accretionary structure.The configurations of two experiments were closely scaled to reflection seismic depth sections across the Peruvian margin. A deformable buttress constructed of compacted rock powder is introduced to replicate the basement rock which allows deformation similar to that in the seismic data. With the sandbox models it is possible to verify a proposed accretionary history derived from seismic and borehole data. The models also help in understanding the mechanisms which control the amount of accretion, subduction and underplating as a function of physical properties, boundary conditions and the duration of convergence.

The link between bottom-simulating reflections and methane flux into the gas hydrate stability zone-new evidence from Lima Basin, Peru Margin

Bottom-simulating reflections (BSRs) are probably the most commonly used indicators for gas hydrates in marine sediments. It is now widely accepted that BSRs are primarily caused by free gas beneath gas-hydrate-bearing sediments. However, our insight into BSR formation to date is mostly limited to theoretical studies. Two endmember processes have been suggested to supply free gas for BSR formation: (i) dissociation of gas hydrates and (ii) migration of methane from below. During a recent campaign of the German Research Vessel Sonne off the shore of Peru, we detected BSRs at locations undergoing both tectonic subsidence and non-sedimentation or seafloor erosion. Tectonic subsidence (and additionally perhaps seafloor erosion) causes the base of gas hydrate stability to migrate downward with respect to gas-hydrate-bearing sediments. This process rules out dissociation of gas hydrates as a source of free gas for BSRs at these locations. Instead, free gas at BSRs is predicted to be absorbed into the gas hydrate stability zone. BSRs appear to be confined to locations where the subsurface structure suggests focusing of fluid flow. We investigated the seafloor at one of these locations with a TV sled and observed fields of rounded boulders and slab-like rocks, which we interpreted as authigenic carbonates. Authigenic carbonates are precipitations typically found at cold vents with methane expulsion. We retrieved a small carbonate-cemented sediment sample from the seafloor above a BSR about 20 km away. This supported our interpretation that the observed slabs and boulders were carbonates. All these observations suggest that BSRs in Lima Basin are maintained predominantly by gas that is supplied from below, demonstrating that this endmember process for BSR formation exists in nature. Results from Ocean Drilling Program Leg 112 showed that methane for gas hydrate formation on the Peru lower slope and the methane in hydrocarbon gases on the upper slope is mostly of biogenic origin. The δ13C composition of the recovered carbonate cement was consistent with biologic methane production below the seafloor (although possibly above the BSR). We speculate that the gas for BSR formation in Lima Basin also is mainly biogenic methane. This would suggest the biologic productivity beneath the gas hydrate zone in Lima Basin to be relatively high in order to supply enough methane to maintain BSRs.

Oblique Convergence along the Chilean Margin: Partitioning, Margin-Parallel Faulting and Force Interaction at the Plate Interface

The Chilean fore-arc exhibits margin-parallel strike-slip faulting and associated fore-arc sliver formation. Comparison of the long-term (geologic timescale) and short-term (human timescale) record of margin-parallel faulting along the oblique Chilean subduction margin between 15° S and 46° S reveals significant spatio-temporal heterogeneity. We have reviewed newly compiled data on the geometric, kinematic and mechanical properties and their variation along-strike of the Chilean margin and evaluated their competing influence on fore-arc deformation. Among the parameters considered are the plate kinematics (e.g., convergence obliquity and rate), overriding plate heterogeneities that affect its capability for localizing horizontal shear (e.g., thermal and structural weaknesses) or resistance to block motion (e.g., plate margin curvature) as well as properties governing or indicating force interaction at the plate interface (e.g., trench sediment-fill, geodetic and seismic coupling depth). Most remarkably, the short-term GPS-derived fore-arc velocity field, dominated by elastic loading processes, shows little variation along-strike of the margin, despite the significantly changing conditions (e.g., trench sediment-fill, mass transfer mode at the tip of the overriding plate, plateau or no plateau, and slab-dip variations). Variations in recent and past strike-slip motion do not appear to depend on the rate or obliquity of convergence, nor on the mode of mass transfer at the subduction front. The frictionally coupled area on the plate contact increases southwards and a decreasing taper along the Chilean margin can be reconciled in the framework of taper theory by a southward decrease of the effective coefficient of friction on the plate interface. The development of fore-arc slivers seems to be primarily controlled by mechanisms that cause effective rheological weakening of parts of the upper plate and/or by geometries that hamper margin-parallel sliver motion. While the seaward concave-shaped margin in North Chile hinders the margin-parallel motion of a fore-arc sliver, the present strike-slip activity of the Liquine-Ofqui Fault Zone in southern Chile is likely facilitated by the superposition of two conditions: a shallowly dipping slab and an exceptionally small arc to trench distance.

Thermo-hydraulics of the Peruvian accretionarycomplex at 12°S

Abstract Coupled heat and fluid transport at the Peruvian convergent margin at 12°S wasstudied with finite element modelling. Structural information was available from two seismicreflection lines. Heat production in the oceanic plate, the metamorphic basement, and sedimentswas estimated from literature. Porosity, permeability, and thermal conductivity for the modelswere partly available from Ocean Drilling Program (ODP) Leg 112; otherwise we used empiricalrelations. Our models accounted for a possible permeability anisotropy. The decollement was bestmodelled as a highly permeable zone (10 −13 m 2 ). Permeabilities of thePeruvian accretionary wedge adopted from the model calculations fall within the range of 2 to7×10 −16 m 2 at the ocean bottom to a few 10 −18 m 2 at the base and need to be anisotropic. Fluid expulsion at the sea floor decreases graduallywith distance from the deformation front and is structure controlled. Small scale variations of heatflux reflected by fluctuations of BSR depths across major faults could be modelled assuming highpermeability in the faults which allow for efficient advective transport along those faults. The models were constrained by the thermal gradient obtained from the depth of bottomsimulating reflectors (BSRs) at the lower slope and some conventional measurements. We foundthat significant frictional heating is required to explain the observed strong landward increase ofheat flux. This is consistent with results from sandbox modelling which predict strong basalfriction at this margin. A significantly higher heat source is needed to match the observed thermalgradient in the southern line.