Benchmarking analogue models of brittle thrust wedges

We performed a quantitative comparison of brittle thrust wedge experiments to evaluate the variability among analogue models and to appraise the reproducibility and limits of model interpretation. Fifteen analogue modeling laboratories participated in this benchmark initiative. Each laboratory received a shipment of the same type of quartz and corundum sand and all laboratories adhered to a stringent model building protocol and used the same type of foil to cover base and sidewalls of the sandbox. Sieve structure, sifting height, filling rate, and details on off-scraping of excess sand followed prescribed procedures. Our analogue benchmark shows that even for simple plane-strain experiments with prescribed stringent model construction techniques, quantitative model results show variability, most notably for surface slope, thrust spacing and number of forward and backthrusts. One of the sources of the variability in model results is related to slight variations in how sand is deposited in the sandbox. Small changes in sifting height, sifting rate, and scraping will result in slightly heterogeneous material bulk densities, which will affect the mechanical properties of the sand, and will result in lateral and vertical differences in peak and boundary friction angles, as well as cohesion values once the model is constructed. Initial variations in basal friction are inferred to play the most important role in causing model variability. Our comparison shows that the human factor plays a decisive role, and even when one modeler repeats the same experiment, quantitative model results still show variability. Our observations highlight the limits of up-scaling quantitative analogue model results to nature or for making comparisons with numerical models. The frictional behavior of sand is highly sensitive to small variations in material state or experimental set-up, and hence, it will remain difficult to scale quantitative results such as number of thrusts, thrust spacing, and pop-up width from model to nature.

Lithium isotopic composition of pore fluids and sediments in the Costa Rica subduction zone

Pore fluid and sediment Li concentrations and isotopic ratios provide important insights on the hydrology, sediment contribution to the arc volcanoes and fluid-sediment reactions at the dominantly non-accretionary Costa Rica subduction zone. Ocean Drilling Program Site 1039 in the trench axis provides a reference section of 400 m of the incoming sediments, and Site 1040, situated arcward from the trench, consists of a deformed sedimentary wedge and apron sediments, the décollement, and the partially dewatered underthrust sediment section. At the reference site, pore fluids show important isotopic variations (delta6Li=-21.7 to -37.8 per mil), reflecting the interplay of in situ alteration of volcanic material and ion exchange with clay minerals. In the basal section, a reversal of Li concentration and delta6Li toward seawater values is observed, providing supporting evidence for a lateral seawater flow system in the upper oceanic basement underlying this sediment section. At Site 1040, pore fluid of the lower deformed wedge sediments and within the décollement is enriched in Li and the isotopic compositions are relatively light, suggesting infiltration of a deep-seated fluid. The delta6Li value of -22 per mil of this Li-enriched fluid (261 µM), when compared with the delta6Li value of the subducted sediment section (-11 per mil), suggests that the deep source fluid originates from mineral fluid dehydration and transformation reactions at temperatures of 100 to 150°C, consistent with the temperature range of the up-dip seismogenic zone and of transformation of smectite to illite. The distribution of Li and its isotopes in the underthrust section are similar to those at the reference site, indicating near complete subduction of the incoming sediments and that early dewatering of the underthrust sediments occurs predominantly by lateral flow into the ocean. The hemipelagic clay-rich sediment section of the subducting plate carries most of the Li into this subduction zone, and the pelagic diatomaceous and nannofossil calcareous oozes contain little Li. The Li isotopes of both the clay-rich hemipelagic sediments and of the pelagic oozes are, however, similar, with delta6Li values of -9 to -12 per mil. The observations that (1) the delta6Li values of the underthrust sediments are distinctly lower than that of the mantle, and (2) the lavas of the Costa Rican volcanoes are enriched in Li and 7Li, provide an approximation of the contribution of the subducted sediments to the arc volcanoes. A first order mass balance calculation suggests that approximately half of the Li flux delivered by subducted sediments and altered oceanic crust into the Middle American Trench is recycled to the Costa Rican arc and at most a quarter of sedimentary Li is returned into the ocean through thrust faults, primarily the décollement thrust.

