Accuracy assessment of lidar-derived digital elevation models

Despite the relatively high cost of airborne lidar-derived digital elevation models (DEMs), such products are usually presented without a satisfactory associated estimate of accuracy. For the most part, DEM accuracy estimates are typically provided by comparing lidar heights against a finite sample of check point coordinates from an independent source of higher accuracy, supposing a normal distribution of the derived height differences or errors. This paper proposes a new methodology to assess the vertical accuracy of lidar DEMs using confidence intervals constructed from a finite sample of errors computed at check points. A non-parametric approach has been tested where no particular error distribution is assumed, making the proposed methodology especially applicable to non-normal error distributions of the type usually found in DEMs derived from lidar. The performance of the proposed model was experimentally validated using Monte Carlo simulation on 18 vertical error data-sets. Fifteen of these data-sets were computed from original lidar data provided by the International Society for Photogrammetry and Remote Sensing Working Group III/3, using their respective filtered reference data as ground truth. The three remaining data-sets were provided by the Natural Environment Research Council's Airborne Research and Survey Facility lidar system, together with check points acquired using high precision kinematic GPS. The results proved promising, the proposed models reproducing the statistical behaviour of vertical errors of lidar using a favourable number of check points, even in the cases of data-sets with non-normally distributed residuals. This research can therefore be considered as a potentially important step towards improving the quality control of lidar-derived DEMs.     

The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.     

The age and source of late Hercynian magmatism in the central Alps: evidence from precise U−Pb ages and initial Hf isotopes

This study presents U–Pb ages for zircon, titanite, allanite and epidote, and initial Hf isotopic compositions for zircon of Upper Carboniferous granites, diorites and syenites from the Aar massif, central Alps. The rocks were emplaced during three magmatic pulses after Hercynian collisional tectonics: (A) a shoshonitic-ultrapotassic series at 334±2.5 Ma; (B) scattered diorites and granites at 308–310 Ma; and (C) a high-K cale-alkaline granite batholith at 298±2 Ma. Inheritance of old zircons is negligible among all three groups. The Southern Aar granite, in contrast, is a syn-tectonic, probably ca. 350 Ma old granite that contains large amounts of inherited Precambrian zircons. Alpine metamorphism caused weak lead loss in many analyzed zircon fractions, but left the titanite U–Pb system undisturbed: thorites were almost completely reset by Alpine and recent lead loss. Mineral isochrons defined by titanite, allanite, epidote and apatite yield initial Pb isotopic compositions that are in agreement with the model values of Stacey and Kramers. Initial Hf isotopic compositions range from _Hf=–8 to +3.5. The data follow a trend of increasing _Hf with decreasing age. The _Hf versus element concentration relationships suggest mixing between a mantle and a crustal component. These relationships can be explained in terms of generation of the melts from a subcontinental mantle that had been enriched during subduction events at about 1 Ga and by 300 Ma had developed an isotopic signature distinct from that of MORB-type mantle. Further contamination of the melts occurred during ascent and differentiation in the crust. This late Hercynian magmatism can be related to post-collisional strike-slip tectonics.     

The behavior of trace elements during schwertmannite precipitation and subsequent transformation into goethite and jarosite

Schwertmannite is a ubiquitous mineral formed from acid rock drainage (ARD), and plays a major role in controlling the water chemistry of many acid streams. The formation of schwertmannite was investigated in the acid discharge of the Monte Romero abandoned mine (Iberian Pyrite Belt, SW, Spain). Schwertmannite precipitated from supersaturated solutions mainly owing to the oxidation of Fe(II) to Fe(III) and transformed with time into goethite and jarosite. In a few hours, schwertmannite precipitation removed more than half of the arsenic load from solution, whereas the concentration of divalent trace metals (Zn, Cu, Pb, Cd, Ni, and Co) remained almost unchanged. In the laboratory, natural schwertmannite was kept in contact with its coexisting acid water in a flask with a solid-liquid mass ratio of 1:5 for 353 days. During this time, the pH of the solution dropped from 3.07 to 1.74 and the concentrations of sulfate and Fe increased. During the first 164 days, schwertmannite transformed into goethite plus H₃O-jarosite but, subsequently, goethite was the only mineral to form. Some of the trace elements, such as Al, Cu, Pb, and As were depleted in solution during the first stage as schwertmannite transformed into goethite plus H₃O-jarosite. On the contrary, the transformation of schwertmannite to goethite (with no jarosite) during the second stage released Al, Cu, and As to the solution. Despite the variation in their concentrations in solution, approximately 80% of the total Al and Cu inventories and more than 99% As and Pb remained in the solid phase throughout the entire aging process.     

HP–LT Variscan metamorphism in the Cubito-Moura schists (Ossa-Morena Zone, southern Iberia)

Multi-equilibrium thermobarometry shows that low-grade metapelites (Cubito-Moura schists) from the Ossa–Morena Zone underwent HP–LT metamorphism from 340–370 °C at 1.0–0.9 GPa to 400–450 °C at 0.8–0.7 GPa. These HP–LT equilibriums were reached by parageneses including white K mica, chlorite and chloritoid, which define the earliest schistosity (S₁) in these rocks. The main foliation in the schists is a crenulation cleavage (S₂), which developed during decompression from 0.8–0.7 to 0.4–0.3 GPa at increasing temperatures from 400–450 °C to 440–465 °C. Fe³⁺ in chlorite decreased greatly during prograde metamorphism from molar fractions of 0.4 determined in syn-S₁ chlorites down to 0.1 in syn-S₂ chlorites. These new data add to previous findings of eclogites in the Moura schists indicating that a pile of allochtonous rocks situated next to the Beja-Acebuches oceanic amphibolites underwent HP–LT metamorphism during the Variscan orogeny. To cite this article: G. Booth-Rea et al., C. R. Geoscience 338 (2006).     

Computer simulation studies of Crystallization under confinement conditions

Crystallization under confinement conditions is a very important process in geochemistry and geophysics. Computer simulations of fluids in nanometer scale pore spaces can provide a unique microscopic insight into the structure, dynamics and forces arising from the crystallization process. We discuss in this paper molecular dynamics computer simulations of crystallization in pores of nanometer dimensions. The crystallization pressure due to the freezing of a model of Argon in a nanopore is computed using molecular dynamics simulations. We also investigate the influence of pore geometry in determining the dynamics of confined fluids, as well as mass separation in binary mixtures. It turns out that the pore geometry reveals itself as an important variable, leading to 1) new mechanisms for fast diffusion in confined spaces, and 2) accumulation of solute in specific regions inside the pore.     

U–Pb ages of the Dalsfjord Complex, SW Norway, and their bearing on the correlation of allochthonous crystalline segments of the Scandinavian Caledonides

The Dalsfjord Complex is an allochthonous crystalline unit in the hanging wall of a major extensional detachment of the Caledonian orogen on the western coast of Norway. U-Pb geochronology has been carried out to test a previously proposed correlation between this unit and the Jotun Nappe Complex, a member of the Middle Allochthon now occupying a foreland position to the east. A monzonite sample of the Dalsfjord suite yields a zircon age of 1,634Dž Ma, whereas a crosscutting gabbro was formed at 1,464Lj Ma. Both rocks were strongly overprinted during the Sveconorwegian orogeny, especially the monzonite whose zircons U-Pb data were strongly pulled towards a lower intercept age of 882ᆱ Ma, and whose titanite indicates a two-stage growth history at >960 Ma and <920 Ma. These relationships support the correlation of these units with the Jotun Complex, and to some degree also with the Lindås Nappe in the Bergen Arcs.