Rare earth element contents in deposits and minerals of the ocean floor

Current understanding of rare earth element (REE) geochemistry in the ocean is given in the book. Chemical properties determining REE migration ability in natural processes, sources of REE in the ocean, behavior of REE in river-sea mixing zones, fractionation of dissolved and particulate REE in ocean waters under aerobic and anaerobic conditions, distribution of REE in terrigenous, authigenic, hydrothermal and biogenic sediment components (clay, bone detritus, barite, phillipsite, Fe- and Mn-oxyhydroxides, Fe-Ca hydroxophosphate, diatoms and foraminiferas) are under consideration.

(Appendix) Recent surface areas and historical lengths of glaciers from Patagonia and Tierra del Fuego, South America

Fluctuations in the length of 72 glaciers in the Northern and Southern Patagonia Icefield (NPI and SPI, respectively) and the Cordillera Darwin Icefield (CDI) were estimated between 1945 and 2005. The information obtained from historical maps based on 1945 aerial photographs was compared to ASTER and Landsat satellite images and to information found in the literature. The majority of glaciers have retreated considerably, with maximum values of 12.2 km for Marinelli Glacier in the CDI, 11.6 km for O'Higgins Glacier in the SPI and 5.7 km for San Rafael Glacier in the NPI. Among the 20 glaciers that have retreated the most relative to their size, small (less than 50 km**2) and medium (between 50 and 200 km**2) glaciers are the most affected. However, no direct relation between glacier retreat and size was found for the 72 glaciers studied. The highest percentage retreat in the CDI was by the CDI-03 Glacier (37.9%) and Marinelli Glacier (37.6%). In the SPI, relative retreats were heterogeneous and fluctuated between 27.2% (Amelia Glacier) and 0.4% (Viedma Glacier). In the NPI, relative retreat was very high for Strindberg and Cachet glaciers (35.9% and 27.6%, respectively) but for the remaining glaciers in this icefield it ranged between 11.8% (Piscis Glacier) and 3.6% (San Quintin Glacier). In addition to surface area, the surface slope (calculated on the basis of the DEM SRTM) was also related to the relative retreat and no straightforward relation was found. From a global point of view, we suggest that glacier retreat in the region is controlled firstly by atmospheric warming, as it has been reported in this area. Besides the general increase in temperature observed, no signal of a geographical pattern for the fluctuations in glacier length was found. Consequently, glaciers appear to initially react to local conditions most probably induced by their exposition, geometry and hypsometry. The heterogeneity of rates of retreat suggests that differences in basin geometry, glacier dynamics and response time are key features to explain fluctuations of each glacier.

Temperature observation of CORK installation during ODP Leg 168

Four boreholes, drilled a few tens of meters into igneous basement on the eastern flank of the Juan de Fuca Ridge during ODP Leg 168, were sealed and instrumented for long-term monitoring to observe the hydrologic state of young sediment-sealed oceanic crust. The thermal regime is dominated by the effects of rapid fluid circulation in uppermost igneous basement driven by very small non-hydrostatic pressure gradients. Upper basement temperatures are uniform laterally between pairs of holes over distances of hundreds of meters to kilometers. In the case of two holes drilled into a sediment-buried basement ridge and adjacent valley, basement temperatures differ by less than 2 K despite the 2.2 km lateral separation of the sites and the 2.5:1 contrast in sediment cover thickness. Under conductive conditions, upper basement temperatures would differ by roughly 50 K. By comparison with modeling results, the observed degree of isothermality suggests a fluid flux of at least 10**-6 m/s (30 m/yr), and an effective permeability in the range of 10**-10 to 10**-9 m**2 in the uppermost igneous crust. The pressure difference available to drive this rapid flux between the ridge and valley, estimated by comparing the observed pressures via the isothermal upper basement hydrostat that is inferred to connect the two sites, is small (~2 kPa) and also suggests high permeability. Relative to the hydrostats defined by the local conductive sediment geotherms, substantial super-hydrostatic pressure (+18 kPa) is present within the buried basement ridge, and sub-hydrostatic pressure is present in the adjacent valley (-26 kPa). Such pressure differentials are the direct consequence of the advection-dominated thermal regime and small pressure losses in high-permeability basement, and are available to drive fluid seepage through sediment sections vertically up above and horizontally away from buried ridges, and down above valleys. No constraints are provided by any of the observations on the depth in the crust to which thermally or chemically significant flow might extend, although just as in the overlying sediments, the pattern of deep flow may be affected by the near-isothermal and near-hydrostatic conditions present in the permeable uppermost crustal section.

