Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Activation of electrical triggers of atrial fibrillation in hyperthyroidism.

CONTEXT: A shortening of the atrial refractory period has been considered as the main mechanism for the increased risk of atrial fibrillation in hyperthyroidism. However, other important factors may be involved. OBJECTIVE: Our objective was to determine the activity of abnormal supraventricular electrical depolarizations in response to elevated thyroid hormones in patients without structural heart disease. PATIENTS AND DESIGN: Twenty-eight patients (25 females, three males, mean age 43+/-11 yr) with newly diagnosed and untreated hyperthyroidism were enrolled in a prospective trial after exclusion of heart disease. Patients were followed up for 16 +/- 6 months and studied at baseline and 6 months after normalization of serum TSH levels. MAIN OUTCOME MEASURES: The incidence of abnormal premature supraventricular depolarizations (SVPD) and the number of episodes of supraventricular tachycardia was defined as primary outcome measurements before the start of the study. In addition, heart rate oscillations (turbulence) after premature depolarizations and heart rate variability were compared at baseline and follow-up. RESULTS: SVPDs decreased from 59 +/- 29 to 21 +/- 8 per 24 h (P = 0.003), very early SVPDs (so called P on T) decreased from 36 +/- 24 to 3 +/- 1 per 24 h (P < 0.0001), respectively, and nonsustained supraventricular tachycardias decreased from 22 +/- 11 to 0.5 +/- 0.2 per 24 h (P = 0.01) after normalization of serum thyrotropin levels. The hyperthyroid phase was characterized by an increased heart rate (93 +/- 14 vs. 79 +/- 8 beats/min, P < 0.0001) and a decreased turbulence slope (3.6 vs. 9.2, P = 0.003), consistent with decreased vagal tone. This was confirmed by a significant decrease of heart rate variability. CONCLUSION: Hyperthyroidism is associated with an increased supraventricular ectopic activity in patients with normal hearts. The activation of these arrhythmogenic foci by elevated thyroid hormones may be an important causal link between hyperthyroidism and atrial fibrillation.

The thermal circuits of the Argentera Massif (western Alps, Italy): An example of low-enthalpy geothermal resources controlled by Neogene alpine tectonics

This work presents geochemistry and structural geology data concerning the low enthalpy geothermal circuits of the Argentera crystalline Massif in northwestern Italian Alps. I n this area some thermal springs (50-60 degreesC), located in the small Bagni di Vinadio village, discharge mixtures made up of a Na-Cl end-member and a Na-SO4 component. The latter is also discharged by the thermal springs of Terme di Valdieri located some kilometres apart within the same tectonic complex. Both end-members share the same meteoric origin and the same reservoir temperature, which is close to 150 degreesC. Explanations are thus required to understand how they reach the surface and how waters of the same origin and circulating in similar rocks can attain such different compositions. Sodium-sulphate waters discharged at both sites, likely represent the common interaction product of meteoric waters with the widespread granitic-migmatitic rocks of the Argentera Massif, whereas Na-CI waters originate through leaching of mineralised cataclastic rocks, which are rich in phyllosilicatic minerals and fluid inclusions, both acting as Cl- sources. Due to the relatively low inferred geothermal gradient of the region, -25C/km, meteoric waters have to descend to depths of 5.5-6 km to attain temperatures of similar to 150 degreesC. These relevant depths can be reached by descending meteoric waters, due to the recent extensional stress field, which allows the development of geothermal circulations at greater depths than in other sectors of the Alps by favouring a greater fractures aperture. The ascent of the thermal waters rakes place along brittle shear zones. In both sites, the thermal waters emerge at the bottoms of the valleys, close to either the lateral termination of a brittle shear zone at Terme di Valdieri, or a step-over between two en-echelon brittle shear zones at Bagni di Vinadio. These observations attest to a strong control operated on the location of outlet regions by both brittle tectonics and the minima in hydraulic potential inside the fractured massif.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

Continuous assessment of electrical epileptic activity in acute stroke.

