A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Petrogenetic modelling of strongly residual metapelitic xenoliths within the southern Platreef,Bushveld Complex,South Africa

Xenoliths of quartz‐absent Fe‐rich aluminous metapelite are common within the platinum group element‐rich mafic/ultramafic magmatic rocks of the Platreef. Relative to well‐characterized protoliths, the xenoliths are strongly depleted in K₂O and H₂O, and have lost a substantial amount of melt (>50 vol.%). Mineral equilibria calculations in the NCKFMASHTO system yield results that are consistent with observations in natural samples. Lower‐grade rocks that lack staurolite constrain peak pressures to ～2.5 kbar in the southern Platreef. Smaller xenoliths and the margins of larger xenoliths comprise micro‐diatexite rich in coarse acicular corundum and spinel, which record evidence for the metastable persistence of lower‐grade hydrous phases and rapid melting consequent on a temperature overstep of several hundred degrees following their incorporation in the mafic/ultramafic magmas. In the cores of larger xenoliths, temperatures increased more slowly enabling progressive metamorphism by continuous prograde equilibration and the loss of H₂O by subsolidus dehydration; the H₂O migrated to xenolith margins where it may have promoted increased melting. According to variations in the original compositional layering, layers became aluminosilicate‐ and/or cordierite‐rich, commonly with spinel but only rarely with corundum. The differing mineralogical and microstructural evolution of the xenoliths depends on heating rates (governed by their size and, therefore, proximity to the Platreef magmas) and the pre‐intrusive metamorphic grade of the protoliths. The presence or absence of certain phases, particularly corundum, is strongly influenced by the degree of metastable retention of lower‐grade hydrates in otherwise identical protolith bulk compositions. The preservation of fine‐scale compositional layering that is inferred to be relict bedding in xenolith cores implies that melt loss by compaction was extremely efficient.     

Timing effects on the fragmentation of small scale blocks of granodiorite

A series of small scale tests, simulating multi-hole blasts have been performed to establish the effect of delays on blast fragmentation. The blasts were performed in high quality granodiorite blocks, which were cut from stone prepared by dimensional stone quarry operations. The pattern used was equilateral triangular, with a distance of 10.2 cm between boreholes, which had a diameter of 11 mm, were loaded with detonating cord and the coupling medium was water. The delays used were achieved using different lengths of detonating cord for the cases of delays between 0 and 100 μs between holes and a sequential blasting machine firing seismic detonators for larger delays up to 4 ms. All fragments were collected and screened. The experiments showed that the worst fragmentation was achieved with simultaneous initiation of all charges. Fragmentation improved with the delay time between holes up to 1 ms between holes. If the experiments are scaled up, the results show that in granodiorite, fragmentation optimization requires delays of few milliseconds per metre of burden. The findings, agree with previously published work, involving larger scale experiments and other rock types.     

An absolute paleotemperature record from 10 to 6 Ka inferred from fluid inclusion D/H ratios of a stalagmite from Vancouver Island, British Columbia, Canada

By using continuous helium flow during the crushing of calcite speleothem samples, we are able to recover liberated inclusion waters without isotopic fractionation. A paleotemperature record for the Jacklah Jill Cave locality, Vancouver Island, BC, was obtained from a 30-cm tall stalagmite that grew 10.3-6.3 Ka ago, using δ¹⁸O values of the crushed calcite and of the inclusion water as inferred from its δD. It is found that the locality experienced mean annual temperature variations up to 11 °C over a 4-Ka period in the early Holocene. At the beginning of the period, local temperature quickly increased from a minimum of ∼1 °C to around 10 °C, but this early climate optimum, about 3 °C warmer than today, only lasted for ∼1200 years. About 8.6 Ka ago, temperature had declined to ∼7 °C, approximately the same as the modern cave temperature. Since then, the study area has experienced only minor temperature fluctuations, but there was a brief fall to ∼4 °C at around 7 Ka ago, which might be caused by a short lived expansion of local alpine glaciers. The long-term T-dependence of δD was 1.47‰/°C, identical to the value in modern precipitation.     

The petroleum evaluation of a tectonically complex area: The western margin of the Southeast Basin (France)

