Correcting an acoustic wavefield for elastic effects

Finite-difference simulations are an important tool for studying elastic and acoustic wave propagation, but remain computationally challenging for elastic waves in three dimensions. Computations for acoustic waves are significantly simpler as they require less memory and operations per grid cell, and more significantly can be performed with coarser grids, both in space and time. In this paper, we present a procedure for correcting acoustic simulations for some of the effects of elasticity, at a cost considerably less than full elastic simulations. Two models are considered: the full elastic model and an equivalent acoustic model with the same P velocity and density. In this paper, although the basic theory is presented for anisotropic elasticity, the specific examples are for an isotropic model. The simulations are performed using the finite-difference method, but the basic method could be applied to other numerical techniques. A simulation in the acoustic model is performed and treated as an approximate solution of the wave propagation in the elastic model. As the acoustic solution is known, the error to the elastic wave equations can be calculated. If extra sources equal to this error were introduced into the elastic model, then the acoustic solution would be an exact solution of the elastic wave equations. Instead, the negative of these sources is introduced into a second acoustic simulation that is used to correct the first acoustic simulation. The corrected acoustic simulation contains some of the effects of elasticity without the full cost of an elastic simulation. It does not contain any shear waves, but amplitudes of reflected P waves are approximately corrected. We expect the corrected acoustic solution to be useful in regions of space and time around a P-wave source, but to deteriorate in some regions, for example, wider angles, and later in time, or after propagation through many interfaces. In this paper, we outline the theory of the correction method, and present results for simulations in a 2-D model with a plane interface. Reflections from a plane interface are simple enough that an analytic analysis is possible, and for plane waves, we give the correction to the acoustic reflection and transmission coefficients. Finally, finite-difference calculations for plane waves are used to confirm the analytic results. Results for wave propagation in more complicated, realistic models will be presented elsewher

Illness perceptions and work participation: a systematic review

Self-regulatory processes play an important role in mediating between the disease and the health outcomes, and potentially also work outcomes. This systematic review aims to explore the relationship between illness perceptions and work participation in patients with somatic diseases and complaints. The bibliographic databases Medline, PsycINFO and Embase were searched from inception to March 2008. Included were cross-sectional or longitudinal studies, patients with somatic diseases or complaints, illness perceptions based on at least four dimensions of the common sense model of self-regulation, and work participation. Two longitudinal and two cross-sectional studies selected for this review report statistically significant findings for one or more illness perception dimensions in patients with various complaints and illnesses, although some dimensions are significant in one study but not in another. Overall, non-working patients perceived more serious consequences, expected their illness to last a longer time, and reported more symptoms and more emotional responses as a result of their illness. Alternatively, working patients had a stronger belief in the controllability of their condition and a better understanding of their disease. The limited number of studies in this review suggests that illness perceptions play a role in the work participation of patients with somatic diseases or complaints, although it is not clear how strong this relationship is and which illness perception dimensions are most useful. Identifying individuals with maladaptive illness perceptions and targeting interventions toward changing these perceptions are promising developments in improving work participatio

Subcellular peptide localization in single identified neurons by capillary microsampling mass spectrometry

Single cell mass spectrometry (MS) is uniquely positioned for the sequencing and identification of peptides in rare cells. Small peptides can take on different roles in subcellular compartments. Whereas some peptides serve as neurotransmitters in the cytoplasm, they can also function as transcription factors in the nucleus. Thus, there is a need to analyze the subcellular peptide compositions in identified single cells. Here, we apply capillary microsampling MS with ion mobility separation for the sequencing of peptides in single neurons of the mollusk Lymnaea stagnalis, and the analysis of peptide distributions between the cytoplasm and nucleus of identified single neurons that are known to express cardioactive Phe-Met-Arg-Phe amide-like (FMRFamide-like) neuropeptides. Nuclei and cytoplasm of Type 1 and Type 2 F group (Fgp) neurons were analyzed for neuropeptides cleaved from the protein precursors encoded by alternative splicing products of the FMRFamide gene. Relative abundances of nine neuropeptides were determined in the cytoplasm. The nuclei contained six of these peptides at different abundances. Enabled by its relative enrichment in Fgp neurons, a new 28-residue neuropeptide was sequenced by tandem MS

Evaluation of a novel magneto-optical method for the detection of malaria parasites

Improving the efficiency of malaria diagnosis is one of the main goals of current malaria research. We have recently developed a magneto-optical (MO) method which allows high-sensitivity detection of malaria pigment (hemozoin crystals) in blood via the magnetically induced rotational motion of the hemozoin crystals. Here, we evaluate this MO technique for the detection of Plasmodium falciparum in infected erythrocytes using in-vitro parasite cultures covering the entire intraerythrocytic life cycle. Our novel method detected parasite densities as low as approximately 40 parasites per microliter of blood (0.0008% parasitemia) at the ring stage and less than 10 parasites/microL (0.0002% parasitemia) in the case of the later stages. These limits of detection, corresponding to approximately 20 pg/microL of hemozoin produced by the parasites, exceed that of rapid diagnostic tests and compete with the threshold achievable by light microscopic observation of blood smears. The MO diagnosis requires no special training of the operator or specific reagents for parasite detection, except for an inexpensive lysis solution to release intracellular hemozoin. The devices can be designed to a portable format for clinical and in-field tests. Besides testing its diagnostic performance, we also applied the MO technique to investigate the change in hemozoin concentration during parasite maturation. Our preliminary data indicate that this method may offer an efficient tool to determine the amount of hemozoin produced by the different parasite stages in synchronized cultures. Hence, it could eventually be used for testing the susceptibility of parasites to antimalarial drugs

