Electrochemical Impedance Spectroscopy-A Tutorial

This tutorial provides the theoretical background, the principles, and applications of Electrochemical Impedance Spectroscopy (EIS) in various research and technological sectors. The text has been organized in 17 sections starting with basic knowledge on sinusoidal signals, complex numbers, phasor notation, and transfer functions, continuing with the definition of impedance in electrical circuits, the principles of EIS, the validation of the experimental data, their simulation to equivalent electrical circuits, and ending with practical considerations and selected examples on the utility of EIS to corrosion, energy related applications, and biosensing. A user interactive excel file showing the Nyquist and Bode plots of some model circuits is provided in the Supporting Information. This tutorial aspires to provide the essential background to graduate students working on EIS, as well as to endow the knowledge of senior researchers on various fields where EIS is involved. We also believe that the content of this tutorial will be a useful educational tool for EIS instructors.

Incorporating local adaptation into forecasts of species' distribution and abundance under climate change

Populations of many species are genetically adapted to local historical climate conditions. Yet most forecasts of species' distributions under climate change have ignored local adaptation (LA), which may paint a false picture of how species will respond across their geographic ranges. We review recent studies that have incorporated intraspecific variation, a potential proxy for LA, into distribution forecasts, assess their strengths and weaknesses, and make recommendations for how to improve forecasts in the face of LA. The three methods used so far (species distribution models, response functions, and mechanistic models) reflect a trade-off between data availability and the ability to rigorously demonstrate LA to climate. We identify key considerations for incorporating LA into distribution forecasts that are currently missing from many published studies, including testing the spatial scale and pattern of LA, the confounding effects of LA to nonclimatic or biotic drivers, and the need to incorporate empirically based dispersal or gene flow processes. We suggest approaches to better evaluate these aspects of LA and their effects on species-level forecasts. In particular, we highlight demographic and dynamic evolutionary models as promising approaches to better integrate LA into forecasts, and emphasize the importance of independent model validation. Finally, we urge closer examination of how LA will alter the responses of central vs. marginal populations to allow stronger generalizations about changes in distribution and abundance in the face of LA.

Dynamics of networks of excitatory and inhibitory neurons in response to time-dependent inputs

We investigate the dynamics of recurrent networks of excitatory (E) and inhibitory (I) neurons in the presence of time-dependent inputs. The dynamics is characterized by the network dynamical transfer function, i.e., how the population firing rate is modulated by sinusoidal inputs at arbitrary frequencies. Two types of networks are studied and compared: (i) a Wilson-Cowan type firing rate model; and (ii) a fully connected network of leaky integrate-and-fire (LIF) neurons, in a strong noise regime. We first characterize the region of stability of the "asynchronous state" (a state in which population activity is constant in time when external inputs are constant) in the space of parameters characterizing the connectivity of the network. We then systematically characterize the qualitative behaviors of the dynamical transfer function, as a function of the connectivity. We find that the transfer function can be either low-pass, or with a single or double resonance, depending on the connection strengths and synaptic time constants. Resonances appear when the system is close to Hopf bifurcations, that can be induced by two separate mechanisms: the I-I connectivity and the E-I connectivity. Double resonances can appear when excitatory delays are larger than inhibitory delays, due to the fact that two distinct instabilities exist with a finite gap between the corresponding frequencies. In networks of LIF neurons, changes in external inputs and external noise are shown to be able to change qualitatively the network transfer function. Firing rate models are shown to exhibit the same diversity of transfer functions as the LIF network, provided delays are present. They can also exhibit input-dependent changes of the transfer function, provided a suitable static non-linearity is incorporated.

A novel, non-invasive, online-monitoring, versatile and easy plant-based probe for measuring leaf water status

