Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 2: Characteristics and petrogenesis of the high-Ti latite and high-Ti and low-Ti voluminous quartz latite eruptives

As a result of their relative concentration towards the respective Atlantic margins, the silicic eruptives of the Parana (Brazil)-Etendeka large igneous province are disproportionately abundant in the Etendeka of Namibia. The NW Etendeka silicic units, dated at similar to132 Ma, occupy the upper stratigraphic levels of the volcanic sequences, restricted to the coastal zone, and comprise three latites and five quartz latites (QL). The large-volume Fria QL is the only low-Ti type. Its trace element and isotopic signatures indicate massive crustal input. The remaining NW Etendeka silicic units are enigmatic high-Ti types, geochemically different from low-Ti types. They exhibit chemical affinities with the temporally overlapping Khumib high-Ti basalt (see Ewart et al. Part 1) and high crystallization temperatures (greater than or equal to980 to 1120degreesC) inferred from augite and pigeonite phenocrysts, both consistent with their evolution from a mafic source. Geochemically, the high-Ti units define three groups, thought genetically related. We test whether these represent independent liquid lines of descent from a common high-Ti mafic parent. Although the recognition of latites reduces the apparent silica gap, difficulty is encountered in fractional crystallization models by the large volumes of two QL units. Numerical modelling does, however, support large-scale open-system fractional crystallization, assimilation of silicic to basaltic materials, and magma mixing, but cannot entirely exclude partial melting processes within the temporally active extensional environment. The fractional crystallization and mixing signatures add to the complexity of these enigmatic and controversial silicic magmas. The existence, however, of temporally and spatially overlapping high-Ti basalts is, in our view, not coincidental and the high-Ti character of the silicic magmas ultimately reflects a mantle signature.

Effects Of Partial Inhibition Of Respiratory Complex I On H2o 2 Production By Isolated Brain Mitochondria In Different Respiratory States.

