Ilmenite composition in the Tellnes Fe–Ti deposit,SW Norway: fractional crystallization,postcumulus evolution and ilmenite–zircon relation

Major and trace element XRF and in situ LA-ICP-MS analyses of ilmenite in the Tellnes ilmenite deposit, Rogaland Anorthosite Province, SW Norway, constrains a two stage fractional crystallization model of a ferrodioritic Fe-Ti-P rich melt. Stage 1 is characterized by ilmenite-plagioclase cumulates, partly stored in the lower part of the ore body (Lower Central Zone, LCZ), and stage 2 by ilmenite-plagioclase-orthopyroxene-olivine cumulates (Upper Central Zone, UCZ). The concentration of V and Cr in ilmenite, corrected for the trapped liquid effect, (1) defines the cotectic proportion of ilmenite to be 17.5 wt% during stage 1, and (2) implies an increase of D _V^Ilm during stage 2, most likely related to a shift in fO₂. The proportion of 17.5 wt% is lower than the modal proportion of ilmenite (ca. 50 wt%) in the ore body, implying accumulation of ilmenite and flotation of plagioclase. The fraction of residual liquid left after crystallization of Tellnes cumulates is estimated at 0.6 and the flotation of plagioclase at 26 wt% of the initial melt mass. The increasing content of intercumulus magnetite with stratigraphic height, from 0 to ca. 3 wt%, results from differentiation of the trapped liquid towards magnetite saturation. The MgO content of ilmenite (1.4–4.4 wt%) is much lower than the expected cumulus composition. It shows extensive postcumulus re-equilibration with trapped liquid and ferromagnesian silicates, correlated with distance to the host anorthosite. The Zr content of ilmenite, provided by in situ analyses, is low (<114 ppm) and uncorrelated with stratigraphy or Cr content. The data demonstrate that zircon coronas observed around ilmenite formed by subsolidus exsolution of ZrO₂ from ilmenite. The U-Pb zircon age of 920 ± 3 Ma probably records this exsolution process. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

An analytical solution to multi-component NAPL dissolution equations

This note presents a novel method for determining the changing composition of a multi-component NAPL body dissolving into moving groundwater, and the consequent changes in the aqueous phase solute concentrations in the surrounding pore water. A canonical system of coupled non-linear governing equations is derived which is suitable for representation of both pooled and residual configurations, and this is solved. Whereas previous authors have handled such problems numerically, it is shown that these governing equations succumb to analytical solution. By a suitable substitution, the equations become decoupled, and the problem collapses to a single first-order equation. The final result is expressed implicitly, with time as a function of the number of moles of the least soluble component, m₁. The number of moles of each other component is expressed explicitly in terms of m₁. It is shown that the time-m₁ relationship has a well behaved inverse. An example is given in which the analytic solution is verified against traditional finite difference analysis, and its computational efficiency is shown.     

Field‐based and spectral indicators for soil erosion mapping in semi‐arid mediterranean environments (Coastal Cordillera of central Chile)

The Coastal Cordillera of central Chile is naturally sensitive to soil erosion due to moderate to steep slopes, intense winter rains when the vegetation cover is scarce, and deeply weathered granitic rocks. In 1965, 60 per cent of its surface was moderately to very severely eroded. Today this process is still largely active, but no data are currently available to evaluate the real extent, distribution and severity of soil degradation on a regional scale. This information is vital to support efficient soil conservation plans. A multi‐scale approach was implemented to produce regional land degradation maps based on remote sensing technologies. Fieldwork has shown that the surface colour or ‘redness’ and the density of coarse fragments are pertinent erosion indicators to describe a typical sequence of soil degradation in the context of mediterranean soil developed on granitic materials and micaschists. Field radiometric experiments concluded that both factors influence the reflectance of natural surfaces and can be modelled using radiometric indices accessible from most satellites operating in the optical domain, i.e. redness index and brightness index. Finally the radiometric indices were successfully applied to SPOT images to produce land degradation maps. Only broad classes of erosion status were discriminated and the detection of the degradation processes was only possible when most of the fertile layer had already been removed. This technology provides decision‐making information required to develop regional soil conservation plans and to prioritize actions between catchment areas, especially in vast inter‐tropical regions where spatialized data are not always readily available. Copyright © 2006 John Wiley & Sons, Ltd.     

Practical method for drawing a VGP path

A statistical method is proposed for the smoothing of polar wander paths and for giving their confidence limits by the mean of successive ellipses. The method is fully parametrical and is based on the relations between the inertia matrix and the parameters of the Fisher distribution, from which a bivariate form is deduced. An elementary tensorial calculation gives the parameters of the confidence ellipse around a vectorial weighted mean, for any unimodal set of vectors. This model can also be used for other statistical tests, wherever the rotational symmetry hypothesis is not consistent (i.e. the fold test).     

The product spectra of gravity and barometric pressure in Europe

