A sampling method for detecting buried objects using electromagnetic scattering

We consider a simple (but fully three-dimensional) mathematical model for the electromagnetic exploration of buried, perfect electrically conducting objects within the soil underground. Moving an electric device parallel to the ground at constant height in order to generate a magnetic field, we measure the induced magnetic field within the device, and factor the underlying mathematics into a product of three operations which correspond to the primary excitation, some kind of reflection on the surface of the buried object(s) and the corresponding secondary excitation, respectively. Using this factorization we are able to give a justification of the so-called sampling method from inverse scattering theory for this particular set-up.

Earthquake forecasting and seismic hazard analysis: some insights on the testing phase and the modeling

This thesis is divided in three chapters. In the first chapter we analyse the results of the world forecasting experiment run by the Collaboratory for the Study of Earthquake Predictability (CSEP). We take the opportunity of this experiment to contribute to the definition of a more robust and reliable statistical procedure to evaluate earthquake forecasting models. We first present the models and the target earthquakes to be forecast. Then we explain the consistency and comparison tests that are used in CSEP experiments to evaluate the performance of the models. Introducing a methodology to create ensemble forecasting models, we show that models, when properly combined, are almost always better performing that any single model. In the second chapter we discuss in depth one of the basic features of PSHA: the declustering of the seismicity rates. We first introduce the Cornell-McGuire method for PSHA and we present the different motivations that stand behind the need of declustering seismic catalogs. Using a theorem of the modern probability (Le Cam's theorem) we show that the declustering is not necessary to obtain a Poissonian behaviour of the exceedances that is usually considered fundamental to transform exceedance rates in exceedance probabilities in the PSHA framework. We present a method to correct PSHA for declustering, building a more realistic PSHA. In the last chapter we explore the methods that are commonly used to take into account the epistemic uncertainty in PSHA. The most widely used method is the logic tree that stands at the basis of the most advanced seismic hazard maps. We illustrate the probabilistic structure of the logic tree, and then we show that this structure is not adequate to describe the epistemic uncertainty. We then propose a new probabilistic framework based on the ensemble modelling that properly accounts for epistemic uncertainties in PSHA.

Total internal reflection fluorescence cross-correlation spectroscopy: theory and application for studying boundary slip phenomenon

I present a new experimental method called Total Internal Reflection Fluorescence Cross-Correlation Spectroscopy (TIR-FCCS). It is a method that can probe hydrodynamic flows near solid surfaces, on length scales of tens of nanometres. Fluorescent tracers flowing with the liquid are excited by evanescent light, produced by epi-illumination through the periphery of a high NA oil-immersion objective. Due to the fast decay of the evanescent wave, fluorescence only occurs for tracers in the ~100 nm proximity of the surface, thus resulting in very high normal resolution. The time-resolved fluorescence intensity signals from two laterally shifted (in flow direction) observation volumes, created by two confocal pinholes are independently measured and recorded. The cross-correlation of these signals provides important information for the tracers’ motion and thus their flow velocity. Due to the high sensitivity of the method, fluorescent species with different size, down to single dye molecules can be used as tracers. The aim of my work was to build an experimental setup for TIR-FCCS and use it to experimentally measure the shear rate and slip length of water flowing on hydrophilic and hydrophobic surfaces. However, in order to extract these parameters from the measured correlation curves a quantitative data analysis is needed. This is not straightforward task due to the complexity of the problem, which makes the derivation of analytical expressions for the correlation functions needed to fit the experimental data, impossible. Therefore in order to process and interpret the experimental results I also describe a new numerical method of data analysis of the acquired auto- and cross-correlation curves – Brownian Dynamics techniques are used to produce simulated auto- and cross-correlation functions and to fit the corresponding experimental data. I show how to combine detailed and fairly realistic theoretical modelling of the phenomena with accurate measurements of the correlation functions, in order to establish a fully quantitative method to retrieve the flow properties from the experiments. An importance-sampling Monte Carlo procedure is employed in order to fit the experiments. This provides the optimum parameter values together with their statistical error bars. The approach is well suited for both modern desktop PC machines and massively parallel computers. The latter allows making the data analysis within short computing times. I applied this method to study flow of aqueous electrolyte solution near smooth hydrophilic and hydrophobic surfaces. Generally on hydrophilic surface slip is not expected, while on hydrophobic surface some slippage may exists. Our results show that on both hydrophilic and moderately hydrophobic (contact angle ~85°) surfaces the slip length is ~10-15nm or lower, and within the limitations of the experiments and the model, indistinguishable from zero.

Application of the specular and diffuse reflection analysis for in vitro diagnostics of dental erosion: correlation with enamel softening, roughness, and calcium release

We present assembly and application of an optical reflectometer for the analysis of dental erosion. The erosive procedure involved acid-induced softening and initial substance loss phases, which are considered to be difficult for visual diagnosis in a clinic. Change of the specular reflection signal showed the highest sensitivity for the detection of the early softening phase of erosion among tested methods. The exponential decrease of the specular reflection intensity with erosive duration was compared to the increase of enamel roughness. Surface roughness was measured by optical analysis, and the observed tendency was correlated with scanning electron microscopy images of eroded enamel. A high correlation between specular reflection intensity and measurement of enamel softening (r(2) ? -0.86) as well as calcium release (r(2) ? -0.86) was found during erosion progression. Measurement of diffuse reflection revealed higher tooth-to-tooth deviation in contrast to the analysis of specular reflection intensity and lower correlation with other applied methods (r(2) = 0.42-0.48). The proposed optical method allows simple and fast surface analysis and could be used for further optimization and construction of the first noncontact and cost-effective diagnostic tool for early erosion assessment in vivo.

Learning My Lessons: Research, Rehearsal, and Reflection on the Greats

A successful actor often requires a specific acting method or style to enhance their performance. Through theatrical research, rehearsal and performance, an actor can narrow down their seemingly endless search for the most productive methodology. By researching, studying, and applying the methods of Constantin Stanislavski, Stella Adler,and Tadashi Suzuki to my rehearsal process, I have found my most effective acting style: Stella Adler?s method. I utilize this acting method during the performance period of my early professional acting career. Experimental research for this thesis was completed inthe studio. I applied each of the three aforementioned methods to a dramatic/classical monologue. The results I gathered helped me to decide upon Adler?s methodology to carry with me through my upcoming professional auditions and career. From casting resulting from the auditions, I will employ the methodology to my professional work asan actress. Each acting teacher has provided the performance world with a new way to experience their stage time. The methods are unique and enable the actor to find the most dynamic performance through engaging technical skill.

TUNING, AVO, AND FLAT-SPOT EFFECTS IN A SEISMIC ANALYSIS OF NORTH SEA BLOCK F3

Reflection seismic data from the F3 block in the Dutch North Sea exhibits many large-amplitude reflections at shallow horizons, typically categorized as “brightspots ” (Schroot and Schuttenhelm, 2003), mainly because of their bright appearance. In most cases, these bright reflections show a significant “flatness” contrasting with local structural trends. While flatspots are often easily identified in thick reservoirs, we have often occasionally observed apparent flatspot tuning effects at fluid contacts near reservoir edges and in thin reservoir beds, while only poorly understanding them. We conclude that many of the shallow large-amplitude reflections in block F3 are dominated by flatspots, and we investigate the thin-bed tuning effects that such flatspots cause as they interact with the reflection from the reservoir’s upper boundary. There are two possible effects to be considered: (1) the “wedge-model” tuning effects of the flatspot and overlying brightspots, dimspots, or polarity-reversals; and (2) the stacking effects that result from possible inclusion of post-critical flatspot reflections in these shallow sands. We modeled the effects of these two phenomena for the particular stratigraphic sequence in block F3. Our results suggest that stacking of post-critical flatspot reflections can cause similar large-amplitude but flat reflections, in some cases even causing an interface expected to produce a ‘dimspot’ to appear as a ‘brightspot’. Analysis of NMO stretch and muting shows the likely exclusion of critical offset data in stacked output. If post-critical reflections are included in stacking, unusual results will be observed. In the North Sea case, we conclude the tuning effect was the primary reason causing for the brightness and flatness of these reflections. However, it is still important to note that care should be taken while applying muting on reflections with wide range of incidence angles and the inclusion of critical offset data may cause some spurious features in the stacked section.

USING AUTO- AND CROSS-CORRELATIONS FROM SEISMIC NOISE TO MONITOR VELOCITY CHANGES AT VILLARRICA VOLCANO, CHILE

We used the Green's functions from auto-correlations and cross-correlations of seismic ambient noise to monitor temporal velocity changes in the subsurface at Villarrica volcano in the Southern Andes of Chile. Campaigns were conducted from March to October 2010 and February to April 2011 with 8 broadband and 6 short-period stations, respectively. We prepared the data by removing the instrument response, normalizing with a root-mean-square method, whitening the spectra, and filtering from 1 to 10 Hz. This frequency band was chosen based on the relatively high background noise level in that range. Hour-long auto- and cross-correlations were computed and the Green's functions stacked by day and total time. To track the temporal velocity changes we stretched a 24 hour moving window of correlation functions from 90% to 110% of the original and cross correlated them with the total stack. All of the stations' auto-correlations detected what is interpreted as an increase in velocity in 2010, with an average increase of 0.13%. Cross-correlations from station V01, near the summit, to the other stations show comparable changes that are also interpreted as increases in velocity. We attribute this change to the closing of cracks in the subsurface due either to seasonal snow loading or regional tectonics. In addition to the common increase in velocity across the stations, there are excursions in velocity on the same order lasting several days. Amplitude decreases as the station's distance from the vent increases suggesting these excursions may be attributed to changes within the volcanic edifice. In at least two occurrences the amplitudes at stations V06 and V07, the stations farthest from the vent, are smaller. Similar short temporal excursions were seen in the auto-correlations from 2011, however, there was little to no increase in the overall velocity.

A new optical detection method to assess the erosion inhibition by in vitro salivary pellicle layer

OBJECTIVES Application of the recently developed optical method based on the monitoring of the specular reflection intensity to study the protective potential of the salivary pellicle layer against early enamel erosion. METHODS The erosion progression was compared between two treatment groups: enamel samples coated by the 15 h-in vitro-formed salivary pellicle layer (group P, n=90) and the non-coated enamel surfaces (control group C, n=90). Different severity of the erosive impact was modelled by the enamel incubation in 1% citric acid (pH=3.6) for 2, 4, 8, 10 or 15 min. Erosion quantification was performed by the optical method as well as by the microhardness and calcium release analyses. RESULTS Optical assessment of the erosion progression showed erosion inhibition by the in vitro salivary pellicle in short term acidic treatments (≤ 4 min) which was also confirmed by microhardness measurements proving significantly less (p<0.05) enamel softening in the group P at 2 and 4 min of erosion compared to the group C. SEM images demonstrated less etched enamel interfaces in the group P at short erosion durations as well. CONCLUSIONS Monitoring of the specular reflection intensity can be successfully applied to quantify early erosion progression in comparative studies. In vitro salivary pellicle (2h) provides erosion inhibition but only in short term acidic exposures. CLINICAL SIGNIFICANCE The proposed optical technique is a promising tool for the fast and non-invasive erosion quantification in clinical studies.

Multibeam and seismic survey, physical oceanography, and sea-bottom videos from the Bay of Biscay

We report the northernmost and deepest known occurrence of deep-water pycnodontine oysters, based on two surveys along the French Atlantic continental margin to the La Chapelle continental slope (2006) and the Guilvinec Canyon (2008). The combined use of multibeam bathymetry, seismic profiling, CTD casts and a remotely operated vehicle (ROV) made it possible to describe the physical habitat and to assess the oceanographic control for the recently described species Neopycnodonte zibrowii. These oysters have been observed in vivo in depths from 540 to 846 m, colonizing overhanging banks or escarpments protruding from steep canyon flanks. Especially in the Bay of Biscay, such physical habitats may only be observed within canyons, where they are created by both long-term turbiditic and contouritic processes. Frequent observations of sand ripples on the seabed indicate the presence of a steady, but enhanced bottom current of about 40 cm/s. The occurrence of oysters also coincides with the interface between the Eastern North Atlantic Water and the Mediterranean Outflow Water. A combination of this water mass mixing, internal tide generation and a strong primary surface productivity may generate an enhanced nutrient flux, which is funnelled through the canyon. When the ideal environmental conditions are met, up to 100 individuals per m² may be observed. These deep-water oysters require a vertical habitat, which is often incompatible with the requirements of other sessile organisms, and are only sparsely distributed along the continental margins. The discovery of these giant oyster banks illustrates the rich biodiversity of deep-sea canyons and their underestimation as true ecosystem hotspots.

Seismic velocities and elastic properties of gabbroic rocks from ODP Hole 118-735B

The nearly continuous recovery of 0.5 km of generally fresh, layer 3 gabbroic rocks at Hole 735B, especially near the bottom of the section, presents scientists an unusual opportunity to study the detailed elastic properties of the lower oceanic crust. Extending compressional-wave and density shipboard measurements at room pressure, Vp and Vs were measured at pressures from 20 to 200 MPa using the pulse transmission method. All of the rocks exhibit significant increases in velocity with increasing pressure up to about 150 MPa, a feature attributed to the closing of microcrack porosity. Measured velocities reflect the mineralogical makeup and microstructures acquired during the tectonic history of Hole 735B. Most of the undeformed and unaltered gabbros are approximately 65:35 plagioclase/clinopyroxene rocks plus olivine or oxide minerals, and the observed densities and velocities are fully consistent with the Voigt-Reuss-Hill (VRH) averages of the component minerals and their proportions. Depending on their olivine content, the predominant olivine gabbros at 200 MPa have average Vp = 7.1 ± 0.2 km/s, Vs = 3.9 ± 0.1 km/s, and grain densities of 2.95 ± 0.5 g/cm3. The less abundant iron-titanium (Fe-Ti) oxide gabbros average Vp = 6.75 ± 0.15 km/s, Vs = 3.70 ± 0.1 km/s, and grain densities of 3.22 ± 0.05 g/cm3, reflecting the higher densities and lower velocities of oxide minerals compared to olivine. About 30% of the core is plastically deformed, and the densities and directionally averaged velocities of these shear-zone tectonites are generally consistent with those of the gabbros, their protoliths. Three sets of observations indicate that the shear-zone metagabbros are elastically anisotropic: (1) directional variations in Vp, both vertical and horizontal and with respect to foliation and lineation; (2) discrepancies among Vp values for the horizontal cores and the VRH averages of the component minerals and their mineral proportions, suggesting preferred crystallographic orientations of anisotropic minerals; and (3) variations of Vs of up to 7%, with polarization directions parallel and perpendicular to foliation. Optical inspection of thin sections of the same samples indicates that plagioclase feldspar, clinopyroxene, and amphibole typically display crystallographic-preferred orientations, and this, plus the elastic anisotropy of these minerals, suggests that preferred orientations are responsible for much of the observed anisotropy, particularly at high pressure. Alteration tends to be localized to brittle faults and brecciated zones, and typical alteration minerals are amphibole and secondary plagioclase, which do not significantly change the velocity-density relationships.

Seismic data from Colle Gnifetti glacie saddle around borehole KCI, Swiss Alps

Two seismic surveys were carried out on the high-altitude glacier saddle, Colle Gnifetti, Monte Rosa, Italy/Switzerland. Explosive and vibroseismic sources were tested to explore the best way to generate seismic waves to deduce shallow and intermediate properties (<100 m) of firn and ice. The explosive source (SISSY) excites strong surface and diving waves, degrading data quality for processing; no englacial reflections besides the noisy bed reflector are visible. However, the strong diving waves are analyzed to derive the density distribution of the firn pack, yielding results similar to a nearby ice core. The vibrator source (ElViS), used in both P- and SH-wave modes, produces detectable laterally coherent reflections within the firn and ice column. We compare these with ice-core and radar data. The SH-wave data are particularly useful in providing detailed, high-resolution information on firn and ice stratigraphy. Our analyses demonstrate the potential of seismic methods to determine physical properties of firn and ice, particularly density and potentially also crystal-orientation fabric.

Physical properties, velocity and attenuation characteristics, and mineralogy of gabbros from ODP Hole 176-735B

Laboratory compressional wave (Vp) and shear wave (Vs) velocities were measured as a function of confining pressure for the gabbros from Hole 735B and compared to results from Leg 118. The upper 500 m of the hole has a Vp mean value of 6895 m/s measured at 200 MPa, and at 500 meters below seafloor (mbsf), Vp measurements show a mean value of 7036 m/s. Vs mean values in the same intervals are 3840 m/s and 3857 m/s, respectively. The mean Vp and Vs values obtained from log data in the upper 600 m are 6520 and 3518 m/s, respectively. These results show a general increase in velocity with depth and the velocity gradients estimate an upper mantle depth of 3.32 km. This value agrees with previous work based on dredged samples and inversion of rare element concentrations in basalts dredged from the conjugate site to the north of the Atlantis Bank. Laboratory measurements show Vp anisotropy ranging between 0.4% and 8.8%, with the majority of the samples having values less than 3.8%. Measurements of velocity anisotropy seem to be associated with zones of high crystal-plastic deformation with predominant preferred mineral orientations of plagioclase, amphiboles, and pyroxenes. These findings are consistent with results on gabbros from the Hess Deep area and suggest that plastic deformation may play an important role in the seismic properties of the lower oceanic crust. In contrast to ophiolite studies, many of the olivine gabbros show a small degree of anisotropy. Log derived Vs anisotropy shows an average of 5.8% for the upper 600 m of Hole 735B and tends to decrease with depth where the overburden pressure and the age of the crustal section suggests closure of cracks and infilling of fractures by alteration minerals. Overall the results indicate that the average shear wave splitting in Hole 735B might be influenced by preferred structural orientations and the average value of shear wave splitting may not be a maximum because structural dips are <90°. The maximum fast-wave orientation values could be influenced by structural features striking slightly oblique to this orientation or by near-field stress concentrations. However, flexural wave dispersion analyses have not been performed to confirm this hypothesis or to indicate to what extent the near-field stresses may be influencing shear wave propagation. Acoustic impedance contrasts calculated from laboratory and logging data were used to generate synthetic seismograms that aid in the interpretation of reflection profiles. Several prominent reflections produced by these calculations suggest that Fe-Ti oxides and shear zones may contribute to the reflective nature of the lower oceanic crust. Laboratory velocity attenuation (Q) measurements from below 500 m have a mean value of 35.1, which is consistent with previous vertical seismic profile (VSP) and laboratory measurements on the upper 500 m.

Ocean bottom seismometer sgy-files of refraction seismic profiles from Professor Logachev cruise TTR16/3

Over the last decade pockmarks have proven to be important seabed features that provide information about fluid flow on continental margins. Their formation and dynamics are still poorly constrained due to the lack of proper three dimensional imaging of their internal structure. Numerous fluid escape features provide evidence for an active fluid-flow system on the Norwegian margin, specifically in the Nyegga region. In June-July 2006 a high-resolution seismic experiment using Ocean Bottom Seismometers (OBS) was carried out to investigate the detailed 3D structure of a pockmark named G11 in the region. An array of 14 OBS was deployed across the pockmark with 1 m location accuracy. Shots fired from surface towed mini GI guns were also recorded on a near surface hydrophone streamer. Several reflectors of high amplitude and reverse polarity are observed on the profiles indicating the presence of gas. Gas hydrates were recovered with gravity cores from less than a meter below the seafloor during the cruise. Indications of gas at shallow depths in the hydrate stability field show that methane is able to escape through the water-saturated sediments in the chimney without being entirely converted into gas hydrate. An initial 2D raytraced forward model of some of the P wave data along a line running NE-SW across the G11 pockmark shows, a gradual increase in velocity between the seafloor and a gas charged zone lying at ~300 m depth below the seabed. The traveltime fit is improved if the pockmark is underlain by velocities higher than in the surrounding layer corresponding to a pipe which ascends from the gas zone, to where it terminates in the pockmark as seen in the reflection profiles. This could be due to the presence of hydrates or carbonates within the sediments.

Dynamic Response of Wind Turbines to Theoretical 3D Seismic Motions Taking into Account the Rotational Component

We study the dynamic response of a wind turbine structure subjected to theoretical seismic motions, taking into account the rotational component of ground shaking. Models are generated for a shallow moderate crustal earthquake in the Madrid Region (Spain). Synthetic translational and rotational time histories are computed using the Discrete Wavenumber Method, assuming a point source and a horizontal layered earth structure. These are used to analyze the dynamic response of a wind turbine, represented by a simple finite element model. Von Mises stress values at different heights of the tower are used to study the dynamical structural response to a set of synthetic ground motion time histories

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography. The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.