Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain

The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups

Comparison of electrohysterogram signal measured by surface electrodes with different designs: A computational study with dipole band and abdomen models

Non-invasive measurement of uterine activity using electrohysterogram (EHG) surface electrodes has been attempted to monitor uterine contraction. This study aimed to computationally compare the performance of acquiring EHG signals using monopolar electrode and three types of Laplacian concentric ring electrodes (bipolar, quasi-bipolar and tri-polar). With the implementation of dipole band model and abdomen model, the performances of four electrodes in terms of the local sensitivity were quantifed by potential attenuation. Furthermore, the efects of fat and muscle thickness on potential attenuation were evaluated using the bipolar and tri-polar electrodes with diferent radius. The results showed that all the four types of electrodes detected the simulated EHG signals with consistency. That the bipolar and tri-polar electrodes had greater attenuations than the others, and the shorter distance between the origin and location of dipole band at 20dB attenuation, indicating that they had relatively better local sensitivity. In addition, ANOVA analysis showed that, for all the electrodes with diferent outer ring radius, the efects of fat and muscle on potential attenuation were signifcant (all p<0.01). It is therefore concluded that the bipolar and tri-polar electrodes had higher local sensitivity than the others, indicating that they can be applied to detect EHG efectively

Sensitivity of MEG and EEG to Source Orientation

An important difference between magnetoencephalography (MEG) and electroencephalography (EEG) is that MEG is insensitive to radially oriented sources. We quantified computationally the dependency of MEG and EEG on the source orientation using a forward model with realistic tissue boundaries. Similar to the simpler case of a spherical head model, in which MEG cannot see radial sources at all, for most cortical locations there was a source orientation to which MEG was insensitive. The median value for the ratio of the signal magnitude for the source orientation of the lowest and the highest sensitivity was 0.06 for MEG and 0.63 for EEG. The difference in the sensitivity to the source orientation is expected to contribute to systematic differences in the signal-to-noise ratio between MEG and EEG.National Institutes of Health (U.S.) (Grant NS057500)National Institutes of Health (U.S.) (Grant NS037462)National Institutes of Health (U.S.) (Grant HD040712)National Center for Research Resources (U.S.) (P41RR14075)Mind Research Networ

Cardiac anisotropy in boundary-element models for the electrocardiogram

The boundary-element method (BEM) is widely used for electrocardiogram (ECG) simulation. Its major disadvantage is its perceived inability to deal with the anisotropic electric conductivity of the myocardial interstitium, which led researchers to represent only intracellular anisotropy or neglect anisotropy altogether. We computed ECGs with a BEM model based on dipole sources that accounted for a “compound” anisotropy ratio. The ECGs were compared with those computed by a finite-difference model, in which intracellular and interstitial anisotropy could be represented without compromise. For a given set of conductivities, we always found a compound anisotropy value that led to acceptable differences between BEM and finite-difference results. In contrast, a fully isotropic model produced unacceptably large differences. A model that accounted only for intracellular anisotropy showed intermediate performance. We conclude that using a compound anisotropy ratio allows BEM-based ECG models to more accurately represent both anisotropies

Spatiotemporally resolved imaging of streamer discharges in air generated in a wire-cylinder reactor with (sub)nanosecond voltage pulses

We use (sub)nanosecond high-voltage pulses to generate streamers in atmospheric-pressure air in a wire-cylinder reactor. We study the effect of reactor length, pulse duration, pulse amplitude, pulse polarity, and pulse rise time on the streamer development, specifically on the streamer distribution in the reactor to relate it to plasma-processing results. We use ICCD imaging with a fully automated setup that can image the streamers in the entire corona-plasma reactor. From the images, we calculate streamer lengths and velocities. We also develop a circuit simulation model of the reactor to support the analysis of the streamer development. The results show how the propagation of the high-voltage pulse through the reactor determines the streamer development. As the pulse travels through the reactor, it generates streamers and attenuates and disperses. At the end of the reactor, it reflects and adds to itself. The local voltage on the wire together with the voltage rise time determine the streamer velocities, and the pulse duration the consequent maximal streamer length

Facilitating arrhythmia simulation: the method of quantitative cellular automata modeling and parallel running