A total of 111 000 hourly values of gravity and barometric pressure from stations in Europe is analysed. The data consist of two sets of records from Brussels, an early set of 36 000 h length and a more recent set of 21 000 h length, a set of records from Bad Homburg of 24 000 h length, and a set of records from Strasbourg of 30 000 h length. All of the gravity measurements were made with similar superconducting instruments and the pressure data were recorded simultaneously at each superconducting gravimeter site. The four sets of records have different time bases, and to bring out common features and suppress individual station systematic errors, the product spectrum is introduced. Spectral density estimates are first computed for a common spectral window for each record, and the product spectrum is formed by multiplying individual spectral estimates across records. The cumulative distribution function is found for the product spectrum and confidence intervals are calculated from it by iteration. The product spectrum in gravity reveals a triplet of resonances in the subtidal band which are shown by an automated computer search to be uniquely associated with the translational modes of the solid inner core. The product spectrum in barometric pressure clearly reveals the first 10 solar heating tides in the atmosphere, but otherwise does not show common features with the gravity product spectrum. In particular, the triplet of resonances in the subtidal band of the gravity product spectrum do not show up in the product spectrum of barometric pressure, climinating the atmosphere as their source.     

A CME mass distribution derived from SOLWIND coronagraph observations

Using estimates of the masses of nearly 1000 CMEs observed by SOLWIND from Howardet al. (1985), we re-plot the numbers of CMEs as a function of CME mass on a log-linear plot. The plot is significant in that it shows a linear trend over more than a decade of CME masses. The plot indicates a simple form for the distribution of the CME masses and allows an easy determination of the total mass ejected into the solar wind in the form of CMEs. We find that approximately 16% of the solar wind at solar maximum can be comprised of CME mass. There is no indication that the numbers of low-mass CMEs increase in number above the trend set by the more massive ones. Specifically, there is no increase in the numbers of small CMEs such that the whole of the solar wind can be comprised of them.     

Brest sea level record: a time series construction back to the early eighteenth century

The completeness and the accuracy of the Brest sea level time series dating from 1807 make it suitable for long-term sea level trend studies. New data sets were recently discovered in the form of handwritten tabulations, including several decades of the eighteenth century. Sea level observations have been made in Brest since 1679. This paper presents the historical data sets which have been assembled so far. These data sets span approximately 300 years and together constitute the longest, near-continuous set of sea level information in France. However, an important question arises: Can we relate the past and the present-day records? We partially provide an answer to this question by analysing the documents of several historical libraries with the tidal data using a ‘data archaeology’ approach advocated by Woodworth (Geophys Res Lett 26:1589–1592, 1999b). A second question arises concerning the accuracy of such records. Careful editing was undertaken by examining the residuals between tidal predictions and observations. It proved useful to remove the worst effects of timing errors, in particular the sundial correction to be applied prior to August 1, 1714. A refined correction based on sundial literature [Savoie, La gnomique, Editions Les Belles Lettres, Paris, 2001] is proposed, which eliminates the systematic offsets seen in the discrepancies in timing of the sea level measurements. The tidal analysis has also shown that shallow-water tidal harmonics at Brest causes a systematic difference of 0.023 m between mean sea level (MSL) and mean tide level (MTL). Thus, MTL should not be mixed with the time series of MSL because of this systematic offset. The study of the trends in MTL and MSL however indicates that MTL can be used as a proxy for MSL. Three linear trend periods are distinguished in the Brest MTL time series over the period 1807–2004. Our results support the recent findings of Holgate and Woodworth (Geophys Res Lett) of an enhanced coastal sea level rise during the last decade compared to the global estimations of about 1.8 mm/year over longer periods (Douglas, J Geophys Res 96:6981–6992, 1991). The onset of the relatively large global sea level trends observed in the twentieth century is an important question in the science of climate change. Our findings point out to an ‘inflexion point’ at around 1890, which is remarkably close to that in 1880 found in the Liverpool record by Woodworth (Geophys Res Lett 26:1589–1592, 1999b).     

Bubble-induced porewater mixing: A 3-D model for deep porewater irrigation

Porewater data from vent sites of the northeastern shelf off Sakhalin Island, Sea of Okhotsk, exhibit bottom-water concentrations down to a sediment depth of up to 300 cm. Below this depth, solute concentrations rapidly change due to methanogenesis and anaerobic methane oxidation (AMO). The profile shapes suggest an irrigation-like process that mixes on a meter scale. At these sites active gas emanation into the overlying water column and near-surface gas hydrates are commonly observed. We propose that methane gas bubbles rise through the soft surface sediments and cause mixing of the porewater.Mathematically, the bubble-induced irrigation can be described by eddy diffusion enhancing the diffusive transport of solutes by several orders of magnitude. A 3-D numerical transport-reaction model was developed to investigate the parameters defining the mixing process, such as bubble rise velocity, tube size, tube distribution in the sediment, and ebullition frequency.Model consistency with the field data requires eddy diffusivities ?1 × 10⁵ cm²/a, tube densities of >4 tubes/m² (equivalent to a tube spacing of <40 cm), active gas seepage for more than a few weeks or months, and moderate to low diagenetic reaction rates of solutes. The corresponding methane gas fluxes that are predicted from the results of the model realizations range from 1 × 10³-5 × 10⁵ L/(m² a). Due to bubble mixing, solute fluxes in these sediments are increased by a factor of 3 and the maximum AMO rate by a factor of 7.