Reproducibility of straylight measurement by C-Quant for assessment of retinal straylight using the compensation comparison method.

BACKGROUND: Straylight gives the appearance of a veil of light thrown over a person's retinal image when there is a strong light source present. We examined the reproducibility of the measurements by C-Quant, and assessed its correlation to characteristics of the eye and subjects' age. PARTICIPANTS AND METHODS: Five repeated straylight measurements were taken using the dominant eye of 45 healthy subjects (age 21-59) with a BCVA of 20/20: 14 emmetropic, 16 myopic, eight hyperopic and seven with astigmatism. We assessed the extent of reproducibility of straylight measures using the intraclass correlation coefficient. RESULTS: The mean straylight value of all measurements was 1.01 (SD 0.23, median 0.97, interquartile range 0.85-1.1). Per 10 years of age, straylight increased in average by 0.10 (95%CI 0.04 to 0.16, p < 0.01]. We found no independent association of refraction (range -5.25 dpt to +2 dpt) on straylight values (0.001; 95%CI -0.022 to 0.024, p = 0.92). Compared to emmetropic subjects, myopia reduced straylight (-.011; -0.024 to 0.02, p = 0.11), whereas higher straylight values (0.09; -0.01 to 0.20, p = 0.09) were observed in subjects with blue irises as compared to dark-colored irises when correcting for age. The intraclass correlation coefficient (ICC) of repeated measurements was 0.83 (95%CI 0.76 to 0.90). CONCLUSIONS: Our study showed that straylight measurements with the C-Quant had a high reproducibility, i.e. a lack of large intra-observer variability, making it appropriate to be applied in long-term follow-up studies assessing the long-term effect of surgical procedures on the quality of vision.

Jurassic sedimentary record and tectonic evolution of the northwestern corner of Africa