BACKGROUND: Many arrhythmias are triggered by abnormal electrical activity at the ionic channel and cell level, and then evolve spatio-temporally within the heart. To understand arrhythmias better and to diagnose them more precisely by their ECG waveforms, a whole-heart model is required to explore the association between the massively parallel activities at the channel/cell level and the integrative electrophysiological phenomena at organ level. METHODS: We have developed a method to build large-scale electrophysiological models by using extended cellular automata, and to run such models on a cluster of shared memory machines. We describe here the method, including the extension of a language-based cellular automaton to implement quantitative computing, the building of a whole-heart model with Visible Human Project data, the parallelization of the model on a cluster of shared memory computers with OpenMP and MPI hybrid programming, and a simulation algorithm that links cellular activity with the ECG. RESULTS: We demonstrate that electrical activities at channel, cell, and organ levels can be traced and captured conveniently in our extended cellular automaton system. Examples of some ECG waveforms simulated with a 2-D slice are given to support the ECG simulation algorithm. A performance evaluation of the 3-D model on a four-node cluster is also given. CONCLUSIONS: Quantitative multicellular modeling with extended cellular automata is a highly efficient and widely applicable method to weave experimental data at different levels into computational models. This process can be used to investigate complex and collective biological activities that can be described neither by their governing differentiation equations nor by discrete parallel computation. Transparent cluster computing is a convenient and effective method to make time-consuming simulation feasible. Arrhythmias, as a typical case, can be effectively simulated with the methods described

Intraoperative electrocorticography using high-frequency oscillations or spikes to tailor epilepsy surgery in the Netherlands (the HFO trial): a randomised, single-blind, adaptive non-inferiority trial

Background Intraoperative electrocorticography is used to tailor epilepsy surgery by analysing interictal spikes or spike patterns that can delineate epileptogenic tissue. High-frequency oscillations (HFOs) on intraoperative electrocorticography have been proposed as a new biomarker of epileptogenic tissue, with higher specificity than spikes. We prospectively tested the non-inferiority of HFO-guided tailoring of epilepsy surgery to spike-guided tailoring on seizure freedom at 1 year.Methods The HFO trial was a randomised, single-blind, adaptive non-inferiority trial at an epilepsy surgery centre (UMC Utrecht) in the Netherlands. We recruited children and adults (no age limits) who had been referred for intraoperative electrocorticography-tailored epilepsy surgery. Participants were randomly allocated (1:1) to either HFO-guided or spike-guided tailoring, using an online randomisation scheme with permuted blocks generated by an independent data manager, stratified by epilepsy type. Treatment allocation was masked to participants and clinicians who documented seizure outcome, but not to the study team or neurosurgeon. Ictiform spike patterns were always considered in surgical decision making. The primary endpoint was seizure outcome after 1 year (dichotomised as seizure freedom [defined as Engel 1A-11 vs seizure recurrence [Engel 1C-4]). We predefined a non-inferiority margin of 10% risk difference. Analysis was by intention to treat, with prespecified subgroup analyses by epilepsy type and for confounders. This completed trial is registered with the Dutch Trial Register, Toetsingonline ABR.NL44527.041.13, and ClinicalTrials.gov, NCT02207673.Findings Between Oct 10, 2014, and Jan 31,2020,78 individuals were enrolled to the study and randomly assigned (39 to HFO-guided tailoring and 39 to spike-guided tailoring). There was no loss to follow-up. Seizure freedom at 1 year occurred in 26 (67%) of 39 participants in the HFO-guided group and 35 (90%) of 39 in the spike-guided group (risk difference -23.5%, 90% CI -39.1 to -7.9; for the 48 patients with temporal lobe epilepsy, the risk difference was -25.5%, -45.1 to -6.0, and for the 30 patients with extratemporal lobe epilepsy it was -20.3%, -46.0 to 5.4). Pathology associated with poor prognosis was identified as a confounding factor, with an adjusted risk difference of-7.9% (90% CI -20.7 to 4.9; adjusted risk difference -12.5%, -31.0 to 5.9, for temporal lobe epilepsy and 5.8%, -7.7 to 19.5, for extratemporal lobe epilepsy). We recorded eight serious adverse events (five in the HFO-guided group and three in the spike-guided group) requiring hospitalisation. No patients died.Interpretation HFO-guided tailoring of epilepsy surgery was not non-inferior to spike-guided tailoring on intraoperative electrocorticography. After adjustment for confounders, HFOs show non-inferiority in extratemporal lobe epilepsy. This trial challenges the clinical value of HFOs as an epilepsy biomarker, especially in temporal lobe epilepsy. Further research is needed to establish whether HFO-guided intraoperative electrocorticography holds promise in extratemporal lobe epilepsy. Copyright (C) 2022 The Author(s). Published by Elsevier Ltd

Basisonderzoek Aardgasverbruik Kleinverbruikers 1998 : daling huishoudelijk gasverbruik zet zich voort

In June 1998 the 20th title survey has been carried out. Every year a representative random sample survey is held amongst family households in the Netherlands to gain insight into the natural gas consumption of the households and the factors which influence this consumption.