The Southeast Basin of France is the thickest onshore French sedimentary basin which contains locally as much as 10 km of Mesozoic-Cenozoic sediment. Basin development occurred in several stages between late Carboniferous and late Cretaceous times. Partial tectonic inversion took place during two compressive events, the so-called ‘Pyrenean’ and ‘Alpine’ phases of late Cretaceous-early Tertiary and late Tertiary ages respectively. They are separated by an intervening stretching event of Oligocene age, which further south resulted in the opening of the western Mediterranean oceanic basin. As a result of this complex tectonic history, structural traps were difficult to image on the seismic data shot during the first phase of exploration prior to 1980. Oil and gas natural seeps, and shows in several wells, indicate that some petroleum systems are, or have been active, at least in some places.The present erosional western margin of the basin is more or less superimposed on the initial Triassic-Jurassic margin. Margin subsidence and Tertiary inversion are discussed using regional sections on which the polyphase history of the entire basin is well shown. These sections are located on three major segments where the Mesozoic margin is either partly preserved (Ardèche), or has been partly inverted in late Tertiary times (Vercors-Chartreuse), or has been completely inverted in early Tertiary times (Corbières-eastern Pyrenees). 1-D ‘Genex’ basin modelling on the Ardèche segment, and 2-D ‘Thrustpack’ structural-maturity modelling in the Vercors-Chartreuse segment are used to further assess the remaining petroleum plays.     

Isotope geochemistry and Re–Os geochronology of the Yanjiagou Mo deposit in the central North China Craton

The Yanjiagou deposit, located in the central North China Craton (NCC), is a newly found porphyry‐type Mo deposit. The Mo mineralization here is spatially associated with the Mapeng batholith. In this study, we identify four stages of ore formation in this deposit: pyrite phyllic stage (I), quartz–pyrite stage (II), quartz–pyrite–molybdenite stage (III), which is the main mineralization stage, and quartz–carbonate stage (IV). We present sulphur and lead isotope data on pyrite, and rhenium and osmium isotopes of molybdenite from the porphyry deposit and evaluate the timing and origin of ore formation. The δ³⁴S values of the pyrite range from ‐1.1‰ to −0.6‰, with an average of −0.875‰, suggesting origin from a mixture of magmatic/mantle sources and the basement rocks. The Pb isotope compositions of the pyrite show a range of 16.369 to 17.079 for ²⁰⁶Pb/²⁰⁴Pb, 15.201 to 15.355 for ²⁰⁷Pb/²⁰⁴Pb, and 36.696 to 37.380 for ²⁰⁸Pb/²⁰⁴Pb, indicating that the ore‐forming materials were derived from a mixture of lower crust (or basement rocks) and mantle. Rhenium contents in molybdenite samples from the main ore stage are between 74.73 to 254.43 ppm, with an average of 147.9 ppm, indicating a mixed crustal‐mantle source for the metal. Eight molybdenite separates yield model ages ranging from 124.17 to 130.80 Ma and a mean model age of 128.46 Ma. An isochron age of 126.7 ± 1.1 Ma (MSWD = 2.1, initial ¹⁸⁷Os = 0.0032 ± 0.0012 ppb) is computed, which reveals a close link between the Mo mineralization and the magmatism that generated the Mapeng batholith. The age is close to the zircon U–Pb age of ca. 130 Ma from the batholith reported in a recent study. The age is also consistent with the timing of mineralization in the Fuping ore cluster in the central NCC, as well as the peak time of lithosphere thinning and destruction of the NCC. We evaluate the spatio‐temporal distribution of the Mo deposits in the NCC and identify three important molybdenum provinces along the northern and southern margins of the craton formed during three distinct episodes: Middle to Late Triassic (240–220 Ma), Early Jurassic (190–175 Ma), and Late Jurassic to Early Cretaceous (150–125 Ma). The third period is considered to mark the most important metallogenic event, coinciding with the peak of lithosphere thinning and craton destruction in the NCC. Copyright © 2014 John Wiley & Sons, Ltd.     

Ion microprobe analysis of oxygen isotope ratios in granulite facies magnetites: diffusive exchange as a guide to cooling history

