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.

Pore pressure measurements on IODP Exp308 sites

Overpressures measured with pore pressure penetrometers during Integrated Ocean Drilling Program (IODP) Expedition 308 reach 70% and 60% of the hydrostatic effective stress (View the MathML source) in the first 200 meters below sea floor (mbsf) at Sites U1322 and U1324, respectively, in the deepwater Gulf of Mexico, offshore Louisiana. High overpressures are present within low permeability mudstones where there have been multiple, very large, submarine landslides during the Pleistocene. Beneath 200 mbsf at Site U1324, pore pressures drop significantly: there are no submarine landslides in this mixture of mudstone, siltstone, and sandstone. The penetrometer measurements did not reach the in situ pressure at the end of the deployment. We used a soil model to determine that an extrapolation approach based on the inverse of square route of time (View the MathML source) requires much less decay time to achieve a desirable accuracy than an inverse time (1/t) extrapolation. Expedition 308 examined how rapid and asymmetric sedimentation above a permeable aquifer drives lateral fluid flow, extreme pore pressures, and submarine landslides. We interpret that the high overpressures observed are driven by rapid sedimentation of low permeability material from the ancestral Mississippi River. Reduced overpressure at depth at Site U1324 suggests lateral flow (drainage) whereas high overpressure at Site U1322 requires inflow from below: lateral flow in the underlying permeable aquifer provides one mechanism for these observations. High overpressure near the seafloor reduces slope stability and provides a mechanism for the large submarine landslides and low regional gradient (2°) offshore from the Mississippi delta.

Trace element geochemistry of zircons frim in situ ocean lithoshere

We characterize the textural and geochemical features of ocean crustal zircon recovered from plagiogranite, evolved gabbro, and metamorphosed ultramafic host-rocks collected along present-day slow and ultraslow spreading mid-ocean ridges (MORs). The geochemistry of 267 zircon grains was measured by sensitive high-resolution ion microprobe-reverse geometry at the USGS-Stanford Ion Microprobe facility. Three types of zircon are recognized based on texture and geochemistry. Most ocean crustal zircons resemble young magmatic zircon from other crustal settings, occurring as pristine, colorless euhedral (Type 1) or subhedral to anhedral (Type 2) grains. In these grains, Hf and most trace elements vary systematically with Ti, typically becoming enriched with falling Ti-in-zircon temperature. Ti-in-zircon temperatures range from 1,040 to 660°C (corrected for a TiO2 ~ 0.7, a SiO2 ~ 1.0, pressure ~ 2 kbar); intra-sample variation is typically ~60-15°C. Decreasing Ti correlates with enrichment in Hf to ~2 wt%, while additional Hf-enrichment occurs at relatively constant temperature. Trends between Ti and U, Y, REE, and Eu/Eu* exhibit a similar inflection, which may denote the onset of eutectic crystallization; the inflection is well-defined by zircons from plagiogranite and implies solidus temperatures of ~680-740°C. A third type of zircon is defined as being porous and colored with chaotic CL zoning, and occurs in ~25% of rock samples studied. These features, along with high measured La, Cl, S, Ca, and Fe, and low (Sm/La)N ratios are suggestive of interaction with aqueous fluids. Non-porous, luminescent CL overgrowth rims on porous grains record uniform temperatures averaging 615 ± 26°C (2SD, n = 7), implying zircon formation below the wet-granite solidus and under water-saturated conditions. Zircon geochemistry reflects, in part, source region; elevated HREE coupled with low U concentrations allow effective discrimination of ~80% of zircon formed at modern MORs from zircon in continental crust. The geochemistry and textural observations reported here serve as an important database for comparison with detrital, xenocrystic, and metamorphosed mafic rock-hosted zircon populations to evaluate provenance.

Present day and paleo-geographic coordinates, thicknesses and isopach area of the Early Eocene +19 ash layer (Table 1)