A high-precision pressure probe is described which allows non-invasive online-monitoring of the water relations of intact leaves. Real-time recording of the leaf water status occurred by data transfer to an Internet server. The leaf patch clamp pressure probe measures the attenuated pressure, P(p), of a leaf patch in response to a constant clamp pressure, P(clamp). P(p) is sensed by a miniaturized silicone pressure sensor integrated into the device. The magnitude of P(p) is dictated by the transfer function of the leaf, T(f), which is a function of leaf patch volume and ultimately of cell turgor pressure, P(c), as shown theoretically. The power function T(f)=f(P(c)) theoretically derived was experimentally confirmed by concomitant P(p) and P(c) measurements on intact leaflets of the liana Tetrastigma voinierianum under greenhouse conditions. Simultaneous P(p) recordings on leaflets up to 10 m height above ground demonstrated that changes in T(f) induced by P(c) changes due to changes of microclimate and/or of the irrigation regime were sensitively reflected in corresponding changes of P(p). Analysis of the data show that transpirational water loss during the morning hours was associated with a transient rise in turgor pressure gradients within the leaflets. Subsequent recovery of turgescence during the afternoon was much faster than the preceding transpiration-induced water loss if the plants were well irrigated. Our data show the enormous potential of the leaf patch clamp pressure probe for leaf water studies including unravelling of the hydraulic communication between neighbouring leaves and over long distances within tall plants (trees).

An animal-to-human scaling law for blast-induced traumatic brain injury risk assessment

Despite recent efforts to understand blast effects on the human brain, there are still no widely accepted injury criteria for humans. Recent animal studies have resulted in important advances in the understanding of brain injury due to intense dynamic loads. However, the applicability of animal brain injury results to humans remains uncertain. Here, we use advanced computational models to derive a scaling law relating blast wave intensity to the mechanical response of brain tissue across species. Detailed simulations of blast effects on the brain are conducted for different mammals using image-based biofidelic models. The intensity of the stress waves computed for different external blast conditions is compared across species. It is found that mass scaling, which successfully estimates blast tolerance of the thorax, fails to capture the brain mechanical response to blast across mammals. Instead, we show that an appropriate scaling variable must account for the mass of protective tissues relative to the brain, as well as their acoustic impedance. Peak stresses transmitted to the brain tissue by the blast are then shown to be a power function of the scaling parameter for a range of blast conditions relevant to TBI. In particular, it is found that human brain vulnerability to blast is higher than for any other mammalian species, which is in distinct contrast to previously proposed scaling laws based on body or brain mass. An application of the scaling law to recent experiments on rabbits furnishes the first physics-based injury estimate for blast-induced TBI in humans.

Predicting direction detection thresholds for arbitrary translational acceleration profiles in the horizontal plane

In previous research, direction detection thresholds have been measured and successfully modeled by exposing participants to sinusoidal acceleration profiles of different durations. In this paper, we present measurements that reveal differences in thresholds depending not only on the duration of the profile, but also on the actual time course of the acceleration. The measurements are further explained by a model based on a transfer function, which is able to predict direction detection thresholds for all types of acceleration profiles. In order to quantify a participant's ability to detect the direction of motion in the horizontal plane, a four-alternative forced-choice task was implemented. Three types of acceleration profiles (sinusoidal, trapezoidal and triangular) were tested for three different durations (1.5, 2.36 and 5.86 s). To the best of our knowledge, this is the first study which varies both quantities (profile and duration) in a systematic way within a single experiment. The lowest thresholds were found for trapezoidal profiles and the highest for triangular profiles. Simulations for frequencies lower than the ones actually measured predict a change from this behavior: Sinusoidal profiles are predicted to yield the highest thresholds at low frequencies. This qualitative prediction is only possible with a model that is able to predict thresholds for different types of acceleration profiles. Our modeling approach represents an important advancement, because it allows for a more general and accurate description of perceptual thresholds for simple and complex translational motions.

Open-loop dynamic and static characteristics of the carotid sinus baroreflex in rats with chronic heart failure after myocardial infarction

We estimated open-loop dynamic characteristics of the carotid sinus baroreflex in normal control rats and chronic heart failure (CHF) rats after myocardial infarction. First, the neural arc transfer function from carotid sinus pressure to splanchnic sympathetic nerve activity (SNA) and its corresponding step response were examined. Although the steady-state response was attenuated in CHF, the negative peak response and the time to peak did not change significantly, suggesting preserved neural arc dynamic characteristics. Next, the peripheral arc transfer function from SNA to arterial pressure (AP) and its corresponding step response were examined. The steady-state response and the initial slope were reduced in CHF, suggesting impaired end-organ responses. In a simulation study based on the dynamic and static characteristics, the percent recovery of AP was reduced progressively as the size of disturbance increased in CHF, suggesting that a reserve for AP buffering is lost in CHF despite relatively maintained baseline AP.

