Strontium isotope ratios of planktonic foraminifera from sediment samples acroos the K/T boundary

Seawater 87Sr/86Sr values increase abruptly by 28 * 10**-6 across the Cretaceous/Tertiary boundary (KTB). This small, but rapid shift is superimposed on the larger scale structure of the seawater Sr isotope curve. The time scale of radiogenic Sr addition appears to be too rapid to reconcile with sources associated with volcanism, and we show that the amount of Sr required to produce even this small increase is too large to be derived from: (1) a KT bolide of the size constrained by the Ir anomaly, (2) continental crust ejecta from the impact of such a bolide, (3) soot from global wildfires initiated by an impact, or (4) any combination of these sources. The probable source of the radiogenic Sr is enhanced continental weathering, but the high rate of increase appears to rule out processes such as sea level regression, glaciation or tectonism. A plausible mechanism for rapid addition of radiogenic Sr to the oceans is enhanced weathering associated with globally distributed acid rain (pH c. 1) which is a proposed by-product of a bolide impact (Prinn and Fegley, 1987, doi:10.1016/0012-821X(87)90046-X).

Landfast sea-ice and platelet-layer thickness and conductivity of Atka Bay, Antarctica, December 2012

Ice shelves strongly impact coastal Antarctic sea-ice and the associated ecosystem through the formation of a sub-sea-ice platelet layer. Although progress has been made in determining and understanding its spatio-temporal variability based on point measurements, an investigation of this phenomenon on a larger scale remains a challenge due to logistical constraints and a lack of suitable methodology. In this study, we applied a laterally-constrained Marquardt-Levenberg inversion to a unique multi-frequency electromagnetic (EM) induction sounding dataset obtained on the landfast sea ice of Atka Bay, eastern Weddell Sea, in 2012. In addition to consistent fast-ice thickness and -conductivities along > 100 km transects; we present the first comprehensive, high resolution platelet-layer thickness and -conductivity dataset recorded on Antarctic sea ice. The reliability of the algorithm was confirmed by using synthetic data, and the inverted platelet-layer thicknesses agreed within the data uncertainty to drill-hole measurements. Ice-volume fractions were calculated from platelet-layer conductivities, revealing that an older and thicker platelet layer is denser and more compacted than a loosely attached, young platelet layer. The overall platelet-layer volume below Atka Bay fast ice suggests that the contribution of ocean/ice-shelf interaction to sea-ice volume in this region is even higher than previously thought. This study also implies that multi-frequency EM induction sounding is an effective approach in determining platelet layer volume on a larger scale than previously feasible. When applied to airborne multi-frequency EM, this method could provide a step towards an Antarctic-wide quantification of ocean/ice-shelf interaction.

Magnetic susceptibility of deep-sea cores from the Scotia Sea

We studied two deep-sea cores from the Scotia Sea to reconstruct past atmospheric circulation in the southern hemisphere and to resolve a long-standing debate on the interpretation of magnetic susceptibility (MS) records in Southern Ocean (SO) sediment. High-sedimentation sites MD07-3134 (0.2 - 1.2 m/kyr) and MD07-3133 (0.3 - 2 m/kyr) cover the last 92.5 kyr and 36 kyr, respectively. Both exhibit a one-to-one coupling of the MS and Ca2+ signal to the non-sea salt (nss) Ca2+ signal of the EDML ice core, clearly identifying atmospheric circulation as means of distribution. Comparison of additional proxies also excludes major influence by volcanic sources, sea-ice, icebergs, or oceanic current transport. The close resemblance of the dust proxies over the last glacial cycle, in turn, allows for the establishment of an age model of unprecedented resolution and precision for SO deep-sea sediment because atmospheric transport involves no major leads or lags. This is of particular importance because MS is routinely measured on deep-sea cores in the SO but the sediments usually lack biogenic carbonate and therefore had only limited stratigraphic control so far. Southern South America (SSA) is the likely source of eolian material because Site MD07-3133, located closer to the continent, has slightly higher MS values than Site MD07-3134, and also the MS record of Patagonian Site SALSA shows comparable variability. Patagonia was the dust source for both the Scotia Sea and East Antarctica. Dust fluxes were several times higher during glacial times, when atmospheric circulation was either stronger or shifted in latitude, sea level was lowered, shelf surfaces were exposed, and environmental conditions in SSA were dominated by glaciers and extended outwash plains. Hence, MS records of SO deep-sea sediment are reliable tracers of atmospheric circulation, allowing for chronologically-constrained reconstructions of the circum Antarctic paleoclimate history.