23 layers of altered volcanic ash (bentonites) originating from the North Atlantic Igneous Province have been recorded in early Eocene deposits of the Austrian Alps, about 1,900 km away from the source area. The Austrian bentonites are distal equivalents of the ''main ash-phase'' in Denmark and the North Sea basin. We have calculated the total eruption volume of this series as 21,000 km**3, which occurred in 600,000 years. The most powerful single eruption of this series took place 54.0 million years ago (Ma) and ejected ca. 1,200 km**3 of ash material, which makes it one of the largest basaltic pyroclastic eruptions in geological history. The clustering of eruptions must have significantly affected the incoming solar radiation in the early Eocene by the continuous production of stratospheric dust and aerosol clouds. This hypothesis is corroborated by oxygen isotope values, which indicate a global decrease of sea surface temperatures between 1 and 2 C during this major phase of explosive volcanism.

Processed 2 minutes-averaged continuous VM-ADCP (vessel-mounted Acoustic Doppler Current Profiler) profiles during POLARSTERN cruise ANT-XXVIII/3

This data set was obtained during the R. V. POLARSTERN cruise ANT-XXVIII/3. Current velocities were measured nearly continuously when outside territorial waters along the ship's track with a vessel-mounted TRD Instruments' 153.6-kHz Ocean Surveyor ADCP. The transducers were located 11 m below the water line and were protected against ice floes by an acoustically transparent plastic window. The current measurements were made using a pulse of 2s and vertical bin length of 4 m. The ship's velocity was calculated from position fixes obtained by the Global Positioning System (GPS). Heading, roll and pitch data from the ship's gyro platforms and the navigation data were used to convert the ADCP velocities into earth coordinates. Accuracy of the ADCP velocities mainly depends on the quality of the position fixes and the ship's heading data. Further errors stem from a misalignment of the transducer with the ship's centerline. The ADCP data were processed using the Ocean Surveyor Sputum Interpreter (OSSI) software developed by GEOMAR Helmholtz-Zentrum für Ozeanforschung Kiel. The averaging interval was set to 120 seconds. The reference layer was set to bins 5 to 16 avoiding near surface effects and biases near bin 1. Sampling interval setting: 2s; Number of bins: 80; Bin length: 4m; Pulse length: 4m; Blank beyond transmit length: 4m. Data processing setting: Top reference bin: 5; Bottom reference bin: 16; Average: 120s; Misalignment amplitude: 1.0276 +/- 0.1611, phase: 0.8100 +/- 0.7190. The precision for single ping and 4m cell size reported by TRDI is 0.30m/s. Resulting from the single ping precision and the number of pings (most of the time 36) during 120seconds the velocity accuracy is nearly 0.05m/s. (Velocity accuracy = single ping precision divided by square root of the number of pings).

Boron isotope and B/Ca ratios of planktonic foraminifera

Here a new analytical methodology is described for measuring the isotopic composition of boron in foraminifera using multicollector inductively coupled plasma mass spectrometry (MC-ICPMS). This new approach is fast (~10 samples analysed in duplicate per analytical session) and accurate (to better than 0.25 per mil at 95% confidence) with acceptable sample size requirements (1-3 mg of carbonate). A core top calibration of several common planktic and two benthic species from geographically widespread localities shows a very close agreement between the isotopic composition measured by MC-ICPMS and the isotopic composition of B(OH)-4 in seawater (as predicted using the recently measured isotopic equilibrium factor of 1.0272) at the depth of habitat. A down core and core top investigation of boron concentration (B/Ca ratio) shows that the partition coefficient is influenced by [CO2-3] complicating the application of this proxy. Nevertheless, it is demonstrated that these two proxies can be used to fully constrain the carbonate system of surface water in the Caribbean Sea (ODP Site 999A) over the last 130 kyr. This reconstruction shows that during much of the Holocene and the last interglacial period surface water at Site 999A was in equilibrium with the atmosphere with respect to CO2. During the intervening colder periods although the surface water pCO2 was lower than the Holocene, it was a minor to significant source of CO2 to the atmosphere possibly due to either an expansion of the eastern equatorial Atlantic upwelling zone, or a more local expansion of coastal upwelling in the southern Caribbean. Such reorganisation of the oceanic carbonate system in favour of a larger source of CO2 to the atmosphere from the equatorial ocean may require mechanisms responsible for lowering atmospheric CO2 during glacial periods to be more efficient than previously supposed.

