SEARCH · Search NASA

Results for “surface‐wave dispersion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗

DOE OSTI2025

Local-Regional Surface Wave Dispersion Analysis using Ambient Noise Data of the Southern Nevada National Security Site Area

This report details seismic ambient noise analysis to improve seismic velocity estimates of the southeastern Nevada National Security Site (NNSS). We compare two different methods for estimating surface wave dispersion curves from ambient noise cross-correlations: frequency time analysis (FTAN) and Aki’s cross-spectral method (XSpec). We find that XSpec performs better for our local-regional dataset and frequency content. Using phase velocity estimates from XSpec, we build a preliminary phase velocity dispersion dataset, which contains data for 1,054 station pairs and 8,905 discrete phase velocity measurements for periods between 0.5 and 13 s. This phase velocity dataset will be incorporated into an updated local-regional P and S wave velocity model of the southeastern NNSS in the future and is expected to improve upon shallow velocity estimates.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Excitation of surface waves in 3D ion beam neutralization

Neutralization of beams with 2D and 3D geometries by the electrons emitted from an external source is studied using particle-in-cell simulations. Our work reveals that the high-energy electrons excite Trivelpiece–Gould (TG) surface waves in the beams with 3D axisymmetric geometries. These high-energy electrons are generated because of a large amplitude electrostatic solitary wave (ESW) that forms near the electron source and has an electric potential amplitude more than three times the electron thermal energy. Here, we also find that surface wave excitation only happens when the beam radius is large enough at the ion source to attract enough electrons that could form the large amplitude ESW. A comparison of the 3D TG surface wave dispersion relation with an expression for 2D surface waves reveals that they become excited in 3D axisymmetric but not in 2D planar beam because of a higher phase speed requirement in the latter case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

CANVAS: An Adjoint Waveform Tomography Model of California and Nevada