Benthic foraminiferal abundances at the Cretaceous/Paleogene boundary of ODP Hole 117-1049C at Blake Nose, Northwest Atlantic (Appendix A)

Sediments recovered at lower bathyal ODP Site 1049 on Blake Nose (Northwestern Atlantic) offer an opportunity to study environmental changes at the Cretaceous/Paleogene (K/P) boundary relatively close to the Chicxulub impact structure on the Yucatan peninsula, Mexico. In Hole 1049C, the boundary is located at the base of a 9-cm-thick layer with abundant spherules, considered to be impact ejecta. Uppermost Maastrichtian oozes below, and lowermost Danian pelagic oozes above the spherulebed contain well-preserved bathyal benthic foraminifera. The spherule-bed itself, in contrast, contains a mixture of shallow (neritic) and deeper (bathyal) species, and specimens vary strongly in preservation. This assemblage was probably formed by reworking and down-slope transport triggered by the K/P impact. Across the spherule-bed (i.e., the K/P boundary) only ~7% of benthic foraminiferal species became extinct, similar to the low extinction rates of benthic foraminifera worldwide. Quantitative analysis of benthic foraminiferal assemblages and morphogroups in the >63-µm size fraction indicates a relatively eutrophic, stable environment during the latest Maastrichtian, interrupted by a sudden decrease in the food supply to the benthos at the K/P boundary and a decrease in diversity of the faunas, followed by a stepped recovery during the earliest Danian. The recovery was probably linked to the gradual recovery of surface-dwelling primary producers.

(Tables 1, 3) Component composition of organic particles in samples from ODP Site 131-808

Analysis of the palynofacies and miospore thermal alteration indices (TAI) of sediments from ODP Site 808 in the Nankai Trough was undertaken to determine (1) the source, depositional environment, and diagenesis of organic matter in the accreted sediments, and (2) the thermal structure and history of the prism and its relationship to fluid flow. Using the Hartax classification system, two palynofacies were recognized in the sedimentary sequence. Facies 1 occurs within the upper 600 m of trench-wedge turbidites (sedimentation rate > 1 km/m.y.) and contains >50% inertite particles. The rest of the assemblage is dominated by well-preserved phytoclasts and contains small amounts of poorly preserved phytoclasts and well-preserved scleratoclasts. Facies 2 occurs within the Shikoku Basin hemipelagites (600-1300 m below seafloor; sedimentation rate <150 m/m.y.) and contains over two-thirds inertite particles. The rest of the assemblage is dominated by poorly preserved phytoclasts. Miospores and marine phytoplankton compose only a small percentage of both palynofacies. Degraded organic matter is most noticeable in Facies 2, whereas its presence in Facies 1 is overshadowed by the high influx of well-preserved primary organic matter. Most of the degraded organic matter and inertite is interpreted to be reworked. Some of the degraded organic matter may be primary, and may have experienced more biodegradation and thermal alteration in Facies 2 than in Facies 1. TAI values indicate an immature stage of organic maturation (< 2) down to about 900 mbsf. Below this, samples show an increase with depth to a mature stage, reaching peak levels of about 3 just above basement. Samples from within the thrust fault and decollement zones do not show levels of maturity significantly greater than those of surrounding samples, leaving uncertain whether hot fluids have migrated along these fault boundaries in the past.

Estimated chronologic ages of benthic microfossil events in ODP Site 105-647 (Table 1)