Electrical resistivity, formation factors, and sound velocity at DSDP Holes 75-530A and 75-530B

From 0 to 277 m at Site 530 are found Holocene to Miocene diatom ooze, nannofossil ooze, marl, clay, and debrisflow deposits; from 277 to 467 m are Miocene to Oligocene mud; from 467 to 1103 m are Eocene to late Albian Cenomanian interbedded mudstone, marlstone, chalk, clastic limestone, sandstone, and black shale in the lower portion; from 1103 to 1121 m are basalts. In the interval from 0 to 467 m, in Holocene to Oligocene pelagic oozes, marl, clay, debris flows, and mud, velocities are 1.5 to 1.8 km/s; below 200 m velocities increase irregularly with increasing depth. From 0 to 100 m, in Holocene to Pleistocene diatom and nannofossil oozes (excluding debris flows), velocities are approximately equivalent to that of the interstitial seawater, and thus acoustic reflections in the upper 100 m are primarily caused by variations in density and porosity. Below 100 or 200 m, acoustic reflections are caused by variations in both velocity and density. From 100 to 467 m, in Miocene-Oligocene nannofossil ooze, clay, marl, debris flows, and mud, acoustic anisotropy irregularly increases to 10%, with 2 to 5% being typical. From 467 to 1103 m in Paleocene to late Albian Cenomanian interbedded mudstone, marlstone, chalk, clastic limestone, and black shale in the lower portion of the hole, velocities range from 1.6 to 5.48 km/s, and acoustic anisotropies are as great as 47% (1.0 km/s) faster horizontally. Mudstone and uncemented sandstone have anisotropies which irregularly increase with increasing depth from 5 to 10% (0.2 km/s). Calcareous mudstones have the greatest anisotropies, typically 35% (0.6 km/s). Below 1103 m, basalt velocities ranged from 4.68 to 4.98 km/s. A typical value is about 4.8 km/s. In situ velocities are calculated from velocity data obtained in the laboratory. These are corrected for in situ temperature, hydrostatic pressure, and porosity rebound (expansion when the overburden pressure is released). These corrections do not include rigidity variations caused by overburden pressures. These corrections affect semiconsolidated sedimentary rocks the most (up to 0.25 km/s faster). These laboratory velocities appear to be greater than the velocities from the sonic log. Reflection coefficients derived from the laboratory data, in general, agree with the major features on the seismic profiles. These indicate more potential reflectors than indicated from the reflection coefficients derived using the Gearhart-Owen Sonic Log from 625 to 940 m, because the Sonic Log data average thin beds. Porosity-density data versus depth for mud, mudstone, and pelagic oozes agree with data for similar sediments as summarized in Hamilton (1976). At depths of about 400 m and about 850 m are zones of relatively higher porosity mudstones, which may suggest anomalously high pore pressure; however, they are more probably caused by variations in grain-size distribution and lithology. Electrical resistivity (horizontal) from 625 to 950 m ranged from about 1.0 to 4.0 ohm-m, in Maestrichtian to Santonian- Coniacian mudstone, marlstone, chalk, clastic limestone, and sandstone. An interstitial-water resistivity curve did not indicate any unexpected lithology or unusual fluid or gas in the pores of the rock. These logs were above the black shale beds. From 0 to 100 m at Sites 530 and 532, the vane shear strength on undisturbed samples of Holocene-Pleistocene diatom and nannofossil ooze uniformly increases from about 80 g/cm**2 to about 800 g/cm**2. From 100 to 300 m, vane shear strength of Pleistocene-Miocene nannofossil ooze, clay, and marl are irregular versus depth with a range of 500 to 2300 g/cm**2; and at Site 532 the vane shear strength appears to decrease irregularly and slightly with increasing depth (gassy zone). Vane shear strength values of gassy samples may not be valid, for the samples may be disturbed as gas evolves, and the sediments may not be gassy at in situ depths.

Boron contents and boron, carbon, and oxygen isotope data for ODP Sites 160-970 and 160-971