ERA-interim reanalysis debiased at FLUXNET sites

(preliminary) Exchanges of carbon, water and energy between the land surface and the atmosphere are monitored by eddy covariance technique at the ecosystem level. Currently, the FLUXNET database contains more than 500 sites registered and up to 250 of them sharing data (Free Fair Use dataset). Many modelling groups use the FLUXNET dataset for evaluating ecosystem model's performances but it requires uninterrupted time series for the meteorological variables used as input. Because original in-situ data often contain gaps, from very short (few hours) up to relatively long (some months), we develop a new and robust method for filling the gaps in meteorological data measured at site level. Our approach has the benefit of making use of continuous data available globally (ERA-interim) and high temporal resolution spanning from 1989 to today. These data are however not measured at site level and for this reason a method to downscale and correct the ERA-interim data is needed. We apply this method on the level 4 data (L4) from the LaThuile collection, freely available after registration under a Fair-Use policy. The performances of the developed method vary across sites and are also function of the meteorological variable. On average overall sites, the bias correction leads to cancel from 10% to 36% of the initial mismatch between in-situ and ERA-interim data, depending of the meteorological variable considered. In comparison to the internal variability of the in-situ data, the root mean square error (RMSE) between the in-situ data and the un-biased ERA-I data remains relatively large (on average overall sites, from 27% to 76% of the standard deviation of in-situ data, depending of the meteorological variable considered). The performance of the method remains low for the Wind Speed field, in particular regarding its capacity to conserve a standard deviation similar to the one measured at FLUXNET stations.

Impact studies from different factors on root development from Stevia rebaudiana bertoni

The demand for an alternative and a high potency sweetener to substitute sugar increases year in year out, more so as a high percentage of the world population becomes increasingly diabetic. The alternative natural sweetener at hand has been Stevia rebaudiana Bertoni, a plant species, native to Paraguay and a member of the family compositae. Stevia is usually propagated by stem cuttings due to low percentage (10 %) seed germination, thus limiting large scales cultivation. To cultivate this crop en mass therefore, there is need to evolve efficient rooting techniques. Influences of irradiation from light, and hormones on rooting have been reported. The rooting efficacy in stem cuttings of this crop under varying light wavelengths, dark and hormone factors was investigated. Evaluated parameters include- (i) day of root emergent, (ii) percentage of rooted cuttings, (iii) average number, (iv) length and (v) width, of roots. Analysis of variance at p<.05 revealed that the number, length and width, of roots differed significantly in each case at p<0.000. Light irradiation was highly effective and a necessary factor for rooting in stems cuttings of this crop. The red light-IBA combined factors served best in stem micro-cutting practice and facilitation of effective mass cultivation in stevia crop.

(Table 3) Air temperatures, and snow and ice characteristics of fresh and brackish water lakes in the Northwest Territories, Canada