OBJECTIVE: To determine the incidence and risk factors of electrical seizures and other electrical epileptic activity using continuous EEG (cEEG) in patients with acute stroke. METHODS: One hundred consecutive patients with acute stroke admitted to our stroke unit underwent cEEG using 10 electrodes. In addition to electrical seizures, repetitive focal sharp waves (RSHWs), repetitive focal spikes (RSPs), and periodic lateralized epileptic discharges (PLEDs) were recorded. RESULTS: In the 100 patients, cEEG was recorded for a mean duration of 17 hours 34 minutes (range 1 hour 12 minutes to 37 hours 10 minutes). Epileptic activity occurred in 17 patients and consisted of RSHWs in seven, RSPs in seven, and PLEDs in three. Electrical seizures occurred in two patients. On univariate Cox regression analysis, predictors for electrical epileptic activity were stroke severity (high score on the National Institutes of Health Stroke Scale) (hazard ratio [HR] 1.12; p = 0.002), cortical involvement (HR 5.71; p = 0.021), and thrombolysis (HR 3.27; p = 0.040). Age, sex, stroke type, use of EEG-modifying medication, and cardiovascular risk factors were not predictors of electrical epileptic activity. On multivariate analysis, stroke severity was the only independent predictor (HR 1.09; p = 0.016). CONCLUSION: In patients with acute stroke, electrical epileptic activity occurs more frequently than previously suspected.

Use of anisotropic modelling in electrical impedance tomography: description of method and preliminary assessment of utility in imaging brain function in the adult human head.

Electrical Impedance Tomography (EIT) is an imaging method which enables a volume conductivity map of a subject to be produced from multiple impedance measurements. It has the potential to become a portable non-invasive imaging technique of particular use in imaging brain function. Accurate numerical forward models may be used to improve image reconstruction but, until now, have employed an assumption of isotropic tissue conductivity. This may be expected to introduce inaccuracy, as body tissues, especially those such as white matter and the skull in head imaging, are highly anisotropic. The purpose of this study was, for the first time, to develop a method for incorporating anisotropy in a forward numerical model for EIT of the head and assess the resulting improvement in image quality in the case of linear reconstruction of one example of the human head. A realistic Finite Element Model (FEM) of an adult human head with segments for the scalp, skull, CSF, and brain was produced from a structural MRI. Anisotropy of the brain was estimated from a diffusion tensor-MRI of the same subject and anisotropy of the skull was approximated from the structural information. A method for incorporation of anisotropy in the forward model and its use in image reconstruction was produced. The improvement in reconstructed image quality was assessed in computer simulation by producing forward data, and then linear reconstruction using a sensitivity matrix approach. The mean boundary data difference between anisotropic and isotropic forward models for a reference conductivity was 50%. Use of the correct anisotropic FEM in image reconstruction, as opposed to an isotropic one, corrected an error of 24 mm in imaging a 10% conductivity decrease located in the hippocampus, improved localisation for conductivity changes deep in the brain and due to epilepsy by 4-17 mm, and, overall, led to a substantial improvement on image quality. This suggests that incorporation of anisotropy in numerical models used for image reconstruction is likely to improve EIT image quality.

Electrical source imaging for presurgical focus localization in epilepsy patients with normal MRI.

PURPOSE: Patients with magnetic resonance (MR)-negative focal epilepsy (MRN-E) have less favorable surgical outcomes (between 40% and 70%) compared to those in whom an MRI lesion guides the site of surgical intervention (60-90%). Patients with extratemporal MRN-E have the worst outcome (around 50% chance of seizure freedom). We studied whether electroencephalography (EEG) source imaging (ESI) of interictal epileptic activity can contribute to the identification of the epileptic focus in patients with normal MRI. METHODS: We carried out ESI in 10 operated patients with nonlesional MRI and a postsurgical follow-up of at least 1 year. Five of the 10 patients had extratemporal lobe epilepsy. Evaluation comprised surface and intracranial EEG monitoring of ictal and interictal events, structural MRI, [(18)F]fluorodeoxyglucose positron emission tomography (FDG-PET), ictal and interictal perfusion single photon emission computed tomography (SPECT) scans. Eight of the 10 patients also underwent intracranial monitoring. RESULTS: ESI correctly localized the epileptic focus within the resection margins in 8 of 10 patients, 9 of whom experienced favorable postsurgical outcomes. DISCUSSION: The results highlight the diagnostic value of ESI and encourage broadening its application to patients with MRN-E. If the surface EEG contains fairly localized spikes, ESI contributes to the presurgical decision process.

Spatial distribution of motor units recruited during electrical stimulation of the quadriceps muscle versus the femoral nerve.

