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.

Downhole variations in grain sizes at DSDP/ODP Hole 504B

The maximum grain sizes of plagioclase and magnetite in the groundmass of the sheeted dike complex drilled at Hole 504B have been measured. Downhole variations through a 440-m-long section show a crude zig-zag pattern consisting of a gradual decrease or increase followed by an abrupt jump. The gradual decrease or increase in grain size extends over many lithologic units, and hence, does not reflect variations in grain size within a single dike. Such a zig-zag pattern is well explained by grain-size variations through multiple dikes. By using the observed inclination of sheeted dikes of 81° ± 2.5°, thickness of the multiple dikes varies from 0.7 to 8.5 m and averages to 4 ± 1 m. The average thickness of individual dikes forming multiple dikes is 0.8 m. We expect such multiple dikes to be formed during rifting events beneath mid-oceanic spreading ridges. If the average expansion at rifting episodes is twice as wide as the average width of the multiple dike units, the full spreading rate of 7.2 cm/yr of Cocos Ridge gives 112 ± 33 yr for a time interval of the rifting. A simple one-dimensional conductive cooling model is applied to solidification of multiple dikes. Numerical simulations show that the grain-size variations observed through the drill hole are more consistent with a model where a new injection of a dike occurs periodically with a constant time interval rather than one where the next dike intrudes just after the solidification of the previous one. Grain-size variations within simple dikes from Iritono, Japan, and those for Makaopuhi lava lake, Hawaii, show that square root of crystallization time is linearly correlated with the logarithm of plagioclase size. By using an empirically derived relationship between these two variables, the variations of plagioclase size through Hole 504B are directly compared with the calculated times for crystallization. Each rifting episode at the Costa Rica Rift lasts for several years, and periodic injection of a new dike occurs into the center of a previously solidified multiple dike at time intervals varying from 1 to 12 months.

(Appendix) Petrographic information, trace element and Sr isotope ratios of Pito Deep samples

Fluid flow through the axial hydrothermal system at fast spreading ridges is investigated using the Sr-isotopic composition of upper crustal samples recovered from a tectonic window at Pito Deep (NE Easter microplate). Samples from the sheeted dike complex collected away from macroscopic evidence of channelized fluid flow, such as faults and centimeter-scale hydrothermal veins, show a range of 87Sr/86Sr from 0.7025 to 0.7030 averaging 0.70276 relative to a protolith with 87Sr/86Sr of ~0.7024. There is no systematic variation in 87Sr/86Sr with depth in the sheeted dike complex. Comparison of these new data with the two other localities that similar data sets exist for (ODP Hole 504B and the Hess Deep tectonic window) reveals that the extent of Sr-isotope exchange is similar in all of these locations. Models that assume that fluid-rock reaction occurs during one-dimensional (recharge) flow lead to significant decreases in the predicted extent of isotopic modification of the rock with depth in the crust. These model results show systematic misfits when compared with the data that can only be avoided if the fluid flow is assumed to be focused in isolated channels with very slow fluid-rock exchange. In this scenario the fluid at the base of the crust is little modified in 87Sr/86Sr from seawater and thus unlike vent fluids. Additionally, this model predicts that some rocks should show no change from the fresh-rock 87Sr/86Sr, but this is not observed. Alternatively, models in which fluid-rock reaction occurs during upflow (discharge) as well as downflow, or in which fluids are recirculated within the hydrothermal system, can reproduce the observed lack of variation in 87Sr/86Sr with depth in the crust. Minimum time-integrated fluid fluxes, calculated from mass balance, are between 1.5 and 2.6 * 10**6 kg/m**2 for all areas studied to date. However, new evidence from both the rocks and a compilation of vent fluid compositions demonstrates that some Sr is leached from the crust. Because this leaching lowers the fluid 87Sr/86Sr without changing the rock 87Sr/86Sr, these mass balance models must underestimate the time-integrated fluid flux. Additionally, these values do not account for fluid flow that is channelized within the crust.

Noble gases sampled from three cruises in the eastern tropical Atlantic in September 2005 and June/July 2006

