Dynamics and kinematics of eruptive activity at Fuego volcano, Guatemala 2005-2009

Volcanoes are the surficial expressions of complex pathways that vent magma and gasses generated deep in the Earth. Geophysical data record at least the partial history of magma and gas movement in the conduit and venting to the atmosphere. This work focuses on developing a more comprehensive understanding of explosive degassing at Fuego volcano, Guatemala through observations and analysis of geophysical data collected in 2005 – 2009. A pattern of eruptive activity was observed during 2005 – 2007 and quantified with seismic and infrasound, satellite thermal and gas measurements, and lava flow lengths. Eruptive styles are related to variable magma flux and accumulation of gas. Explosive degassing was recorded on broadband seismic and infrasound sensors in 2008 and 2009. Explosion energy partitioning between the ground and the atmosphere shows an increase in acoustic energy from 2008 to 2009, indicating a shift toward increased gas pressure in the conduit. Very-long-period (VLP) seismic signals are associated with the strongest explosions recorded in 2009 and waveform modeling in the 10 – 30 s band produces a best-fit source location 300 m west and 300 m below the summit crater. The calculated moment tensor indicates a volumetric source, which is modeled as a dike feeding a SW-dipping (35°) sill. The sill is the dominant component and its projection to the surface nearly intersects the summit crater. The deformation history of the sill is interpreted as: 1) an initial inflation due to pressurization, followed by 2) a rapid deflation as overpressure is explosively release, and finally 3) a reinflation as fresh magma flows into the sill and degasses. Tilt signals are derived from the horizontal components of the seismometer and show repetitive inflation deflation cycles with a 20 minute period coincident with strong explosions. These cycles represent the pressurization of the shallow conduit and explosive venting of overpressure that develops beneath a partially crystallized plug of magma. The energy released during the strong explosions has allowed for imaging of Fuego’s shallow conduit, which appears to have migrated west of the summit crater. In summary, Fuego is becoming more gas charged and its summit centered vent is shifting to the west - serious hazard consequences are likely.

Optimal current control for sub-unity power factor correction

Switching mode power supplies (SMPS) are subject to low power factor and high harmonic distortions. Active power-factor correction (APFC) is a technique to improve the power factor and to reduce the harmonic distortion of SMPSs. However, this technique results in double frequency output voltage variation which can be reduced by using a large output capacitance. Using large capacitors increases the cost and size of the converter. Furthermore, the capacitors are subject to frequent failures mainly caused by evaporation of the electrolytic solution which reduce the converter reliability. This thesis presents an optimal control method for the input current of a boost converter to reduce the size of the output capacitor. The optimum current waveform as a function of weighing factor is found by using the Euler Lagrange equation. A set of simulations are performed to determine the ideal weighing which gives the lowest possible output voltage variation as the converter still meets the IEC-61000-3-2 class-A harmonics requirements with a power factor of 0.8 or higher. The proposed method is verified by the experimental work. A boost converter is designed and it is run for different power levels, 100 W, 200 W and 400 W. The desired output voltage ripple is 10 V peak to peak for the output voltage of 200 Vdc. This ripple value corresponds to a ± 2.5% output voltage ripple. The experimental and the simulation results are found to be quite matching. A significant reduction in capacitor size, as high as 50%, is accomplished by using the proposed method.

Production of liquid core-polymer shell microcapsules