A three-dimensional seismic tomographic study of the gas hydrate stability zone, offshore Vancouver Island

Methane hydrate bottom-simulating reflectors (BSRs) are widespread on the northern Cascadia margin offshore Vancouver Island. We conducted a three-dimensional tomographic seismic study of the hydrate stability zone in an area around Ocean Drilling Program Site 889 using two deployments of five ocean bottom hydrophones and air gun shots along a series of closely spaced profiles in various orientations. Further constraints on reflector geometry come from coincident single-channel reflection profiles. Travel times of reflected and refracted phases were inverted with a regularized three-dimensional inversion using perturbation ray tracing through smooth isotropic media for the forward step. The seismic data allow us to constrain the velocity structure in a ∼6 km2 area around the drill site. Mean velocities range from 1.50 km s−1 at the seabed to 1.84 km s−1 at the BSR, and velocities at Site 889 match well those measured using a vertical seismic profile. At equivalent depths below the seafloor, velocities vary laterally by typically ∼0.15 km s−1. Close to the seafloor, velocities may be controlled primarily by lithology, but close to the BSR we infer hydrate contents of up to 15% of the pore space from effective medium modeling. The mean hydrate saturation in the well-constrained volume of the velocity model is estimated to be 2.2%. There is no correlation between the seismic velocity above the BSR and the reflection coefficient at the BSR, so the latter is likely controlled primarily by the distribution of free gas beneath the hydrate stability zone

Three-dimensional tomographic inversion of combined reflection and refraction seismic traveltime data

A tomographic inversion method is presented for the determination of 3-D velocity and interface structure from a wide range of body-wave seismic traveltime data types. It is applicable to refraction, wide-angle reflection, normal-incidence and multichannel seismic data, and is best suited to a combination of these that provides good independent constraints on seismic velocities and interface depths. The inversion process seeks a layer-interface minimum-structure model that is able to explain the given data satisfactorily by inverting to minimize data misfit and model roughness norms simultaneously. This regularized inversion, and the use of smooth functions to describe velocities and depths, allows the highly non-linear tomographic problem to be approximated as a series of linear steps. The inversion process begins by optimizing the fit to the data of a highly-smoothed initial model. In each subsequent step, structure is allowed to develop in the model with successively greater detail evolving until a satisfactory fit to the data is obtained. Parameter uncertainties for the final model are then estimated using an a posteriori covariance matrix analysis. Smooth layer-interface models are parametrized using regular grids of velocity and depth nodes from which spline-interpolated interface surfaces and velocity fields are defined. Forward modelling is achieved using ray perturbation theory and a two-point ray tracing method that is optimized for a large number of closely-spaced shot or receiver points. The method may be used to generate 1- and 2-D models (from, for example vertical seismic profile data or 2-D surveys) in which the 3-D geometry of a survey is correctly accounted for. The ability of the method to resolve typical target structures is tested in a synthetic salt dome inversion. From a set of noisy traveltime data, the model converges quickly to a well-resolved final model from different starting models. The application of this method to real data is demonstrated with a combined 3-D inversion of refraction and reflection data which provide P-wave velocity constraints on the methane hydrate stability zone in the Cascadia Margin offshore Vancouver Island

A three-dimensional seismic tomographic study of the gas hydrate stability zone, offshore Vancouver Island

Methane hydrate bottom-simulating reflectors (BSRs) are widespread on the northern Cascadia margin offshore Vancouver Island. We conducted a three-dimensional tomographic seismic study of the hydrate stability zone in an area around Ocean Drilling Program Site 889 using two deployments of five ocean bottom hydrophones and air gun shots along a series of closely spaced profiles in various orientations. Further constraints on reflector geometry come from coincident single-channel reflection profiles. Travel times of reflected and refracted phases were inverted with a regularized three-dimensional inversion using perturbation ray tracing through smooth isotropic media for the forward step. The seismic data allow us to constrain the velocity structure in a ?6 km2 area around the drill site. Mean velocities range from 1.50 km s?1 at the seabed to 1.84 km s?1 at the BSR, and velocities at Site 889 match well those measured using a vertical seismic profile. At equivalent depths below the seafloor, velocities vary laterally by typically ?0.15 km s?1. Close to the seafloor, velocities may be controlled primarily by lithology, but close to the BSR we infer hydrate contents of up to 15% of the pore space from effective medium modeling. The mean hydrate saturation in the well-constrained volume of the velocity model is estimated to be 2.2%. There is no correlation between the seismic velocity above the BSR and the reflection coefficient at the BSR, so the latter is likely controlled primarily by the distribution of free gas beneath the hydrate stability zone