Changes in vascular properties, not ventricular properties, predominantly contribute to baroreflex regulation of arterial pressure

Baroreflex modulates both the ventricular and vascular properties and stabilizes arterial pressure (AP). However, how changes in those mechanical properties quantitatively impact the dynamic AP regulation remains unknown. We developed a framework of circulatory equilibrium, in which both venous return and cardiac output are expressed as functions of left ventricular (LV) end-systolic elastance (Ees), heart rate (HR), systemic vascular resistance (R), and stressed blood volume (V). We investigated the contribution of each mechanical property using the framework of circulatory equilibrium. In six anesthetized dogs, we vascularly isolated carotid sinuses and randomly changed carotid sinus pressure (CSP), while measuring the LV Ees, aortic flow, right and left atrial pressure, and AP for at least 60 min. We estimated transfer functions from CSP to Ees, HR, R, and V in each dog. We then predicted these parameters in response to changes in CSP from the transfer functions using a data set not used for identifying transfer functions and predicted changes in AP using the equilibrium framework. Predicted APs matched reasonably well with those measured (r2=0.85-0.96, P<0.001). Sensitivity analyses indicated that Ees and HR (ventricular properties) accounted for 14±4 and 4±2%, respectively, whereas R and V (vascular properties) accounted for 32±4 and 39±4%, respectively, of baroreflex-induced AP regulation. We concluded that baroreflex-induced dynamic AP changes can be accurately predicted by the transfer functions from CSP to mechanical properties using our framework of circulatory equilibrium. Changes in the vascular properties, not the ventricular properties, predominantly determine baroreflex-induced AP regulation.

Compression on trigger points in the leg muscle increases parasympathetic nervous activity based on heart rate variability

Massotherapy, the therapeutic use of massage, is used to treat various chronic pain syndromes. One type of massotherapy, pressure stimulus applied over trigger points (TPs), is reported to have excellent therapeutic effects. Its effect is possibly mediated through changes in the autonomic nervous system although little research has been conducted to assess autonomic activity during TP compression. We have investigated how compression applied over TPs affects the autonomic nervous system. Six healthy young adult females whose daily working routine was carried out predominantly in a standing position were enrolled in the study cohort. After a day's work, the subjects were asked to rest supine, and electrocardiograms (ECGs), instantaneous lung volume (ILV) and systolic and diastolic blood pressures (SBP, DBP) were measured before and after pressure application over the TPs in those lower limb muscles where the subjects felt muscle fatigue or discomfort. The subjects were also asked to coordinate breathing with the beeping sounds. The therapeutic effects of TP compression were assessed by a subjective fatigue scale. Parasympathetic nervous activity was also assessed by spectral analysis of heart rate (HR) variability. The transfer function from ILV to HR was evaluated using linear analysis. The results indicated that TP compression (1) decreased HR, SBP and DBP, (2) increased parasympathetic activity, (3) increased the gain from ILV to HR, and (4) improved the fatigue scores. These findings suggest that an increase in parasympathetic nervous activity after the TP compression induced a reduction of fatigue. The therapeutic mechanisms of TP compression to enhance parasympathetic nervous system are discussed.

Oxygen Level and LFP in Task-Positive and Task-Negative Areas: Bridging BOLD fMRI and Electrophysiology

The human default mode network (DMN) shows decreased blood oxygen level dependent (BOLD) signals in response to a wide range of attention-demanding tasks. Our understanding of the specifics regarding the neural activity underlying these "task-negative" BOLD responses remains incomplete. We paired oxygen polarography, an electrode-based oxygen measurement technique, with standard electrophysiological recording to assess the relationship of oxygen and neural activity in task-negative posterior cingulate cortex (PCC), a hub of the DMN, and visually responsive task-positive area V3 in the awake macaque. In response to engaging visual stimulation, oxygen, LFP power, and multi-unit activity in PCC showed transient activation followed by sustained suppression. In V3, oxygen, LFP power, and multi-unit activity showed an initial phasic response to the stimulus followed by sustained activation. Oxygen responses were correlated with LFP power in both areas, although the apparent hemodynamic coupling between oxygen level and electrophysiology differed across areas. Our results suggest that oxygen responses reflect changes in LFP power and multi-unit activity and that either the coupling of neural activity to blood flow and metabolism differs between PCC and V3 or computing a linear transformation from a single LFP band to oxygen level does not capture the true physiological process.