The history of the opening seaway from the westernmost Tethys to the Central Atlantic is traced by the analysis of the sedimentary facies development in the external Rif basin of Northern Morocco and the geological and seismic data from the Moroccan Atlantic continental margin. In the Rif basin, after the early Sinemurian, sedimentary facies dated by ammonites, foraminifers and brachiopods, indicate a progression of rapid subsidence resulting from extensional tectonic (tilted blocks, escarpment fault breccias, neptunian dykes etc.) from the N and NE to the S and SW. From the Toarcian to the Bajocian, deltas progress from the W and SW into the `'Rides sud-rifaines'' realm. From the late Bathonian to the Oxfordian, deep-sea fans develop in the external Rif. During the same period, deltaic sediments fill in the Middle Atlas basin of Eastern Morocco and progress into the external Rif. The top of the Jurassic is characterised by carbonate deposits. At the northwestern corner of Africa, the subsidence of the sedimentary basins by rifting is initiated in the late Triassic; however, at the Mazagan transect of the Atlantic continental margin, the tectonic pattern characteristic of a passive continental margin appears clearly only in the early Jurassic. At the foot of the Mazagan escarpment, the sedimentary record shows a foundering of the first bloc during early to middle Lias. A thermal uplift phase is indicated by emersion of the African margin shoulder in late Liassic, and thermal relaxation starts in the middle Jurassic. The morphology of this transect, compared with the conjugate side of the American continent is most easily explained by the uniform sense simple shear model.

Three-dimensional sequence stratigraphy and subtle stratigraphic traps associated with systems tracts: West Cameron region, offshore Louisiana, Gulf of Mexico

Three-dimensional sequence stratigraphy is a potent exploration and development tool for the discovery of subtle stratigraphic traps. Reservoir morphology, heterogeneity and subtle stratigraphic trapping mechanisms can be better understood through systematic horizontal identification of sedimentary facies of systems tracts provided by three-dimensional attribute maps used as an important complement to the sequential analysis on the two-dimensional seismic lines and the well log data. On new prospects as well as on already-producing fields, the additional input of sequential analysis on three-dimensional data enables the identification, location and precise delimitation of new potentially productive zones. The first part of this paper presents four typical horizontal seismic facies assigned to the successive systems tracts of a third- or fourth-order sequence deposited in inner to outer neritic conditions on a elastic shelf. The construction of this synthetic representative sequence is based on the observed reproducibility of the horizontal seismic facies response to cyclic eustatic events on more than 35 sequences registered in the Gulf coast Plio-Pleistocene and Late Miocene, offshore Louisiana in the West Cameron region of the Gulf of Mexico. The second part shows how three-dimensional sequence stratigraphy can contribute in localizing and understanding sedimentary facies associated with productive zones. A case study in the early Middle Miocene Cibicides opima sands shows multiple stacked gas accumulations in the top slope fan, prograding wedge and basal transgressive systems tract of the third-order sequence between SB15.5 and SB 13.8 Ma.

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.

Global mass wasting during the Middle Ordovician: Meteoritic trigger or plate-tectonic environment ?

Mass wasting at continental margins on a global scale during the Middle Ordovician has recently been related to high meteorite influx. Although a high meteorite influx during the Ordovician should not be neglected, we challenge the idea that mass wasting was mainly produced by meteorite impacts over a period of almost 10 Ma. Having strong arguments against the impact-related hypothesis, we propose an alternative explanation, which is based on a re-evaluation of the mass wasting sites, considering their plate-tectonic distribution and the global sea level curve. A striking and important feature is the distribution of most of the mass wasting sites along continental margins characterised by periods of magmatism, terrane accretion and continental or back-arc rifting, respectively, related to subduction of oceanic lithosphere. Such processes are commonly connected with seismic activity causing earthquakes, which can cause downslope movement of sediment and rock. Considering all that, it seems more likely that most of this mass wasting was triggered by earthquakes related to plate-tectonic processes, which caused destabilisation of continental margins resulting in megabreccias and debris flows. Moreover, the period of mass wasting coincides with sea level drops during global sea level lowstand. In some cases, sea level drops can release pore-water overpressure reducing sediment strength and hence promoting instability of sediment at continental margins. Reduced pore-water overpressure can also destabilise gas hydrate-bearing sediment, causing slope failure, and thus resulting in submarine mass wasting. Overall, the global mass wasting during the Middle Ordovician does not need meteoritic trigger. (C) 2010 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

Geochemistry and 40Ar/39Ar geochronology of pseudotachylyte associated with UHP whiteschists from the Dora Maira massif, Italy

In situ UV-Iaser ablation Ar-40/(39) Ar geochronological and geochemical data, together with rock and mineral compositional data, have been determined from pseudotachylyte and surrounding mylonitic gneiss associated with the UHP whiteschists of the Dora Maira Massif, Italy. Several generations of fresh pseudotachylyte occur as irregular veins up to a few cur thick both parallel and at high angles to the foliation. Whole rock XRF data collected from representative lithologies of mylonitic gneiss are uniformly consistent with a mildly alkalic granitic protolith. Minimal compositional variation is observed between the pseudotachylyte and its surrounding mylonitic gneiss. The pseudotachylyte contains newly crystallized grains of biotite and K-feldspar in a matrix of glass with partially fused grains of quartz, zircon, apatite, and titanite. Electron microprobe analyses of the glass show significant compositional variation that is probably strongly influenced by micrometer-scale changes in mineralogy. UV-Iaser ablation ICP-MS traverses across the mylonitic gneiss-pseudotachylyte contact are consistent with cataclastic communition of REE carriers such as epidote, monazite, allanite, zircon, and apatite before melting as an efficient mechanism of REE homogenization in the pseudotachylyte. The 40Ar/39Ar data from one band of pseudotachylyte indicate formation at 20.1 +/- 0.5 Ma, when the mylonitic gneisses were already in a near surface position. The variable effects of top-to-the-west shear deformation within outcrops of the coesite-bearing unit are reflected in localized zones of protomylonite, cataclasite, ultracataclasite, and pseudotachylyte. Preservation of several generations of pseudotachylyte suggests that seismic events may have played a significant role in triggering late unroofing of the UHP rocks. It is speculated that deeper crustal seismic events potentially played a role in the unroofing of the UHP rocks at earlier stages in their exhumation history. (c) 2005 Elsevier B.V. All rights reserved.

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantan that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a similar to 50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a similar to 20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantan and illustrate the potential of ERT surveying for paleoseismic studies.

Frequency-dependent attenuation as a potential indicator of oil saturation

At seismic frequencies, wave-induced fluid flow is a major cause of P-wave attenuation in partially saturated porous rocks. Attenuation is of great importance for the oil industry in the interpretation of seismic field data. Here, the effects on P-wave attenuation resulting from changes in oil saturation are studied for media with coexisting water, oil, and gas. For that, creep experiments are numerically simulated by solving Biot's equations for consolidation of poroelastic media with the finite-element method. The experiments yield time-dependent stress?strain relations that are used to calculate the complex P-wave modulus from which frequency-dependent P-wave attenuation is determined. The models are layered media with periodically alternating triplets of layers. Models consisting of triplets of layers having randomly varying layer thicknesses are also considered. The layers in each triplet are fully saturated with water, oil, and gas. The layer saturated with water has lower porosity and permeability than the layers saturated with oil and gas. These models represent hydrocarbon reservoirs in which water is the wetting fluid preferentially saturating regions of lower porosity. The results from the numerical experiments showed that increasing oil saturation, connected to a decrease in gas saturation, resulted in a significant increase of attenuation at low frequencies (lower than 2 Hz). Furthermore, replacing the oil with water resulted in a distinguishable behavior of the frequency-dependent attenuation. These results imply that, according to the physical mechanism of wave-induced fluid flow, frequency-dependent attenuation in media saturated with water, oil, and gas is a potential indicator of oil saturation.

Fluid evolution at the Variscan front in the vicinity of the Aachen thrust

Quartz-carbonate-chlorite veins were studied in borehole samples of the RWTH-1 well in Aachen. Veins formed in Devonian rocks in the footwall of the Aachen thrust during Variscan deformation and associated fluid flow. Primary fluid inclusions indicate subsolvus unmixing of a homogenous H(2)O-CO(2)-CH(4)-(N(2))-Na-(K)-Cl fluid into a H(2)O-Na-(K)-Cl solution and a vapour-rich CO(2)-(H(2)O, CH(4), N(2)) fluid. The aqueous end-member composition resembles that of metamorphic fluids of the Variscan front zone with salinities ranging from 4 to 7% NaCl equiv. and maximum homogenisation temperatures of close to 400A degrees C. Pressure estimates indicate a burial depth between 4,500 and 8,000 m at geothermal gradients between 50 and 75A degrees C/26 MPa, but pressure decrease to sublithostatic conditions is also indicated, probably as a consequence of fracture opening during episodic seismic activity. A second fluid system, mainly preserved in pseudo-secondary and secondary fluid inclusions, is characterised by fluid temperatures between 200 and 250A degrees C and salinities of < 5% NaCl equiv. Bulk stable isotope analyses of fluids released from vein quartz, calcite, and dolomite by decrepitation yielded delta D(H2O) values from -89 to -113 aEuro degrees, delta(13)C(CH4) from -26.9 to -28.9aEuro degrees (VPDB) and delta(13)C(CO2) from -12.8 to -23.3aEuro degrees (VPDB). The low delta D and delta(13)C range of the fluids is considered to be due to interaction with cracked hydrocarbons. The second fluid influx caused partial isotope exchange and disequilibrium. It is envisaged that an initial short lived flux of hot metamorphic fluids expelled from the epizonal metamorphic domains of the Stavelot-Venn massif. The metamorphic fluid was focused along major thrust faults of the Variscan front zone such as the Aachen thrust. A second fluid influx was introduced from formation waters in the footwall of the Aachen thrust as a consequence of progressive deformation. Mixing of the cooler and lower salinity formation water with the hot metamorphic fluid during episodic fluid trapping resulted in an evolving range of physicochemical fluid inclusion characteristics.

Upper Triassic to Cretaceous radiolaria from Nicaragua and northern Costa Rica - The Mesquito composite oceanic terrane

We propose a new terrane subdivision of Nicaragua and Northern Costa Rica, based on Upper Triassic to Upper Cretaceous radiolarian biochronology of ribbon radiolarites, the newly studied Siuna Serpentinite Mélange, and published 40Ar/39Ar dating and geochemistry of mafic and ultramafic igneous rock units of the area. The new Mesquito Composite Oceanic Terrane (MCOT) comprises the southern half of the Chortis Block, that was assumed to be a continental fragment of N-America. The MCOT is defined by 4 corner localities characterized by ultramafic and mafic oceanic rocks and radiolarites of Late Triassic, Jurassic and Early Cretaceous age: 1. The Siuna Serpentinite Mélange (NE-Nicaragua), 2. The El Castillo Mélange (Nicaragua/Costa Rica border), 3.The Santa Elena Ultramafics (N-Costa Rica) and, 4. DSDP Legs 67/84. 1. The Siuna Serpentinite Mélange contains, high pressure metamorphic mafics and Middle Jurassic (Bajocian-Bathonian) radiolarites in original, sedimentary contact with arc-metandesites. The Siuna Mélange also contains Upper Jurassic black detrital chert formed in a marginal (fore-arc?) basin shortly before subduction. A phengite 40Ar/39Ar -cooling age dates the exhumation of the high pressure rocks as 139 Ma (earliest Cretaceous). 2. The El Castillo Mélange comprises a radiolarite block tectonically embedded in serpentinite that yielded a diverse Rhaetian (latest Triassic) radiolarian assemblage, the oldest fossils recovered so far from S-Central America. 3. The Santa Elena Ultramafics of N-Costa Rica together with the serpentinite outcrops near El Castillo (2) in Southern Nicaragua, are the southernmost outcrops of the MCOT. The Santa Elena Unit (3) itself is still undated, but it is thrust onto the middle Cretaceous Santa Rosa Accretionary Complex (SRAC), that contains Lower to Upper Jurassic, highly deformed radiolarite blocks, probably reworked from the MCOT, which was the upper plate with respect to the SRAC. 4. Serpentinites, metagabbros and basalts have long been known from DSDP Leg 67/84 (3), drilled off Guatemala in the Nicaragua-Guatemala forearc basement. They have been restudied and reveal 40Ar/39Ar dated Upper Triassic to middle Cretaceous enriched Ocean Island Basalts and Jurassic to Lower Cretaceous depleted Island arc rocks of probable Pacific origin. The area between localities 1-4 is largely covered by Tertiary to Recent arcs, but we suspect that its basement is made of oceanic/accreted terranes. Earthquake seismic studies indicate an ill-defined, shallow Moho in this area. The MCOT covers most of Nicaragua and could extend to Guatemala to the W and form the Lower (southern) Nicaragua Rise to the NE. Some basement complexes of Jamaica, Hispaniola and Puerto Rico may also belong to the MCOT. The Nicoya Complex s. str. has been regarded as an example of Caribbean crust and the Caribbean Large Igneous Province (CLIP). However, 40Ar/39Ar - dates on basalts and intrusives indicate ages as old as Early Cretaceous. Highly deformed Jurassic and Lower Cretaceous radiolarites occur as blocks within younger intrusives and basalts. Our interpretation is that radiolarites became first accreted to the MCOT, then became reworked into the Nicoya Plateau in Late Cretaceous times. This implies that the Nicoya Plateau formed along the Pacific edge of the MCOT, independent form the CLIP and most probably unrelated with he Galapagos hotspot. No Jurassic radiolarite, no older sediment age than Coniacian-Santonian, and no older 40Ar/39Ar age than 95 Ma is known from S-Central America between SE of Nicoya and Colombia. For us this area represents the trailing edge of the CLIP s. str.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.

A new method to estimate the infilling of alluvial sediment of glacial valleys using a sloping local base level

A new method is used to estimate the volumes of sediments of glacial valleys. This method is based on the concept of sloping local base level and requires only a digital terrain model and the limits of the alluvial valleys as input data. The bedrock surface of the glacial valley is estimated by a progressive excavation of the digital elevation model (DEM) of the filled valley area. This is performed using an iterative routine that replaces the altitude of a point of the DEM by the mean value of its neighbors minus a fixed value. The result is a curved surface, quadratic in 2D. The bedrock surface of the Rhone Valley in Switzerland was estimated by this method using the free digital terrain model Shuttle Radar Topography Mission (SRTM) (~92 m resolution). The results obtained are in good agreement with the previous estimations based on seismic profiles and gravimetric modeling, with the exceptions of some particular locations. The results from the present method and those from the seismic interpretation are slightly different from the results of the gravimetric data. This discrepancy may result from the presence of large buried landslides in the bottom of the Rhone Valley.