Fission track analyses of ODP Hole 129-800A samples

Four models of fission track annealing in apatite are compared with measured fission track lengths in samples from Site 800 in the East Mariana Basin, Ocean Drilling Program Leg 129, given an independently determined temperature history. The temperature history of Site 800 was calculated using a one-dimensional, compactive, conductive heat flow model assuming two end-member thermal cases: one for cooling of Jurassic ocean crust that has experienced no subsequent heating, and one for cooling of Cretaceous ocean crust. Because the samples analyzed were only shallowly buried and because the tectonic history of the area since sample deposition is simple, resolution of the temperature history is high. The maximum temperature experienced by the sampled bed is between 16°-21°C and occurs at 96 Ma; temperatures since the Cretaceous have dropped in spite of continued pelagic sediment deposition because heat flow has continued to decay exponentially and bottom-water temperatures have dropped. Fission tracks observed within apatite grains from the sampled bed are 14.6 +/- 0.1 µm (1 sigma) long. Given the proposed temperature history of the samples, one unpublished and three published models of fission track annealing predict mean track lengths from 14.8 to 15.9 µm. These models require temperatures as much as 40°C higher than the calculated paleotemperature maximum of the sampled bed to produce the same degree of track annealing. Measured and predicted values are different because annealing models are based on extrapolation of high temperature laboratory data to geologic times. The model that makes the closest prediction is based on the greatest number of experiments performed at low temperature and on an apatite having composition closest to that of the core samples.

Dust flux in the Arabian Sea

Eolian grain size and flux were measured on samples from 11 Arabian Sea sediment traps deployed 200-1250 km offshore. The timing of increased grain size is coincident with the onset of strong summer monsoon winds and dust storm activity over the Arabian Peninsula and Middle East. Data spanning a full annual cycle show that eolian grain size is highly correlated with barometric pressure (r=-0.91) and wind speed (r=0.84), enabling calibration of the downcore record in terms of these primary meteorological variables. Eolian flux is highly correlated with organic carbon flux (r=0.80); both increase 6-8 weeks after the grain size increase and summer monsoon onset. This lag, and the low correlation between eolian grain size and eolian flux (r=0.36), likely result from the differential sinking rates of large and small dust particles in the surface waters as well as biological scavenging associated with monsoon-induced productivity.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Borehole temperature, submarine permafrost methane concentrations and pore water chemistry measured on borehole BK2, central Laptev Sea shelf

Submarine permafrost degradation has been invoked as a cause for recent observations of methane emissions from the seabed to the water column and atmosphere of the East Siberian shelf. Sediment drilled 52 m down from the sea ice in Buor Khaya Bay, central Laptev Sea revealed unfrozen sediment overlying ice-bonded permafrost. Methane concentrations in the overlying unfrozen sediment were low (mean 20 µM) but higher in the underlying ice-bonded submarine permafrost (mean 380 µM). In contrast, sulfate concentrations were substantially higher in the unfrozen sediment (mean 2.5 mM) than in the underlying submarine permafrost (mean 0.1 mM). Using deduced permafrost degradation rates, we calculate potential mean methane efflux from degrading permafrost of 120 mg/m**2 per year at this site. However, a drop of methane concentrations from 190 µM to 19 µM and a concomitant increase of methane d13C from -63 per mil to -35 per mil directly above the ice-bonded permafrost suggest that methane is effectively oxidized within the overlying unfrozen sediment before it reaches the water column. High rates of methane ebullition into the water column observed elsewhere are thus unlikely to have ice-bonded permafrost as their source.

Chronology and petrology in MIS4-2 of DSDP Site 94-609

Ice-rafting evidence for a '1500-year cycle' sparked considerable debate on millennial-scale climate change and the role of solar variability. Here, we reinterpret the last 70,000 years of the subpolar North Atlantic record, focusing on classic DSDP Site 609, in the context of newly available raw data, the latest radiocarbon calibration (Marine09) and ice core chronology (GICC05), and a wider range of statistical methodologies. A ~1500-year oscillation is primarily limited to the short glacial Stage 4, the age of which is derived solely from an ice flow model (ss09sea), subject to uncertainty, and offset most from the original chronology. Results from the most well-dated, younger interval suggest that the original 1500 ± 500 year cycle may actually be an admixture of the ~1000 and ~2000 cycles that are observed within the Holocene at multiple locations. In Holocene sections these variations are coherent with 14C and 10Be estimates of solar variability. Our new results suggest that the '1500-year cycle' may be a transient phenomenon whose origin could be due, for example, to ice sheet boundary conditions for the interval in which it is observed. We therefore question whether it is necessary to invoke such exotic explanations as heterodyne frequencies or combination tones to explain a phenomenon of such fleeting occurrence that is potentially an artifact of arithmetic averaging.

Age model of ODP Site 162-983

At ODP Site 983, relative geomagnetic paleointensity and planktic and benthic delta18O records have been acquired for the last 350 kyr. The mean sedimentation rate in this interval is 11.3 cm/kyr. Magnetic properties and hysteresis ratios indicate that pseudo-single domain magnetite is the remanence carrier. Volume susceptibility (kappa), anhysteretic (ARM) and isothermal (IRM) remanence values vary by a factor of 3-4, well within the criteria usually cited for paleointensity studies. Natural remanent magnetization (NRM) is normalized by ARM and IRM to acquire the paleointensity proxy. Arithmetic means of NRM/ARM and NRM/IRM, calculated for five demagnetization steps in the 25-45 mT range, constitute the relative paleointensity estimates. Some paleointensity lows (particularly those at ~40, ~120 and ~188 ka) are associated with directional excursions of the field, especially the event at ~188 ka (referred to here as the Iceland Basin Event) that constitutes a short-lived polarity reversal. For the last 200 kyr, the records can be correlated with other high-resolution paleointensity records such as those from the Labrador Sea, Mediterranean/Somali Basin and Sulu Sea, implying that the millennial scale features are globally synchronous. A labeling system for paleointensity features is proposed that ties prominent highs and lows to oxygen isotope stages.