Products of two mud volcanoes from the distal part of the Mediterranean Ridge accretionary complex have been investigated regarding their B, C, and O stable isotope signatures. The mud breccias have been divided into mud matrix, lithified clasts, biogenic deposits, and authigenic cements and crusts related to fluid flow and cementation. Isotope geochemistry is used to evaluate the depth of mobilization of each phase in the subduction zone. B contents and isotope ratios of the mud and mud clasts show a general trend of B enrichment and decreasing d11B values with increasing consolidation (i.e., depth). However, the majority of the clast and matrix samples relate to moderate depths of mobilization within the wedge (1-2 km below seafloor). The carbonate cements of most of these clasts as well as the authigenic crusts, however, provide evidence for a deep fluid influence, probably associated with the décollement at 5-6 km depth. This interpretation is supported by d13C ratios of the crust, which indicate precipitation of C from thermogenic methane, and by the d11B ratios of pore-water samples of mud-breccia drill cores. Clams (Vesicomya sp.) living adjacent to fluid vents have d11B and d18O values corresponding to brines known in the area, which acted as the parent solution for shell precipitation. Such brines are most likely Miocene pore waters trapped at deep levels within the backstop to the accretionary prism, probably prior to desiccation of the Mediterranean in the Messinian (6-5 Ma). Combining all results, deep fluid circulation and expulsion are identified as the main processes triggering mud liquefaction and extrusion, whereas brines contribute only locally. Given the high B contents, mud extrusion has to be considered a major backflux mechanism of B into the hydrosphere.

Geochemistry of DSDP Hole 86-577B samples

THE chemical, mineralogical and isotopic characteristics of deposits at the Cretaceous/Tertiary (K/T) boundary are suggestive of a large impact event, the prime candidate (Sharpton et al., 1992, doi:10.1038/359819a0) being the Chicxulub crater in Yucatan, Mexico. Spinel-bearing spherules, which may be associated with such impacts, have been reported (Smit and Romein, 1985, doi:10.1016/0012-821X(85)90019-6) at several K/T boundary sites worldwide, but their origin is still uncertain. We have examined the spinel-bearing material recovered from K/T boundary deposits at site 577 in the Pacific Ocean (Heath et al., 1985, doi:10.2973/dsdp.proc.86.104.1985) and find two distinct populations of particles: spherules with dendritic spinel textures dispersed throughout the grains and irregularly shaped fragments with spinels essentially confined to the rim. The morphology and composition of the particles are characteristic of melted and partially melted meteoritic ablation debris, but their location is difficult to reconcile with an impact on the Yucatan peninsula, some 10,000 km away. We suggest instead that the spinel-bearing particles at site 577 are derived from the impact of a 2-km asteroid in the Pacific Ocean, and that several accretionary events of this type are required to explain the global distribution of spinel-bearing spherules at the K/T boundary.

Geochemistry of vein rocks of ODP Hole 118-735B

Rock samples from Hole 735B, Southwest Indian Ridge, were examined to determine the principal vein-related types of alteration that occurred, the nature of fluids that were present, and the temperatures and pressures of these fluids. Samples studied included veined metagabbro, veined mylonitic metagabbro, felsic trondhjemite, and late-stage leucocratic diopside-bearing veins. The methods used were standard petrographic analysis, mineral chemical analysis by electron microprobe, fluid inclusion petrography and analysis by heating/freezing techniques and laser Raman microspectroscopy, and oxygen isotopic analyses of mineral separates. Alteration in lithologic Units I and II (above the level of Core 118-735B-3OR; approximately 140 meters below the seafloor) is dominated by hydration by seawater-derived fluids at high temperature, up to about 700°C, and low water/rock ratio, during and immediately after pervasive ductile deformation. Below Core 118-735B-30R, pervasive deformation is less common, and brittle veining and brecciation are the major alteration styles. Leucocratic centimeter-scale veins, often containing diopside and plagioclase, were produced by interaction of hot (about 500°C) seawater-derived fluid and gabbro. The water/rock ratio was locally high at the veins and breccia zones, but the integrated water/rock ratio for the lower part of the hole is probably low. Accessory hydrous magmatic or deuteric phases formed from magmatic volatiles in some gabbro and in trondhjemite. Most subsequent alteration was affected by fluids that were seawater-derived, based on isotopic and chemical analyses of minerals and analyses of fluid inclusions. Many early-generation fluid inclusions, associated with high-temperature veining, contain appreciable methane as well as saline water. The source of methane is unclear, but it may have formed as seawater was reduced during low water/rock interaction with ultramafic upper mantle or ultramafic and mafic layer 3. Temperatures of alteration were calculated on the basis of coexisting mineral chemistry and isotopic values. Hydrothermal metamorphism commenced at about 720°C and continued to about 550°C. Leucocratic veining took place at about 500°C. Alteration within brecciated horizons was also at about 500° to less than 400°C, and the trondhjemite was altered at about 550° to below 490°C. Pressures calculated from a diopside-bearing vein, based on a combination of fluid inclusion and isotopic analysis, were 90 to 100 MPa. This pressure places the sample, from Core 118-735B-70R in Unit V, at about 2 km below the seafloor.