Abstract We present the California‐Nevada Adjoint Simulations (CANVAS) model, an adjoint waveform tomography model of the crust and uppermost mantle of the states of California and Nevada. We used WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) as the starting model and iteratively decreased the minimum period of CANVAS from 30 to 12 s. CANVAS was iterated in two distinct stages: the first stage with source mechanisms from the Global Centroid Moment Tensor (GCMT) catalog and the second stage with inverted moment tensors (MT) using the CANV_WUS model (Doody et al., 2023, https://doi.org/10.1029/2023jb026463 ). We show that updating the MTs with 3D Green's functions improved waveform fits and azimuthal coverage of windowed data used to calculate the gradients. As for the model itself, we improved waveform fits over WUS256, particularly in the dispersed surface waves. CANVAS resolved tectonic features seen in other models and accurately defined the depth to basement of major basins, including the Central Valley and the Ventura Basin. We propose CANVAS as a starting model for crustal tomography models on smaller scales.

58 GEOSCIENCES↗

DOE OSTI2025

Seismic imaging of the Ecuadorian forearc and arc from joint ambient noise, local, and teleseismic tomography: catching the Nazca slab in the act of flattening

SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.

58 GEOSCIENCES↗

DOE OSTI2022

Unified Wind-Wave Growth and Spectrum Functions for All Water Depths: Field Observations and Model Results

Abstract Wind-wave development is governed by the fetch- or duration-limited growth principle that is expressed as a pair of similarity functions relating the dimensionless elevation variance (wave energy) and spectral peak frequency to fetch or duration. Combining the pair of similarity functions, the fetch or duration variable can be removed to form a dimensionless function of elevation variance and spectral peak frequency, which is interpreted as the wave energy evolution with wave age. The relationship is initially developed for quasi-neural stability and quasi-steady wind forcing conditions. Further analyses show that the same fetch, duration, and wave-age similarity functions are applicable to unsteady wind forcing conditions, including rapidly accelerating and decelerating mountain gap wind episodes and tropical cyclone (TC) wind fields. Here it is shown that with the dimensionless frequency converted to dimensionless wavenumber using the surface wave dispersion relationship, the same similarity function is applicable in all water depths. Field data collected in shallow to deep waters and mild to TC wind conditions and synthetic data generated by spectrum model computations are assembled to illustrate the applicability. For the simulation work, the finite-depth wind-wave spectrum model and its shoaling function are formulated for variable spectral slopes. Given wind speed, wave age, and water depth, the measured and spectrum-computed significant wave heights and the associated growth parameters are in good agreement in forcing conditions from mild to TC winds and in all depths from deep ocean to shallow lake. Significance Statement This paper presents a growth function and spectrum model to describe wind-wave development in all water depths. Their applicability covers a wide range of wind forcing conditions including steady, accelerating, decelerating, and tropical cyclone events. Support for the unified spectrum model and growth function is presented with field observations and numerical computations.

Hwang, Paul A.↗

DOE OSTI2021

Report on the LLNL Global Full-waveform Inversion Workflow and Progress

LLNL has developed the SPiRaL global seismic tomography model based upon millions of body wave travel times and surface wave dispersion curves in the period range of 25-200 seconds (Simmons et al., 2021). The motivation to construct the SPiRaL model (which stands for S-, P-, Rayleigh, and Love waves) is based upon the guidance provided by a 3-D model workshop held in Berkeley in 2007 and documented in a paper by Zucca et al. (2009). It was concluded in Zucca et al. (2009) that a “practical and obtainable” model is one that predicts body wave travel times for event location while also predicting surface wave velocities for moment tensor estimation. We have largely followed that guidance which has led to the construction of SPiRaL. Looking forward, we now seek to improve the global model through full waveform inversion (FWI). This necessitates the development/customization of an efficient and flexible modeling workflow to automate the process using LLNL’s supercomputing platforms. This report briefly outlines the work that being performed at LLNL to develop the necessary tools and framework to update the SPiRaL global model with full waveform information.

58 GEOSCIENCES↗

DOE OSTI2021

A New Approach to Predict Hydrogeological Parameters Using Shear Waves from the Multichannel Analysis of Surface Waves Method

For near-surface contaminant characterization, the accurate prediction of hydrogeological parameters in anisotropic and heterogeneous environments has been a challenge since the last decades. However, recent advances in near-surface geophysics have facilitated the use of geophysical data for hydrogeological characterization in the last few years. A pseudo 3-D high resolution P-wave shallow seismic reflection survey was performed at the P Reactor Area, Savannah River Site, South Carolina in order to delineate and predict migration pathways of a large contaminant plume including trichloroethylene. This contaminant plume originates from the northwest section of the reactor facility that is located within the Upper Atlantic Coastal Plain. The data were collected with 40 Hz geophones, an accelerated weight-drop as seismic source and 1 m receiver spacing with near- and far-offsets of 0.5 and 119.5 m, respectively. In such areas with near-surface contaminants, a detailed subsurface characterization of the vadose zone hydraulic parameters is very important. Indeed, an inexpensive method of deriving such parameters by the use of seismic reflection surveys is beneficial, and our approach uses the relationship between seismic velocity and hydrogeological parameters together with empirical observations relating porosity to permeability and hydraulic conductivity. Shear wave velocity ( V s ) profiles were estimated from surface wave dispersion analysis of the seismic reflection data and were subsequently used to derive hydraulic parameters such as porosity, permeability, and hydraulic conductivity. Additional geophysical data including core samples, vertical seismic profiling, surface electrical resistivity tomography, natural gamma and electrical resistivity logs allowed for a robust assessment of the validity and geological significance of the estimated V s and hydrogeological models. The results demonstrate the usefulness of this approach for the upper 15 m of shallow unconsolidated sediments even though the survey design parameters were not optimal for surface wave analysis due to the higher than desired frequency geophones.

Engineering↗

DOE OSTI2023

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Pinzón, Juan I. (ORCID:0000000338078068)↗

DOE OSTI2021

Data–Driven Velocity Model Evaluation Using K–Means Clustering

In this work, we develop a data-driven clustering method to evaluate a velocity model using surface wave velocity dispersion. This is done by first computing theoretical dispersion curves for 1-D velocity profiles of all the grid locations and then splitting the resulting dispersion curves into a certain number of groups via the K-means clustering. The observed dispersion curves are also clustered following the same procedure and the velocity model is assessed by comparing the spatial patterns obtained for the observed and synthetic data sets. The method is applied to evaluate two community velocity models in southern California, CVM-S4.26 and CVM-H15.1, using phase velocity maps derived for 3–16 s Rayleigh waves. We found a good correlation in the spatial distribution of clusters between the result of CVM-S4.26 and that of the observed data, suggesting that the CVM-S4.26 fits the observed dispersion maps better than the CVM-H15.1 in terms of features extracted from the clustering analysis.

58 GEOSCIENCES↗

DOE OSTI2022

WUS256: An Adjoint Waveform Tomography Model of the Crust and Upper Mantle of the Western United States for Improved Waveform Simulations

Abstract We report a new model (WUS256) of radially anisotropic seismic wavespeeds of the crust and upper mantle of the western United States (WUS) obtained from adjoint waveform tomography for the purpose of improving synthetic waveform fits to observed data. WUS256 is based on inversion of over 94,000 waveforms from 72 earthquakes recorded by nearly 3,400 stations. We started with the SPiRaL global model (Simmons et al., 2021, https://doi.org/10.1093/gji/ggab277 ) and waveforms in the period band of 50–120 s. We followed a conservative multiscale inversion approach with eight stages and 256 total inversion iterations which enabled monotonic misfit reduction to 20‐s minimum‐period waves. WUS256 relied on time‐frequency (TF) phase misfits and a trust region limited memory Broyden–Fletcher–Goldfarb–Shanno (L‐BFGS) optimization. Hessian‐vector products were used to qualitatively assess model resolution. Results indicate that WUS256 has good coverage of the continental regions to depths of about 150 km and is able to resolve features on lateral scales of about 200 km. We quantify waveform fits by the reduction in TF and normalized amplitude difference misfits between WUS256 and the SPiRaL starting model. WUS256 significantly improves waveform fits with misfit reduction 64% for both inversion and validation data sets compared to the SPiRaL starting model and shows even better fits compared to other models. Waveform fits illustrate that WUS256 reproduces body‐waves, fundamental mode surface waves as well as late arriving dispersed and/or scattered short period surface waves. The improvement in waveform fit indicates that WUS256 can be used to reproduce path effects on regional complete waveforms and moment tensor inversions.

58 GEOSCIENCES↗

DOE OSTI2022

Probabilistic inversion of circular phase spectra: application to two-station phase-velocity dispersion estimation in western Canada

SUMMARY Periodic directional and temporal measurements are common in seismology, and necessitate specific statistical analyses that are appropriate for circular quantities. In this work, we explore the use of a von Mises distribution as a representation of errors on circular seismological observations. Specifically, we automate the estimation of surface-wave phase-velocity dispersion for the teleseismic two-station method, which generally suffers from a 2π phase ambiguity. The use of Bayesian inverse techniques, which aim to rigorously quantify model parameter uncertainty, have become widespread throughout seismology over the last decade. Here, we apply Bayesian inversion to measurements of surface-wave phase spectra in order to estimate 1-D, path-averaged Earth structure between station pairs. The dispersion curve and associated uncertainties are additional results of the inversion, which can then be used as input for subsequent analyses (e.g. tomography). We demonstrate this technique through application to surface-wave recordings from long-running seismic stations throughout western Canada. Our results for over 10 000 station pairs reveal first-order tectonic features consistent with previous studies, which provides confidence in our approach as well as an initial step towards resolving a full 3-D seismic velocity model for the region. This work also presents a foundation for the inversion of surface-wave phase spectra to estimate 3-D Earth structure directly. Finally, the ideas presented in this work are not limited to the inversion of surface-wave phase spectra, but can also be considered for Bayesian geophysical inversion of any circular quantities.

Gosselin, Jeremy M. (ORCID:0000000203754102)↗

DOE OSTI2023

High-Resolution Near-Surface Imaging at the Basin Scale Using Dark Fiber and Distributed Acoustic Sensing: Toward Site Effect Estimation in Urban Environments

Near-surface seismic structure, particularly the shear wave velocity (V s ), can strongly affect local site response, and should be accurately estimated for ground motion prediction during seismic hazard assessment. The Imperial Valley (California), occupying the southern end of the Salton Trough, is a seismically active basin with thick surficial lacustrine sedimentary deposits. In this study, we utilize ambient noise records and local earthquake events for high-resolution near-surface characterization and site effect estimation with an unlit fiber-optic telecommunication infrastructure (dark fiber) in Imperial Valley by using the distributed acoustic sensing (DAS) technique. We apply ambient noise interferometry to retrieve coherent surface waves from DAS records, and evaluate performances of three different surface wave methods on DAS ambient noise dispersion imaging. We develop a quality control workflow to improve the dispersion measurement of noisy portions of the DAS data set by using a data selection strategy. Using the joint inversion of both the fundamental mode and higher overtones of Rayleigh waves, a high resolution two-dimensional (2D) V s structure down to 70 m depth is obtained. We successfully achieve an improved V s 30 (the time-averaged shear-wave velocity in the top 30 m) model with higher spatial-resolution and reliability compared to the existing community model for the area. We also explore the potential for utilizing DAS earthquake events for site amplification estimation. The preliminary results reveal a clear anti-correlation between the approximated site response and the V s 30 profile. In conclusion, our results indicate the potential utility of DAS deployed on dark fiber for near-surface characterization in appropriate contexts.

58 GEOSCIENCES↗

DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗

DOE OSTI2025

A curvature-based mechanism for the spontaneous emission of electromagnetic radiation at curved plasma–vacuum interfaces

In this work, we derive from first principles the dispersion relation for electromagnetic radiation that is spontaneously emitted at a 2D curved plasma–vacuum interface in the absence of both background and external electric and magnetic fields. This plasma geometry is motivated by plasmas generated through hypervelocity impact projectiles—an area that is of importance to satellite safety in the aerospace community, and may be of importance for the developing study of projectile impact-driven inertial confinement fusion. We find spontaneous radiation occurs when thermal surface waves propagate along the edge of a circular plasma, with a frequency lower than the plasma frequency and spatial decay proportional to 1/r. This is distinct from more well-known mechanisms of radiation, in which the radiation frequency is greater than or equal to the plasma frequency. This theory is validated against 2D, electromagnetic particle-in-cell simulations and shows reasonable agreement. Finally, we discuss the effects of plasma geometry and attribute the radiation generation mechanism to the acceleration of charged surface waves as they traverse azimuthally along a curved path.

Lau, Raymond (ORCID:0000000299318101)↗

DOE OSTI2021

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗

DOE OSTI2022

Real-time laser ultrasonic monitoring of laser-induced thermal processes

Abstract Intra- and inter-layer integrity of components fabricated with advanced manufacturing techniques, such as laser powder bed fusion, is dependent upon rapid heating, melting, and solidification processes. There is a need for new techniques to provide in situ feedback of these processes. Here a laser-based ultrasonic technique to probe thermal effects induced by a high-power continuous wave laser in titanium samples is described. Numerical simulations were performed to show that, for a spatially uniform heating beam, laser-induced surface acoustic waves are strongly influenced by surface heating conditions, are dispersive in the case of rapid heating, and that an abrupt velocity reduction happens upon the onset of surface melting. Furthermore, laser-based ultrasound experimental results which monitor the transient change of surface wave travel time associated with high power laser surface heating are provided. A pulsed laser is used to generate high frequency surface acoustic waves that propagate through the laser-heated region and are detected using a photorefractive crystal-based interferometer. Qualitative agreement is observed between theory and experiment with both showing a rapid reduction in the surface wave velocity at the onset of illumination and further decrease in surface wave velocity associated with melting. It is demonstrated that changes in the surface wave velocity can be used to track local heating and detect the onset of surface melting in real time.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Saturation and decay of an induced scattering instability driven by static field asymmetries in a non-neutral plasma

Experimental data on the saturation and decay of a previously described induced scattering instability are presented. The instability is driven by static background field asymmetries and initially produces a single growing mode (ω 1 , k 1 ). As the instability develops, several harmonically related (ω n =nω 1 , k n =nk 1 ) eigenmodes of the plasma column grow up, resulting in a spatially localized moving potential perturbation with δϕ/T ≈ 20%. After a period of rapid growth, the mode amplitudes saturate and vary slowly for roughly 10 5 ω 1 −1 . This quasi-equilibrium ends with an unusual change in the eigenmode phases, wherein δϕ changes from positive to negative. Plasma heating during this period increases the damping of all modes and rapidly quenches the instability. This instability is of particular interest because the coupling with the field asymmetry produces an enhanced level of radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