Mg/Ca ratios in planktonic foraminifera

The primary Mg/Ca ratio of foraminiferal shells is a potentially valuable paleoproxy for sea surface temperature (SST) reconstructions. However, the reliable extraction of this ratio from sedimentary calcite assumes that we can overcome artifacts related to foraminiferal ecology and partial dissolution, as well as contamination by secondary calcite and clay. The standard batch method for Mg/Ca analysis involves cracking, sonicating, and rinsing the tests to remove clay, followed by chemical cleaning, and finally acid-digestion and single-point measurement. This laborious procedure often results in substantial loss of sample (typically 30-60%). We find that even the earliest steps of this procedure can fractionate Mg from Ca, thus biasing the result toward a more variable and often anomalously low Mg/Ca ratio. Moreover, the more rigorous the cleaning, the more calcite is lost, and the more likely it becomes that any residual clay that has not been removed by physical cleaning will increase the ratio. These potentially significant sources of error can be overcome with a flow-through (FT) sequential leaching method that makes time- and labor-intensive pretreatments unnecessary. When combined with time-resolved analysis (FT-TRA) flow-through, performed with a gradually increasing and highly regulated acid strength, produces continuous records of Mg, Sr, Al, and Ca concentrations in the leachate sorted by dissolution susceptibility of the reacting material. Flow-through separates secondary calcite from less susceptible biogenic calcite and clay, and further resolves the biogenic component into primary and more resistant fractions. FT-TRA reliably separates secondary calcite (which is not representative of original life habitats) from the more resistant biogenic calcite (the desired signal) and clay (a contaminant of high Mg/Ca, which also contains Al), and further resolves the biogenic component into primary and more resistant fractions that may reflect habitat or other changes during ontogeny. We find that the most susceptible fraction of biogenic calcite in surface dwelling foraminifera gives the most accurate value for SST and therefore best represents primary calcite. Sequential dissolution curves can be used to correct the primary Mg/Ca ratio for clay, if necessary. However, the temporal separation of calcite from clay in FT-TRA is so complete that this correction is typically <=2%, even in clay-rich sediments. Unlike hands-on batch methods, that are difficult to reproduce exactly, flow-through lends itself to automation, providing precise replication of treatment for every sample. Our automated flow-through system can process 22 samples, two system blanks, and 48 mixed standards in <12 hours of unattended operation. FT-TRA thus represents a faster, cheaper, and better way to determine Mg/Ca ratios in foraminiferal calcite.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Sea surafe temperture calculation for the Arabian Sea

Sea surface temperature (SST) and seawater d18O (d18Ow) were reconstructed in a suite of sediment cores from throughout the Arabian Sea for four distinct time intervals (0 ka, 8 ka, 15 ka, and 20 ka) with the aim of understanding the history of the Indian Monsoon and the climate of the Arabian Sea region. This was accomplished through the use of paired Mg/Ca and d18O measurements of the planktonic foraminifer Globigerinoides ruber. By analyzing basin-wide changes and changes in cross-basinal gradients, we assess both monsoonal and regional-scale climate changes. SST was colder than present for the majority of sites within all three paleotime slices. Furthermore, both the Indian Monsoon and the regional Arabian Sea mean climate have varied substantially over the past 20 kyr. The 20 ka and 15 ka time slices exhibit average negative temperature anomalies of 2.5°-3.5°C attributable, in part, to the influences of glacial atmospheric CO2 concentrations and large continental ice sheets. The elimination of the cross-basinal SST gradient during these two time slices likely reflects a decrease in summer monsoon and an increase in winter monsoon strength. Changes in d18Ow that are smaller than the d18O signal due to global ice volume reflect decreased evaporation and increased winter monsoon mixing. SSTs throughout the Arabian Sea were still cooler than present by an average of 1.4°C in the 8 ka time slice. These cool SSTs, along with lower d18Ow throughout the basin, are attributed to stronger than modern summer and winter monsoons and increased runoff and precipitation. The results of this study underscore the importance of taking a spatial approach to the reconstruction of processes such as monsoon upwelling.

Phosphorus components, CaCO3 and biogenic opal in sediments from ODP Leg 167 sites

The ocean history of reactive phosphorus (P) (i.e., dissolved P available to fuel oceanic primary productivity) is of interest because of the role of P as a biolimiting nutrient, and knowledge of P burial in marine sediments is key to testing hypotheses about temporal changes in P input or output fluxes. Our understanding of the history of the P cycle over the Cenozoic has increased substantially with temporal records of reactive P mass accumulation rates from open-ocean Pacific and Atlantic equatorial sites. However, questions about the relative importance of nutrient burial in ocean-margin sediments relative to burial in open-ocean sediments and about the extent of P remobilization in organic-rich, reducing environments characteristic of margin sediments remain unresolved. Nutrient burial in oceanic boundary current systems has been suggested to have a controlling role in oceanic nutrient budgets in certain time intervals (Vincent and Berger, 1985, doi:10.1029/GM032p0455), with higher sediment accumulation rates balancing the limited spatial extent of these sediments. Some investigators suggest that remobilization of P from reducing sediments in margin settings is a significant positive feedback to primary productivity (e.g., Van Cappellan and Ingall, 1994, doi:10.1029/94PA01455), whereas other results indicate that both P uptake and P release may occur in these settings depending on the balance of organic carbon and iron supply to the sediments and on the oxygenation of bottom waters (McManus et al., 1997, doi:10.1016/S0016-7037(97)00138-5). It is important to quantitatively understand the geochemistry of reactive P in margin sediments, where productivity and delivery of organic-rich material to the sediments in relatively shallow-water settings is often sufficient to promote anoxia in interstitial waters. To address these questions, we determined the P concentrations and geochemistry in sediment samples from eight sites drilled during Ocean Drilling Program (ODP) Leg 167, California margin (Sites 1010-1012, 1014, 1016-1017, and 1021-1022). These results are the first records of reactive P concentrations on long time scales-required for the calculation of P accumulation rates-for sediments from a highly productive eastern boundary current setting. In addition, we determined calcium carbonate contents and biogenic silica concentrations to define the environments of sedimentary production, burial, and diagenesis.

