Thermodynamic controls on element partitioning between titanomagnetite and andesitic–dacitic silicate melts

Titanomagnetite–melt partitioning of Mg, Mn, Al, Ti, Sc, V, Co, Ni, Cu, Zn, Ga, Zr, Nb, Mo, Hf and Ta was investigated experimentally as a function of oxygen fugacity (fO₂) and temperature (T) in an andesitic–dacitic bulk-chemical compositional range. In these bulk systems, at constant T, there are strong increases in the titanomagnetite–melt partitioning of the divalent cations (Mg²⁺, Mn²⁺, Co²⁺, Ni²⁺, Zn²⁺) and Cu²⁺/Cu⁺ with increasing fO₂ between 0.2 and 3.7 log units above the fayalite–magnetite–quartz buffer. This is attributed to a coupling between magnetite crystallisation and melt composition. Although melt structure has been invoked to explain the patterns of mineral–melt partitioning of divalent cations, a more rigorous justification of magnetite–melt partitioning can be derived from thermodynamic principles, which accounts for much of the supposed influence ascribed to melt structure. The presence of magnetite-rich spinel in equilibrium with melt over a range of fO₂ implies a reciprocal relationship between a(Fe²⁺O) and a(Fe³⁺O_1.5) in the melt. We show that this relationship accounts for the observed dependence of titanomagnetite–melt partitioning of divalent cations with fO₂ in magnetite-rich spinel. As a result of this, titanomagnetite–melt partitioning of divalent cations is indirectly sensitive to changes in fO₂ in silicic, but less so in mafic bulk systems.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

Areal extent of the 1988 flood in Bangladesh: How much did the satellite imagery show?

This paper presents the methods and results of visual interpretation of satellite imagery for estimating areal extent of the 1988 flood in Bangladesh. The main limitation of this method was that cloud-free images were not available for the peak flood period in the first week of September. Relatively cloud-free images on three selected dates in the last three weeks of September showed that areas of inundation on those dates ranged from 44000 to 60000 km² ± 5% (31 to 42% of Bangladesh). The Government of Bangladesh official estimates of the areal extent of flooding were tentative in nature, since these were based on a choropleth map of relative proportions (percentages) of flooding in different parts of the country. The official estimate of the maximum extent of flooding was 82000 km² (57% of Bangladesh), i.e. 15 to 26% larger than the area shown on satellite imagery. The actual extent of flooding was certainly larger than the estimates from satellite imagery but, perhaps, smaller than the official estimates, since a reduction of 15 to 26% of flood area in one to three weeks seemed unlikely. An alternative method of mapping flood-affected areas by using newspaper-interpreted data was attempted, but the method had limited values because of reporting bias.     

Acid mine drainage at Cerro Rico de Potosí I: unabated high-strength discharges reflect a five century legacy of mining

Intensive mining and processing of Ag, Sn, Pb and Zn ores have occurred in various locations within and around the city of Potosí, Bolivia since 1545. Surface and subsurface waters, stream sediments and soils are contaminated with various ecotoxic metals in the headwaters of the economically vital, yet highly impacted, upper Rio Pilcomayo watershed. Previous studies have documented downstream trace metal contamination, however, not addressed specific sources. The AMD discharges identified in this study help link downstream pollution to primary origins. The majority of AMD would be considered high-strength due to metal concentrations and acidity orders of magnitude greater than typical AMD. Discharges from both operating and abandoned portals as well as tailings-related deposits displayed a high degree of heterogeneity with total metal concentrations ranging from 0.11 to 7480, <0.022 to 889, <0.0006 to 65.3, <0.001 to 310, 0.12 to 72,100, 0.3 to 402, <0.012 to 34.8, and 0.24 to 19,600 mg/L of Al, As, Cd, Cu, Fe, Mn, Pb and Zn, respectively. Net acidity and pH ranged from −10 to 246,000 mg/L as CaCO₃ equivalent and 0.90–6.94 standard units, respectively. Data were gathered during two sampling events centered around the most extreme periods of the dry and wet seasons of one water-year. Loadings to local streams were marginally greater for most metals in the wet season. If observed loadings are historically representative, Cerro Rico AMD has contributed thousands of tonnes of ecotoxic metals to the upper Rio Pilcomayo over the last five centuries. Metal and hydrogen ion concentrations in the majority of AMD sampled were several orders of magnitude above discharge limits set by the Bolivian government, yet no action has historically or contemporarily been taken.     

Occurrence and Distribution of Organic Contaminants in the Aquatic System in Berlin.Part II: Substituted Phenols in Berlin Surface Water