The algorithms designed to estimate snow water equivalent (SWE) using passive microwave measurements falter in lake-rich high-latitude environments due to the emission properties of ice covered lakes on low frequency measurements. Microwave emission models have been used to simulate brightness temperatures (Tbs) for snowpack characteristics in terrestrial environments but cannot be applied to snow on lakes because of the differing subsurface emissivities and scattering matrices present in ice. This paper examines the performance of a modified version of the Helsinki University of Technology (HUT) snow emission model that incorporates microwave emission from lake ice and sub-ice water. Inputs to the HUT model include measurements collected over brackish and freshwater lakes north of Inuvik, Northwest Territories, Canada in April 2008, consisting of snowpack (depth, density, and snow water equivalent) and lake ice (thickness and ice type). Coincident airborne radiometer measurements at a resolution of 80x100 m were used as ground-truth to evaluate the simulations. The results indicate that subsurface media are simulated best when utilizing a modeled effective grain size and a 1 mm RMS surface roughness at the ice/water interface compared to using measured grain size and a flat Fresnel reflective surface as input. Simulations at 37 GHz (vertical polarization) produce the best results compared to airborne Tbs, with a Root Mean Square Error (RMSE) of 6.2 K and 7.9 K, as well as Mean Bias Errors (MBEs) of -8.4 K and -8.8 K for brackish and freshwater sites respectively. Freshwater simulations at 6.9 and 19 GHz H exhibited low RMSE (10.53 and 6.15 K respectively) and MBE (-5.37 and 8.36 K respectively) but did not accurately simulate Tb variability (R= -0.15 and 0.01 respectively). Over brackish water, 6.9 GHz simulations had poor agreement with airborne Tbs, while 19 GHz V exhibited a low RMSE (6.15 K), MBE (-4.52 K) and improved relative agreement to airborne measurements (R = 0.47). Salinity considerations reduced 6.9 GHz errors substantially, with a drop in RMSE from 51.48 K and 57.18 K for H and V polarizations respectively, to 26.2 K and 31.6 K, although Tb variability was not well simulated. With best results at 37 GHz, HUT simulations exhibit the potential to track Tb evolution, and therefore SWE through the winter season.

Mineraly of oceanic crust from the Kane Transform, Mid-Atlantic Ridge

Gabbroic cumulates drilled south of the Kane Transform Fault on the slow-spread Mid-Atlantic Ridge preserve up to three discrete magnetization components. Here we use absolute age constraints derived from the paleomagnetic data to develop a model for the magmatic construction of this section of the lower oceanic crust. By comparing the paleomagnetic data with mineral compositions, and based on thermal models of local reheating, we infer that magmas that began crystallizing in the upper mantle intruded into the lower oceanic crust and formed meter-scale sills. Some of these magmas were crystal-laden and the subsequent expulsion of interstitial liquid from them produced '"cumulus" sills. These small-scale magmatic injections took place over at least 210 000 years and at distances of ~3 km from the ridge axis and may have formed much of the lower crust. This model explains many of the complexities described in this area and can be used to help understand the general formation of oceanic crust at slow-spread ridges.

GPS raw data (control points and ground control points) from the Liguria Region, Borghetto Santo Spirito, Italy

Monitoring the impact of sea storms on coastal areas is fundamental to study beach evolution and the vulnerability of low-lying coasts to erosion and flooding. Modelling wave runup on a beach is possible, but it requires accurate topographic data and model tuning, that can be done comparing observed and modeled runup. In this study we collected aerial photos using an Unmanned Aerial Vehicle after two different swells on the same study area. We merged the point cloud obtained with photogrammetry with multibeam data, in order to obtain a complete beach topography. Then, on each set of rectified and georeferenced UAV orthophotos, we identified the maximum wave runup for both events recognizing the wet area left by the waves. We then used our topography and numerical models to simulate the wave runup and compare the model results to observed values during the two events. Our results highlight the potential of the methodology presented, which integrates UAV platforms, photogrammetry and Geographic Information Systems to provide faster and cheaper information on beach topography and geomorphology compared with traditional techniques without losing in accuracy. We use the results obtained from this technique as a topographic base for a model that calculates runup for the two swells. The observed and modeled runups are consistent, and open new directions for future research.

Diatom abundance in surface sediments of the Southern Ocean