Polylactide (PLA) is a biodegradable polymer that has been used in particle form for drug release, due to its biocompatibility, tailorable degradation kinetics, and desirable mechanical properties. Active pharmaceutical ingredients (APIs) may be either dissolved or encapsulated within these biomaterials to create micro- or nanoparticles. Delivery of an AIP within fine particles may overcome solubility or stability issues that can result in early elimination or degradation of the AIP in a hostile biological environment. Furthermore, it is a promising method for controlling the rate of drug delivery and dosage. The goal of this project is to develop a simple and cost-effective device that allows us to produce monodisperse micro- and nanocapsules with controllable size and adjustable sheath thickness on demand. To achieve this goal, we have studied the dual-capillary electrospray and pulsed electrospray. Dual-capillary electrospray has received considerable attention in recent years due to its ability to create core-shell structures in a single-step. However, it also increases the difficulty of controlling the inner and outer particle morphology, since two simultaneous flows are required. Conventional electrospraying has been mainly conducted using direct-current (DC) voltage with little control over anything but the electrical potential. In contrast, control over the input voltage waveform (i.e. pulsing) in electrospraying offers greater control over the process variables. Poly(L-lactic acid) (PLLA) microspheres and microcapsules were successfully fabricated via pulsed-DC electrospray and dual-capillary electrospray, respectively. Core shell combinations produced include: Water/PLLA, PLLA/polyethylene glycol (PEG), and oleic Acid/PLLA. In this study, we designed a novel high-voltage pulse forming network and a set of new designs for coaxial electrospray nozzles. We also investigated the effect of the pulsed voltage characteristics (e.g. pulse frequency, pulse amplitude and pulse width) on the particle’s size and uniformity. We found that pulse frequency, pulse amplitude, pulse width, and the combinations of these factors had a statistically significant effect on the particle’s size. In addition, factors such as polymer concentration, solvent type, feed flow rate, collection method, temperature, and humidity can significantly affect the size and shape of the particles formed.

Parametric estimation of pulse arrival time: a robust approach to pulse wave velocity

Pulse wave velocity (PWV) is a surrogate of arterial stiffness and represents a non-invasive marker of cardiovascular risk. The non-invasive measurement of PWV requires tracking the arrival time of pressure pulses recorded in vivo, commonly referred to as pulse arrival time (PAT). In the state of the art, PAT is estimated by identifying a characteristic point of the pressure pulse waveform. This paper demonstrates that for ambulatory scenarios, where signal-to-noise ratios are below 10 dB, the performance in terms of repeatability of PAT measurements through characteristic points identification degrades drastically. Hence, we introduce a novel family of PAT estimators based on the parametric modeling of the anacrotic phase of a pressure pulse. In particular, we propose a parametric PAT estimator (TANH) that depicts high correlation with the Complior(R) characteristic point D1 (CC = 0.99), increases noise robustness and reduces by a five-fold factor the number of heartbeats required to obtain reliable PAT measurements.

Application of Signal Advance Technology to Electrophysiology

Medical instrumentation used in diagnosis and treatment relies on the accurate detection and processing of various physiological events and signals. While signal detection technology has improved greatly in recent years, there remain inherent delays in signal detection/ processing. These delays may have significant negative clinical consequences during various pathophysiological events. Reducing or eliminating such delays would increase the ability to provide successful early intervention in certain disorders thereby increasing the efficacy of treatment. In recent years, a physical phenomenon referred to as Negative Group Delay (NGD), demonstrated in simple electronic circuits, has been shown to temporally advance the detection of analog waveforms. Specifically, the output is temporally advanced relative to the input, as the time delay through the circuit is negative. The circuit output precedes the complete detection of the input signal. This process is referred to as signal advance (SA) detection. An SA circuit model incorporating NGD was designed, developed and tested. It imparts a constant temporal signal advance over a pre-specified spectral range in which the output is almost identical to the input signal (i.e., it has minimal distortion). Certain human patho-electrophysiological events are good candidates for the application of temporally-advanced waveform detection. SA technology has potential in early arrhythmia and epileptic seizure detection and intervention. Demonstrating reliable and consistent temporally advanced detection of electrophysiological waveforms may enable intervention with a pathological event (much) earlier than previously possible. SA detection could also be used to improve the performance of neural computer interfaces, neurotherapy applications, radiation therapy and imaging. In this study, the performance of a single-stage SA circuit model on a variety of constructed input signals, and human ECGs is investigated. The data obtained is used to quantify and characterize the temporal advances and circuit gain, as well as distortions in the output waveforms relative to their inputs. This project combines elements of physics, engineering, signal processing, statistics and electrophysiology. Its success has important consequences for the development of novel interventional methodologies in cardiology and neurophysiology as well as significant potential in a broader range of both biomedical and non-biomedical areas of application.

Reply to the editor