Upwelling velocities w in the equatorial band are too small to be directly observed. Here, we apply a recently proposed indirect method, using the observed helium isotope (3He or 4He) disequilibria in the mixed layer. The helium data were sampled from three cruises in the eastern tropical Atlantic in September 2005 and June/July 2006. A one-dimensional two-box model was applied, where the helium air-sea gas exchange is balanced by upwelling from 3He-rich water below the mixed layer and by vertical mixing. The mixing coefficients Kv were estimated from microstructure measurements, and on two of the cruises, Kv exceeded 1 x 10**-4 m**2/s, making the vertical mixing term of the same order of magnitude as the gas exchange and the upwelling term. In total, helium disequilibrium was observed on 54 stations. Of the calculated upwelling velocities, 48% were smaller than 1.0 x 10**-5 m/s, 19% were between 1.0 and 2.0 x 10**-5 m/s, 22% were between 2.0 and 4.0 x 10**-5 m/s, and on 11% of upwelling velocities exceeded this limit. The highest upwelling velocities were found in late June 2006. Meridional upwelling distribution indicated an equatorial asymmetry with higher vertical velocities between the equator and 1° to 2° south compared to north of the equator, particularly at 10°W. Associated heat flux into the mixed layer could be as high as 138 W/m**2, but this depends strongly on the chosen depths where the upwelled water comes from. By combining upwelling velocities with sea surface temperature and productivity distributions, a mean monthly equatorial upwelling rate of 19 Sv was estimated for June 2006 and a biweekly mean of 24 Sv was estimated for September 2005.

Helicopter-borne sea ice thickness measurements in the Storfjord, Spitsbergen from flight HEM_STJ06_01

Results from electromagnetic induction surveys of sea-ice thickness in Storfjorden, Svalbard, reveal large interannual ice-thickness variations in a region which is typically characterized by a reoccurring polynya. The surveys were performed in March 2003, May 2006 and March 2007 with helicopter- and ship-based sensors. The thickness distributions are influenced by sea-ice and atmospheric boundary conditions 2 months prior to the surveys, which are assessed with synthetic aperture radar (SAR) images, regional QuikSCAT backscatter maps and wind information from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset. Locally formed thin ice from the Storfjorden polynya was frequently observed in 2003 and 2007 (mean thickness 0.55 and 0.37 m, respectively) because these years were characterized by prevailing northeasterly winds. In contrast, the entire fjord was covered with thick external sea ice in 2006 (mean thickness 2.21 m), when ice from the Barents Sea was driven into the fjord by predominantly southerly winds. The modal thickness of this external ice in 2006 increased from 1.2 m in the northern fjord to 2.4 m in the southern fjord, indicating stronger deformation in the southern part. This dynamically thickened ice was even thicker than multi-year ice advected from the central Arctic Ocean in 2003 (mean thickness 1.83 m). The thermodynamic ice thickness of fast ice as boundary condition is investigated with a one-dimensional sea-ice growth model (1DICE) forced with meteorological data from the weather station at the island of Hopen, southeast of Storfjorden. The model results are in good agreement with the modal thicknesses of fast-ice measurements in all years.

Cenozoic global ice volume and temperature simulations over the past 40 million years

Variations in global ice volume and temperature over the Cenozoic era have been investigated with a set of one-dimensional (1-D) ice-sheet models. Simulations include three ice sheets representing glaciation in the Northern Hemisphere, i.e. in Eurasia, North America and Greenland, and two separate ice sheets for Antarctic glaciation. The continental mean Northern Hemisphere surface-air temperature has been derived through an inverse procedure from observed benthic d18O records. These data have yielded a mutually consistent and continuous record of temperature, global ice volume and benthic d18O over the past 35 Ma. The simple 1-D model shows good agreement with a comprehensive 3-D ice-sheet model for the past 3 Ma. On average, differences are only 1.0°C for temperature and 6.2 m for sea level. Most notably, over the 35 Ma period, the reconstructed ice volume-temperature sensitivity shows a transition from a climate controlled by Southern Hemisphere ice sheets to one controlled by Northern Hemisphere ice sheets. Although the transient behaviour is important, equilibrium experiments show that the relationship between temperature and sea level is linear and symmetric, providing limited evidence for hysteresis. Furthermore, the results show a good comparison with other simulations of Antarctic ice volume and observed sea level.