Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Composition of hydrothermal manganese oxide deposits from Galapagos mounds, DSDP Leg 70, Hole 509B

During DSDP Leg 70, a 1.60 m thick manganese oxide layer was sampled in hole 509B. This deposit is formed of alternating layers of hard plates of pure todorokite, about 2 mm thick, and of a more powdery material deeply impregnated with manganese oxide, about 3 mm thick. A SEM study of the plates and the associated powder shows that the powdery material is a transformation of a pre-existing sediment, while the plates are a direct precipitation from a hydrothermal solution. The uranium series disequilibrium method was used to determine the ages of the plates. They are found to be in good chronological sequence and in accordance with the sedimentation rate of the area (4.9 cm/10^3 years) which implies that they have been formed at the sediment-seawater interface during a pulsed injection of hydrothermal solution. The powder presents systematically an "older age" which is explained by a slowing down of the injection while the normal sediment settles; the older age is due to the 230Th excess of the sediment.

Mineralogy, chemistry and stable isotope composition of hydrothermal altered sheeted dikes of ODP Hole 111-504B

During ODP Leg 111 Hole 504B was extended 212 m deeper into the sheeted dikes of oceanic Layer 2, for a total penetration of 1288 m within basement. Study of the mineralogy, chemistry, and stable isotopic compositions of the rocks recovered on Leg 111 has confirmed and extended the previous model for hydrothermal alteration at the site: axial greenschist hydrothermal metamorphism was followed by seawater recharge and subsequent off-axis alteration. The dikes are depleted in 18O (mean delta18O = +5.1 ? +/- 0.6 ?) relative to fresh mid-ocean ridge basalt. Oxygen isotopic data on whole rocks and isolated secondary minerals indicate temperatures during axial metamorphism of 250°-350°C and water/rock ratios about one. Increasing amounts of actinolite with depth in the dike section, however, suggest that temperatures increased downward in the dikes. Pyrite + pyrrhotite + chalcopyrite + magnetite was the stable sulfide + oxide mineral assemblage during axial alteration, but these minerals partly re-equilibrated later at temperatures less than 200°C. The dikes sampled on Leg 111 contain an average of 500 ppm sulfur, slightly lower than igneous values. The delta34S values of sulfide average 0?, which indicates the presence of basaltic sulfide and incorporation of little or no seawater-derived sulfide into the rocks. These data are consistent with models for the presence of rock-dominated sulfur in deep hydrothermal fluids. The presence of anhydrite at 1176 m within basement indicates that unaltered seawater can penetrate to significant depths in the crust during recharge.

Composition and age of ores and bottom sediments accumulated within and near the Rainbow hydrothermal field

Results of direct geological and geochemical observations of the modern Rainbow hydrothermal field (Mid-Atlantic Ridge, 36°14'N; 33°54'W) carried out from the deep-sea manned Mir submersibles during Cruises 41 and 42 of the R/V Akademik Mstislav Keldysh in 1998-1999 and data of laboratory studies of collected samples are under consideration in the paper. The field lacks neovolcanic rocks and the axial part of the rift is filled in with a serpentinite protrusion. In this field there occur metalliferous sediments, as well as active and relict sulfide edifices composed of sulfide minerals; pyrrhotite, chalcopyrite, isocubanite, sphalerite, marcasite, pyrite, bornite, chalcosine, digenite, magnetite, anhydrite, rare troilite, wurtzite, millerite, and pentlandite have been determined. Sulfide ores are characterized by concentric-zoned textures. During in situ measurements during 35 minutes temperature of hydrothermal fluids was varying within a range from 250 to 350°C. Calculated chemical and isotopic composition of hydrothermal fluid shows elevated concentrations of Cl, Ni, Co, CH4, and H2. Values of d34S of H2S range from +2.4 to +3.1 per mil, of d13C of CH4 from -15.2 to -11.2 per mil, and d13C of CO2 from +1.0 to -4.0 per mil. Fluid inclusions are homogenized at temperatures from 140 to 360°C, whereas salinity of the fluid varies from 4.2 to 8.5 wt %. d34S values of sulfides range from +1.3 to +12.5 per mil. 3He/4He ratio in mineral-forming fluid contained in the fluid inclusions from sulfides of the Rainbow field varies from 0.00000374 to 0.0000101. It is shown that hydrothermal activity in the area continues approximately during 100 ka. It is assumed that the fluid and sulfide edifices contain components from the upper mantle. A hypothesis of phase separation of a supercritical fluid that results in formation of brines is proposed. Hydrothermal activity is related to the tectonic, not volcanic, phase of the Mid-Atlantic Ridge evolution.

Methane, helium isotope and hydrogen data from the Nibelungen hydrothermal area, Mid Atlantic Ridge

Hydrothermal emission of mantle helium appears to be directly related to magma production rate, but other processes can generate methane and hydrogen on mid-ocean ridges. In an on-going effort to characterize these processes in the South Atlantic, the flux and distribution of these gases were investigated in the vicinity of a powerful black smoker recently discovered at 8°17.9' S, 13°30.4' W. The vent lies on the shoulder of an oblique offset in the Mid-Atlantic Ridge and discharges high concentrations of methane and hydrogen. Measurements during expeditions in 2004 and 2006 show that the ratio of CH4 to 3He in the neutrally buoyant plume is quite high, 4 x 10**8. The CTD stations were accompanied by velocity measurements with lowered acoustic Doppler current profilers (LADCP), and from these data we estimate the methane transport to have been 0.5 mol/sec in a WSW-trending plume that seems to develop during the ebb tidal phase. This transport is an order of magnitude greater than the source of CH4 calculated from its concentration in the vent fluid and the rise height of the plume. From this range of methane fluxes, the source of 3He is estimated to be between 0.14 and 1.2 nmol/sec. In either case, the 3He source is significantly lower than expected from the spreading rate of the Mid-Atlantic Ridge. From the inventory of methane in the rift valley adjacent to the vent, it appears that the average specific rate of oxidation is 2.6 to 23/yr, corresponding to a turnover time between 140 and 16 days. Vertical profiles of methane in the surrounding region often exhibited Gaussian-like distributions, and the variances appear to increase with distance from the vent. Using a Gaussian plume model, we obtained a range of vertical eddy diffusivities between 0.009 and 0.08 m2m2/sec. These high values may be due to tidally driven internal waves across the promontory on which the vent is located.

(Table 3) Relative content of n-alkane groups in samples from hydrothermal fields of the Mid-Atlantic Ridge

Optimum conditions were selected for chromatographic separation of model mixtures of C12-C40 n-alkanes. For one of samples of hydrothermal deposits extraction conditions of hydrocarbons were studied and a sample preparation procedure was selected. The procedure was proposed to determine n-alkanes in samples of hydrothermal deposits by means of gas chromatography - mass spectrometry (GC-MS). Detection limit for n-alkanes was 3x10**-9 to 10**-8% depending on components. On the basis of the proposed procedure composition of n-alkanes was studied in samples of hydrothermal deposits collected at the Mid-Atlantic Ridge (Broken Spur, Lost City, and Rainbow hydrothermal fields). Analyses showed that samples contained C14-C35 n-alkanes. Concentrations of the n-alkanes were rather low and varied from 0.002 to 0.038 µg/g. Hypotheses concerning genesis of identified n-alkanes were offered.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.