In a surface water screening, 30 representative surface water samples collected from rivers, canals, and lakes in Berlin were investigated for the presence of 22 substituted phenols. The phenols selected include the 11 phenols considered as “priority pollutants” by the US Environmental Protection Agency (US-EPA). Surface water samples were extracted applying solid-phase extraction with styrenedivinylbenzene adsorbent. The recoveries, determined in spiking experiments, were between 80 % and 103 %. After derivatization with N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBSTFA) the samples were analyzed by capillary gas chromatography-mass spectrometry (GC-MS). Phenol, cresols, 2-ethylphenol, 2-chlorophenol, 4-chloro-3-methylphenol, pentachlorophenol, 2-nitrophenol, and 4-nitrophenol were detected in the surface water samples at concentrations between 0.02 μg/L and 7.8 μg/L, respectively. The distribution of these residues in the Berlin surface waters showed that the phenolic residues, with the exception of pentachlorophenol and 2-ethylphenol, do not originate primarily from municipal sewage treatment plants discharges. Some of the phenols are formed naturally or occur as ubiquitous anthropogenic contaminants in the aquatic system.     

Diffusive equilibrium in thin films provides evidence of suppression of hyporheic exchange and large‐scale nitrate transformation in a groundwater‐fed river

The hyporheic zone of riverbed sediments has the potential to attenuate nitrate from upwelling, polluted groundwater. However, the coarse‐scale (5–10 cm) measurement of nitrogen biogeochemistry in the hyporheic zone can often mask fine‐scale (<1 cm) biogeochemical patterns, especially in near‐surface sediments, leading to incomplete or inaccurate representation of the capacity of the hyporheic zone to transform upwelling NO₃^?. In this study, we utilised diffusive equilibrium in thin‐films samplers to capture high resolution (cm‐scale) vertical concentration profiles of NO₃^?, SO₄^2?, Fe and Mn in the upper 15 cm of armoured and permeable riverbed sediments. The goal was to test whether nitrate attenuation was occurring in a sub‐reach characterised by strong vertical (upwelling) water fluxes. The vertical concentration profiles obtained from diffusive equilibrium in thin‐films samplers indicate considerable cm‐scale variability in NO₃^? (4.4 ± 2.9 mg N/L), SO₄^2? (9.9 ± 3.1 mg/l) and dissolved Fe (1.6 ± 2.1 mg/l) and Mn (0.2 ± 0.2 mg/l). However, the overall trend suggests the absence of substantial net chemical transformations and surface‐subsurface water mixing in the shallow sediments of our sub‐reach under baseflow conditions. The significance of this is that upwelling NO₃^?‐rich groundwater does not appear to be attenuated in the riverbed sediments at <15 cm depth as might occur where hyporheic exchange flows deliver organic matter to the sediments for metabolic processes. It would appear that the chemical patterns observed in the shallow sediments of our sub‐reach are not controlled exclusively by redox processes and/or hyporheic exchange flows. Deeper‐seated groundwater fluxes and hydro‐stratigraphy may be additional important drivers of chemical patterns in the shallow sediments of our study sub‐reach. © 2015 The Authors. Hydrological Processes Published by John Wiley & Sons Ltd.     

Long-range weather forecasts through numerical and empirical methods

There has been much improvement in numerical weather prediction since L.F. Richardson (1922, Weather Prediction by Numerical Process, Cambridge University Press, Cambridge, p. 236) wrote his famous book. NWP has primarily been successful in improving day-by-day forecasts starting from an observed detailed Initial Condition (IC) out to about a week. The purpose of this paper is to discuss first the state of the art in long-range NWP by presenting results of a new large numerical experiment (named DERF90; from Dynamical Extended Range Forecasting in 1990 out to 90 days) conducted at the National Meteorological Center (NMC) during the summer and autumn of 1990 (Section 2). One hundred and twenty eight 90-day global forecasts were made from successive daily initial conditions (IC), thus giving us ample opportunity to assess skill of forecasts at lead times beyond 1 week.We then move on to define the notion of a limit of predictability (LOP), and following a procedure by Lorenz (1982), give a numerical estimate of the LOP using the DERF90 data set. We then produce a list of reasons, as to why this estimate (LOP = 18 days) should not be taken too literally. In particular, we argue that the LOP varies as a function of the flow itself, and it would be (much) larger if we had, as we will ultimately, a coupled ocean-atmosphere model for making long-lead forecasts. Last, but not least, we present results of empirical forecasts that point to modest but significant skill well beyond the traditional LOP (a few weeks).A specific recent example of empirical forecasting is discussed. Through Canonical Correlation Analysis (CCA), experimental forecasts are being made for the United States surface temperatures at lead times of several seasons. While modest, the skill is significant in that it defies the existence or a 3-week LOP, and so demonstrates the potential for model improvements.     

Comparison of the LEW88516 and ALHA77005 martian meteorites: Similar but distinct