Concentrations of tektites and chemical elements in sediments from ODP Holes 121-758B and 124-769A

Recently published studies of Ocean Drilling Project (ODP) cores from near southeast Asia revealed microtektite contents much higher than those in previously studied cores, suggesting that Ir contents might be enhanced in the tektite-bearing horizons. We determined a positive Ir anomaly in ODP core 758B from the Ninetyeast Ridge, eastern Indian Ocean; the peak Ir concentration of 190 pg/ g was ~2X the continuum level. The net Ir fluence is 1.8+/-0.5 ng/cm**2 over the depth interval from 10.87 to 11.32 m; a small additional peak also associated with microtektites contributes another 0.5 ng Ir/cm**2. Concentrations of Ir in core 769A show more scatter, but a small Ir enhancement is associated with the peak microtektite abundance; our best estimate of the poorly constrained fluence is 1.3+/-0.7 ng/cm**2. Data on deep-sea cores show that the microtektite fluence falls exponentially away from southeast Asia, the fluence dropping a factor of 2 in ~400 km. In southeast Asia the trend merges with a roughly estimated mass fluence of ~1.1 g/cm**2 inferred from evidence of a melt sheet in northeast Thailand. Integration of the inferred distribution yields a total mass of Australasian tektites of 3.2x10**16 g, much higher than previous estimates. Assuming a similar fallout distribution for the impactor and a chondritic composition allows us to calculate its mass to be 1.5x10**15 g, about 3 orders of magnitude smaller than the minimum mass of the impactor responsible for the extinctions at the end of the Cretaceous.

Uranium concentrations and isotope ratios Bahamas seawater and ODP holes pore-water

The geometry, timing, and rate of fluid-flow through carbonate margins and platforms is not well constrained. In this study, we use U concentrations and isotope ratios measured on small volumes of pore-water from Bahamas slope sediment, coupled with existing chlorinity data, to place constraints on the fluid-flow in this region and, by implication, other carbonate platforms. These data also allow an assessment of the behaviour of U isotopes in an unusually well constrained water-rock system. We report pore-water U concentrations which are controlled by dissolution of high-U organic material at shallow depths in the sediment and by reduction of U to its insoluble 4+ state at greater depths. The dominant process influencing pore-water (234U/238U) is alpha recoil. In Holocene sediments, the increase of pore-water (234U/238U) due to recoil provides an estimate of the horizontal flow rate of 11 cm/year, but with considerable uncertainty. At depths in the sediment where conditions are reducing, features in the U concentration and (234U/238U) profiles are offset from one another which constrains the effective diffusivity for U in these sediments to be c. 1-2 * 10**-8 cm**2/s. At depths between the Holocene and these reducing sediments, pore-water (234U/238U) values are unusually low due to a recent increase in the dissolution rate of grain surfaces. This suggests a strengthening of fluid flow, probably due to the flooding of the banks at the last deglaciation and the re-initiation of thermally-driven venting of fluid on the bank top and accompanying recharge on the slopes. Interpretation of existing chlorinity data, in the light of this change in flow rate, constrain the recent horizontal flow rate to be 10.6 ( 3.4) cm/year. Estimates of flow rate from (234U/238U) and Cl[-] are therefore in agreement and suggest flow rates close to those predicted by thermally-driven models of fluid flow. This agreement supports the idea that flow within the Bahamas Banks is mostly thermally driven and suggests that flow rates on the order of 10 cm/year are typical for carbonate platforms where such flow occurs.

U-Th dating of marine sediments from the Bahamas slope