Benthic foraminifers were examined from the Paleogene of Ocean Drilling Program (ODP) Site 647 and Deep Sea Drilling Program (DSDP) Site 112 in the southern Labrador Sea. The Paleogene sequence of the deep Labrador Sea can be subdivided into seven assemblages, based on the ranges and relative abundance of characteristic taxa. The first occurrences (FOs) and last occurrences (LOs) of important benthic taxa are calibrated to a standard biochronology, by interpolating from our age model for Site 647. The biostratigraphy of Site 647 is used to improve the age estimates of Site112 cores. Fifteen microfossil events in Site 647 also are found in the sedimentary wedge along the Labrador Margin. A comparison of the probabilistic microfossil sequence from the Labrador Margin with that at Site 647 yields four isochronous benthic foraminifer LOs. Two new species are described from Sites 647 and 112: Hyperammina kenmilleri, Kaminski n.sp., and Ammodiscus nagyi Kaminski n.sp. Significant faunal turnovers are observed at the Ypresian/Lutetian and Eocene/Oligocene boundaries. The Ypresian/Lutetian boundary is characterized by a Glomospira-facies and is attributed to a rise in the CCD (carbonate compensation depth) associated with the NP14 lowstand in sea level. The Eocene/Oligocene boundary is delimited by the LO of Spiroplectammina spectabilis and Reticulophragmium amplectens. The change from an Eocene agglutinated assemblageto a predominantly calcareous assemblage in the early Oligocene took place gradually, over a period of about 4 Ma, but the rate of change accelerated near the boundary. This faunal turnover is attributed to changes in the preservationof agglutinated foraminifers, as delicate species disappeared first. Increasingly poorer preservation of agglutinated foraminifers in the late Eocene to earliest Oligocene reflects the first appearance of cool, nutrient-poor deep water in the southern Labrador Sea. The approximately coeval disappearance of agglutinated assemblages along the Labrador Margin was caused by a regional trend from slope to shelf environments, accentuated by the 'mid'-Oligocene lowstand in sea level.

Radiocarbon ages, dose rates, water isotopic composition and hydrochemistry of samples from Komakuk Beach and Herschel Island

Terrestrial permafrost archives along the Yukon Coastal Plain (northwest Canada) have recorded landscape development and environmental change since the Late Wisconsinan at the interface of unglaciated Beringia (i.e. Komakuk Beach) and the northwestern limit of the Laurentide Ice Sheet (i.e. Herschel Island). The objective of this paper is to compare the late glacial and Holocene landscape development on both sides of the former ice margin based on permafrost sequences and ground ice. Analyses at these sites involved a multi-proxy approach including: sedimentology, cryostratigraphy, palaeoecology of ostracods, stable water isotopes in ground ice, hydrochemistry, and AMS radiocarbon and infrared stimulated luminescence (IRSL) dating. AMS and IRSL age determinations yielded full glacial ages at Komakuk Beach that is the northeastern limit of ice-free Beringia. Herschel Island to the east marks the Late Wisconsinan limit of the northwest Laurentide Ice Sheet and is composed of ice-thrust sediments containing plant detritus as young as 16.2 cal ka BP that might provide a maximum age on ice arrival. Late Wisconsinan ice wedges with sediment-rich fillings on Herschel Island are depleted in heavy oxygen isotopes (mean d18O of -29.1 per mil); this, together with low d-excess values, indicates colder-than-modern winter temperatures and probably reduced snow depths. Grain-size distribution and fossil ostracod assemblages indicate that deglaciation of the Herschel Island ice-thrust moraine was accompanied by alluvial, proluvial, and eolian sedimentation on the adjacent unglaciated Yukon Coastal Plain until ~11 cal ka BP during a period of low glacio-eustatic sea level. The late glacial-Holocene transition was marked by higher-than-modern summer temperatures leading to permafrost degradation that began no later than 11.2 cal ka BP and caused a regional thaw unconformity. Cryostructures and ice wedges were truncated while organic matter was incorporated and soluble ions were leached in the thaw zone. Thermokarst activity led to the formation of ice-wedge casts and deposition of thermokarst lake sediments. These were subsequently covered by rapidly accumulating peat during the early Holocene Thermal Maximum. A rising permafrost table, reduced peat accumulation, and extensive ice-wedge growth resulted from climate cooling starting in the middle Holocene until the late 20th century. The reconstruction of palaeolandscape dynamics on the Yukon Coastal Plain and the eastern Beringian edge contributes to unraveling the linkages between ice sheet, ocean, and permafrost that have existed since the Late Wisconsinan.

