Randomized Two-Process Wait-Free Test-and-Set

We present the first explicit, and currently simplest, randomized algorithm for 2-process wait-free test-and-set. It is implemented with two 4-valued single writer single reader atomic variables. A test-and-set takes at most 11 expected elementary steps, while a reset takes exactly 1 elementary step. Based on a finite-state analysis, the proofs of correctness and expected length are compressed into one table.Comment: 9 pages, 4 figures, LaTeX source; Submitte

Full Waveform Inversion for Time-Distance Helioseismology

Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions, and likely for supergranular flows as well. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the sub-surface structure of large-amplitude perturbations. Defining the penalty functional as the $L_2$ norm of wave travel-time deviations, we compute the the total misfit gradient of this functional with respect to the relevant model parameters %(only sound speed in this case) at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or quasi-Newton limited-memory BFGS. Performing nonlinear iterations requires privileging pixels (such as those in the near-field of the scatterer), a practice not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a 2-D inverse problem with only a sound-speed perturbation.Comment: 24 pages, 10 figures, to appear in Ap

Performance of Two 18-Story Steel Moment-Frame Buildings in Southern California During Two Large Simulated San Andreas Earthquakes

Using state-of-the-art computational tools in seismology and structural engineering, validated using data from the Mw=6.7 January 1994 Northridge earthquake, we determine the damage to two 18-story steel moment-frame buildings, one existing and one new, located in southern California due to ground motions from two hypothetical magnitude 7.9 earthquakes on the San Andreas Fault. The new building has the same configuration as the existing building but has been redesigned to current building code standards. Two cases are considered: rupture initiating at Parkfield and propagating from north to south, and rupture propagating from south to north and terminating at Parkfield. Severe damage occurs in these buildings at many locations in the region in the north-to-south rupture scenario. Peak velocities of 1 m.s−1 and 2 m.s−1 occur in the Los Angeles Basin and San Fernando Valley, respectively, while the corresponding peak displacements are about 1 m and 2 m, respectively. Peak interstory drifts in the two buildings exceed 0.10 and 0.06 in many areas of the San Fernando Valley and the Los Angeles Basin, respectively. The redesigned building performs significantly better than the existing building; however, its improved design based on the 1997 Uniform Building Code is still not adequate to prevent serious damage. The results from the south-to-north scenario are not as alarming, although damage is serious enough to cause significant business interruption and compromise life safety

Algorithmic Statistics

While Kolmogorov complexity is the accepted absolute measure of information content of an individual finite object, a similarly absolute notion is needed for the relation between an individual data sample and an individual model summarizing the information in the data, for example, a finite set (or probability distribution) where the data sample typically came from. The statistical theory based on such relations between individual objects can be called algorithmic statistics, in contrast to classical statistical theory that deals with relations between probabilistic ensembles. We develop the algorithmic theory of statistic, sufficient statistic, and minimal sufficient statistic. This theory is based on two-part codes consisting of the code for the statistic (the model summarizing the regularity, the meaningful information, in the data) and the model-to-data code. In contrast to the situation in probabilistic statistical theory, the algorithmic relation of (minimal) sufficiency is an absolute relation between the individual model and the individual data sample. We distinguish implicit and explicit descriptions of the models. We give characterizations of algorithmic (Kolmogorov) minimal sufficient statistic for all data samples for both description modes--in the explicit mode under some constraints. We also strengthen and elaborate earlier results on the ``Kolmogorov structure function'' and ``absolutely non-stochastic objects''--those rare objects for which the simplest models that summarize their relevant information (minimal sufficient statistics) are at least as complex as the objects themselves. We demonstrate a close relation between the probabilistic notions and the algorithmic ones.Comment: LaTeX, 22 pages, 1 figure, with correction to the published journal versio

An investigation of the hydration properties of chemically vapour-deposited silicon dioxide films by means of ellipsometry

Time-dependent ellipsometric measurements were made of the water absoprtion in SiO2 films grown by chemical vapour deposition at low temperatures. These films were exposed to water vapour at different pressurures. We determined in the ellipsometric volume percentage of water in the oxide film from the changes in the ellipsometric parameters due to water absorption. We also obtained a diffusion coefficient of water of the order of 1.3 x 10-13 cm2 s-1. From the reversible changes in the ellipsometric parameters due to dehydration we conclude that the water is physically adsorbed

Waveform modeling of the slab beneath Japan

The tomographic P wave model for the Japan subduction zone derived by Zhao et al. (1994) has two very striking features: a slab about 90 km thick with P wave velocities 3–6% higher than the surrounding mantle and a mantle wedge with −6% low-velocity anomalies. We study three-component seismograms from more than 600 Hi-net stations produced by two earthquakes which occurred in the downgoing Pacific Plate at depths greater than 400 km. We simulate body wave propagation in the three-dimensional (3-D) P wave model using 2-D finite difference (FDM) and 3-D spectral element (SEM) methods. As measured by cross correlation between synthetics and data, the P wave model typically explains about half of the traveltime anomaly and some of the waveform complexity but fails to predict the extended SH wave train. In this study we take advantage of the densely distributed Hi-net stations and use 2-D FDM modeling to simulate the P-SV and SH waveforms. Our 2-D model suggests that a thin, elongated low-velocity zone exists atop the slab, extending down to a depth of 300 km with an S wave velocity reduction of 14% if a thickness of 20 km is assumed. Further, 3-D SEM simulations confirm that this model explains a strong secondary arrival which cannot easily be imaged with standard tomographic techniques. The low-velocity layer could explain the relatively weak coupling associated with most subduction zones at shallow depths (<50 km), generally involving abundant volcanic activity and silent earthquakes, and it may also help to further our understanding of the water-related phase transition of ultramafic rocks, and the nature of seismicity at intermediate depths (~70–300 km)

Tsunami generation by ocean floor rupture front propagation: Hamiltonian description

The Hamiltonian method is applied to the problem of tsunami generation caused by a propagating rupture front and deformation of the ocean floor. The method establishes an alternative framework for analyzing the tsunami generation process and produces analytical expressions for the power and directivity of tsunami radiation (in the far-field) for two illustrative cases, with constant and gradually varying speeds of rupture front propagation

Finite-frequency sensitivity of body waves to anisotropy based upon adjoint methods

We investigate the sensitivity of finite-frequency body-wave observables to mantle anisotropy based upon kernels calculated by combining adjoint methods and spectral-element modelling of seismic wave propagation. Anisotropy is described by 21 density-normalized elastic parameters naturally involved in asymptotic wave propagation in weakly anisotropic media. In a 1-D reference model, body-wave sensitivity to anisotropy is characterized by ‘banana–doughnut’ kernels which exhibit large, path-dependent variations and even sign changes. P-wave travel-times appear much more sensitive to certain azimuthally anisotropic parameters than to the usual isotropic parameters, suggesting that isotropic P-wave tomography could be significantly biased by coherent anisotropic structures, such as slabs. Because of shear-wave splitting, the common cross-correlation travel-time anomaly is not an appropriate observable for S waves propagating in anisotropic media. We propose two new observables for shear waves. The first observable is a generalized cross-correlation travel-time anomaly, and the second a generalized ‘splitting intensity’. Like P waves, S waves analysed based upon these observables are generally sensitive to a large number of the 21 anisotropic parameters and show significant path-dependent variations. The specific path-geometry of SKS waves results in favourable properties for imaging based upon the splitting intensity, because it is sensitive to a smaller number of anisotropic parameters, and the region which is sampled is mainly limited to the upper mantle beneath the receiver