INTRODUCTION: In this study we investigated differences in the spatial recruitment of motor units (MUs) in the quadriceps when electrical stimulation is applied over the quadriceps belly versus the femoral nerve. METHODS: M-waves and mechanical twitches were evoked using over-the-quadriceps and femoral nerve stimulation of gradually increasing intensity from 22 young, healthy subjects. Spatial recruitment was investigated using recruitment curves of M-waves recorded from the vastus medialis (VM) and vastus lateralis (VL) and of twitches recorded from the quadriceps. RESULTS: At maximal stimulation intensity (Imax), no differences were found between nerve and over-the-quadriceps stimulation. At submaximal intensities, VL M-wave amplitude was higher for over-the-quadriceps stimulation at 40% Imax, and peak twitch force was greater for nerve stimulation at 60% and 80% Imax. CONCLUSIONS: For the VM, MU spatial recruitment during nerve and over-the-quadriceps stimulation of increasing intensity occurred in a similar manner, whereas significant differences were observed for the VL. Muscle Nerve, 2013.

Emergence of oriented circuits driven by synaptic pruning associated with spike-timing-dependent plasticity (STDP)

ABSTRACT: Massive synaptic pruning following over-growth is a general feature of mammalian brain maturation. Pruning starts near time of birth and is completed by time of sexual maturation. Trigger signals able to induce synaptic pruning could be related to dynamic functions that depend on the timing of action potentials. Spike-timing-dependent synaptic plasticity (STDP) is a change in the synaptic strength based on the ordering of pre- and postsynaptic spikes. The relation between synaptic efficacy and synaptic pruning suggests that the weak synapses may be modified and removed through competitive "learning" rules. This plasticity rule might produce the strengthening of the connections among neurons that belong to cell assemblies characterized by recurrent patterns of firing. Conversely, the connections that are not recurrently activated might decrease in efficiency and eventually be eliminated. The main goal of our study is to determine whether or not, and under which conditions, such cell assemblies may emerge out of a locally connected random network of integrate-and-fire units distributed on a 2D lattice receiving background noise and content-related input organized in both temporal and spatial dimensions. The originality of our study stands on the relatively large size of the network, 10,000 units, the duration of the experiment, 10E6 time units (one time unit corresponding to the duration of a spike), and the application of an original bio-inspired STDP modification rule compatible with hardware implementation. A first batch of experiments was performed to test that the randomly generated connectivity and the STDP-driven pruning did not show any spurious bias in absence of stimulation. Among other things, a scale factor was approximated to compensate for the network size on the ac¬tivity. Networks were then stimulated with the spatiotemporal patterns. The analysis of the connections remaining at the end of the simulations, as well as the analysis of the time series resulting from the interconnected units activity, suggest that feed-forward circuits emerge from the initially randomly connected networks by pruning. RESUME: L'élagage massif des synapses après une croissance excessive est une phase normale de la ma¬turation du cerveau des mammifères. L'élagage commence peu avant la naissance et est complété avant l'âge de la maturité sexuelle. Les facteurs déclenchants capables d'induire l'élagage des synapses pourraient être liés à des processus dynamiques qui dépendent de la temporalité rela¬tive des potentiels d'actions. La plasticité synaptique à modulation temporelle relative (STDP) correspond à un changement de la force synaptique basé sur l'ordre des décharges pré- et post- synaptiques. La relation entre l'efficacité synaptique et l'élagage des synapses suggère que les synapses les plus faibles pourraient être modifiées et retirées au moyen d'une règle "d'appren¬tissage" faisant intervenir une compétition. Cette règle de plasticité pourrait produire le ren¬forcement des connexions parmi les neurones qui appartiennent à une assemblée de cellules caractérisée par des motifs de décharge récurrents. A l'inverse, les connexions qui ne sont pas activées de façon récurrente pourraient voir leur efficacité diminuée et être finalement éliminées. Le but principal de notre travail est de déterminer s'il serait possible, et dans quelles conditions, que de telles assemblées de cellules émergent d'un réseau d'unités integrate-and¬-fire connectées aléatoirement et distribuées à la surface d'une grille bidimensionnelle recevant à la fois du bruit et des entrées organisées dans les dimensions temporelle et spatiale. L'originalité de notre étude tient dans la taille relativement grande du réseau, 10'000 unités, dans la durée des simulations, 1 million d'unités de temps (une unité de temps correspondant à une milliseconde), et dans l'utilisation d'une règle STDP originale compatible avec une implémentation matérielle. Une première série d'expériences a été effectuée pour tester que la connectivité produite aléatoirement et que l'élagage dirigé par STDP ne produisaient pas de biais en absence de stimu¬lation extérieure. Entre autres choses, un facteur d'échelle a pu être approximé pour compenser l'effet de la variation de la taille du réseau sur son activité. Les réseaux ont ensuite été stimulés avec des motifs spatiotemporels. L'analyse des connexions se maintenant à la fin des simulations, ainsi que l'analyse des séries temporelles résultantes de l'activité des neurones, suggèrent que des circuits feed-forward émergent par l'élagage des réseaux initialement connectés au hasard.