We appreciate the comments and concerns expressed by Arakawa and colleagues regarding our article, titled “Pulsatile control of rotary blood pumps: Does the modulation waveform matter?”1 Unfortunately, we have to disagree with Arakawa and colleagues. As is obvious from the title of our article, it investigates the effect of different waveforms on the heart–device interaction. In contrast to the authors' claim, this is the first article in the literature that uses basic waveforms (sine, triangle, saw tooth, and rectangular) with different phase shifts to examines their impact on left ventricular unloading. The previous publications2, 3 and 4 just varied the pump speed during systole and diastole, which was first reported by Bearnson and associates5 in 1996, and studied its effect on aortic pressure, coronary flow, and end-diastolic volume. We should mention that dp/dtmax is a load-sensitive parameter of contractility and not representative for the degree of unloading. Moreover, none of the aforementioned reports has studied mechanical unloading and in particular the stroke work of the left ventricle. Our method is unique because we do not just alternate between high and low speed but have accurate control of the waveform because of the direct drive system of Levitronix Technologies LLC (Waltham, Mass) and a custom-developed pump controller. Without referring, Arakawa and associates state “several previous studies have already reported the coronary flow diminishes as the left ventricular assist device support increases.” It should be noted that all the waveforms used in our study have 2000 rpm average value with 1000 rpm amplitude, which is not an excessive speed for the CentriMag rotary pump (Levitronix) to collapse the ventricle and diminish the coronary flow. We agree with Arakawa and coworkers that there is a need for a heart failure model to come to more relevant results with respect to clinical expectations. However, we have explored many existing models, including species and breeds that have a native proneness to cardiomyopathy, but all of them differ from the genetic presentation in humans. We certainly do not believe that the use of microembolization, in which the coronary circulation is impaired by the injection of microspheres, would form a good model from which to draw conclusions about coronary flow change under different loading conditions. A model would be needed in which either an infarct is created to mimic ischemic heart failure or the coronary circulation remains untouched to simulate, for instance, dilated cardiomyopathy. Furthermore, in discussion we clearly mention that “lack of heart failure is a major limitation of our study.” We also believe that unloading is not the only factor of the cardiac functional recovery, and an excessive unloading of the left ventricle might lead to cardiac tissue atrophy. Therefore, in our article we mention that control of the level of cardiac unloading by assist devices has been suggested as a mechanical tool to promote recovery, and more studies are required to find better strategies for the speed modulation of rotary pumps and to achieve an optimal heart load control to enhance myocardial recovery. Finally, there are many publications about pulsing rotary blood pumps and it was impossible to include them all. We preferred to reference some of the earlier basic works such as an original research by Bearnson and coworkers5 and another article published by our group,6 which is more relevant.

Motile response patterns and ultrastructural observations of the isolated outer hair cell

The tonotopic organization of the mammalian cochlea is accompanied by structural gradients which include the somatic lengths of outer hair cells (OHCs). These receptors rest upon the vibrating portion of the basilar membrane and have been reported to exhibit motile responses following chemical and electrical stimulation. These movements were examined in detail in this dissertation. It was found that isolated OHCs cultured in vitro respond to chemical depolarization with slow tonic movements, and to electrical waveforms with bi-directional, frequency following movements extending from DC to at least 10 kHz.^ Slow contractions were also elicited following electrical stimulation, bath incubation in carbachol (a cholinergic agonist), and increases in extracellular K+ concentration as little as 50 mM.^ Isolated OHCs display anatomical features which are remarkable when contrasted with those prepared from intact receptor organs. A complex structure located between the cuticular plate and the nuclear membrane was consistently observed and was examined by serial cross-sections which revealed a network of non-membrane bound densities. This corresponded to a granular complex seen at the light microscope level. The complex was composed of dense regions of organelles, striated structures embedded within the core, and a circumferential network of microtubules residing in the peri-nuclear portion of the cell. In cells which had lost their nuclear attachment to the terminal synaptic body, the granular complex could be made to contract without effecting any change in cellular length, implying that the complex may be the driving force behind certain aspects of the motile response.^ Most cells displayed movements which revealed asymmetries analogous to those reported for OHC receptor potentials in vivo. The contraction phase (for longer cells) was shown to have a small time constant (approximately 400 microseconds) and saturated with limited displacements. The expansion phase had time constants as large as 1.3 milliseconds but yielded displacements as much as 60 percent larger than those seen for contractions.^ Additional waveform characteristics seen in the in vivo response could be emulated either by biasing the cell's resting length with either direct current, triggering contractions via large electrical displacements, or incubation with depolarizing compounds.^ Alternatively, short (20-30 um) cells revealed more linear response characteristics to the probe stimulus. Partial saturation was achieved and revealed a DC component which was opposite in polarity to that seen in longer cells. (Abstract shortened with permission of author.) ^