The aim of this work was to characterize the effects of partial inhibition of respiratory complex I by rotenone on H2O2 production by isolated rat brain mitochondria in different respiratory states. Flow cytometric analysis of membrane potential in isolated mitochondria indicated that rotenone leads to uniform respiratory inhibition when added to a suspension of mitochondria. When mitochondria were incubated in the presence of a low concentration of rotenone (10 nm) and NADH-linked substrates, oxygen consumption was reduced from 45.9 ± 1.0 to 26.4 ± 2.6 nmol O2 mg(-1) min(-1) and from 7.8 ± 0.3 to 6.3 ± 0.3 nmol O2 mg(-1) min(-1) in respiratory states 3 (ADP-stimulated respiration) and 4 (resting respiration), respectively. Under these conditions, mitochondrial H2O2 production was stimulated from 12.2 ± 1.1 to 21.0 ± 1.2 pmol H2O2 mg(-1) min(-1) and 56.5 ± 4.7 to 95.0 ± 11.1 pmol H2O2 mg(-1) min(-1) in respiratory states 3 and 4, respectively. Similar results were observed when comparing mitochondrial preparations enriched with synaptic or nonsynaptic mitochondria or when 1-methyl-4-phenylpyridinium ion (MPP(+)) was used as a respiratory complex I inhibitor. Rotenone-stimulated H2O2 production in respiratory states 3 and 4 was associated with a high reduction state of endogenous nicotinamide nucleotides. In succinate-supported mitochondrial respiration, where most of the mitochondrial H2O2 production relies on electron backflow from complex II to complex I, low rotenone concentrations inhibited H2O2 production. Rotenone had no effect on mitochondrial elimination of micromolar concentrations of H2O2. The present results support the conclusion that partial complex I inhibition may result in mitochondrial energy crisis and oxidative stress, the former being predominant under oxidative phosphorylation and the latter under resting respiration conditions.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Nickel complexes with tris(2-aminoethyl)amine (tren): [Ni-3(tren)(4)(H2O)(2)][Cr(ox)(3)](2)center dot 6H(2)O (ox = oxalate), {[Ni-2(tren)(3)][ClO4](4)center dot H2O}(n), and [Ni-2(tren)(2)(aepd)][ClO4](4)center dot 2H(2)O (aepd = N-(2-aminoethyl)pyrrolidi

Reaction of K-3[Cr(ox)(3)] (ox = oxalate) with nickel(II) and tris(2-aminoethyl)amine (tren) in aqueous solution resulted in isolation of the bimetallic assembly [Ni-3(tren)(4)(H2O)(2)][Cr(ox)(3)](2). 6H(2)O. The polymeric complex {[Ni-2(tren)(3)][ClO4](4). H2O}(n) has been prepared by reaction of nickel(II) perchlorate and tren in aqueous solution. From the same reaction mixture the complex [Ni-2(tren)(2)(aepd)][ClO4](4). 2H(2)O (aepd = N-(2-aminoethyl)pyrrolidine-3,4-diamine), in which a bridging tren ligand contains a carbon-carbon bond between two arms forming a substituted pyrrolidine, has been isolated. The complexes have been characterized by X-ray crystallography. The magnetic susceptibility (300-4.2 K) and magnetization data (2, 4 K, H = 0-5 T) for {[Ni-2(tren)(3)][ClO4](4). H2O}(n) (300 K , 4.23 mu(B)) exhibit evidence of weak antiferromagnetic coupling and zero field splitting (2J = -1.8 cm(-1); \ D\ = 2 cm(-1)) at low temperature. For [Ni-3(tren)(4)(H2O)(2)][Cr(ox)(3)](2). 6H(2)O the susceptibility data at 300 K are indicative of uncoupled nickel(II) and chromium(III) sites with zero-field splitting and intramolecular antiferromagnetic coupling predicted at low temperature.

Stratigraphy, facies architecture and tectonic implications of the Upper Devonian to Lower Carboniferous Campwyn Volcanics of the northern New England Fold Belt

Upper Devonian to Lower Carboniferous strata of the Campwyn Volcanics of east central Queensland preserve a substantial sequence of first-cycle volcaniclastic sedimentary and coeval volcanic rocks that record prolonged volcanic activity along the northern New England Fold Belt. The style and scale of volcanism varied with time, producing an Upper Devonian sequence of mafic volcano-sedimentary rocks overlain by a rhyolitic ignimbrite-dominated sequence that passes upward into a Lower Carboniferous limestone-bearing sedimentary sequence. We define two facies associations for the Campwyn Volcanics. A lower facies association is dominated by mafic volcanic-derived sedimentary breccias with subordinate primary mafic volcanic rocks comprising predominantly hyaloclastite and peperite. Sedimentary breccias record episodic and high energy, subaqueous depositional events with clastic material sourced from a mafic lava-dominated terrain. Some breccias contain a high proportion of attenuated dense, glassy mafic juvenile clasts, suggesting a syn-eruptive origin. The lower facies association coarsens upwards from a lithic sand-dominated sequence through a thick interval of pebble- to boulder-grade polymict volcaniclastic breccias, culminating in facies that demonstrate subaerial exposure. The silicic upper facies association marks a significant change in eruptive style, magma composition and the nature of eruptive sources, as well as the widespread development of subaerial depositional conditions. Crystal-rich, high-grade, low- to high-silica rhyolite ignimbrites dominate the base of this facies association. Biostratigraphic age controls indicate that the ignimbrite-bearing sequences are Famennian to lower-mid Tournaisian in age. The ignimbrites represent extra-caldera facies with individual units up to 40 m thick and mostly lacking coarse lithic breccias. Thick deposits of pyroclastic material interbedded with fine-grained siliceous sandstone and mudstone (locally radiolarian-bearing) were deposited from pyroclastic flows that crossed palaeoshorelines or represent syn-eruptive, resedimented pyroclastic material. Some block-bearing lithic-pumice-crystal breccias may also reflect more proximal subaqueous silicic explosive eruptions. Crystal-lithic sandstones interbedded with, and overlying the ignimbrites, contain abundant detrital volcanic quartz and feldspar derived from the pyroclastic deposits. Limestone is common in the upper part of the upper facies association, and several beds are oolitic (cf. Rockhampton Group of the Yarrol terrane). Overall, the upper facies association fines upward and is transgressive, recording a return to shallow-marine conditions. Palaeocurrent data from all stratigraphic levels in the Campwyn Volcanics indicate that the regional sediment-dispersal direction was to the northwest, and opposed to the generally accepted notion of easterly sediment dispersal from a volcanic arc source. The silicic upper facies association correlates in age and lithology to Early Carboniferous silicic volcanism in the Drummond (Cycle 1) and Burdekin Basins, Connors Arch, and in the Yarrol terranes of eastern Queensland. The widespread development of silicic volcanism in the Early Carboniferous indicates that silicic (rift-related) magmatism was not restricted to the Drummond Basin, but was part of a more substantial silicic igneous province.

Low-pressure, water-assisted anatexis of basic dykes in a contact metamorphic aureole, Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite-gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative delta O-18 whole-rock values (+0.2 to -3.4 parts per thousand), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest delta O-18 values (-2.8, -4.2, - 4.4 and -7.6 parts per thousand, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic delta O-18 values up to the anatectic zone, where the values are lower and more heterogeneous. Low delta O-18 values, decreasing towards the intrusion, can be ascribed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation.

Magma generation at the easternmost section of the Hellenic arc: Hf, Nd, Pb and Sr isotope geochemistry of Nisyros and Yali volcanoes (Greece)

Geochemical and petrographical studies of lavas and ignimbrites from the Quaternary Nisyros-Yali volcanic system in the easternmost part of the Hellenic arc (Greece) reveal insight into magma generating processes. A compositional gap between 61 and 68 wt.% SiO2 is recognized that coincides with the stratigraphic distinction between pre-caldera and postcaldera volcanic units. Trace element systematics support the subdivision of Nisyros and Yali volcanic units into two distinct suites of rocks. The variation of Nd and Hf present day isotope data and the fact that they are distinct from the isotope compositions of MORB rule out an origin by pure differentiation and require assimilation of a crustal component. Lead isotope ratios of Nisyros and Yali volcanic rocks support mixing of mantle material with a lower crust equivalent. However, Sr-87/Sr-86 ratios of 0.7036-0.7048 are incompatible with a simple binary mixing scenario and give low depleted mantle extraction ages (< 0.1 Ga), in contrast with Pb model ages of 0.3 Ga and Hf and Nd model ages of ca. 0.8 Ga. The budget of fluid-mobile elements Sr and Pb is likely to be dominated by abundant hydrous fluids characterised by mantle-like Sr isotope ratios. Late stage fluids probably were enriched in CO2, needed to explain the high Th concentrations. The occurrence of hydrated minerals (e.g., amphibole) in the first post-caldera unit with the lowermost Sr-87/Sr-86 ratio of 0.7036 +/- 2 can be interpreted as the result of the increased water activity in the source. The presence of two different plagioclase phenocryst generations in the first lava subsequent to the caldera-causing event is indicative for a longer storage time of this magma at a shallower level. A model capable of explaining these observations involves three evolutionary stages. First stage, assimilation of lower crustal material by a primitive magma of mantle origin (as modelled by Nd-Hf isotope systematics). This stage ended by an interruption in replenishment that led to an increase of crystallization and, hence, an increase in viscosity, suppressing eruption. During this time gap, differentiation by fractional crystallization led to enrichment of incompatible species, especially aqueous fluids, to silica depolymerisation and to a decrease in viscosity, finally enabling eruption again in the third stage. (c) 2005 Elsevier B.V. All rights reserved.

Oxidation of methane at the CH4/H2O-(CO2) transition zone in the external part of the Central Alps, Switzerland: Evidence from stable isotope investigations

With the aim of understanding the mechanisms that control the metamorphic transition from the CH4- to the H2O-(CO2)-dominated fluid zone in the Helvetic domain of the Central Alps of Switzerland, fluid inclusions in quartz, illite ``crystallinity'' index, vitrinite reflectance, and the stable isotope compositions of vein and whole rock minerals and fluids trapped in quartz were investigated along four cross-sections. Increasing temperature during prograde metamorphism led to the formation of dry gas by hydrocarbon cracking in the CH4-zone. Fluid immiscibility in the H2O-CH4-(CO2)-NaCl system resulted in cogenetic, CH4- and H2O-dominated fluid inclusions. In the CH4-zone, fluids were trapped at temperatures <= 270 +/- 5 degrees C. The end of the CH4-zone is markedby a sudden increase of CO2 content in the gas phase of fluid inclusions. At temperatures > 270 +/- 5 degrees C, in the H2O-zone, the total amount of volatiles within the fluid decreased below 1 mol% with no immiscibility. This resulted m total homogenization temperatures of H2O-(CO2-CH4)-NaCl inclusions below 180 degrees C. Hydrogen isotope compositions of methane in fluid inclusion have delta D values of less than -100 parts per thousand in the CH4-zone, typical for an origin through cracking of higher hydrocarbons, but where the methane has not equilibrated with the pore water. delta D values of fluid inclusion water are around -40 parts per thousand., in isotopic equilibrium with phyllosilicates of the whole rocks. Within the CH4 to H2O(CO2) transition zone, delta D(H2O) values in fluid inclusions decrease to -130 parts per thousand interpreted to reflect the contribution of deuterium depleted water from methane oxidation. In the H2O-zone, delta D(H2O) values increase again towards an average of -30 parts per thousand which is again consistent with isotopic equilibrium with host-rock phyllosilicates. delta C-13 values of methane in fluid inclusions from the CH4-zone are around -27 parts per thousand in isotopic equilibrium with calcite in veins and whole rocks. The delta C-13(CH4) values decrease to less than -35 parts per thousand at the transition to the H2O-zone and are no longer in equilibrium with the carbonates in the whole rocks. delta C-13 values of CO, are variable but too low to be in equilibrium with the wall rock fluids, compatible with a contribution of CO2 from closed system oxidation of methane. Differences in isotopic composition between host-rock and Alpine fissure carbonate are generally small, suggesting that the amount of CO2 produced by oxidation of methane was small compared to the C-budget in the rocks and local pore fluids were buffered by the wall rocks during precipitation of calcite within the fissures. (c) 2006 Elsevier B.V. All rights reserved.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Clastic and core lava components of a silicic lava dome

The formation of a lava dome involves fractionation of the lava into core and clastic components. We show that for three separate, successive andesitic lava domes that grew at Soufrière Hills volcano, Montserrat, between 1999 and 2007, the volumetric proportion of the lava converted to talus or pyroclastic flow deposits was 50%–90% of the lava extruded. Currently, only 8% of the total magma extruded during the 1995–2007 eruption remains as core lava. The equivalent representation in the geological record will probably be even lower. Most of the lava extruded at the surface flowed no further than 150–300 m from the vent before disaggregation, resulting in a lava core whose shape tends to a cylinder. Moderate to high extrusion rates at the Soufrière Hills domes may have contributed to the large clastic fraction observed. Creating talus dissipates much of the energy that would otherwise be stored in the core lava of domes. The extreme hazards from large pyroclastic flows and blasts posed by wholesale collapse of a lava dome depend largely on the size of the lava core, and hence on the aggregate history of the partitioning process, not on the size of the dome.

Investigation of the prototype silylene reaction, SiH2+H2O(and D2O): time-resolved gas-phase kinetic studies, isotope effects, RRKM calculations, and quantum chemical calculations of the reaction energy surface

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 296 and at 339 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied over the pressure range 10-200 Torr with SF6 as bath gas. The second-order rate constants obtained were pressure dependent, indicating that the reaction is a third-body assisted association process. Rate constants at 339 K were about half those at 296 K. Isotope effects, k(H)/k(D), were small averaging 1.076 0.080, suggesting no involvement of H- (or D-) atom transfer in the rate determining step. RRKM modeling was undertaken based on a transition state appropriate to formation of the expected zwitterionic donoracceptor complex, H2Si...OH2. Because the reaction is close to the low pressure (third order) region, it is difficult to be definitive about the activated complex structure. Various structures were tried, both with and without the incorporation of rotational modes, leading to values for the high-pressure limiting (i.e., true secondorder) rate constant in the range 9.5 x 10(-11) to 5 x 10(-10) cm(3) molecule' s(-1). The RRKM modeling and mechanistic interpretation is supported by ab initio quantum calculations carried out at the G2 and G3 levels. The results are compared and contrasted with the previous studies.

Synthesis, crystal structure, magnetic behavior and thermal property of three polynuclear complexes: [M(dca)(2)(H2O)2](n)center dot(hmt)(n) [M = Mn(II), Co(II)] and [Co(dca)(2)(bpds)](n) [dca, dicyanamide; hint, hexamethylenetetramine; bpds, 4,4 '-bipyridyl disulfide]

Three mu(1.5)-dicyanamide bridged Mn(II) and Co(II) complexes having molecular formula [Mn(dca)(2)(H2O)(2)](n)center dot(hmt)(n) (1), [Co(dca)(2) (H2O)(2)](n)center dot(hmt)(n) (2) and [Co(dca)(2)(bpds)](n) (3) [dca = dicyanamide; hmt = hexamethylenetetramine; bpds = 4,4'-bipyridyl disulfide] have been synthesized and characterized by single crystal X-ray diffraction study, low temperature (300-2 K) magnetic measurement and thermal behavior. The X-ray diffraction analysis of 1 and 2 reveals that they are isostructural, comprising of 1D coordination polymers [M(dca)(2)(H2O)(2)](n) [M = Mn(II), Co(II) for 1 and 2. respectively] with uncoordinated hmt molecules located among the chains. The [M(dca)(2)(H2O)(2)](n) chains and the lattice hint molecules are connected through H-bonds resulting in a 3D supramolecular architecture. The octahedral N4O2 chromophore surrounding the metal ion forms via two trans located water oxygens and four nitrogens from four nitrile dca. Complex 3 is a 1D chain formed by two mu(1.5)-dca and one bridging bpds. The octahedral N-6 coordination sphere surrounding the cobalt ions comprises four nitrogens from dca and two from bpds. Low temperature magnetic study indicates small antiferromagnetic coupling for all the complexes. Best fit parameters for 1: J = -0.17 cm(-1), g = -2.03 with R = 6.1 x 10(-4), for 2, J = -0.50 cm(-1), and for 3, J = -0.95 cm(-1). (c) 2005 Elsevier B.V. All rights reserved.