Middle Miocene isotope stratigraphy and paleoceanographic evolution of the northwest and southwest Australian margins (Wombat Plateau and Great Australian Bight)

The benthic stable isotope record from ODP Site 761 (Wombat Plateau, NW Australia, 2179.3 m water depth) documents complete recovery of the middle Miocene delta13C excursion corresponding to the climatic optimum and subsequent expansion of the East Antarctic Ice Sheet. The six main delta13C maxima of the "Monterey Excursion" between 16.4 and 13.6 Ma and the characteristic stepped increase in delta18O between 14.5 and 13.9 Ma are clearly identified. The sedimentary record of the shallower ODP Sites 1126 and 1134 [Great Australian Bight (GAB), SWAustralia, 783.8 and 701 m water depth, respectively] is truncated by several unconformities. However, a composite benthic stable isotope curve for these sites provides a first middle Miocene bathyal record for southwest Australia. The delta18O and delta13C curves for Sites 1126 and 1134 indicate a cooler, better-ventilated water mass at ~700 m water depth in the Great Australian Bight since approximately 16 Ma. This cooler and younger water mass probably originated from a close southern source. Cooling of the bottom water at ~16 Ma started much earlier than at other sites of equivalent paleodepths in the central and western parts of the Indian Ocean. At Site 761, the delta18O curve shows an excellent match with the global sea level curve between ~11.5 and 15.1 Ma, and thus closely reflects changes in global ice volume. Prior to 15.1 Ma, the mismatch between the delta18O curve and the sea level curve indicates that delta18O fluctuations are mainly due to changes in bottom water temperature.

Evolution of thermokarst in East Siberian ice-rich permafrost

Thermokarst lakes and basins are major components of ice-rich permafrost landscapes in East Siberian coastal lowlands and are regarded as indicators of regional climatic changes. We investigate the temporal and spatial dynamics of a 7.5 km**2, partly drained thermokarst basin (alas) using field investigations, remote sensing, Geographic Information Systems (GIS), and sediment analyses. The evolution of the thermokarst basin proceeded in two phases. The first phase started at the Pleistocene/Holocene transition (13 to 12 ka BP) with the initiation of a primary thermokarst lake on the Ice Complex surface. The lake expanded and persisted throughout the early Holocene before it drained abruptly about 5.7 ka BP, thereby creating a > 20 m deep alas with residual lakes. The second phase (5.7 ka BP to present) is characterized by alternating stages of lower and higher thermokarst intensity within the alas that were mainly controlled by local hydrological and relief conditions and accompanied by permafrost aggradation and degradation. It included diverse concurrent processes like lake expansion and stepwise drainage, polygonal ice-wedge growth, and the formation of drainage channels and a pingo, which occurred in different parts of the alas. This more dynamic thermokarst evolution resulted in a complex modern thermokarst landscape. However, on the regional scale, the changes during the second evolutionary phase after drainage of the initial thermokarst lakes were less intense than the early Holocene extensive thermokarst development in East Siberian coastal lowlands as a result of a significant regional change to warmer and wetter climate conditions.

Diatom biostratigraphy of sediments from ODP Leg 170 sites

The following data paper summarizes diatom biostratigraphic data from sediments drilled in the Costa Rica accretionary wedge during Ocean Drilling Program Leg 170. Quaternary through lower Miocene diatom zones characteristic of the equatorial Pacific region are recognized in the reference section, Site 1039, which was drilled on the downgoing Cocos plate. At Sites 1040-1043, where the recovered silty clay units are primarily wedge and apron sediments that overlie the underthrust sections, diatoms are generally low in abundance, and complete zonation of the cores was not possible above the décollement surface.