Single acquisition electrical property mapping based on relative coil sensitivities: A proof-of-concept demonstration.

PURPOSE: All methods presented to date to map both conductivity and permittivity rely on multiple acquisitions to compute quantitatively the magnitude of radiofrequency transmit fields, B1+. In this work, we propose a method to compute both conductivity and permittivity based solely on relative receive coil sensitivities ( B1-) that can be obtained in one single measurement without the need to neither explicitly perform transmit/receive phase separation nor make assumptions regarding those phases. THEORY AND METHODS: To demonstrate the validity and the noise sensitivity of our method we used electromagnetic finite differences simulations of a 16-channel transceiver array. To experimentally validate our methodology at 7 Tesla, multi compartment phantom data was acquired using a standard 32-channel receive coil system and two-dimensional (2D) and 3D gradient echo acquisition. The reconstructed electric properties were correlated to those measured using dielectric probes. RESULTS: The method was demonstrated both in simulations and in phantom data with correlations to both the modeled and bench measurements being close to identity. The noise properties were modeled and understood. CONCLUSION: The proposed methodology allows to quantitatively determine the electrical properties of a sample using any MR contrast, with the only constraint being the need to have 4 or more receive coils and high SNR. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Comparison of electrical nerve stimulation, electrical muscle stimulation and magnetic nerve stimulation to assess the neuromuscular function of the plantar flexor muscles.

INTRODUCTION: As it might lead to less discomfort, magnetic nerve stimulation (MNS) is increasingly used as an alternative to electrical stimulation methods. Yet, MNS and electrical nerve stimulation (ENS) and electrical muscle stimulation (EMS) have not been formally compared for the evaluation of plantar flexor neuromuscular function. METHODS: We quantified plantar flexor neuromuscular function with ENS, EMS and MNS in 10 volunteers in fresh and fatigued muscles. Central alterations were assessed through changes in voluntary activation level (VAL) and peripheral function through changes in M-wave, twitch and doublet (PS100) amplitudes. Discomfort associated with 100-Hz paired stimuli delivered with each method was evaluated on a 10-cm visual analog scale. RESULTS: VAL, agonist and antagonist M-wave amplitudes and PS100 were similar between the different methods in both fresh and fatigued states. Potentiated peak twitch was lower in EMS compared to ENS, whereas no difference was found between ENS and MNS for any parameter. Discomfort associated with MNS (1.5 ± 1.4 cm) was significantly less compared to ENS (5.5 ± 1.9 cm) and EMS (4.2 ± 2.6 cm) (p < 0.05). CONCLUSION: When PS100 is used to evaluate neuromuscular properties, MNS, EMS and ENS can be used interchangeably for plantar flexor neuromuscular function assessment as they provide similar evaluation of central and peripheral factors in unfatigued and fatigued states. Importantly, electrical current spread to antagonist muscles was similar between the three methods while discomfort from MNS was much less compared to ENS and EMS. MNS may be potentially employed to assess neuromuscular function of plantar flexor muscles in fragile populations.

Extra Forces induced by wide-pulse, high-frequency electrical stimulation: Occurrence, magnitude, variability and underlying mechanisms.

OBJECTIVE: In contrast to conventional (CONV) neuromuscular electrical stimulation (NMES), the use of "wide-pulse, high-frequencies" (WPHF) can generate higher forces than expected by the direct activation of motor axons alone. We aimed at investigating the occurrence, magnitude, variability and underlying neuromuscular mechanisms of these "Extra Forces" (EF). METHODS: Electrically-evoked isometric plantar flexion force was recorded in 42 healthy subjects. Additionally, twitch potentiation, H-reflex and M-wave responses were assessed in 13 participants. CONV (25Hz, 0.05ms) and WPHF (100Hz, 1ms) NMES consisted of five stimulation trains (20s on-90s off). RESULTS: K-means clustering analysis disclosed a responder rate of almost 60%. Within this group of responders, force significantly increased from 4% to 16% of the maximal voluntary contraction force and H-reflexes were depressed after WPHF NMES. In contrast, non-responders showed neither EF nor H-reflex depression. Twitch potentiation and resting EMG data were similar between groups. Interestingly, a large inter- and intrasubject variability of EF was observed. CONCLUSION: The responder percentage was overestimated in previous studies. SIGNIFICANCE: This study proposes a novel methodological framework for unraveling the neurophysiological mechanisms involved in EF and provides further evidence for a central contribution to EF in responders.