Assessment of central blood pressure in patients with type 2 diabetes: a comparison between SphygmoCor and invasively measured values

BACKGROUND

The SphygmoCor is used for noninvasive assessment of ascending aortic blood pressure (BP). However, the validity of the SphygmoCor transfer function has not been tested in an exclusively type 2 diabetic patient sample. Calibration with systolic (SBP) and diastolic (DBP) brachial BP has previously been associated with substantial imprecision of central BP estimates. We hypothesized that different noninvasive calibration strategies might improve the accuracy of the estimated ascending aortic BPs.

METHODS

In 34 patients with type 2 diabetes we estimated ascending aortic SBP and DBP using the SphygmoCor device and compared these data with invasively recorded data. The validity of the transfer function was assessed by calibrating with invasively recorded DBP and mean BP (MBP). The influence of noninvasive calibration strategies was assessed by calibrating with brachial oscillometric SBP+DBP vs. DBP+MBP using a form factor (ff) of 0.33 and 0.40, respectively.

RESULTS

When calibrating with invasive BP, the difference between estimated and invasively measured ascending aortic SBP and DBP was -2.3±5.6/1.0±0.9 mm Hg. When calibrating with oscillometric brachial BPs, the differences were -9.6±8.1/14.1±6.2 mm Hg (calibration with SBP and DBP), -8.3±11.7/13.9±6.1 mm Hg (DBP and MBP; ff = 0.33), and 1.9±12.2/14.1±6.2 mm Hg (DBP and MBP; ff = 0.40), respectively. Calibration with the average of 3 brachial BPs did not improve accuracy.

CONCLUSIONS

The SphygmoCor transfer function seems valid in patients with type 2 diabetes. Noninvasive calibration with DBP and MBP (ff = 0.40) enables accurate estimation of mean ascending aortic SBP at the group level. However, the wide limits of agreement indicate limited accuracy in the individual patient.

CLINICAL TRIALS REGISTRATION

Clinical Trials No. NCT01538290.

Lighting up yeast cell factories by transcription factor-based biosensors

Our ability to rewire cellular metabolism for the sustainable production of chemicals, fuels and therapeutics based on microbial cell factories has advanced rapidly during the last two decades. Especially the speed and precision by which microbial genomes can be engineered now allow for more advanced designs to be implemented and tested. However, compared to the methods developed for engineering cell factories, the methods developed for testing the performance of newly engineered cell factories in high throughput are lagging far behind, which consequently impacts the overall biomanufacturing process. For this purpose, there is a need to develop new techniques for screening and selection of best-performing cell factory designs in multiplex. Here we review the current status of the sourcing, design and engineering of biosensors derived from allosterically regulated transcription factors applied to the biotechnology work-horse budding yeast Saccharomyces cerevisiae. We conclude by providing a perspective on the most important challenges and opportunities lying ahead in order to harness the full potential of biosensor development for increasing both the throughput of cell factory development and robustness of overall bioprocesses.

MRI-based transfer function determination for the assessment of implant safety

PURPOSE

We introduce a new MR-based method to determine the transfer function (TF) for radiofrequency (RF) safety assessment of active implantable medical devices. Transfer functions are implant-specific measures that relate the incident tangential electric field on an (elongated) implant to a scattered electric field at its tip. The proposed method allows for TF determination with a high spatial resolution in relatively fast measurements without requiring dedicated bench setups from MRI images.

THEORY AND METHODS

The principle of reciprocity is used in conjunction with the potential to measure currents with MRI to determine TF. Low-flip angle 3D dual gradient echo MRI data are acquired with an implant as transceive antenna, which requires minimal hardware adaptations. The implant-specific TF is determined from the acquired MRI data, with two different postprocessing methods for comparison.