Estimation of aortic pressure waveforms from 4D phase-contrast MRI

Several approaches for the non-invasive MRI-based measurement of the aortic pressure waveform over the heart cycle have been proposed in the last years. These methods are normally based on time-resolved, two-dimensional phase-contrast sequences with uni-directionally encoded velocities (2D PC-MRI). In contrast, three-dimensional acquisitions with tridirectional velocity encoding (4D PC-MRI) have been shown to be a suitable data source for detailed investigations of blood flow and spatial blood pressure maps. In order to avoid additional MR acquisitions, it would be advantageous if the aortic pressure waveform could also be computed from this particular form of MRI. Therefore, we propose an approach for the computation of the aortic pressure waveform which can be completely performed using 4D PC-MRI. After the application of a segmentation algorithm, the approach automatically computes the aortic pressure waveform without any manual steps. We show that our method agrees well with catheter measurements in an experimental phantom setup and produces physiologically realistic results in three healthy volunteers.

Antireflection coatings optimized for single-cycle THz pulses

We show that a single-layer antireflection coating on a THz source of high refractive index can substantially increase the transmission of emitted THz pulses. Calculations indicate that the optimum coating thickness depends on the exact shape of the generated THz waveform and whether the transmitted waveform is to be optimized for the highest peak (temporal) amplitude, peak spectral amplitude, or pulse energy. We experimentally demonstrate a 15% increase in peak amplitude, a 33% increase in peak spectral amplitude, and a 48% increase in energy for a 100 μm thick fused silica AR coating on a lithium niobate crystal used as THz emitter.

APPLICATION OF THE BRAINSTEM EVOKED RESPONSE (BER) TO NEUROTOXICITY EVALUATION OF THE AUDITORY SYSTEM: ASSESSMENT OF INORGANIC LEAD-INDUCED OTOTOXICITY

An experimental procedure was developed using the Brainstem Evoked Response (BER) electrophysiological technique to assess the effect of neurotoxic substances on the auditory system. The procedure utilizes Sprague-Dawley albino rats who have had dural electrodes implanted in their skulls, allowing neuroelectric evoked potentials to be recorded from their brainstems. Latency and amplitude parameters derived from the evoked potentials help assess the neuroanatomical integrity of the auditory pathway in the brainstem. Moreover, since frequency-specific auditory stimuli are used to evoke the neural responses, additional audiometric information is obtainable. An investigation on non-exposed control animals shows the BER threshold curve obtained by tests at various frequencies very closely approximates that obtained by behavioral audibility tests. Thus, the BER appears to be a valid measure of both functional and neuroanatomical integrity of the afferent auditory neural pathway.^ To determine the usefulness of the BER technique in neurobehavioral toxicology research, a known neurotoxic agent, Pb, was studied. Female Sprague-Dawley rats were dosed for 45 days with low levels of Pb acetate in their drinking water, after which BER recordings were obtained. The Pb dosages were determined from the findings of an earlier pilot study. One group of 6 rats received normal tap water, one group of 7 rats received a solution of 0.1% Pb, and another group of 7 rats received a solution of 0.2% Pb. After 45 days, the three groups exhibited blood Pb levels of 4.5 (+OR-) 0.43 (mu)g/100 ml, 37.8 (+OR-) 4.8 (mu)g/100 ml and 47.3 (+OR-) 2.7 (mu)g/100 ml, respectively.^ The results of the BER recording indicated evoked response waveform latency abnormalities in both the Pb-treated groups when midrange frequency (8 kHz to 32 kHz) stimuli were used. For the most part, waveform amplitudes did not vary significantly from control values. BER recordings obtained after a 30-day recovery period indicated the effects seen in the 0.1% Pb group had disappeared. However, those anomalies exhibited by the 0.2% Pb group either remained or increased in number. This outcome indicates a longer lasting or possibly irreversible effect on the auditory system from the higher dose of Pb. The auditory pathway effect appears to be in the periphery, at the level of the cochlea or the auditory (VIII) nerve. The results of this research indicate the BER technique is a valuable and sensitive indicator of low-level toxic effects on the auditory system.^