Based on the quantitative analysis of diatom assemblages preserved in 274 surface sediment samples recovered in the Pacific, Atlantic and western Indian sectors of the Southern Ocean we have defined a new reference database for quantitative estimation of late-middle Pleistocene Antarctic sea ice fields using the transfer function technique. The Detrended Canonical Analysis (DCA) of the diatom data set points to a unimodal distribution of the diatom assemblages. Canonical Correspondence Analysis (CCA) indicates that winter sea ice (WSI) but also summer sea surface temperature (SSST) represent the most prominent environmental variables that control the spatial species distribution. To test the applicability of transfer functions for sea ice reconstruction in terms of concentration and occurrence probability we applied four different methods, the Imbrie and Kipp Method (IKM), the Modern Analog Technique (MAT), Weighted Averaging (WA), and Weighted Averaging Partial Least Squares (WAPLS), using logarithm-transformed diatom data and satellite-derived (1981-2010) sea ice data as a reference. The best performance for IKM results was obtained using a subset of 172 samples with 28 diatom taxa/taxa groups, quadratic regression and a three-factor model (IKM-D172/28/3q) resulting in root mean square errors of prediction (RMSEP) of 7.27% and 11.4% for WSI and summer sea ice (SSI) concentration, respectively. MAT estimates were calculated with different numbers of analogs (4, 6) using a 274-sample/28-taxa reference data set (MAT-D274/28/4an, -6an) resulting in RMSEP's ranging from 5.52% (4an) to 5.91% (6an) for WSI as well as 8.93% (4an) to 9.05% (6an) for SSI. WA and WAPLS performed less well with the D274 data set, compared to MAT, achieving WSI concentration RMSEP's of 9.91% with WA and 11.29% with WAPLS, recommending the use of IKM and MAT. The application of IKM and MAT to surface sediment data revealed strong relations to the satellite-derived winter and summer sea ice field. Sea ice reconstructions performed on an Atlantic- and a Pacific Southern Ocean sediment core, both documenting sea ice variability over the past 150,000 years (MIS 1 - MIS 6), resulted in similar glacial/interglacial trends of IKM and MAT-based sea-ice estimates. On the average, however, IKM estimates display smaller WSI and slightly higher SSI concentration and probability at lower variability in comparison with MAT. This pattern is a result of different estimation techniques with integration of WSI and SSI signals in one single factor assemblage by applying IKM and selecting specific single samples, thus keeping close to the original diatom database and included variability, by MAT. In contrast to the estimation of WSI, reconstructions of past SSI variability remains weaker. Combined with diatom-based estimates, the abundance and flux pattern of biogenic opal represents an additional indication for the WSI and SSI extent.

Reconstructions of eolian dust accumulation in northwest African margin sediments

Reconstructions of eolian dust accumulation in northwest African margin sediments provide important continuous records of past changes in atmospheric circulation and aridity in the region. Existing records indicate dramatic changes in North African dust emissions over the last 20 ka, but the limited spatial extent of these records and the lack of high-resolution flux data do not allow us to determine whether changes in dust deposition occurred with similar timing, magnitude and abruptness throughout northwest Africa. Here we present new records from a meridional transect of cores stretching from 31°N to 19°N along the northwest African margin. By combining grain size endmember modeling with 230Th-normalized fluxes for the first time, we are able to document spatial and temporal changes in dust deposition under the North African dust plume throughout the last 20 ka. Our results provide quantitative estimates of the magnitude of dust flux changes associated with Heinrich Stadial 1, the Younger Dryas, and the African Humid Period (AHP; ~11.7-5 ka), offering robust targets for model-based estimates of the climatic and biogeochemical impacts of past changes in North African dust emissions. Our data suggest that dust fluxes between 8 and 6 ka were a factor of ~5 lower than average fluxes during the last 2 ka. Using a simple model to estimate the effects of bioturbation on dust input signals, we find that our data are consistent with abrupt, synchronous changes in dust fluxes in all cores at the beginning and end of the AHP. The mean ages of these transitions are 11.8±0.2 ka (1Sigma) and 4.9±0.2 ka, respectively.