RESULTS

TFs of linear and helical implants can be determined accurately (with a Pearson correlation coefficient R ≥ 0.7 between measurements and simulations, and a difference in field at the tip ΔEtip ≤ 19%) from relatively quick (t < 20 minutes) MRI acquisitions with (several) millimeter spatial resolution.

CONCLUSION

Transfer function determination with MRI for RF safety assessment of implantable medical devices is possible. The proposed MR-based method allows for TF determination in more realistic exposure scenarios and solid media. Magn Reson Med 78:2449-2459, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Tube-load model parameter estimation for monitoring arterial hemodynamics

A useful model of the arterial system is the uniform, lossless tube with parametric load. This tube-load model is able to account for wave propagation and reflection (unlike lumped-parameter models such as the Windkessel) while being defined by only a few parameters (unlike comprehensive distributed-parameter models). As a result, the parameters may be readily estimated by accurate fitting of the model to available arterial pressure and flow waveforms so as to permit improved monitoring of arterial hemodynamics. In this paper, we review tube-load model parameter estimation techniques that have appeared in the literature for monitoring wave reflection, large artery compliance, pulse transit time, and central aortic pressure. We begin by motivating the use of the tube-load model for parameter estimation. We then describe the tube-load model, its assumptions and validity, and approaches for estimating its parameters. We next summarize the various techniques and their experimental results while highlighting their advantages over conventional techniques. We conclude the review by suggesting future research directions and describing potential applications.

Effects of autoregulation and CO2 reactivity on cerebral oxygen transport

Both autoregulation and CO(2) reactivity are known to have significant effects on cerebral blood flow and thus on the transport of oxygen through the vasculature. In this paper, a previous model of the autoregulation of blood flow in the cerebral vasculature is expanded to include the dynamic behavior of oxygen transport through binding with hemoglobin. The model is used to predict the transfer functions for both oxyhemoglobin and deoxyhemoglobin in response to fluctuations in arterial blood pressure and arterial CO(2) concentration. It is shown that only six additional nondimensional groups are required in addition to the five that were previously found to characterize the cerebral blood flow response. A resonant frequency in the pressure-oxyhemoglobin transfer function is found to occur in the region of 0.1 Hz, which is a frequency of considerable physiological interest. The model predictions are compared with results from the published literature of phase angle at this frequency, showing that the effects of changes in breathing rate can significantly alter the inferred phase dynamics between blood pressure and hemoglobin. The question of whether dynamic cerebral autoregulation is affected under conditions of stenosis or stroke is then examined.

Experimental investigation of rotational tomography in reconstructed middle ears with clinical implications

A large air-bone-gap after ossiculoplasty may be due to a malpositioned or displaced prosthesis. Rotational tomography (RT) has the potential to provide high-resolution images of implants without artifacts and with less radiation dosage than CT scan. Twenty-seven temporal bone specimens underwent measurements of middle ear transfer function using Laser-Doppler-Vibrometry (LDV) before and after placement of ossicular replacement prostheses (PORPs, TORPs) made of titanium. RT was performed on all specimens. RT allowed 3-dimensional viewing of the temporal bone, accurate localization of implants within the reconstructed middle ear and determination of angles between the inserted prostheses and the tympanic membrane (TM) and/or the malleus handle (MH). Presence or absence of contact between the implant and the TM, malleus or stapes could be clearly visualized. Displaced prostheses were readily identified. The functional LDV-measurements for TORPs showed a trend favoring coupling to the malleus handle, while for PORPs, coupling to the TM was favored. For PORPs, sound transmission was worse with increasing angles between the PORP and stapes superstructure (p<0.05). Following our experimental results RT is an innovative, relevant and useful imaging technique to obtain immediate postoperative feedback after ossicular reconstruction and to precisely determine the position of middle ear implants.

Rapid acquisition of the 3D MRI gradient impulse response function using a simple phantom measurement

PURPOSE

To provide a simple tool for rapid measurement of the 3D gradient modulation transfer function (GMTF) of clinical MRI systems using a phantom. Knowledge of the transfer function is useful for gradient chain characterization, system calibration, and improvement of image reconstruction results.