Abstract By mineral and bulk compositions, the Lewis Cliff (LEW) 88516 meteorite is quite similar to the ALHA77005 martian meteorite. These two meteorites are not paired because their mineral compositions are distinct, they were found 500 km apart in ice fields with different sources for meteorites, and their terrestrial residence ages are different. Minerals in LEW88516 include: olivine, pyroxenes (low- and high-Ca), and maskelynite (after plagioclase); and the minor minerals chromite, whitlockite, ilmenite, and pyrrhotite. Mineral grains in LEW88516 range up to a few mm. Texturally, the meteorite is complex, with regions of olivine and chromite poikilitically enclosed in pyroxene, regions of interstitial basaltic texture, and glass-rich (shock) veinlets. Olivine compositions range from Fo₆₄ to Fo₇₀, (avg. Fo₆₇), more ferroan and with more variation than in ALHA77005 (Fo₆₉ to Fo₇₃). Pyroxene compositions fall between En₇₇Wo₄ and En₆₅Wo₁₅ and in clusters near En₆₃Wo₉ and En₅₃Wo₃₃, on average more magnesian and with more variation than in ALHA77005. Shock features in LEW88516 range from weak deformation through complete melting. Bulk chemical analyses by modal recombination of electron microprobe analyses, instrumental neutron activation, and radiochemical neutron activation confirm that LEW88516 is more closely related to ALHA77005 than to other known martian meteorites. Key element abundance ratios are typical of martian meteorites, as is its non-chondritic rare earth pattern. Differences between the chemical compositions of LEW88516 and ALHA77005 are consistent with slight differences in the proportions of their constituent minerals and not from fundamental petrogenetic differences. Noble gas abundances in LEW88516, like those in ALHA77005, show modest excesses of ⁴⁰Ar and ¹²⁹Xe from trapped (shock-implanted) gas. As with other ALHA77005 and the shergottite martian meteorites (except EETA79001), noble gas isotope abundances in LEW88516 are consistent with exposure to cosmic rays for 2.5–3 Ma. The absence of substantial effects of shielding from cosmic rays suggest LEW88516 spent this time as an object no larger than a few cm in diameter.     

Petrogenesis and paleotectonic history of the Wild Bight Group,an Ordovician rifted island arc in central Newfoundland

The Wild Bight Group (WBG) is a sequence of early and middle Ordovician volcanic, subvolcanic and epiclastic rocks, part of the Dunnage Tectonostratigraphic Zone of the Newfoundland Appalachians. A detailed geochemical and Nd-isotopic study of the volcanic and subvolcanic rocks has been carried out to determine the geochemical characteristics of the rocks, interpret their palcotectonic environments and constrain their petrogenetic history. The lower and central stratigraphic levels of the WBG contain mafic volcanic rocks with island-arc geochemical signatures, including LREE-enriched are tholeiites with _Nd(t) =-0.1 to +2.2 (type A-I), LREE-depleted arc tholeiites with _Nd(t) =+5.6 to +7.1 (type A-II) and an unusual suite of strongly incompatible-element depleted tholeiites in which _Nd(t) ranges from-0.9 to +4.6 and is negatively correlated with¹⁴⁷Sm/¹⁴⁴Nd (type A-III). High-silica, low-K rhyolites occur locally in the central part of the stratigraphy, associated with mafic rocks of arc affinity, and have _Nd(t) =+4.7 to +5.4. The upper stratigraphic levels of the WBG dominantly contain rocks with non-arc geochemical signatures, including alkalic basalts with _Nd(t) =+4.6 to +5.5 (type N-I), strongly LREE- and incompatible element-enriched tholeiites that are transitional between alkalic and non-alkalic rocks with _Nd(t) =+4.4 to +7.0 (type N-II) and rocks with flat to slightly LREE-enriched patterns and _Nd(t) =+5.1 to +7.4 (type N-III). Rocks with non-arc and arc signatures are locally interbedded near the stratigraphic type of the WBG. Nd-isotopic data in the type A-I and A-II rocks are generally compatible with mixing/partial melting models involving depleted mantle, variably contaminated by a subducted crustally-derived sediment. The petrogenesis of type A-III rocks must involve source mixing and multi-stage partial melting, but the details are not clear. The geochemistry and Nd isotope data for types N-I, N-II and N-III rocks are compatible with petrogenetic models involving variable partial melting of a source similar to that postulated for modern oceanic island basalts. Comparison of the WBG with modern analogues suggests a 3-stage developmental model: stage 1) island-arc volcanism (eruption of type mafic volcancs); stage 2) arc-rifting (continued eruption of type A-I, A-I, eruption of types A-II and A-III mafic volcanics and high-silica, low-K rhyolites); and stage 3) back-arc basin volcanism (continued minor eruption of type A-I basalts, eruption of types N-I, N-II, N-III basalts). Stages 1 and 2 volcanism involved partial melting of subduction contaminated mantle, while stage 3 volcanism utilized depleted-mantle sources not affected by the subducting slab. This model provides a basis for interpreting coeval sequences in central Newfoundland and a comparative framework for some early Paleozoic oceanic volcanic sequences elsewhere in the Appalachian orogen.