In order to understand the driving forces for Pleistocene climate change more fully we need to compare the timing of climate events with their possible forcing. In contrast to the last interglacial (marine isotope stage (MIS) 5) the timing of the penultimate interglacial (MIS 7) is poorly constrained. This study constrains its timing and structure by precise U-Th dating of high-resolution delta18O records from aragonite-rich Bahamian slope sediments of ODP Leg 166 (Sites 1008 and 1009). The major glacial-interglacial cycles in delta18O are distinct within these cores and some MIS 7 substages can be identified. These sediments are well suited for U-Th dating because they have uranium concentrations of up to 12 ppm and very low initial 230Th contributions with most samples showing 230Th/232Th activity ratio of >75. U and Th concentrations and isotope ratios were measured by thermal ionisation mass spectrometry and multiple collector inductively coupled plasma mass spectrometry, with the latter providing dramatically better precision. Twenty-nine of the 41 samples measured have a delta234U value close to modern seawater suggesting that they have experienced little diagenesis. Ages from 27 of the 41 samples were deemed reliable on the basis of both their U and their Th isotope ratios. Ages generally increase with depth, although we see a repeated section of stratigraphy in one core. Extrapolation of constant sedimentation rate through each substage suggests that the peak of MIS 7e lasted from ~237 to 228 ka and that 7c began at 215 ka. This timing is consistent with existing low precision radiometric dates from speleothem deposits. The beginning of both these substages appears to be slightly later than in orbitally tuned timescales. The end of MIS 7 is complex, but also appears to be somewhat later than is suggested by orbitally tuned timescales, although this event is not particularly well defined in these cores.

Geochemical composition and radiocarbon ages of cores from Lago Fugnano, Tierra del Fuego

Recent advances in the chronology and the palaeoclimatic understanding of Antarctic ice core records point towards a larger heterogeneity of latitudinal climate fluctuations than previously thought. Thus, realistic palaeoclimate reconstructions rely in the development of a tight array of well-constrained records with a dense latitudinal coverage. Climatic records from southernmost South America are critical cornerstones to link these Antarctic palaeoclimatic archives with their South American counterparts. At 54° S on the Island of Tierra del Fuego, Lago Fagnano is located in one of the most substantially and extensively glaciated regions of southernmost South America during the Late Pleistocene. This elongated lake is the largest (~110km long) and non-ice covered lake at high southern latitudes. A multi-proxy study of selected cores allows the characterisation of a Holocene sedimentary record. Detailed petrophysical, sedimentological and geochemical studies of a complete lacustrine laminated sequence reveal variations in major and trace elements, as well as organic content, suggesting high variability in environmental conditions. Comparison of these results with other regional records allows the identification of major known late Holocene climatic intervals and the proposal for a time for the onset of the Southern Westerlies in Tierra del Fuego. These results improve our understanding of the forcing mechanisms behind climate change in southernmost Patagonia.

Miocene to Pleistocene diatom biostratigraphy with a focus on Thalassiosira in sediments of ODP Leg 186 sites

Late Neogene biostratigraphy of diatoms has been investigated from two sites occupied during Ocean Drilling Program (ODP) Leg 186 off the coast of northeast Japan. A unique aspect of ODP Leg 186 was the installation of two permanent borehole geophysical observatories at the deep-sea terrace along the Japan Trench. The Neogene subsidence history of the forearc was documented from both Sites 1150 and 1151, and Quaternary to middle Miocene (16 Ma) sediments represent a nearly continuous stratigraphic sequence including numerous ash records, especially during the past 9 m.y. Diatoms are found in most samples in variable abundance and in a moderately well preserved state throughout the sequence. The assemblages are characterized consistently by age-diagnostic species of Denticulopsis and Neodenticula found in regions of high surface water productivity typical of middle to high latitudes. The Neogene North Pacific diatom zonation divides the Miocene to Quaternary sequences fundamentally well, except that the latest Miocene through early Pliocene Thalassiosira oestrupii Subzone is not applicable. Miocene and late Pliocene through Pleistocene diatom datum levels that have been proven to be of great stratigraphic utility in the North Pacific Ocean appear to be nearly isochronous within the level of resolution constrained by core catcher sample spacing. The taxonomy and stratigraphy of previously described species determined to be useful across the Miocene/Pliocene boundary have been investigated on the basis of the evolutionary changes within the Thalassiosira trifulta group. The biostratigraphically important forms belonging to the genus Thalassiosira have been illustrated with scanning electron micrographs.