METHODS

Starting from the well-established thin slice method used for phantom-based measurement of the 1D GMTF, we add phase encoding to partition the thin slices into voxels that act as localized field probes. From the signal phase evolution measured at the 3D voxel positions, the GMTF can be derived for cross and higher order spatial terms represented by spherical harmonics up to 3rd order.

RESULTS

Using spherical phantoms, 16 GMTFs representing all terms up to 3rd order harmonics can be determined in a scan time of <2 min. A large voxel volume of >1 mL yields high SNR, enabling signal acquisition using the system's body coil. The method is applied for improving system calibration and for characterizing the effect of additional hardware in the bore.

CONCLUSION

The presented method seems well-suited for rapid measurement of the GMTF of a clinical system, as it delivers high-quality results in a short scan time.

Measurement of 100 nm and 60 nm Particle Standards by Differential Mobility Analysis

The peak particle size and expanded uncertainties (95 % confidence interval) for two new particle calibration standards are measured as 101.8 nm ± 1.1 nm and 60.39 nm ± 0.63 nm. The particle samples are polystyrene spheres suspended in filtered, deionized water at a mass fraction of about 0.5 %. The size distribution measurements of aerosolized particles are made using a differential mobility analyzer (DMA) system calibrated using SRM(®) 1963 (100.7 nm polystyrene spheres). An electrospray aerosol generator was used for generating the 60 nm aerosol to almost eliminate the generation of multiply charged dimers and trimers and to minimize the effect of non-volatile contaminants increasing the particle size. The testing for the homogeneity of the samples and for the presence of multimers using dynamic light scattering is described. The use of the transfer function integral in the calibration of the DMA is shown to reduce the uncertainty in the measurement of the peak particle size compared to the approach based on the peak in the concentration vs. voltage distribution. A modified aerosol/sheath inlet, recirculating sheath flow, a high ratio of sheath flow to the aerosol flow, and accurate pressure, temperature, and voltage measurements have increased the resolution and accuracy of the measurements. A significant consideration in the uncertainty analysis was the correlation between the slip correction of the calibration particle and the measured particle. Including the correlation reduced the expanded uncertainty from approximately 1.8 % of the particle size to about 1.0 %. The effect of non-volatile contaminants in the polystyrene suspensions on the peak particle size and the uncertainty in the size is determined. The full size distributions for both the 60 nm and 100 nm spheres are tabulated and selected mean sizes including the number mean diameter and the dynamic light scattering mean diameter are computed. The use of these particles for calibrating DMAs and for making deposition standards to be used with surface scanning inspection systems is discussed.

A Refresher of the Original Bloch's Law Paper (Bloch, July 1885)

In 1885, Adolphe-Moïse Bloch asked the following simple question “Is there a law describing the relationship between the duration of a light and its perceived intensity?” Based on a series of experiments using a Foucault regulator and a candle, Bloch concluded that “when the lighting duration varies from 0.00173 to 0.0518 seconds (…) the [visible] light is markedly in inverse proportion to its duration”—his famous law. As this law pertains to the more general and hotly debated question of accumulation of sensory information over time, it is timely to offer the public a full translation of Bloch’s original paper (from French) and to present it within the context of contemporary research.

MRI-based, wireless determination of the transfer function of a linear implant: Introduction of the transfer matrix

PURPOSE

We introduce the transfer matrix (TM) that makes MR-based wireless determination of transfer functions (TFs) possible. TFs are implant specific measures for RF-safety assessment of linear implants. The TF relates an incident tangential electric field on an implant to a scattered electric field at its tip that generally governs local heating. The TM extends this concept and relates an incident tangential electric field to a current distribution in the implant therewith characterizing the RF response along the entire implant. The TM is exploited to measure TFs with MRI without hardware alterations.

THEORY AND METHODS

A model of rightward and leftward propagating attenuated waves undergoing multiple reflections is used to derive an analytical expression for the TM. This allows parameterization of the TM of generic implants, e.g., (partially) insulated single wires, in a homogeneous medium in a few unknowns that simultaneously describe the TF. These unknowns can be determined with MRI making it possible to measure the TM and, therefore, also the TF.