Oxygen-isotope composition, temperature and depth of formation of chert at DSDP Leg 62 Holes

We measured oxygen-isotope compositions of 16 siliceous rocks from Deep Sea Drilling Project Sites 463, 464, 465, and 466 (Leg 62). Samples are from deposits that range in age from about 40 to 103 m.y. and that occur at sub-bottom depths of 9 to 461 meters. Mean d18O values range from 28.4 to 36.8 per mil and 36.0 ± 0.3 per mil for quartz-rich and opal-CTrich rocks, respectively. d18O values in chert decrease with increasing sub-bottom depth; the slope of the d18O/depth curve is less steep for Site 464 than for the other sites which indicates that chert at Site 464 formed at higher temperatures than chert at Sites 463, 465, and 466. Temperatures of formation of cherts were 7 to 42°C, using the silica-water fractionation factor of Knauth and Epstein (1976), or 19 to 56°C, using the equation of Clayton et al. (1972). Temperatures in the sediment where the cherts now occur are lower than their isotopically determined temperatures of formation, which means that the cherts record an earlier history when temperatures in the sediment section were greater. Estimated sediment temperatures when the cherts formed are comparable to, but generally slightly lower than, those calculated from Knauth and Epstein's equation. The isotopic composition of cherts is more closely related to environment of formation (diagenetic environment) or paleogeothermal gradients, than to paleoclimates (bottom-water temperatures). Opal-CT-rich rocks may better record paleo-bottom-water temperature. In Leg 62 cherts, better crystallinity of quartz corresponds to lower d18O values; this implies progressively higher temperatures of equilibration between quartz and water during maturation of quartz. The interrelationship of d18O and crystallinity is noted also in continental-margin deposits such as the Monterey Formation - but for higher temperatures. The apparent temperature difference between open-ocean and continental-margin deposits can be explained by the dominant control of temperature on silica transformation in the rapidly deposited continental-margin deposits, whereas time, as well as temperature, has a strong influence on the transformations in open-ocean deposits. Comparisons between the chemistry and d18O values of cherts reveal two apparent trends: both boron and SiO2 increase as d18O increases. However, the correspondence between SiO2 and d18O is only apparent, because the two cherts lowest in SiO2 are also the most deeply buried, so the trend actually reflects depth of burial. The correspondence between boron and d18O supports the conclusion that boron is incorporated in the quartz crystal structure during precipitation

Homogenized radiosonde record at station Ny-Ålesund, Spitsbergen, 1993-2014

Radiosonde measurements obtained at the Arctic site Ny-Ålesund (78.9° N, 11.9° E), Svalbard, from 1993 to 2014 have been homogenized accounting for instrumentation discontinuities by correcting known errors in the manufacturer provided profiles. From the homogenized data record, the first Ny-Ålesund upper-air climatology of wind, temperature and humidity is presented, forming the background for the analysis of changes during the 22-year period. Particularly during the winter season, a strong increase in atmospheric temperature and humidity is observed, with a significant warming of the free troposphere in January and February up to 3 K per decade. This winter warming is even more pronounced in the boundary layer below 1 km, presumably amplified by mesoscale processes including e.g. orographic effects or the boundary layer capping inversion. Though the largest contribution to the increasing atmospheric water vapour column in winter originates from the lowermost 2 km, no increase in the contribution by specific humidity inversions is detected. Instead, we find an increase in the humidity content of the large scale background humidity profiles. At the same time, the tropospheric flow in winter is found to occur less frequent from northerly directions and to the same amount more frequent from the South. We conclude that changes in the atmospheric circulation lead to an enhanced advection of warm and moist air from lower latitudes to the Svalbard region in the winter season, causing the warming and moistening of the atmospheric column above Ny-Ålesund.

Ice and snow characteristics in Prydz Bay in 2006