Ion microprobe analysis of magnetites from the Adirondack Mountains, NY, yields oxygen isotope ratios with spatial resolution of 2–8 m and precision in the range of 1 (1 sigma). These analyses represent 11 orders of magnitude reduction in sample size compared to conventional analyses on this material and they are the first report of routinely reproducible precision in the 1 per mil range for analysis of ¹⁸O at this scale. High precision micro-analyses of this sort will permit wide-ranging new applications in stable isotope geochemistry. The analyzed magnetites form nearly spherical grains in a calcite matrix with diopside and monticellite. Textures are characteristic of granulite facies marbles and show no evidence for retrograde recrystallization of magnetite. Magnetites are near to Fe₃O₄ in composition, and optically and chemically homogeneous. A combination of ion probe plus conventional BrF₅ analysis shows that individual grains are homogeneous with ¹⁸O=8.9±1 SMOW from the core to near the rim of 0.1–1.2 mm diameter grains. Depth profiling into crystal growth faces of magnetites shows that rims are 9 depleted in ¹⁸O. These low ¹⁸O values increase in smooth gradients across the outer 10 m of magnetite rims in contact with calcite. These are the sharpest intracrystalline gradients measured to date in geological materials. This discovery is confirmed by bulk analysis of 150–350 m diameter magnetites which average 1.2 lower in ¹⁸O than coarse magnetites due to low ¹⁸O rims. Conventional analysis of coexisting calcite yields °¹⁸O=18.19, suggesting that bulk ¹⁸O (Cc-Mt)=9.3 and yielding an apparent equilibration temperature of 525° C, over 200° C below the temperature of regional metamorphism. Consideration of experimental diffusion data and grain size distribution for magnetite and calcite suggests two contrasting cooling histories. The data for oxygen in calcite under hydrothermal conditions at high P(H₂O) indicates that diffusion is faster in magnetite and modelling of the low ¹⁸O rims on magnetite would suggest that the Adirondacks experienced slow cooling after Grenville metamorphism, followed by a brief period of rapid cooling, possibly related to uplift. Conversely, the data for calcite at low P(H₂O) show slower oxygen diffusion than in magnetite. Modelling based on these data is consistent with geochronology that shows slow cooling through the blocking temperature of both minerals, suggesting that the low ¹⁸O rims form by exchange with late, low temperature fluids similar to those that infiltrated the rock to serpentinize monticellite and which infiltrated adjacent anorthosite to form late calcite veinlets. In either case, the ion microprobe results indicate that two distinct events are recorded in the post-metamorphic exchange history of these magnetites. Recognition of these events is only possible through microanalysis and has important implications for geothermometry.     

Petrology of ultramafic and mafic xenoliths in picrite of Kahoolawe Island,Hawaii

A picrite lava (22 wt% MgO; 35 vol.% ol) along the western shore of the1.3–1.4 Ma Kahoolawe tholeiitic shield, Hawaii, contains small xenoliths of harzburgite, lherzolite, norite, and wehrlite. The various rock types have textures where either orthopyroxene, clinopyroxene, or plagioclase is in a poikilitic relationship with olivine. The Mg#s of the olivine, orthopyroxene, and clinopyroxene in this xenolith suite range between 86 and 82; spinel Mg#s range from 60 to 49, and plagioclase is An_75–80. A ⁸⁷Sr/⁸⁶Sr ratio for one ol-norite xenolith is 0.70444. In comparison, the host picrite has olivine phenocrysts with an average Mg# of 86.2 (range 87.5–84.5), and a whole-rock ⁸⁷Sr/⁸⁶Sr ratio of 0.70426. Textural and isotopic information together with mineral compositions indicate that the xenoliths are related to Kahoolawe tholeiitic magmatism, but are not crystallization products of the magma represented by their host picrite. Rather, the xenoliths are crystalline products of earlier primitive liquids (FeO/MgO ranging 1 to 1.3) at 5–9 kbar in the cumulate environment of a magma reservoir or conduit system. The presence of ultramafic xenoliths in picrite but not in typical Kahoolawe tholeiitic lava (6–9 wt% MgO) is consistent with replenishment of reservoirs by dense Mg-rich magma emplaced beneath resident, less dense tholeiitic magma. Mg-rich magmas have proximity to reservoir cumulate zones and are therefore more likely than fractionated residual liquids to entrain fragments of cumulate rock.     

Combining European Earth Observation products with Dynamic Global Vegetation Models for estimating Essential Biodiversity Variables

ABSTRACT

Global, fast and accessible monitoring of biodiversity is one of the main pillars of the efforts undertaken in order to revert it loss. The Group on Earth Observations Biodiversity Observation Network (GEO-BON) provided an expert-based definition of the biological properties that should be monitored, the Essential Biodiversity Variables (EBVs). Initiatives to provide indicators for EBVs rely on global, freely available remote sensing (RS) products in combination with empirical models and field data, and are invaluable for decision making. In this study, we provide alternatives for the expansion and improvement of the EBV indicators, by suggesting current and future data from the European Space Agencýs COPERNICUS and explore the potential of RS-integrated Dynamic Global Vegetation Models (DGVMs) for the estimation of EBVs. Our review found that mainly due to the inclusion of the Sentinel constellation, Copernicus products have similar or superior potential for EBV indicator estimation in relation to their NASA counterparts. DGVMs simulate the ecosystem level EBVs (ecosystem function and structure), and when integrated with remote sensing data have great potential to not only offer improved estimation of current states but to provide projection of ecosystem impacts. We suggest that focus on producing EBV relevant outputs should be a priority within the research community, to support biodiversity preservation efforts.