Geochemistry of Pliocene-Pleistocene volcanic sands of ODP Site 136-842

Several distinct, thin (2-7 cm), volcanic sand layers ("ashes") were recovered in the upper portions of Holes 842A and 842B. These holes were drilled 320 km west of the island of Hawaii on the outer side of the arch that surrounds the southern end of the Hawaiian chain. These layers are Pliocene to Pleistocene in age, graded, and contain fresh glass and mineral fragments (mainly olivine, plagioclase, and clinopyroxene) and tests of Pleistocene to Eocene radiolarians. The glass fragments are weakly vesicular and blocky to platy in shape. The glass and olivine fragments from individual layers have large ranges in composition (i.e, larger than expected for a single eruption). These features are inconsistent with an explosive eruption origin for the sands. The only other viable mechanism for transporting these sands hundreds of kilometers from their probable source, the Hawaiian Islands, is turbidity currents. These currents were probably related to several of the giant debris slides that were identified from Gloria sidescan images around the islands. These currents would have run over the ~500-m-high Hawaiian Arch on their way to Site 842. This indicates that the turbidity currents were at least 325 m thick. Paleomagnetic and biostratigraphic data allow the ages of the sands to be constrained and, thus, related to particular Hawaiian debris flows. These correlations were checked by comparing the compositions of the glasses from the sands with those of glasses and rocks from islands with debris flows directed toward Site 842. Good correlations were found for the 110-ka slide from Mauna Loa and the ~1.4-Ma slide from Lanai. The correlation with Kauai is poor, probably because the data base for that volcano is small. The low to moderate sulfur content of the sand glasses indicates that they were derived from moderately to strongly degassed lavas (shallow marine or subaerially erupted), which correlates well with the location of the landslide scars on the flanks of the Hawaiian volcanoes. The glass sands may have been formed by brecciation during the landslide events or spallation and granulation as lava erupted into shallow water.

Geochemistry of boninite and harzburgite of ODP Leg 125 basement samples

Holes drilled into the volcanic and ultrabasic basement of the Izu-Ogasawara and Mariana forearc terranes during Leg 125 provide data on some of the earliest lithosphere created after the start of Eocene subduction in the Western Pacific. The volcanic basement contains three boninite series and one tholeiite series. (1) Eocene low-Ca boninite and low-Ca bronzite andesite pillow lavas and dikes dominate the lowermost part of the deep crustal section through the outer-arc high at Site 786. (2) Eocene intermediate-Ca boninite and its fractionation products (bronzite andesite, andesite, dacite, and rhyolite) make up the main part of the boninitic edifice at Site 786. (3) Early Oligocene intermediate-Ca to high-Ca boninite sills or dikes intrude the edifice and perhaps feed an uppermost breccia unit at Site 786. (4) Eocene or Early Oligocene tholeiitic andesite, dacite, and rhyolite form the uppermost part of the outer-arc high at Site 782. All four groups can be explained by remelting above a subduction zone of oceanic mantle lithosphere that has been depleted by its previous episode of partial melting at an ocean ridge. We estimate that the average boninite source had lost 10-15 wt% of melt at the ridge before undergoing further melting (5-10%) shortly after subduction started. The composition of the harzburgite (<2% clinopyroxene, Fo content of about 92%) indicates that it underwent a total of about 25% melting with respect to a fertile MORB mantle. The low concentration of Nb in the boninite indicates that the oceanic lithosphere prior to subduction was not enriched by any asthenospheric (OIB) component. The subduction component is characterized by (1) high Zr and Hf contents relative to Sm, Ti, Y, and middle-heavy REE, (2) light REE-enrichment, (3) low contents of Nb and Ta relative to Th, Rb, or La, (4) high contents of Na and Al, and (5) Pb isotopes on the Northern Hemisphere Reference Line. This component is unlike any subduction component from active arc volcanoes in the Izu-Mariana region or elsewhere. Modeling suggests that these characteristics fit a trondhjemitic melt from slab fusion in amphibolite facies. The resulting metasomatized mantle may have contained about 0.15 wt% water. The overall melting regime is constrained by experimental data to shallow depths and high temperatures (1250? C and 1.5 kb for an average boninite) of boninite segregation. We thus envisage that boninites were generated by decompression melting of a diapir of metasomatized residual MORB mantle leaving the harzburgites as the uppermost, most depleted residue from this second stage of melting. Thermal constraints require that both subducted lithosphere and overlying oceanic lithosphere of the mantle wedge be very young at the time of boninite genesis. This conclusion is consistent with models in which an active transform fault offsetting two ridge axes is placed under compression or transpression following the Eocene plate reorganization in the Pacific. Comparison between Leg 125 boninites and boninites and related rocks elsewhere in the Western Pacific highlights large regional differences in petrogenesis in terms of mantle mineralogy, degree of partial melting, composition of subduction components, and the nature of pre-subduction lithosphere. It is likely that, on a regional scale, the initiation of subduction involved subducted crust and lithospheric mantle wedge of a range of ages and compositions, as might be expected in this type of tectonic setting.