RESULTS

The TM is able to predict an induced current due to an incident electric field and can be accurately parameterized with a limited number of unknowns. Using this description the TF is determined accurately (with a Pearson correlation coefficient R ≥ 0.9 between measurements and simulations) from MRI acquisitions.

CONCLUSION

The TM enables measuring of TFs with MRI of the tested generic implant models. The MR-based method does not need hardware alterations and is wireless hence making TF determination in more realistic scenarios conceivable.

Aortic pulse wave velocity and reflecting distance estimation from peripheral waveforms in humans: detection of age- and exercise training-related differences

We hypothesized that demographic/anthropometric parameters can be used to estimate effective reflecting distance (EfRD), required to derive aortic pulse wave velocity (APWV), a prognostic marker of cardiovascular risk, from peripheral waveforms and that such estimates can discriminate differences in APWV and EfRD with aging and habitual endurance exercise in healthy adults. Ascending aortic pressure waveforms were derived from peripheral waveforms (brachial artery pressure, n = 25; and finger volume pulse, n = 15) via a transfer function and then used to determine the time delay between forward- and backward-traveling waves (Δtf-b). True EfRDs were computed as directly measured carotid-femoral pulse wave velocity (CFPWV) × 1/2Δtf-b and then used in regression analysis to establish an equation for EfRD based on demographic/anthropometric data (EfRD = 0.173·age + 0.661·BMI + 34.548 cm, where BMI is body mass index). We found good agreement between true and estimated APWV (Pearson's R² = 0.43; intraclass correlation = 0.64; both P < 0.05) and EfRD (R² = 0.24; intraclass correlation = 0.40; both P < 0.05). In young sedentary (22 ± 2 years, n = 6), older sedentary (62 ± 1 years, n = 24), and older endurance-trained (61 ± 2 years, n = 14) subjects, EfRD (from demographic/anthropometric parameters), APWV, and 1/2Δtf-b (from brachial artery pressure waveforms) were 52.0 ± 0.5, 61.8 ± 0.4, and 60.6 ± 0.5 cm; 6.4 ± 0.3, 9.6 ± 0.2, and 8.1 ± 0.2 m/s; and 82 ± 3, 65 ± 1 and 76 ± 2 ms (all P < 0.05), respectively. Our results demonstrate that APWV derived from peripheral waveforms using age and BMI to estimate EfRD correlates with CFPWV in healthy adults. This method can reliably detect the distal shift of the reflecting site with age and the increase in APWV with sedentary aging that is attenuated with habitual endurance exercise.

Tabla: A Proof-of-Concept Auscultatory Percussion Device for Low-Cost Pneumonia Detection

Pneumonia causes the deaths of over a million people worldwide each year, with most occurring in countries with limited access to expensive but effective diagnostic methods, e.g., chest X-rays. Physical examination, the other major established method of diagnosis, suffers from several drawbacks, most notably low accuracy and high interobserver error. We sought to address this diagnostic gap by developing a proof-of-concept non-invasive device to identify the accumulation of fluid in the lungs (consolidation) characteristic of pneumonia. This device, named Tabla after the percussive instrument of the same name, utilizes the technique of auscultatory percussion; a percussive input sound is sent through the chest and recorded with a digital stethoscope for analysis. Tabla analyzes differences in sound transmission through the chest at audible frequencies as a marker for lung consolidation. This paper presents preliminary data from five pneumonia patients and eight healthy subjects. We demonstrate 92.3% accuracy in distinguishing between healthy subjects and patients with pneumonia after data analysis with a K-nearest neighbors algorithm. This prototype device is low cost and simple to implement and may offer a rapid and inexpensive method for pneumonia diagnosis appropriate for general use and in areas with limited medical infrastructure.

Improved pulse wave velocity estimation using an arterial tube-load model

Pulse wave velocity (PWV) is the most important index of arterial stiffness. It is conventionally estimated by noninvasively measuring central and peripheral blood pressure (BP) and/or velocity (BV) waveforms and then detecting the foot-to-foot time delay between the waveforms wherein wave reflection is presumed absent. We developed techniques for improved estimation of PWV from the same waveforms. The techniques effectively estimate PWV from the entire waveforms, rather than just their feet, by mathematically eliminating the reflected wave via an arterial tube-load model. In this way, the techniques may be more robust to artifact while revealing the true PWV in absence of wave reflection. We applied the techniques to estimate aortic PWV from simultaneously and sequentially measured central and peripheral BP waveforms and simultaneously measured central BV and peripheral BP waveforms from 17 anesthetized animals during diverse interventions that perturbed BP widely. Since BP is the major acute determinant of aortic PWV, especially under anesthesia wherein vasomotor tone changes are minimal, we evaluated the techniques in terms of the ability of their PWV estimates to track the acute BP changes in each subject. Overall, the PWV estimates of the techniques tracked the BP changes better than those of the conventional technique (e.g., diastolic BP root-mean-squared errors of 3.4 versus 5.2 mmHg for the simultaneous BP waveforms and 7.0 versus 12.2 mmHg for the BV and BP waveforms (p <; 0.02)). With further testing, the arterial tube-load model-based PWV estimation techniques may afford more accurate arterial stiffness monitoring in hypertensive and other patients.

Testate amoeba records indicate regional 20th-century lowering of water tables in ombrotrophic peatlands in central-northern Alberta, Canada

Testate amoebae are abundant in the surface layers of northern peatlands. Analysis of their fossilized shell (test) assemblages allows for reconstructions of local water-table depths (WTD). We have reconstructed WTD dynamics for five peat cores from peatlands ranging in distance from the Athabasca bituminous sands (ABS) region in western Canada. Amoeba assemblages were combined with plant macrofossil records, acid-insoluble ash (AIA) fluxes and instrumental climate data to identify drivers for environmental change. Two functional traits of testate amoebae, mixotrophy and the tendency to integrate xenogenic mineral matter in test construction, were quantified to infer possible effects of AIA flux on testate amoeba presence. Age-depth models showed the cores each covered at least the last ~315 years, with some spanning the last millennium. Testate amoeba assemblages were likely affected by permafrost development in two of the peatlands, yet the most important shift in assemblages was detected after 1960 CE. This shift represents a significant apparent lowering of water tables in four out of five cores, with a mean drop of ~15 cm. Over the last 50 years, assemblages shifted towards more xerophilous taxa, a trend which was best explained by increasing Sphagnum s. Acutifolia and, to a lesser extent, mean summer temperature. This trend was most evident in the two cores from the sites located farthest away from the ABS region. AIA flux variations did not show a clear effect on mineral-agglutinating taxa, nor on S. s. Acutifolia presence. We therefore suggest the drying trend was forced by the establishment of S. s. Acutifolia, driven by enhanced productivity following regional warming. Such recent apparent drying of peatlands, which may only be reconstructed by appropriate indicators combined with high chronological control, may affect vulnerability to future burning and promote emissions of CO₂ .

Evaluation of an independent linear model for acoustic noise on a conventional MRI scanner and implications for acoustic noise reduction

PURPOSE

To evaluate an independent linear model for gradient acoustic noise on a conventional MRI scanner, and to explore implications for acoustic noise reduction in routine imaging.

METHODS

Acoustic noise generated from each physical gradient axis was modeled as the prescribed gradient waveform passed through a linear time-invariant system. Homogeneity and superposition properties were experimentally determined. We also developed a new method to correct relative time shifts between the measured impulse responses for different physical gradient axes. Model accuracy was determined by comparing predicted and measured sound using normalized energy difference. Transfer functions were also measured in subjects with different body habitus and at multiple microphone locations.

RESULTS

Both superposition and homogeneity held for each physical gradient axis with errors less than 3%. When all gradients were on simultaneous sound prediction, error was reduced from 32% to 4% after time-shift correction. Transfer functions also showed high sensitivity to body habitus and microphone location.

CONCLUSION

The independent linear model predicts MRI acoustic noise with less than 4% error. Acoustic transfer functions are highly sensitive to body habitus and position within the bore, making it challenging to produce a general approach to acoustic noise reduction based on avoiding system resonance peaks.