The sandstone-hosted Beverley uranium Deposit, Lake Frome Basin, South Australia: Mineralogy, geochemistry, and a time-constrained model for its genesis

The sandstone-hosted Beverley uranium deposit is located in terrestrial sediments in the Lake Frome basin in the North Flinders Ranges, South Australia. The deposit is 13 km from the U-rich Mesoproterozoic basement of the Mount Painter inlier, which is being uplifted 100 to 200 m above the basin by neotectonic activity that probably initiated in the early Pliocene. The mineralization was deposited mainly in organic matter-poor Miocene lacustrine sands and partly in the underlying reductive strata comprising organic matter-rich clays and silts. The bulk of the mineralization consists of coffinite and/or uraninite nodules, growing around Co-rich pyrite with an S isotope composition (delta S-34 = 1.0 +/- 0.3 parts per thousand), suggestive of an early diagenetic lacustrine origin. In contrast, authigenic sulfides in the bulk of the sediments have a negative S isotope signature (delta S-34 ranges from -26.2 to -35.5 parts per thousand), indicative of an origin via bacterially mediated sulfate reduction. Minor amounts of Zn-bearing native copper and native lead also support the presence of specific, reducing microenvironments in the ore zone. Small amounts of carnotite are associated with the coffinite ore and also occur beneath a paleosoil horizon overlying the uranium deposit. Provenance studies suggest that the host Miocene sediments were derived from the reworking of Early Cretaceous glacial or glaciolacustrine sediments ultimately derived from Paleozoic terranes in eastern Australia. In contrast, the overlying Pliocene strata were in part derived from the Mesoproterozoic basement inlier. Mass-balance and geochemical data confirm that granites of the Mount Painter domain were the ultimate source of U and BEE at Beverley. U-Pb dating of coffinite and carnotite suggest that the U mineralization is Pliocene (6.7-3.4 Ma). The suitability of the Beverley deposit for efficient mining via in situ leaching, and hence its economic value, are determined by the nature of the hosting sand unit, which provides the permeability and low reactivity required for high fluid flow and low chemical consumption. These favorable sedimentologic and geometrical features result from a complex conjunction of factors, including deposition in lacustrine shore environment, reworking of angular sands of glacial origin, deep Pliocene weathering, and proximity to an active fault exposing extremely U rich rocks.

3-D density structure and geological evolution of Stromboli volcano (Aeolian Islands, Italy) inferred from land-based and sea-surface gravity data

We present the first density model of Stromboli volcano (Aeolian Islands, Italy) obtained by simultaneously inverting land-based (543) and sea-surface (327) relative gravity data. Modern positioning technology, a 1 x 1 m digital elevation model, and a 15 x 15 m bathymetric model made it possible to obtain a detailed 3-D density model through an iteratively reweighted smoothness-constrained least-squares inversion that explained the land-based gravity data to 0.09 mGal and the sea-surface data to 5 mGal. Our inverse formulation avoids introducing any assumptions about density magnitudes. At 125 m depth from the land surface, the inferred mean density of the island is 2380 kg m(-3), with corresponding 2.5 and 97.5 percentiles of 2200 and 2530 kg m-3. This density range covers the rock densities of new and previously published samples of Paleostromboli I, Vancori, Neostromboli and San Bartolo lava flows. High-density anomalies in the central and southern part of the island can be related to two main degassing faults crossing the island (N41 and NM) that are interpreted as preferential regions of dyke intrusions. In addition, two low-density anomalies are found in the northeastern part and in the summit area of the island. These anomalies seem to be geographically related with past paroxysmal explosive phreato-magmatic events that have played important roles in the evolution of Stromboli Island by forming the Scari caldera and the Neostromboli crater, respectively. (C) 2014 Elsevier B.V. All rights reserved.

Evaluation of the reconstruction limits of a frequency-independent crosshole georadar waveform inversion scheme in the presence of dispersion

Waveform tomographic imaging of crosshole georadar data is a powerful method to investigate the shallow subsurface because of its ability to provide images of pertinent petrophysical parameters with extremely high spatial resolution. All current crosshole georadar waveform inversion strategies are based on the assumption of frequency-independent electromagnetic constitutive parameters. However, in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behavior. In this paper, we evaluate synthetically the reconstruction limits of a recently published crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. Our results indicate that, when combined with a source wavelet estimation procedure that provides a means of partially accounting for the frequency-dependent effects through an "effective" wavelet, the inversion algorithm performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Western Préalpes Médianes Romandes: Timing and structure. A review.

Between the original position and their present day location as klippen, the Prealpes Medianes underwent a complex history of paleotectonics and alpine tectonics. Due to the opening of the Piemont ocean the Brianconnais sedimentation realm of the Prealpes Medianes evolved as a rim basin of the northern passive margin during Jurassic to Eocene times. Different paleotectonic features (normal faults, synsedimentary growth structures, inversion structures) developed and were active above a basal detachment in evaporitic layers. The tectonic movements were a consequence of thermal events in the crust. Isolated from the Iberic continent at the end of the Late Cretaceous, the Brianconnais exotic terrain was incorporated into the accretionary prism of the closing Piemont ocean and the incipient alpine orogeny during the Lutetian-Bartonian. The Prealpes Medianes were detached from their homeland during the Bartonian-Priabonian and were transported onto the foreland. The tectonic style is one of a thin-skinned foreland fold and thrust belt. Fault associated fold development above a main decollement, together with internal deformation, represent the Prealpes Medianes main structural features. The very low-grade metamorphic conditions have their origin in the heat flux induced by tectonic burial by overriding nappes in the accretionary prism. After having been transported on top of the developing Helvetic nappes the Prealpes were emplaced in their present day position in front of the Alpine mountain belt during Oligocene times. Post-emplacement and out of sequence thrusting, possibly younger than Oligocene, is observed and can be related to thrusting in the sedimentary substratum and the basement.