Fractional abundance and indices for brGDGTs and crenarchaeol in dust samples, surface water samples and surface sediment samples

Branched glycerol dialkyl glycerol tetraethers (brGDGTs) are membrane lipids produced by soil bacteria and occur in near coastal marine sediments as a result of soil organic matter input. Their abundance relative to marine-derived crenarchaeol, quantified in the BIT index, generally decreases offshore. However, in distal marine sediments, low relative amounts of brGDGTs can often still be observed. Sedimentary in situ production as well as dust input have been suggested as potential, though as yet not well constrained, sources. In this study brGDGT distributions in dust were examined and compared with those in distal marine sediments. Dust was sampled along the equatorial West African coast and brGDGTs were detected in most of the samples, albeit in low abundance. Their degree of methylation and cyclisation, expressed in the MBT' (methylation index of branched tetraethers) and DC (degree of cyclisation) indices, respectively, were comparable with those for African soils, their presumed source. Comparison of DC index values for brGDGTS in global soils, Congo deep-sea river fan sediments and dust with those of distal marine sediments clearly showed, however, that distal marine sediments had significantly higher values. This distinctive distribution is suggestive of sedimentary in situ production as a source of brGDGTs in marine sediments, rather than dust input. The presence of in situ produced brGDGTs in marine sediments means that caution should be exercised when applying the MBT'-CBT palaeothermometer to sediments with low BIT index values, i.e. < 0.1, based on our dataset.

Geochemical composition of surface sediments in the Trondheimsfjord, central Norway

High sedimentation rates in fjords provide excellent possibilities for high resolution sedimentary and geochemical records over the Holocene. As a baseline for an improved interpretation of geochemical data from fjord sediment cores, this study aims to investigate the inorganic/organic geochemistry of surface sediments and to identify geochemical proxies for terrestrial input and river discharge in the Trondheimsfjord, central Norway. Sixty evenly distributed surface sediment samples were analysed for their elemental composition, total organic carbon (Corg), nitrogen (Norg) and organic carbon stable isotopes (d13Corg), bulk mineral composition and grain size distribution. Our results indicate carbonate marine productivity to be the main CaCO3 source. Also, a strong decreasing gradient of marine-derived organic matter from the entrance towards the fjord inner part is consistent with modern primary production data. We show that the origin of the organic matter as well as the distribution of CaCO3 in Trondheimsfjord sediments can be used as a proxy for the variable inflow of Atlantic water and changes in river runoff. Furthermore, the comparison of grain size independent Al-based trace element ratios with geochemical analysis from terrigenous sediments and bedrocks provides evidence that the distribution of K/Al, Ni/Al and K/Ni in the fjord sediments reflect regional sources of K and Ni in the northern and southern drainage basin of the Trondheimsfjord. Applying these findings to temporally well-constrained sediment records will provide important insights into both the palaeoenvironmental changes of the hinterland and the palaeoceanographic modifications in the Norwegian Sea as response to rapid climate changes and associated feedback mechanisms.

Geochemistry of oceanic flood basalts from the Caribbean plate

