HARPA : a versatile three-dimensional Hamiltonian ray-tracing program for acoustic waves in the atmosphere above irregular terrain /

Mode of access: Internet.

Surface waves in the Beaufort Sea

Thesis (Master's)--University of Washington, 2016-06

Surface quality and surface waves on subwavelength-structured silver films

We analyze the physical-chemical surface properties of single-slit, single-groove subwavelength-structured silver films with high-resolution transmission electron microscopy and calculate exact solutions to Maxwell’s equations corresponding to recent far-field interferometry experiments using these structures. Contrary to a recent suggestion the surface analysis shows that the silver films are free of detectable contaminants. The finite-difference time-domain calculations, in excellent agreement with experiment, show a rapid fringe amplitude decrease in the near zone (slit-groove distance out to 3–4 wavelengths). Extrapolation to slit-groove distances beyond the near zone shows that the surface wave evolves to the expected bound surface plasmon polariton (SPP). Fourier analysis of these results indicates the presence of a distribution of transient, evanescent modes around the SPP that dephase and dissipate as the surface wave evolves from the near to the far zone.

The optical response of nanostructured surfaces and the composite diffracted evanescent wave model

Investigations of the optical response of subwavelength-structure arrays milled into thin metal films have revealed surprising phenomena, including reports of unexpectedly high transmission of light. Many studies have interpreted the optical coupling to the surface in terms of the resonant excitation of surface plasmon polaritons (SPPs), but other approaches involving composite diffraction of surface evanescent waves (CDEW) have also been proposed. Here we present a series of measurements on very simple one-dimensional subwavelength structures to test the key properties of the surface waves, and compare them to the CDEW and SPP models. We find that the optical response of the silver metal surface proceeds in two steps: a diffractive perturbation in the immediate vicinity (2–3 mu m) of the structure, followed by excitation of a persistent surface wave that propagates over tens of micrometres. The measured wavelength and phase of this persistent wave are significantly shifted from those expected for resonance excitation of a conventional SPP on a pure silver surface.

The Hurst exponent as an indicator of the behaviour of a model monopile in an ocean wave testing basin

With the importance of renewable energy well-established worldwide, and targets of such energy quantified in many cases, there exists a considerable interest in the assessment of wind and wave devices. While the individual components of these devices are often relatively well understood and the aspects of energy generation well researched, there seems to be a gap in the understanding of these devices as a whole and especially in the field of their dynamic responses under operational conditions. The mathematical modelling and estimation of their dynamic responses are more evolved but research directed towards testing of these devices still requires significant attention. Model-free indicators of the dynamic responses of these devices are important since it reflects the as-deployed behaviour of the devices when the exposure conditions are scaled reasonably correctly, along with the structural dimensions. This paper demonstrates how the Hurst exponent of the dynamic responses of a monopile exposed to different exposure conditions in an ocean wave basin can be used as a model-free indicator of various responses. The scaled model is exposed to Froude scaled waves and tested under different exposure conditions. The analysis and interpretation is carried out in a model-free and output-only environment, with only some preliminary ideas regarding the input of the system. The analysis indicates how the Hurst exponent can be an interesting descriptor to compare and contrast various scenarios of dynamic response conditions.

Aged embankment imaging and assessment using surface waves

Rapid, non-intrusive surface wave surveys provide depth profiles from which ground models can be generated for use in earthwork condition assessment. Stiffness throughout earthworks controls the behaviour under static and dynamic loads, and characterising heterogeneity is of interest in relation to the stability of engineered backfill and life-cycle deterioration in aged utility and transportation infrastructure. Continuous surface wave methods were used to identify interfaces between fine- and coarse-grained fill in an end-tipped embankment along the Great Central Railway in Nottinghamshire, UK. Multichannel analysis of surface wave (MASW) methods were used to characterise subsurface voiding in a canal embankment along the Knottingley and Goole canal near Eggborough, Yorkshire. MASW methods are currently being used to study extreme weather impacts on the stability of a highplasticity clay embankment along the Gloucestershire–Warwickshire railway near Laverton. Optimal results were obtained using equipment capable of generating and detecting over wide frequency ranges.

Airborne remote sensing of ocean wave directional wavenumber spectra in the marginal ice zone

Interactions between surface waves and sea ice are thought to be an important, but poorly understood, physical process in the atmosphere-ice-ocean system. In this work, airborne scanning lidar was used to observe ocean waves propagating into the marginal ice zone (MIZ). These represent the first direct spatial measurements of the surface wave field in the polar MIZ. Data were compared against two attenuation models, one based on viscous dissipation and one based on scattering. Both models were capable of reproducing the measured wave energy. The observed wavenumber dependence of attenuation was found to be consistent with viscous processes, while the spectral spreading of higher wavenumbers suggested a scattering mechanism. Both models reproduced a change in peak direction due to preferential directional filtering. Floe sizes were recorded using co-located visible imagery, and their distribution was found to be consistent with ice breakup by the wave field.

Nonlocal magnetorotational instability

An analytical theory of the nonlocal magnetorotational instability (MRI) is developed for the simplest astrophysical plasma model. It is assumed that the rotation frequency profile has a steplike character, so that there are two regions in which it has constant different values, separated by a narrow transition layer. The surface wave approach is employed to investigate the MRI in this configuration. It is shown that the main regularities of the nonlocal MRI are similar to those of the local instability and that driving the nonaxisymmetric MRI is less effective than the axisymmetric one, also for the case of the nonlocal instability. The existence of nonlocal instabilities in nonmagnetized plasma is predicted. (c) 2008 American Institute of Physics.

New methods for the assessment of in vitro and in vivo stress cracking in biomedical polyurethanes

This article describes a new test method for the assessment of the severity of environmental stress cracking of biomedical polyurethanes in a manner that minimizes the degree of subjectivity involved. The effect of applied strain and acetone pre-treatment on degradation of Pellethane 2363 80A and Pellethane 2363 55D polyurethanes under in vitro and in vivo conditions is studied. The results are presented using a magnification-weighted image rating system that allows the semi-quantitative rating of degradation based on distribution and severity of surface damage. Devices for applying controlled strain to both flat sheet and tubing samples are described. The new rating system consistently discriminated between. the effects of acetone pre-treatments, strain and exposure times in both in vitro and in vivo experiments. As expected, P80A underwent considerable stress cracking compared with P55D. P80A produced similar stress crack ratings in both in vivo and in vitro experiments, however P55D performed worse under in vitro conditions compared with in vivo. This result indicated that care must be taken when interpreting in vitro results in the absence of in vivo data. (C) 2001 Elsevier Science Ltd. All rights reserved.

Seismic imaging of the western Iberian crust using ambient noise: Boundaries and internal structure of the Iberian Massif

We present new Rayleigh-wave dispersion maps of the western Iberian Peninsula for periods between 8 and 30 s, obtained from correlations of seismic ambient noise, following the recent increase in seismic broadband network density in Portugal and Spain. Group velocities have been computed for each station pair using the empirical Green's functions generated by cross-correlating one-day-length seismic ambient-noise records. The resulting high-path density allows us to obtain lateral variations of the group velocities as a function of period in cells of 0.5 degrees x 0.5 degrees with an unprecedented resolution. As a result we were able to address some of the unknowns regarding the lithospheric structure beneath SW Iberia. The dispersion maps allow the imaging of the major structural units, namely the Iberian Massif, and the Lusitanian and Algarve Meso-Cenozoic basins. The Cadiz Gulf/Gibraltar Strait area corresponds to a strong low-velocity anomaly, which can be followed to the largest period inverted, although slightly shifted to the east at longer periods. Within the Iberian Massif, second-order perturbations in the group velocities are consistent with the transitions between tectonic units composing the massif. (C) 2013 Elsevier B.V. All rights reserved.

Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Ambient seismic noise monitoring of a clay landslide: Toward failure prediction

Given that clay-rich landslides may become mobilized, leading to rapid mass movements (earthflows and debris flows), they pose critical problems in risk management worldwide. The most widely proposed mechanism leading to such flow-like movements is the increase in water pore pressure in the sliding mass, generating partial or complete liquefaction. This solid-to-liquid transition results in a dramatic reduction of mechanical rigidity in the liquefied zones, which could be detected by monitoring shear wave velocity variations. With this purpose in mind, the ambient seismic noise correlation technique has been applied to measure the variation in the seismic surface wave velocity in the Pont Bourquin landslide (Swiss Alps). This small but active composite earthslide-earthflow was equipped with continuously recording seismic sensors during spring and summer 2010. An earthslide of a few thousand cubic meters was triggered in mid-August 2010, after a rainy period. This article shows that the seismic velocity of the sliding material, measured from daily noise correlograms, decreased continuously and rapidly for several days prior to the catastrophic event. From a spectral analysis of the velocity decrease, it was possible to determine the location of the change at the base of the sliding layer. These results demonstrate that ambient seismic noise can be used to detect rigidity variations before failure and could potentially be used to predict landslides.

Assessment of Nondestructive Testing Technologies for Quality Control/Quality Assurance of Asphalt Mixtures, TR-653, 2015

Asphalt pavements suffer various failures due to insufficient quality within their design lives. The American Association of State Highway and Transportation Officials (AASHTO) Mechanistic-Empirical Pavement Design Guide (MEPDG) has been proposed to improve pavement quality through quantitative performance prediction. Evaluation of the actual performance (quality) of pavements requires in situ nondestructive testing (NDT) techniques that can accurately measure the most critical, objective, and sensitive properties of pavement systems. The purpose of this study is to assess existing as well as promising new NDT technologies for quality control/quality assurance (QC/QA) of asphalt mixtures. Specifically, this study examined field measurements of density via the PaveTracker electromagnetic gage, shear-wave velocity via surface-wave testing methods, and dynamic stiffness via the Humboldt GeoGauge for five representative paving projects covering a range of mixes and traffic loads. The in situ tests were compared against laboratory measurements of core density and dynamic modulus. The in situ PaveTracker density had a low correlation with laboratory density and was not sensitive to variations in temperature or asphalt mix type. The in situ shear-wave velocity measured by surface-wave methods was most sensitive to variations in temperature and asphalt mix type. The in situ density and in situ shear-wave velocity were combined to calculate an in situ dynamic modulus, which is a performance-based quality measurement. The in situ GeoGauge stiffness measured on hot asphalt mixtures several hours after paving had a high correlation with the in situ dynamic modulus and the laboratory density, whereas the stiffness measurement of asphalt mixtures cooled with dry ice or at ambient temperature one or more days after paving had a very low correlation with the other measurements. To transform the in situ moduli from surface-wave testing into quantitative quality measurements, a QC/QA procedure was developed to first correct the in situ moduli measured at different field temperatures to the moduli at a common reference temperature based on master curves from laboratory dynamic modulus tests. The corrected in situ moduli can then be compared against the design moduli for an assessment of the actual pavement performance. A preliminary study of microelectromechanical systems- (MEMS)-based sensors for QC/QA and health monitoring of asphalt pavements was also performed.

Reflective Crack Mitigation Guide for Flexible Pavements, TR-641, 2015

Reflective cracks form in pavements when hot-mix asphalt (HMA) overlays are placed over jointed and/or severely cracked rigid and flexible pavements. In the first part of the research, survival analysis was conducted to identify the most appropriate rehabilitation method for composite pavements and to evaluate the influence of different factors on reflective crack development. Four rehabilitation methods, including mill and fill, overlay, heater scarification (SCR), and rubblization, were analyzed using three performance indicators: reflective cracking, international roughness index (IRI), and pavement condition index (PCI). It was found that rubblization can significantly retard reflective cracking development compared to the other three methods. No significant difference for PCI was seen among the four rehabilitation methods. Heater scarification showed the lowest survival probability for both reflective cracking and IRI, while an overlay resulted in the poorest overall pavement condition based on PCI. In addition, traffic level was found not to be a significant factor for reflective cracking development. An increase in overlay thickness can significantly delay the propagation of reflective cracking for all four treatments. Soil types in rubblization pavement sites were assessed, and no close relationship was found between rubblized pavement performance and subgrade soil condition. In the second part of the research, the study objective was to evaluate the modulus and performance of four reflective cracking treatments: full rubblization, modified rubblization, crack and seat, and rock interlayer. A total of 16 pavement sites were tested by the surface wave method (SWM), and in the first four sites both falling weight deflectometer (FWD) and SWM were conducted for a preliminary analysis. The SWM gave close concrete layer moduli compared to the FWD moduli on a conventional composite pavement. However, the SWM provided higher moduli for the rubblized concrete layer. After the preliminary analysis, another 12 pavement sites were tested by the SWM. The results showed that the crack and seat method provided the highest moduli, followed by the modified rubblization method. The full rubblization and the rock interlayer methods gave similar, but lower, moduli. Pavement performance surveys were also conducted during the field study. In general, none of the pavement sites had rutting problems. The conventional composite pavement site had the largest amount of reflective cracking. A moderate amount of reflective cracking was observed for the two pavement sites with full rubblization. Pavements with the rock interlayer and modified rubblization treatments had much less reflective cracking. It is recommended that use of the modified rubblization and rock interlayer treatments for reflective cracking mitigation are best.

Evolution of superficial lake water temperature profile under diurnal radiative forcing

In lentic water bodies, such as lakes, the water temperature near the surface typically increases during the day, and decreases during the night as a consequence of the diurnal radiative forcing (solar and infrared radiation). These temperature variations penetrate vertically into the water, transported mainly by heat conduction enhanced by eddy diffusion, which may vary due to atmospheric conditions, surface wave breaking, and internal dynamics of the water body. These two processes can be described in terms of an effective thermal diffusivity, which can be experimentally estimated. However, the transparency of the water (depending on turbidity) also allows solar radiation to penetrate below the surface into the water body, where it is locally absorbed (either by the water or by the deployed sensors). This process makes the estimation of effective thermal diffusivity from experimental water temperature profiles more difficult. In this study, we analyze water temperature profiles in a lake with the aim of showing that assessment of the role played by radiative forcing is necessary to estimate the effective thermal diffusivity. To this end we investigate diurnal water temperature fluctuations with depth. We try to quantify the effect of locally absorbed radiation and assess the impact of atmospheric conditions (wind speed, net radiation) on the estimation of the thermal diffusivity. The whole analysis is based on the results of fiber optic distributed temperature sensing, which allows unprecedented high spatial resolution measurements (∼4 mm) of the temperature profile in the water and near the water surface.