Inversion prepared coronary MR angiography: direct visualization of coronary blood flow.

PURPOSE: Visualization of coronary blood flow by means of a slice-selective inversion pre-pulse in concert with bright-blood coronary MRA. MATERIALS AND METHODS: Coronary magnetic resonance angiography (MRA) of the right coronary artery (RCA) was performed in eight healthy adult subjects on a 1.5 Tesla MR system (Gyroscan ACS-NT, Philips Medical Systems, Best, NL) using a free-breathing navigator-gated and cardiac-triggered 3D steady-state free-precession (SSFP) sequence with radial k-space sampling. Imaging was performed with and without a slice-selective inversion pre-pulse, which was positioned along the main axis of the coronary artery but perpendicular to the imaging volume. Objective image quality parameters such as SNR, CNR, maximal visible vessel length, and vessel border definition were analyzed. RESULTS: In contrast to conventional bright-blood 3D coronary MRA, the selective inversion pre-pulse provided a direct measure of coronary blood flow. In addition, CNR between the RCA and right ventricular blood pool was increased and the vessels had a tendency towards better delineation. Blood SNR and CNR between right coronary blood and epicardial fat were comparable in both sequences. CONCLUSION: The combination of a free-breathing navigator-gated and cardiac-triggered 3D SSFP sequence with a slice-selective inversion pre-pulse allows for direct and directional visualization of coronary blood flow with the additional benefit of improved contrast between coronary and right ventricular blood pool.

La nappe du Lebendun entre Alte Kaserne et le Val Cairasca (massif du Simplon): nouvelles observations et interprétations

The lithostratigraphic description of the covers of three Lower Penninic nappes (Monte Leone, Lebendun and Antigorio) allows the comparison of their sedimentary content and their thickness. It has been established that the Lebendun nappe is formed by an ante-Triassic paragneissic core (Valgrande gneiss), and a Mesozoic sedimentary cover in reversed position. The cover series shows a continuous detritic sedimentation, off which the material comes from a continental erosion related to the early Lias rifting phase of the Alpine Tethys. The erosion has reached the basement, resedimented as pebbles and sandstones. This can be observed in both Lebendun and Antigorio covers. The definition of a unit named <<serie intermediaire>> between the Lebendun and the Antigorio covers has important palinspastic implications for both nappes. The unit is composed of a banded marble, a garnet bearing gneiss and a calcschist with great blocks. The comparison between the thickness of Antigorio and Lebendun covers suggests a shoulder position for Antigorio. and a proximal rift basin position tor Lebendun. The general thickness decrease of the series towards the SW points to a NE origin for the Lebendun clastics, taking into account the increase of tectonic deformation in the region trending from east to west. The detritic sedimentation ends with the basin drowning during the Malm, represented by a pure marble sealing the erosive disconformity of the Antigorio cover, and the clastic deposits of Lebendun. Three hypotheses are proposed for the calcschists age and attribution of the <<serie intermediaire>>: A: they belong entirely or partially to the Lebendun cover and correspond to a conglomeratic deposit of Cretaceous-Tertiary Niesen flysch type, of proximal facies. The tectonic limit could be situated in the middle of the calcschists at the level of the huge blocks encountered. B: they belong to Antigorio and correspond to an upper Lias-Dogger synrift deposit, then the marble is liassic. C: they belong to Antigorio and have been deposited following the Lebendun basin inversion (Cretaceous-Tertiary). that generates Tertiary wildflysch deposits, coming from the South for the ultrahelvetic and from the North for the Niesen.

Free-breathing renal magnetic resonance angiography with steady-state free-precession and slab-selective spin inversion combined with radial k-space sampling and water-selective excitation.

The impact of radial k-space sampling and water-selective excitation on a novel navigator-gated cardiac-triggered slab-selective inversion prepared 3D steady-state free-precession (SSFP) renal MR angiography (MRA) sequence was investigated. Renal MRA was performed on a 1.5-T MR system using three inversion prepared SSFP approaches: Cartesian (TR/TE: 5.7/2.8 ms, FA: 85 degrees), radial (TR/TE: 5.5/2.7 ms, FA: 85 degrees) SSFP, and radial SSFP combined with water-selective excitation (TR/TE: 9.9/4.9 ms, FA: 85 degrees). Radial data acquisition lead to significantly reduced motion artifacts (P < 0.05). SNR and CNR were best using Cartesian SSFP (P < 0.05). Vessel sharpness and vessel length were comparable in all sequences. The addition of a water-selective excitation could not improve image quality. In conclusion, radial k-space sampling reduces motion artifacts significantly in slab-selective inversion prepared renal MRA, while SNR and CNR are decreased. The addition of water-selective excitation could not improve the lower CNR in radial scanning.

Régime métropolitain, migrations résidentielles et recompositions urbaines : l'hypothèse de l'inversion du modèle de mobilité

Phénomène resté assez marginal jusqu'aux années quatre-vingt, la gentrification est devenue, dans certaines villes européennes, un principe de production de l'espace urbain. La ville que désertaient les classes moyennes parce qu'elle était associée à une densité excessive et à des nuisances de toutes sortes redevient ainsi le lieu d'une valorisation immobilière renouvelée, mais aussi d'un entre-soi sélectif, du spectacle et du plaisir répondant aux nouveaux styles d'habiter des nouveaux cadres de l'urbain globalisé. Qu'en est-il en Suisse ? L'article présente quelques résultats issus de deux recherches menées à l'échelle du système urbain et de l'agglomération lausannoise. L'observation des dynamiques à l'oeuvre ne permet pas de corroborer la thèse du « retour en ville ».

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.