The thick oceanic crust of the Caribbean plate appears to be the tectonized remnant of an eastern Pacific oceanic plateau that has been inserted between North and South America. The emplacement of the plateau into its present position has resulted in the obduction and exposure of its margins, providing an opportunity to study the age relations, internal structure and compositional features of the plateau. We present the results of 40Ar-39Ar radiometric dating, major-, trace-element, and isotopic compositions of basalts from some of the exposed sections as well as drill core basalt samples from Leg 15 of the Deep Sea Drilling Project. Five widely spaced, margin sections yielded ages ranging from 91 to 88 Ma. Less well-constrained radiometric ages from the drill cores, combined with the biostratigraphic age of surrounding sediments indicate a minimum crystallization age of ~90 Ma in the Venezuelan Basin. The synchroneity of ages across the region is consistent with a flood basalt origin for the bulk of the Caribbean plateau i.e., large volume, rapidly erupted, regionally extensive volcanism.. The ages and compositions are also consistent with plate reconstructions that place the Caribbean plateau in the vicinity of the Galápagos hotspot at its inception. The trace-element and isotopic compositions of the ~90 Ma rocks indicate a depleted mantle and an enriched, plume-like mantle were involved in melting to varying degrees across the plateau. Within the same region, a volumetrically secondary, but widespread magmatic event occurred at 76 Ma, as is evident in Curacao, western Colombia, Haiti, and at DSDP Site 152/ODP Site 1001 near the Hess Escarpment. Limited trace-element data indicate that this phase of magmatism was generally more depleted than the first. We speculate that magmatism may have resulted from upwelling of mantle, still hot from the 90 Ma event, during lithospheric extension attending gravitational collapse of the plateau, andror tectonic emplacement of the plateau between North and South America. Still younger volcanics are found in the Dominican Republic (69 Ma) and the Quepos Peninsula of Costa Rica (63 Ma). The latter occurrence conceivably formed over the Galápagos hotspot and subsequently accreted to the western edge of the plateau during subduction of the Farallon plate.

Specific microbial biomarker signatures from Congo fan sediments of ODP site 175-1075, Congo soils, wetland and estuary sediments

Methane (CH4) is a strong greenhouse gas known to have perturbed global climate in the past, especially when released in large quantities over short time periods from continental or marine sources. It is therefore crucial to understand and, if possible, quantify the individual and combined response of these variable methane sources to natural climate variability. However, past changes in the stability of greenhouse gas reservoirs remain uncertain and poorly constrained by geological evidence. Here, we present a record from the Congo fan of a highly specific bacteriohopanepolyol (BHP) biomarker for aerobic methane oxidation (AMO), 35-aminobacteriohopane-30,31,32,33,34-pentol (aminopentol), that identifies discrete periods of increased AMO as far back as 1.2 Ma. Fluctuations in the concentration of aminopentol, and other 35-aminoBHPs, follow a pattern that correlates with late Quaternary glacial-interglacial climate cycles, with highest concentrations during warm periods. We discuss possible sources of aminopentol, and the methane consumed by the precursor methanotrophs, within the context of the Congo River setting, including supply of methane oxidation markers from terrestrial watersheds and/or marine sources (gas hydrate and/or deep subsurface gas reservoir). Compound-specific carbon isotope values of -30 per mil to -40 per mil for BHPs in ODP 1075 and strong similarities between the BHP signature of the core and surface sediments from the Congo estuary and floodplain wetlands from the interior of the Congo River Basin, support a methanotrophic and likely terrigenous origin of the 35-aminoBHPs found in the fan sediments. This new evidence supports a causal connection between marine sediment BHP records of tropical deep sea fans and wetland settings in the feeding river catchments, and thus tropical continental hydrology. Further research is needed to better constrain the different sources and pathways of methane emission. However, this study identifies the large potential of aminoBHPs, in particular aminopentol, to trace and, once better calibrated and understood, quantify past methane sources and fluxes from terrestrial and potentially also marine sources.

Geochemistry of oceanic crust beneath Costa Rica

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge (Feigenson et al., 2004, doi:10.1029/2003GC000621; Herrstom et al., 1995, doi:10.1130/0091-7613(1995)023<0617:VILCAW>2.3.CO;2; Abratis and Woerner, 2001) or eroded fore-arc complexes (Goss and Kay, 2006, doi:10.1029/2005GC001163) but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (~63-190 mm/yr) and are comparable to the magnitude of subducting Cocos plate motion (approx85 mm/yr). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.