Chemistry and boron isotope ratios of Izu arc front volcanic rocks

The process of fluid release from the subducting slab beneath the Izu arc volcanic front (Izu VF) was examined by measuring B concentrations and B isotope ratios in the Neogene fallout tephra (ODP Site 782A). Both were measured by secondary ion mass spectrometry, in a subset of matrix glasses and glassy plagioclase-hosted melt inclusions selected from material previously analyzed for major and trace elements (glasses) and radiogenic isotopes (Sr, Nd, Pb; bulk tephra). These tephra glasses have high B abundances (~10-60 ppm) and heavy delta11B values (+4.5? to +12.0?), extending the previously reported range for Izu VF rocks (delta11B, +7.0? to +7.3?). The glasses show striking negative correlations of delta11B with large ion lithophile element (LILE)/Nb ratios. These correlations cannot be explained by mixing two separate slab fluids, originating from the subducting sediment and the subducting basaltic crust, respectively (model A). Two alternative models (models B and C) are proposed. Model B proposes that the inverse correlations are inherited from altered oceanic crust (AOC), which shows a systematic decrease of B and LILE with increasing depth (from basaltic layer 2A to layer 3), paralleled by an increase in delta11B (from ~ +1? to +10? to +24?). In this model, the contribution of sedimentary B is insignificant (<4% of B in the Izu VF rocks). Model C explains the correlation as a mixture of a low-delta11B (~ +1?) 'composite' slab fluid (a mixture of metasediment- and metabasalt-derived fluids) with a metasomatized mantle wedge containing elevated B (~1-2 ppm) and heavy delta11B (~ +14?). The mantle wedge was likely metasomatized by 11B-rich fluids beneath the outer forearc, and subsequently down dragged to arc front depths by the descending slab. Pb-B isotope systematics indicate that, at arc front depths, ~ 53% of the B in the Izu VF is derived from the wedge. This implies that the heavy delta11B values of Izu VF rocks are largely a result of fluid fractionation, and do not reflect variations in slab source provenance (i.e. subducting sediment vs. basaltic crust). Since the B content of the peridotite at the outer forearc (7-58 ppm B, mean 24 +/- 16 ppm) is much higher than beneath the arc front (~1-2 ppm B), the hydrated mantle wedge must have released a B-rich fluid on its downward path. This 'wedge flux' can explain (1) the across-arc decrease in B and delta11B (e.g. Izu, Kuriles), without requiring a progressive decrease in fluid flux from the subducting slab, and (2) the thermal structure of volcanic arcs, as reflected in the B and delta11B variations of volcanic arc rocks.

Grain size composition of sediment cores from the Weddell Sea, Antarctica

Sedimentary processes in the southeastern Weddell Sea are influenced by glacial-interglacial ice-shelf dynamics and the cyclonic circulation of the Weddell Gyre, which affects all water masses down to the sea floor. Significantly increased sedimentation rates occur during glacial stages, when ice sheets advance to the shelf edge and trigger gravitational sediment transport to the deep sea. Downslope transport on the Crary Fan and off Dronning Maud and Coats Land is channelized into three huge channel systems, which originate on the eastern-, the central and the western Crary Fan. They gradually turn from a northerly direction eastward until they follow a course parallel to the continental slope. All channels show strongly asymmetric cross sections with well-developed levees on their northwestern sides, forming wedge-shaped sediment bodies. They level off very gently. Levees on the southeastern sides are small, if present at all. This characteristic morphology likely results from the process of combined turbidite-contourite deposition. Strong thermohaline currents of the Weddell Gyre entrain particles from turbidity-current suspensions, which flow down the channels, and carry them westward out of the channel where they settle on a surface gently dipping away from the channel. These sediments are intercalated with overbank deposits of high-energy and high-volume turbidity currents, which preferentially flood the left of the channels (looking downchannel) as a result of Coriolis force. In the distal setting of the easternmost channel-levee complex, where thermohaline currents are directed northeastward as a result of a recirculation of water masses from the Enderby Basin, the setting and the internal structures of a wedge-shaped sediment body indicate a contourite drift rather than a channel levee. Dating of the sediments reveals that the levees in their present form started to develop with a late Miocene cooling event, which caused an expansion of the East Antarctic Ice Sheet and an invigoration of thermohaline current activity.