Magnetic susceptibility anisotropy in deformed sediments of ODP Leg 110

The anisotropy of magnetic susceptibility documents the generation of tectonically produced fabrics in sediments that macroscopically show no evidence of this disruption. The fabric observed in initial accretion is largely produced by overprinting of the original sedimentary susceptibility anisotropy by an E-W horizontal tectonic shortening and vertical extension. The response of the sediments to stress during initial accretion is variable, particularly near the sediment surface, and appears to reflect the inhomogeneous distribution of strain rate in the overthrust sequence. The susceptibility anisotropy of sediments possessing scaly fabric is consistent with the strong orientation of Phyllosilicates seen in thin section, producing a Kmin normal to the scalyness. The slope sediments deposited on the accreted sequence are also affected by tectonic shortening. The accreted sequences at Sites 673 and 674 show a complex history of fabric modification, with previous tectonic fabrics overprinted by later fabric modifications, pointing to continued tectonic shortening during the accretion process. The form of the susceptibility anisotropy axes at Sites 673 and 674 is consistent with NESW shortening, probably reflected in the NW-SE surface expression of the out-of-sequence thrusts. The susceptibility anisotropy appears to document a downhole change in the trend of shortening from E to W at the surface to more NESW at depth, probably as a result of the obliquely trending basement ridge, the Tiburon Rise.

Organic geochemistry of sediments at DSDP Legs 56-57 Holes

The quantity, type, and maturity of the organic matter of Quaternary and Tertiary sediments from the Japan Trench (DSDP Leg 56, Sites 434 and 436; and Leg 57, Sites 438, 439 and 440) were determined. The hydrocarbons in lipid extracts were analyzed by capillary- column gas chromatography and combined gas chromatography/ mass spectrometry. Kerogen concentrates were investigated by microscopy, and vitrinite-reflectance values were determined. Measured organic-carbon values were in the range of 0.13 to 1.00 per cent. Extract yields, however, were extremely low. Normalized to organic carbon, total extracts ranged from 4.1 to 15.7 mg/g Corg. Gas chromatography of non-aromatic hydrocarbons showed that all sediments, except one Oligocene sample, contained very immature, mainly terrigenous organic material. This was indicated by n-alkane maxima at C29 and C31 and high odd-carbon-number predominances. Unsaturated steroid hydrocarbons were found to be major cyclic compounds in lower- and middle-Miocene samples from the upper inner trench slope (Sites 438 and 439). Perylene was the dominating aromatic hydrocarbon in all but the Oligocene sample. Microscopy showed kerogens rich in terrigenous organic particles, with a major portion of recycled vitrinite. Nevertheless, almost all the liptinite particles appeared to be primary. This is a paradox, as the bulk of the samples were composed of hemipelagic mineral matter with a major siliceous biogenic (planktonic) component. A trend of reduced size and increased roundness can be seen for the vitrinite/ inertinite particles from west to east (from upper inner trench slope to outer trench slope). All sediments but one are relatively immature, with mean huminite-reflectance values (Ro)in the range of 0.30 to 0.45 per cent. The oldest and deepest sediment investigated, an Oligocene sandstone from Site 439, yielded a mean vitrinitereflectance value of 0.74 per cent and a mature n-alkane distribution. This sample may indicate a geothermal event in late Oligocene time. It failed to affect the overlying lower Miocene and may have been caused by an intrusion. Boulders of acidic igneous rocks in the Oligocene can be interpreted as witnesses of nearby volcanic activity accompanied by intrusions.

Geochemistry and grain-size of sediment core PC29A, Quitralco Fjord, Chilean Patagonia

The climate of Chilean Patagonia is strongly influenced by the southern westerlies, which control the amount and latitudinal distribution of precipitation in the southern Andes. In austral summer, the Southern Westerly Wind Belt (SWWB) is restricted to the high latitudes. It expands northward in winter, which results in a strong precipitation seasonality between 35 and 45°S. Here, we present a new precipitation seasonality proxy record from Quitralco fjord (46°S), where relatively small latitudinal shifts in the SWWB result in large changes in precipitation seasonality. Our 1400 yr record is based on sedimentological and geochemical data obtained on a sediment core collected in front of a small river that drains the Patagonian Andes, which makes this site particularly sensitive to changes in river discharge. Our results show Fe/Al and Ti/Al values that are low between 600 and 1200 CE, increasing at 1200-1500 CE, and high between 1500 and 1950 CE. The increasing Fe/Al and Ti/Al values reflect a decrease in mean sediment grain-size from 30 to 20 µm, which is interpreted as a decrease in seasonal floods resulting from an equatorward shift of the SWWB. Our results suggest that, compared to present-day conditions, the SWWB was located in a more poleward position before 1200 CE. It gradually shifted towards the equator in 1200-1500 CE, where it remained in a sustained position until 1950 CE. The comparison of our record with published regional sea surface temperature (SST) reconstructions for the late Holocene shows that equatorward shifts in the SWWB are systematically coeval with decreasing SSTs and vice versa, which resembles fluctuations over glacial-interglacial timescales. We argue that the synchronicity between SST and SWWB changes during the last 1400 years represents the response of the SWWB to temperature changes in the Southern Hemisphere.