Thermal diagenesis of clay minerals within volcanogenic material from ODP Site 135-841

Clay minerals are examined in detail in the sediment from the Tonga Trench margin at Site 841 (Leg 135 ODP). The changes in amount and nature of secondary clays with depth provide an alternative explanation for the intensive alteration of volcanogenic material at convergent margins. A characteristic distribution of clay minerals with depth shows four distinct zones unexplainable by simple burial diagenesis processes. These are named the upper, reactive, lower and rhyolitic zones. The reactive zone is intercalated with numerous sills and is characterized by the dominant iron-rich clays such as saponite, corrensite and chlorite associated with analcime. The occurrence of such iron-rich clays, mostly associated with a large amount of analcime, yields chemical and mineralogical evidence for thermal diagenesis. The required heat for the diagenetic process was transferred from recently intruded basaltic andesite sills. In the vicinity of these intrusions, the iron-rich clay minerals may have formed at temperatures up to 200°C. A zoning with respect to clay and zeolite minerals indicates that the influence of the palaeoheat flow decreased with the distance from the intrusion. The formation of interlayered I/S, illite, kaolinite and aluminous chlorite, which are recognized as major secondary minerals within the rhyolitic complex, was mainly controlled by both early diagenesis at moderately elevated temperatures, and since the Eocene by burial diagenesis at low temperatures. The occurrence of a steam zone in an early stage of the intrusion is restricted to Miocene tuffs and has overprinted the early alteration of the volcanogenic material within the tuffs and has changed the originally pristine composition of the pore fluids.

Flow rates, hydraulic conductivities and permeabilities of sediments from ODP Holes 190-1173A and 190-1174B

Understanding the role of fluids in active accretionary prisms requires quantitative knowledge of parameters such as permeability. We report here the results of permeability tests on four samples from Ocean Drilling Program Leg 190 at the Nankai Trough accretionary prism-two from Site 1173 and two from Site 1174. Volcanic ash is present in one of the samples; otherwise, the material is hemipelagic mud. A constant-rate-of-flow technique was used at various effective pressures and rates of flow. The permeability of the four samples ranges between 10**-15 and 10**-18 m**2, with the ash-bearing sample showing the highest values.

Debris flow magnitudes: a global catalog

The episodic occurrence of debris flow events in response to stochastic precipitation and wildfire events makes hazard prediction challenging. Previous work has shown that frequency-magnitude distributions of non-fire-related debris flows follow a power law, but less is known about the distribution of post-fire debris flows. As a first step in parameterizing hazard models, we use frequency-magnitude distributions and cumulative distribution functions to compare volumes of post-fire debris flows to non-fire-related debris flows. Due to the large number of events required to parameterize frequency-magnitude distributions, and the relatively small number of post-fire event magnitudes recorded in the literature, we collected data on 73 recent post-fire events in the field. The resulting catalog of 988 debris flow events is presented as an appendix to this article. We found that the empirical cumulative distribution function of post-fire debris flow volumes is composed of smaller events than that of non-fire-related debris flows. In addition, the slope of the frequency-magnitude distribution of post-fire debris flows is steeper than that of non-fire-related debris flows, evidence that differences in the post-fire environment tend to produce a higher proportion of small events. We propose two possible explanations: 1) post-fire events occur on shorter return intervals than debris flows in similar basins that do not experience fire, causing their distribution to shift toward smaller events due to limitations in sediment supply, or 2) fire causes changes in resisting and driving forces on a package of sediment, such that a smaller perturbation of the system is required in order for a debris flow to occur, resulting in smaller event volumes.

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.

Isotope analysis and heat flow measurements of Maria S. Merian cruise MSM21/4 on the western Svalbard margin

Methane hydrate is an ice-like substance that is stable at high-pressure and low temperature in continental margin sediments. Since the discovery of a large number of gas flares at the landward termination of the gas hydrate stability zone off Svalbard, there has been concern that warming bottom waters have started to dissociate large amounts of gas hydrate and that the resulting methane release may possibly accelerate global warming. Here, we can corroborate that hydrates play a role in the observed seepage of gas, but we present evidence that seepage off Svalbard has been ongoing for at least three thousand years and that seasonal fluctuations of 1-2°C in the bottom-water temperature cause periodic gas hydrate formation and dissociation, which focus seepage at the observed sites.

Table 2: Stream flow measurement at the gauge station Vernagtbach

From 1974 to 1982 repeated tracer tests using fluorescent dyes were carried out in the highly glaciated drainage basin of Vernagtbach. These tests enabled the quantitative determination of the runoff in the forefield of the Vernagtferner, the calculation of travel times of the stream water and estimations of the relative contributions to the entire runoff originating from individual streams. In addition, tracer tests were carried out in the firn area of the glacier resulting in data concerning the storage and travel time of meltwater inside the glacier.