Clay mineralogy of ODP Site 131-808 sediments

Drill core recovered at Ocean Drilling Program Site 808 (Leg 131) proves that the wedge of trench sediment within the central region of the Nankai Trough comprises approximately 600 m of hemipelagic mud, sandy turbidites, and silty turbidites. The stratigraphic succession thickens and coarsens upward, with hemipelagic muds and volcanic-ash layers of the Shikoku Basin overlain by silty and sandy trench-wedge deposits. Past investigations of clay mineralogy and sand petrography within this region have led to the hypothesis that most of the detritus in the Nankai Trough was derived from the Izu-Honshu collision zone and transported southwestward via axial turbidity currents. Shipboard analyses of paleocurrent indicators, on the other hand, show that most of the ripple cross-laminae within silty turbidites of the outer marginal trench-wedge facies are inclined to the north and northwest; thus, many of the turbidity currents reflected off the seaward slope of the trench rather than moving straight down the trench axis. Shore-based analyses of detrital clay minerals demonstrate that the hemipelagic muds and matrix materials within sandy and silty turbidites are all enriched in illite; chlorite is the second-most abundant clay mineral, followed by smectite. In general, the relative mineral percentages change relatively little as a function of depth, and the hemipelagic clay-mineral population is virtually identical to the turbidite-matrix population. Comparisons between different size fractions (<2 µm and 2-6 µm) show modest amounts of mineral partitioning, with chlorite content increasing in the coarser fraction and smectite increasing in the finer fraction. Values of illite crystallinity index are consistent with conditions of advanced anchimetamorphism and epimetamorphism within the source region. Of the three mica polytypes detected, the 2M1 variety dominates over the 1M and 1Md polytypes; these data are consistent with values of illite crystallinity. Measurements of mica bo lattice spacing show that the detrital illite particles were eroded from a zone of intermediate-pressure metamorphism. Collectively, these data provide an excellent match with the lithologic and metamorphic character of the Izu-Honshu collision zone. Data from Leg 131, therefore, confirm the earlier interpretations of detrital provenance. The regional pattern of sediment dispersal is dominated by a combination of southwest-directed axial turbidity currents, radial expansion of the axial flows, oblique movement of suspended clouds onto and beyond the seaward slope of the Nankai Trough, and flow reflection back toward the trench axis.

Compilation of Mg/Ca ratios of Neogloboquadrina pachyderma (sinistral) and SST estimates from the northern North Atlantic

A map of estimated calcification temperatures of the planktic foraminifer Neogloboquadrina pachyderma sinistral (TNps) for the Nordic Seas and the northern North Atlantic for the Last Glacial Maximum was produced from oxygen isotopes with support of Mg/Ca ratios. To arrive at the reconstruction, several constraints concerning the plausible salinity and ?18O-fields were employed. The reconstruction indicates inflow of temperate waters in a wedge along the eastern border of the Nordic Seas and at least seasonally ice-free waters. The reconstruction from oxygen isotopes shows similarities with Mg/Ca based paleotemperatures in the southern and southeastern sector, while unrealistically high Mg/Ca values in the central Nordic Seas prevent the application of the method in this area. The oxygen isotope based reconstruction shows some agreement with temperature reconstructions based on the modern analogue technique, but with somewhat lower temperatures and a stronger internal gradient inside the Nordic Seas. All told, our results suggest a much more ice-free and dynamic high latitude ocean than the CLIMAP reconstruction.