Digital surface model, hillshade and Lambertian reflectance model of the rock glaciers Oelgrube and Aeusseres Hochebenkar (Oetztal Alps, Tyrol, Austria)

With full-waveform (FWF) lidar systems becoming increasingly available from different commercial manufacturers, the possibility for extracting physical parameters of the scanned surfaces in an area-wide sense, as addendum to their geometric representation, has risen as well. The mentioned FWF systems digitize the temporal profiles of the transmitted laser pulse and of its backscattered echoes, allowing for a reliable determination of the target distance to the instrument and of physical target quantities by means of radiometric calibration, one of such quantities being the diffuse Lambertian reflectance. The delineation of glaciers is a time-consuming task, commonly performed manually by experts and involving field trips as well as image interpretation of orthophotos, digital terrain models and shaded reliefs. In this study, the diffuse Lambertian reflectance was compared to the glacier outlines mapped by experts. We start the presentation with the workflow for analysis of FWF data, their direct georeferencing and the calculation of the diffuse Lambertian reflectance by radiometric calibration; this workflow is illustrated for a large FWF lidar campaign in the Ötztal Alps (Tyrol, Austria), operated with an Optech ALTM 3100 system. The geometric performance of the presented procedure was evaluated by means of a relative and an absolute accuracy assessment using strip differences and orthophotos, resp. The diffuse Lambertian reflectance was evaluated at two rock glaciers within the mentioned lidar campaign. This feature showed good performance for the delineation of the rock glacier boundaries, especially at their lower parts.

Table 5 Comparison of surface elevations by Seasat and optical survey along EGIG line

Surface elevation maps of the southern half of the Greenland subcontinent are produced from radar altimeter data acquired by the Seasat satellite. A summary of the processing procedure and examples of return waveform data are given. The elevation data are used to generate a regular grid which is then computer contoured to provide an elevation contour map. Ancillary maps show the statistical quality of the elevation data and various characteristics of the surface. The elevation map is used to define ice flow directions and delineate the major drainage basins. Regular maps of the Jakobshavns Glacier drainage basin and the ice divide in the vicinity of Crete Station are presented. Altimeter derived elevations are compared with elevations measured both by satellite geoceivers and optical surveying.

Physical properties and sedimentology on core GeoB1510-1

Detailed information about the sediment properties and microstructure can be provided through the analysis of digital ultrasonic P wave seismograms recorded automatically during full waveform core logging. The physical parameter which predominantly affects the elastic wave propagation in water-saturated sediments is the P wave attenuation coefficient. The related sedimentological parameter is the grain size distribution. A set of high-resolution ultrasonic transmission seismograms (ca. 50-500 kHz), which indicate downcore variations in the grain size by their signal shape and frequency content, are presented. Layers of coarse-grained foraminiferal ooze can be identified by highly attenuated P waves, whereas almost unattenuated waves are recorded in fine-grained areas of nannofossil ooze. Color-encoded pixel graphics of the seismograms and instantaneous frequencies present full waveform images of the lithology and attenuation. A modified spectral difference method is introduced to determine the attenuation coefficient and its power law a = kfn. Applied to synthetic seismograms derived using a "constant Q" model, even low attenuation coefficients can be quantified. A downcore analysis gives an attenuation log which ranges from ca. 700 dB/m at 400 kHz and a power of n = 1-2 in coarse-grained sands to few decibels per meter and n ? 0.5 in fine-grained clays. A least squares fit of a second degree polynomial describes the mutual relationship between the mean grain size and the attenuation coefficient. When it is used to predict the mean grain size, an almost perfect coincidence with the values derived from sedimentological measurements is achieved.