Under the Chinese National Antarctic Research Expedition program in 2006, the annual thermal mass balance of landfast ice in the vicinity of Zhongshan Station, Prydz Bay, east Antarctica, was investigated. Sea ice formed from mid-February onward, and maximum ice thickness occurred in late November. Snow cover remained thin, and blowing snow caused frequent redistribution of the snow. The vertical ice salinity showed a 'question-mark-shaped' profile for most of the ice growth season, which only turned into an 'I-shaped' profile after the onset of ice melt. The oceanic heat flux as estimated from a flux balance at ice-ocean interface using internal ice temperatures decreased from 11.8 (±3.5) W/m**2 in April to an annual minimum of 1.9 (±2.4) W/m**2 in September. It remained low through late November, in mid-December it increased sharply to about 20.0 W/m**2. Simulations applying the modified versions of Stefan's law, taking account the oceanic heat flux and ice-atmosphere coupling, compare well with observed ice growth. There was no obvious seasonal cycle for the thermal conductivity of snow cover, which was also derived from internal ice temperatures. Its annual mean was 0.20 (±0.04) W/m/°C.

Pore-water chemistry of ODP sites 111-677 and 111-678

Sites 677 and 678 were drilled on ODP Leg 111 to test hypotheses about the nature and pattern of hydrothermal circulation on a mid-ocean ridge flank. Together with earlier results from DSDP Site 501/504 and several heatflow and piston coring surveys covering a 100-km**2 area surrounding the three drill sites, they confirm that hydrothermal circulation persists in this 5.9-m.y.-old crust, both in basement and through the overlying sediments (Langseth et al., 1988, doi:10.2973/odp.proc.ir.111.102.1988). Profiles of sediment pore-water composition with depth at the three drill sites show both vertical and horizontal gradients. The shapes of the profiles and their variation from one site to another result from a combination of vertical and horizontal diffusion, convection, and reaction in the sediments and basement. Chemical species that are highly reactive in the siliceous-calcareous biogenic sediments include bicarbonate (alkalinity), ammonium, sulfate, manganese, calcium, strontium, lithium, silica, and possibly potassium. Reactions include bacterial sulfate reduction, mobilization of Mn2+, precipitation of CaCO3, and recrystallization of calcareous and siliceous oozes to chalk, limestone, and chert. Species with profiles more affected by reaction in basaltic basement than in the sediments include Mg, Ca, Na, K, and oxygen isotopes. Reaction in basement at 60?C and at higher temperatures has produced a highly altered basement formation water that is uniform in composition over distances of several kilometers. As inferred from the composition of the basal sediment pore water at the three sites, this uniformity extends from up flow zone to downflow zone in basement and the sediments. It exists in spite of large variations in heat flow and depth to basement, apparently as a result of homogenization by hydrothermal circulation in basement. Profiles for chlorinity, Na, Mg, and other species in the sediment pore waters confirm that Site 678, drilled on a localized heatflow high identified by Langseth et al. (1988), is a site of long-lived upwelling of warm water from basement through the sediments at velocities of 1 to 2 mm/yr. The upflow through the anomalously thin sediments is apparently localized above an uplifted fault block in basement. This site and other similar sites in the survey area give rise to lateral diffusion and possibly flow through the sediments, which produces lateral gradients in sediment pore-water composition at sites such as 501/504. The complementary pore-water profiles at the low-heatflow Site 677 2 km to the south indicate that downflow is occurring through the sediments there, at comparable rates of 1 to 2 mm/yr.

Topographic model of Greenland (elevation lines in PanMap and ASCII format)

In this paper, a new high-resolution elevation model of Greenland, including the ice sheet as well as the ice free regions, is presented. It is the first published full coverage model, computed with an average resolution of 2 km and providing an unprecedented degree of detail. The topography is modeled from a wide selection of data sources, including satellite radar altimetry from Geosat and ERS 1, airborne radar altimetry and airborne laser altimetry over the ice sheet, and photogrammetric and manual map scannings in the ice free region. The ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations. The mean accuracy of the ice sheet elevations is estimated to be 12-13 m, and it is found that on surfaces of a slope between 0.2° and 0.8°, corresponding to approximately 50% of the ice sheet, the model presents a 40% improvement over models based on satellite altimetry alone. On coastal bedrock, the model is compared with stereo triangulated reference points, and it is found that the model accuracy is of the order of 25-35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere.