Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

MRD's approach: Ideally suited to support Future Earth research agenda in mountains

Both Future Earth and Mountain Research and Development (MRD) aim to support production and dissemination of knowledge for sustainable development. As shown in Future Earth’s Strategic Research Agenda 2014, the global research community has begun to acknowledge its societal role and the need for a new type of research in which scientists link disciplines and coproduce transformation knowledge with stakeholders. Future Earth has defined three research themes that conceptualize the issues to be dealt with at the same time as the way in which this should be done. In many ways, MRD’s policy has made the journal a forerunner of Future Earth’s stipulated “step-change in research”. Indeed, MRD’s section policies aim to support similar contents and ways of producing these forms of knowledge. MRD publishes “systems knowledge” in its MountainResearch section, “target knowledge” in its MountainAgenda section, and “transformation knowledge” in its MountainDevelopment section. Each of these sections has dedicated review criteria to assess and enhance the quality of the knowledge presented in the papers. In this poster, we provide examples from each of the three sections of what the knowledge types look like, how they are assessed, and how they contribute to the three Future Earth themes.

Variations of clinopyroxene/melt element partitioning during assimilation of olivine/peridotite by low-Mg diorite magma

The effects of crystal chemistry and melt composition on the control of clinopyroxene/melt element partitioning (D) during the assimilation of olivine/peridotite by felsic magma have been investigated in Mesozoic high-Mg diorites from North China. The assimilation resulted in significant increase of Mg, Cr and Ni and only slight (< 30%) decrease of incompatible elements of the magma, and the compositional variations have been mirrored by the normally and reversely zoned clinopyroxene microphenocrysts formed at the early stage of the magma evolution. The Mg# [100 × Mg / (Mg + Fe)] values of the reversely zoned clinopyroxenes increase from 65 to 75 in the core to 85–90 in the high-Mg midsection, and reduce back to 73–79 at the rim. Trace element profiles across all these clinopyroxene domains have been measured by LA-ICP-MS. The melt trace element composition has been constrained from bulk rock analyses of the fine-grained low- and high-Mg diorites. Clinopyroxene/melt partition coefficients for rare earth elements (REE) and Y in the high-Mg group zonings (Mg# > 73–79, DDy < 1.2) are positively correlated with tetrahedral IVAl and increase by a factor of 3–4 as tetrahedral IVAl increases from 0.01 to 0.1 per formula unit (pfu). These systematic variations are interpreted to be controlled by the clinopyroxene composition. In contrast, partition coefficients for low-Mg group zonings (Mg# < 75–79, DDy > 1.2) are elevated by up to an order of magnitude (for REE and Y) or more (for Zr and Hf) at similar IVAl, indicating dominant control of melt composition/structure. DZr and DHf show a larger sensitivity to the compositional change of crystal and melt than DREE. DTi values for the low- and high-Mg zonings show a uniform dependence on IVAl. DSr and DLi are insensitive to the compositional change of clinopyroxene and melt, resulting in Sr depletions in the clinopyroxene zonings with elevated REE without crystallization of plagioclase. Our observations show that crystal chemistry and melt composition/structure may alternatively control clinopyroxene/melt partitioning during the assimilation of peridotite by felsic magma, and may be useful for deciphering clinopyroxene compositions and related crust–mantle processes.

Impact of the arc length on GNSS analysis results

Homogeneously reprocessed combined GPS/GLONASS 1- and 3-day solutions from 1994 to 2013, generated by the Center for Orbit Determination in Europe (CODE) in the frame of the second reprocessing campaign REPRO-2 of the International GNSS Service, as well as GPS- and GLONASS-only 1- and 3-day solutions for the years 2009 to 2011 are analyzed to assess the impact of the arc length on the estimated Earth Orientation Parameters (EOP, namely polar motion and length of day), on the geocenter, and on the orbits. The conventional CODE 3-day solutions assume continuity of orbits, polar motion components, and of other parameters at the day boundaries. An experimental 3-day solution, which assumes continuity of the orbits, but independence from day to day for all other parameters, as well as a non-overlapping 3-day solution, is included into our analysis. The time series of EOPs, geocenter coordinates, and orbit misclosures, are analyzed. The long-arc solutions were found to be superior to the 1-day solutions: the RMS values of EOP and geocenter series are typically reduced between 10 and 40 %, except for the polar motion rates, where RMS reductions by factors of 2–3 with respect to the 1-day solutions are achieved for the overlapping and the non-overlapping 3-day solutions. In the low-frequency part of the spectrum, the reduction is even more important. The better performance of the orbits of 3-day solutions with respect to 1-day solutions is also confirmed by the validation with satellite laser ranging.

Origin and history of ureilitic material in the solar system: The view from asteroid 2008 TC3 and the Almahata Sitta meteorite

Asteroid 2008 TC3 (approximately 4m diameter) was tracked and studied in space for approximately 19h before it impacted Earth's atmosphere, shattering at 44-36km altitude. The recovered samples (>680 individual rocks) comprise the meteorite Almahata Sitta (AhS). Approximately 50-70% of these are ureilites (ultramafic achondrites). The rest are chondrites, mainly enstatite, ordinary, and Rumuruti types. The goal of this work is to understand how fragments of so many different types of parent bodies became mixed in the same asteroid. Almahata Sitta has been classified as a polymict ureilite with an anomalously high component of foreign clasts. However, we calculate that the mass of fallen material was 0.1% of the pre-atmospheric mass of the asteroid. Based on published data for the reflectance spectrum of the asteroid and laboratory spectra of the samples, we infer that the lost material was mostly ureilitic. Therefore, 2008 TC3 probably contained only a few percent nonureilitic materials, similar to other polymict ureilites except less well consolidated. From available data for the AhS meteorite fragments, we conclude that 2008 TC3 samples essentially the same range of types of ureilitic and nonureilitic materials as other polymict ureilites. We therefore suggest that the immediate parent of 2008 TC3 was the immediate parent of all ureilitic material sampled on Earth. We trace critical stages in the evolution of that material through solar system history. Based on various types of new modeling and re-evaluation of published data, we propose the following scenario. (1) The ureilite parent body (UPB) accreted 0.5-0.6Ma after formation of calcium-aluminum-rich inclusions (CAI), beyond the ice line (outer asteroid belt). Differentiation began approximately 1Ma after CAI. (2) The UPB was catastrophically disrupted by a major impact approximately 5Ma after CAI, with selective subsets of the fragments reassembling into daughter bodies. (3) Either the UPB (before breakup), or one of its daughters (after breakup), migrated to the inner belt due to scattering by massive embryos. (4) One daughter (after forming in or migrating to the inner belt) became the parent of 2008 TC3. It developed a regolith, mostly 3.8Ga ago. Clasts of enstatite, ordinary, and Rumuruti-type chondrites were implanted by low-velocity collisions. (5) Recently, the daughter was disrupted. Fragments were injected or drifted into Earth-crossing orbits. 2008 TC3 comes from outer layers of regolith, other polymict ureilites from deeper regolith, and main group ureilites from the interior of this body. In contrast to other models that have been proposed, this model invokes a stochastic history to explain the unique diversity of foreign materials in 2008 TC3 and other polymict ureilites.

On the theory of high and low pressure areas: The significance of divergence in pressure areas

This is the edited and translated version of the article by Richard Scherhag “Zur Theorie der Hoch- und Tiefdruckgebiete. Die Bedeutung der Divergenz in Druckfeldern” (On the theory of high and low pressure areas: The significance of divergence in pressure areas), which was published in Meteorologische Zeitschrift 51, 129–138.

Appendix A. Whole-rock major, trace and rare earth element analyses of the Ulubey volcanics, Appendix B. Analytical standards and detection limits for analyses

Collisional and post-collisional volcanic rocks in the Ulubey (Ordu) area at the western edge of the Eastern Pontide Tertiary Volcanic Province (EPTVP) in NE Turkey are divided into four suites; Middle Eocene (49.4-44.6 Ma) aged Andesite-Trachyandesite (AT), Trachyandesite-Trachydacite-Rhyolite (TTR), Trachydacite-Dacite (TD) suites, and Middle Miocene (15.1 Ma) aged Trachybasalt (TB) suite. Local stratigraphy in the Ulubey area starts with shallow marine environment sediments of the Paleocene-Eocene time and then continues extensively with sub-aerial andesitic to rhyolitic and rare basaltic volcanism during Eocene and Miocene time, respectively. Petrographically, the volcanic rocks are composed primarily of andesites/trachyandesites, with minor trachydacites/rhyolites, basalts/trachybasalts and pyroclastics, and show porphyric, hyalo-microlitic porphyric and rarely glomeroporphyric, intersertal, intergranular, fluidal and sieve textures. The Ulubey (Ordu) volcanic rocks indicate magma evolution from tholeiitic-alkaline to calc-alkaline with medium-K contents. Primitive mantle normalized trace element and chondrite normalized rare earth element (REE) patterns show that the volcanic rocks have moderate light rare earth element (LREE)/heavy rare earth element (HREE) ratios relative to E-Type MORB and depletion in Nb, Ta and Ti. High Th/Yb ratios indicate parental magma(s) derived from an enriched source formed by mixing of slab and asthenospheric melts previously modified by fluids and sediments from a subduction zone. All of the volcanic rocks share similar incompatible element ratios (e.g., La/Sm, Zr/Nb, La/Nb) and chondrite-normalized REE patterns, indicating that the basic to acidic rocks originated from the same source. The volcanic rocks were produced by the slab dehydration-induced melting of an existing metasomatized mantle source, and the fluids from the slab dehydration introduced significant large ion lithophile element (LILE) and LREE to the source, masking its inherent HFSE-enriched characteristics. The initial 87Sr/86Sr (0.7044-0.7050) and eNd (-0.3 to +3.4) ratios of the volcanics suggest that they originated from an enriched lithospheric mantle source with low Sm/Nd ratios. Integration of the geochemical, petrological and isotopical with regional and local geological data suggest that the Tertiary volcanic rocks from the Ulubey (Ordu) area were derived from an enriched mantle, which had been previously metasomatized by fluids derived from subducted slab during Eocene to Miocene in collisional and post-collisional extension-related geodynamic setting following Late Mesozoic continental collision between the Eurasian plate and the Tauride-Anatolide platform.

A low-latitude Upper Pleistocene oxygen isotope stack

Below oxygen isotope stage 16, the orbitally derived time-scale developed by Shackleton et al. (1990) from ODP site 677 in the equatorial Pacific differs significantly from previous ones (e.g. Kominz and Pisias, 1979 doi:10.1126/science.204.4389.171; Morley and Hays, 1981 doi:10.1016/0012-821X(81)90034-0, Imbrie et al. 1984), yielding estimated ages for the last Earth magnetic reversals that are 5-7% older than the K/Ar values (Mankinen and Dalrymple, 1979 doi:10.1029/JB084iB02p00615; Berggren et al., 1985; Harland and Armstrong, 1989) but are in good agreement with recent Ar/Ar dating (Baksi et al., 1991; 1992 doi:10.1126/science.256.5055.356; Spell and McDougall, 1992 doi:10.1029/92GL01125). These results suggest that in the lower Brunhes and upper Matuyama chronozones most deep-sea climatic records retrieved so far apparently missed or misinterpreted several oscillations predicted by the astronomical theory of climate. To test this hypothesis, we studied a high-resolution oxygen isotope record from giant piston core MD900963 (Maldives area, tropical Indian Ocean) in which precession-related oscillations in delta18O are particularly well expressed, owing to the superimposition of a local salinity signal on the global ice volume signal (Rostek et al., 1993 doi:10.1038/364319a0). Three additional precession-related cycles are observed in oxygen isotope stages 17 and 18 of core MD900963, compared to the SPECMAP composite curves (Imbrie et al., 1984; Prell et al., 1986 doi:10.1029/PA001i002p00137), and stage 21 clearly presents three precession oscillations, as predicted by Shackleton et al. (1990). The precession peaks found in the delta18O record from core MD900963 are in excellent agreement with climatic oscillations predicted by the astronomical theory of climate. Our delta18O record therefore permits the development of an accurate astronomical time-scale. Based on our age model, the Brunhes-Matuyama reversal is dated at 775 +/- 10 ka, in good agreement with the age estimate of 780 ka obtained by Shackleton et al. (1990) and recent radiochronological Ar/Ar datings on lavas (Baksi et al., 1991; 1992; Spell and McDougall, 1992). We developed a new low-latitude, Upper Pleistocene delta18O reference record by stacking and tuning the delta18O records from core MD900963 and site 677 to orbital forcing functions.

Rare earth elements determined in Antarctic ice by inductively coupled plasma

Inductively coupled plasma mass spectrometry (ICP-MS) is a suitable tool for multi-element analysis at low concentration levels. Rare earth element (REE) determinations in standard reference materials and small volumes of molten ice core samples from Antarctica have been performed with an ICP-time of flight-MS (ICP-TOF-MS) system. Recovery rates for REE in e.g. SPS-SW1 amounted to not, vert, similar ~103%, and the relative standard deviations were 3.4% for replicate analysis at REE concentrations in the lower ng/l range. Analyses of REE concentrations in Antarctic ice core samples showed that the ICP-TOF-MS technique meets the demands of restricted sample mass. The data obtained are in good agreement with ICP-Quadrupole-MS (ICP-Q-MS) and ICP-Sector Field-MS (ICP-SF-MS) results. The ICP-TOF-MS system determines accurately and precisely REE concentrations exceeding 5 ng/l while between 0.5 and 5 ng/l accuracy and precision are element dependent.

Laboratory investigation of flow structure and resistance of high- and low-angle dunes

A prominent control on the flow over subaqueous dunes is the slope of the downstream leeside. While previous work has focused on steep (~30°), asymmetric dunes with permanent flow separation, little is known about dunes with lower lee-slope angles for which flow separation is absent or intermittent. Here, we present a laboratory investigation where we systematically varied the dune lee-slope, holding other geometric parameters and flow hydraulics constant, to explore effects on the turbulent flow field and flow resistance. Three sets of fixed dunes (lee-slopes of 10°, 20° and 30°) were separately installed in a 15 m long and 1 m wide flume and subjected to 0.20 m deep flow. Measurements consisted of high-frequency, vertical profiles collected with a Laser Doppler Velocimeter (LDV). We show that the temporal and spatial occurrence of flow separation decreases with dune lee-slope. Velocity gradients in the dune leeside depict a free shear layer downstream of the 30° dunes and a weaker shear layer closer to the bed for the 20° and 10° dunes. The decrease in velocity gradients leads to lower magnitude of turbulence production for gentle lee-slopes. Aperiodic, strong ejection events dominate the shear layer, but decrease in strength and frequency for low-angle dunes. Flow resistance of dunes decreases with lee-slope; the transition being non-linear. Over the 10°, 20° and 30° dunes, shear stress is 8%, 33% and 90 % greater than a flat bed, respectively. Our results demonstrate that dune lee-slope plays an important, but often ignored role in flow resistance.

Major, trace and rare earth elements and minerals at DSDP Legs 69-70 Holes

We obtained major and trace element data on 113 samples from basalts drilled during DSDP Legs 69 and 70 in the Costa Rica Rift area. The majority have major and trace element characteristics typical of ocean-ridge tholeiities. Most of the basalts are relatively MgO rich (MgO > 8 wt.%) and have Mg values (MgO/MgO + 0.85FeO x 100) of about 53, characteristics that clearly indicate that the various magmas underwent only a small amount of crystal fractionation before being erupted onto the seafloor. According to their normative mineralogies, the rocks are olivine tholeiites. A few samples plot close to the diopside-hypersthene join of the projected basalt tetrahedron. Except for basalts from two thin intervals in Hole 504B, which differ significantly from all the other basalts of the hole, practically no chemical downhole variation could be established. In the two exceptional intervals, both TiO2 and P2O5 contents are markedly enriched among the major oxides. The trace elements in these intervals are distinguished by relatively high contents of magmatophile elements and have flat to enriched chondrite-normalized distribution patterns of light rare earth elements (LREE). Most of the rocks outside these intervals are strongly depleted in large-ionlithophile (LIL) elements and LREE. We offer no satisfactory hypothesis for the origin of these basalts at this time. They might have originated within pockets of mantle materials that were more primitive than the LIL-element-depleted magmas that were the source of the other basalts. A significant change with depth in the type of alteration occurs in the 561 meters of basalt cored in Hole 504B. According to the behavior of such alteration-sensitive species as K2O, H2O-, CO2, S, Tl, and the iron oxidation ratio, the alteration is oxidative in the upper part and nonoxidative or even reducing in the lower part. The oxidative alteration may have resulted from low temperature basalt/seawater interaction, whereas hydrothermal solutions may be responsible for the nonoxidative alteration.

Major oxides, trace elements and rare earth elements of selected basalt samples at DSDP Hole 83-504B

DSDP Hole 504B is the deepest section drilled into oceanic basement, penetrating through a 571.5-m lava pile and a 209-m transition zone of lavas and dikes into 295 m of a sheeted dike complex. To define the basement composition 194 samples of least altered basalts, representing all lithologic units, were analyzed for their major and 26 trace elements. As is evident from the alteration-sensitive indicators H2O+, CO2, S, K, Mn, Zn, Cu, and the iron oxidation ratio, all rocks recovered are chemically altered to some extent. Downhole variation in these parameters enables us to distinguish five depth-related alteration zones that closely correlate with changes in alteration mineralogy. Alteration in the uppermost basement portion is characterized by pronounced K-uptake, sulfur loss, and iron oxidation and clearly demonstrates low-temperature seawater interaction. A very spectacular type of alteration is confined to the depth range from 910 to 1059 m below seafloor (BSF). Rocks from this basement portion exhibit the lowest iron oxidation, the highest H2O+ contents, and a considerable enrichment in Mn, S, Zn, and Cu. At the top of this zone a stockwork-like sulfide mineralization occurs. The chemical data suggest that this basement portion was at one time within a hydrothermal upflow zone. The steep gradient in alteration chemistry above this zone and the ore precipitation are interpreted as the result of mixing of the upflowing hydrothermal fluids with lower-temperature solutions circulating in the lava pile. Despite the chemical alteration the primary composition and variation of the rocks can be reliably established. All data demonstrate that the pillow lavas and the dikes are remarkably uniform and display almost the same range of variation. A general characteristic of the rocks that classify as olivine tholeiites is their high MgO contents (up to 10.5 wt.%) and their low K abundances (-200 ppm). According to their mg-values, which range from 0.60 to 0.74, most basalts appear to have undergone some high-level crystal fractionation. Despite the overall similarity in composition, there are two major basalt groups that have significantly different abundances and ratios of incompatible elements at similar mg-values. The majority of the basalts from the pillow lava and dike sections are chemically closely related, and most probably represent differentiation products of a common parental magma. They are low in Na2O, TiO2, and P2O5, and very low in the more hygromagmaphile elements. Interdigitated with this basalt group is a very rarely occurring basalt that is higher in Na2O, TiO2, P2O5, much less depleted in hygromagmaphile elements, and similar to normal mid-ocean ridge basalt (MORB). The latter is restricted to Lithologic Units 5 and 36 of the pillow lava section and Lithologic Unit 83 of the dike section. The two basalt groups cannot be related by differentiation processes but have to be regarded as products of two different parental magmas. The compositional uniformity of the majority of the basalts suggests that the magma chamber beneath the Costa Rica Rift reached nearly steady-state conditions. However, the presence of lavas and dikes that crystallized from a different parental magma requires the existence of a separate conduit-magma chamber system for these melts. Occasionally mixing between the two magma types appears to have occurred. The chemical characteristics of the two magma types imply some heterogeneity in the mantle source underlying the Costa Rica Rift. The predominant magma type represents an extremely depleted source, whereas the rare magma type presumably originated from regions of less depleted mantle material (relict or affected by metasomatism).

Major and rare earth element composition of ODP Site 191-1179 sediments, northwest Pacific

Depth profiles of major element and rare earth element (REE) abundances in sediment samples (mainly siliceous ooze and clay) recovered from Holes 1179B and 1179C at Site 1179, Ocean Drilling Program Leg 191 (41.4°N, 159.6°E) were determined. The oxidation states of Mn and Ce were determined by X-ray absorption near-edge structure. Some geochemical indicators were tested, including the MnO/TiO2 ratios, a bivariate diagram of La/Ce vs. Al2O3/(Al2O3+Fe2O3), and other discrimination diagrams. The oxidation state of Mn is reduced Mn(II) in the depth profile below 0.60 meters below seafloor (mbsf), which is consistent with relatively low and high abundances of Mn in the sediments and pore waters, respectively. It is possible that the diagenetic effect on the oxidation state and abundance of Mn makes it difficult for the MnO/TiO2 ratio to reflect the depositional environment. The normalized ratio of La and Ce does not change very much with depth, suggesting that the diagenetic effect does not affect the REE signature in the sediments. On the diagram of La/Ce vs. Al2O3/(Al2O3+Fe2O3), the sediments studied here plot at the boundary of the pelagic and continental margin fields. This suggests that continental material has contributed to the sediment to some degree, even though Site 1179 is in a pelagic region of the northwestern Pacific Ocean, >1600 km from Japan.

Laboratory shearing experiments at low sliding velocities of IODP Hole 343-C0019E

The 2011 Tohoku-Oki earthquake demonstrated that the shallowest reaches of plate boundary subduction megathrusts can host substantial coseismic slip that generates large and destructive tsunamis, contrary to the common assumption that the frictional properties of unconsolidated clay-rich sediments at depths less than View the MathML source should inhibit rupture. We report on laboratory shearing experiments at low sliding velocities (View the MathML source) using borehole samples recovered during IODP Expedition 343 (JFAST), spanning the plate-boundary décollement within the region of large coseismic slip during the Tohoku earthquake. We show that at sub-seismic slip rates the fault is weak (sliding friction µs=0.2-0.26), in contrast to the much stronger wall rocks (µs>~0.5). The fault is weak due to elevated smectite clay content and is frictionally similar to a pelagic clay layer of similar composition. The higher cohesion of intact wall rock samples coupled with their higher amorphous silica content suggests that the wall rock is stronger due to diagenetic cementation and low clay content. Our measurements also show that the strongly developed in-situ fabric in the fault zone does not contribute to its frictional weakness, but does lead to a near-cohesionless fault zone, which may facilitate rupture propagation by reducing shear strength and surface energy at the tip of the rupture front. We suggest that the shallow rupture and large coseismic slip during the 2011 Tohoku earthquake was facilitated by a weak and cohesionless fault combined with strong wall rocks that drive localized deformation within a narrow zone.