Autonomous navigation of stratospheric balloons using reinforcement learning.

Efficiently navigating a superpressure balloon in the stratosphere1 requires the integration of a multitude of cues, such as wind speed and solar elevation, and the process is complicated by forecast errors and sparse wind measurements. Coupled with the need to make decisions in real time, these factors rule out the use of conventional control techniques2,3. Here we describe the use of reinforcement learning4,5 to create a high-performing flight controller. Our algorithm uses data augmentation6,7 and a self-correcting design to overcome the key technical challenge of reinforcement learning from imperfect data, which has proved to be a major obstacle to its application to physical systems8. We deployed our controller to station Loon superpressure balloons at multiple locations across the globe, including a 39-day controlled experiment over the Pacific Ocean. Analyses show that the controller outperforms Loon's previous algorithm and is robust to the natural diversity in stratospheric winds. These results demonstrate that reinforcement learning is an effective solution to real-world autonomous control problems in which neither conventional methods nor human intervention suffice, offering clues about what may be needed to create artificially intelligent agents that continuously interact with real, dynamic environments.

Molecular evidence of widespread benzimidazole drug resistance in Ancylostoma caninum from domestic dogs throughout the USA and discovery of a novel ß-tubulin benzimidazole resistance mutation.

Ancylostoma caninum is an important zoonotic gastrointestinal nematode of dogs worldwide and a close relative of human hookworms. We recently reported that racing greyhound dogs in the USA are infected with A. caninum that are commonly resistant to multiple anthelmintics. Benzimidazole resistance in A. caninum in greyhounds was associated with a high frequency of the canonical F167Y(TTC>TAC) isotype-1 ß-tubulin mutation. In this work, we show that benzimidazole resistance is remarkably widespread in A. caninum from domestic dogs across the USA. First, we identified and showed the functional significance of a novel benzimidazole isotype-1 ß-tubulin resistance mutation, Q134H(CAA>CAT). Several benzimidazole resistant A. caninum isolates from greyhounds with a low frequency of the F167Y(TTC>TAC) mutation had a high frequency of a Q134H(CAA>CAT) mutation not previously reported from any eukaryotic pathogen in the field. Structural modeling predicted that the Q134 residue is directly involved in benzimidazole drug binding and that the 134H substitution would significantly reduce binding affinity. Introduction of the Q134H substitution into the C. elegans ß-tubulin gene ben-1, by CRISPR-Cas9 editing, conferred similar levels of resistance as a ben-1 null allele. Deep amplicon sequencing on A. caninum eggs from 685 hookworm positive pet dog fecal samples revealed that both mutations were widespread across the USA, with prevalences of 49.7% (overall mean frequency 54.0%) and 31.1% (overall mean frequency 16.4%) for F167Y(TTC>TAC) and Q134H(CAA>CAT), respectively. Canonical codon 198 and 200 benzimidazole resistance mutations were absent. The F167Y(TTC>TAC) mutation had a significantly higher prevalence and frequency in Western USA than in other regions, which we hypothesize is due to differences in refugia. This work has important implications for companion animal parasite control and the potential emergence of drug resistance in human hookworms.

Restoring the Fluctuation-Dissipation Theorem in Kardar-Parisi-Zhang Universality Class through a New Emergent Fractal Dimension.

The Kardar-Parisi-Zhang (KPZ) equation describes a wide range of growth-like phenomena, with applications in physics, chemistry and biology. There are three central questions in the study of KPZ growth: the determination of height probability distributions; the search for ever more precise universal growth exponents; and the apparent absence of a fluctuation-dissipation theorem (FDT) for spatial dimension d>1. Notably, these questions were answered exactly only for 1+1 dimensions. In this work, we propose a new FDT valid for the KPZ problem in d+1 dimensions. This is achieved by rearranging terms and identifying a new correlated noise which we argue to be characterized by a fractal dimension dn. We present relations between the KPZ exponents and two emergent fractal dimensions, namely df, of the rough interface, and dn. Also, we simulate KPZ growth to obtain values for transient versions of the roughness exponent α, the surface fractal dimension df and, through our relations, the noise fractal dimension dn. Our results indicate that KPZ may have at least two fractal dimensions and that, within this proposal, an FDT is restored. Finally, we provide new insights into the old question about the upper critical dimension of the KPZ universality class.

Mixtures of self-propelled particles interacting with asymmetric obstacles.

In the presence of an obstacle, active particles condensate into a surface "wetting" layer due to persistent motion. If the obstacle is asymmetric, a rectification current arises in addition to wetting. Asymmetric geometries are therefore commonly used to concentrate microorganisms like bacteria and sperms. However, most studies neglect the fact that biological active matter is diverse, composed of individuals with distinct self-propulsions. Using simulations, we study a mixture of "fast" and "slow" active Brownian disks in two dimensions interacting with large half-disk obstacles. With this prototypical obstacle geometry, we analyze how the stationary collective behavior depends on the degree of self-propulsion "diversity," defined as proportional to the difference between the self-propulsion speeds, while keeping the average self-propulsion speed fixed. A wetting layer rich in fast particles arises. The rectification current is amplified by speed diversity due to a superlinear dependence of rectification on self-propulsion speed, which arises from cooperative effects. Thus, the total rectification current cannot be obtained from an effective one-component active fluid with the same average self-propulsion speed, highlighting the importance of considering diversity in active matter.

TREXIO: A file format and library for quantum chemistry.

TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an overview of the TREXIO file format and library. The library consists of a front-end implemented in the C programming language and two different back-ends: a text back-end and a binary back-end utilizing the hierarchical data format version 5 library, which enables fast read and write operations. It is compatible with a variety of platforms and has interfaces for Fortran, Python, and OCaml programming languages. In addition, a suite of tools have been developed to facilitate the use of the TREXIO format and library, including converters for popular quantum chemistry codes and utilities for validating and manipulating data stored in TREXIO files. The simplicity, versatility, and ease of use of TREXIO make it a valuable resource for researchers working with quantum chemistry data.

Oxidative potential of metal-polluted urban dust as a potential environmental stressor for chronic diseases.

Oxidative stress (OS) associated with metals in urban dust has become a public health concern. Chronic diseases linked to general inflammation are particularly affected by OS. This research analyzes the spatial distribution of metals associated with OS, the urban dust´s oxidative potential (OP), and the occurrence of diseases whose treatments are affected by OS. We collected 70 urban dust samples during pre- and post-monsoon seasons to achieve this. We analyzed particle size distribution and morphology by scanning electron microscopy, as well as metal(loid)s by portable X-ray fluorescence, and OP of dust in artificial lysosomal fluid by using an ascorbic acid depletion assay. Our results show that the mean concentration of Fe, Pb, As, Cr, Cu, and V in pre-monsoon was 83,984.6, 98.4, 23.5, 165.8, 301.3, and 141.9 mg kg-1, while during post-monsoon was 50,638.8, 73.9, 16.7, 124.3, 178.9, and 133.5 mg kg-1, respectively. Impoverished areas with the highest presence of cardiovascular, cancer, diabetes, and respiratory diseases coincide with contaminated areas where young adults live. We identified significant differences in the OP between seasons. OP increases during the pre-monsoon (from 7.8 to 237.5 nmol AA min-1) compared to the post-monsoon season (from 1.6 to 163.2 nmol AA min-1). OP values are much higher than measured standards corresponding to contaminated soil and urban particulate matter, which means that additional sources beside metals cause the elevated OP. The results show no risk from chronic exposure to metals; however, our results highlight the importance of studying dust as an environmental factor that may potentially increase oxidative stress.

Ethical issues of organ donation after circulatory death: Considerations for a successful implementation in Chile.

Organ transplantation is a lifesaving procedure for end-organ damage and remains up to today as the most cost-effective alternative to treat these conditions. However, the main limitation to performing organ transplants is the availability of donor organs suitable for transplantation. To increase the donor pool, expanding organ donation from the conventional neurologic determination of death (NDD) to include circulatory determination of death (DCD) has been a well-established method of increasing donors in other countries. In this article, we discuss the clinical and ethical considerations for introducing DCD in Chile. The concepts we have used could very well be translatable to other similar countries which have not implemented this donation system yet. The most relevant issue to date is that DCD needs to alter the care of dying patients to obtain quality donor organs. In some countries, including Chile, there are some cultural barriers regarding withdrawal-of-care. These barriers include confusing withdrawal of care with acceleration of death, which leads to many practitioners refusing to remove artificial life support, and in turn only minimize ventilatory support or switch to a T-tube (without extubation). This cultural barrier could be overcome with careful consideration of the opinions of healthcare workers, family members, community and policy-based stakeholders. We also identified ethical issues related to informed consent of both donor and recipients, among other relevant ethical considerations. In conclusion, DCD donation in Chile can increase organ donation numbers in one of Latin America's countries with the lowest effective donor rate. However, this opportunity must be taken with caution to avoid the opposite effect if this policy is not well implemented, respecting the sound ethical principles mentioned in this paper.

Bidirectional astrocytic GLUT1 activation by elevated extracellular K.

The acute rise in interstitial K+ that accompanies neural activity couples the energy demand of neurons to the metabolism of astrocytes. The effects of elevated K+ on astrocytes include activation of aerobic glycolysis, inhibition of mitochondrial respiration and the release of lactate. Using a genetically encoded FRET glucose sensor and a novel protocol based on 3-O-methylglucose trans-acceleration and numerical simulation of glucose dynamics, we report that extracellular K+ is also a potent and reversible modulator of the astrocytic glucose transporter GLUT1. In cultured mouse astrocytes, the stimulatory effect developed within seconds, engaged both the influx and efflux modes of the transporter, and was detected even at 1 mM incremental K+ . The modulation of GLUT1 explains how astrocytes are able to maintain their glucose pool in the face of strong glycolysis stimulation. We propose that the stimulation of GLUT1 by K+ supports the production of lactate by astrocytes and the timely delivery of glucose to active neurons.

Improved short-term outcomes of kidney transplants in controlled donation after the circulatory determination of death with the use of normothermic regional perfusion.

Normothermic regional perfusion (NRP) allows the in situ perfusion of organs with oxygenated blood in donation after the circulatory determination of death (DCDD). We aimed at evaluating the impact of NRP on the short-term outcomes of kidney transplants in controlled DCDD (cDCDD). This is a multicenter, nationwide, retrospective study comparing cDCDD kidneys obtained with NRP versus the standard rapid recovery (RR) technique. During 2012-2018, 2302 cDCDD adult kidney transplants were performed in Spain using NRP (n = 865) or RR (n = 1437). The study groups differed in donor and recipient age, warm, and cold ischemic time and use of ex situ machine perfusion. Transplants in the NRP group were more frequently performed in high-volume centers (≥90 transplants/year). Through matching by propensity score, two cohorts with a total of 770 patients were obtained. After the matching, no statistically significant differences were observed between the groups in terms of primary nonfunction (p = .261) and mortality at 1 year (p =  .111). However, the RR of kidneys was associated with a significantly increased odds of delayed graft function (OR 1.97 [95% CI 1.43-2.72]; p < .001) and 1-year graft loss (OR 1.77 [95% CI 1.01-3.17]; p = .034). In conclusion, compared with RR, NRP appears to improve the short-term outcomes of cDCDD kidney transplants.

Active mixtures in a narrow channel: motility diversity changes cluster sizes.

The persistent motion of bacteria produces clusters with a stationary cluster size distribution (CSD). Here we develop a minimal model for bacteria in a narrow channel to assess the relative importance of motility diversity (i.e. polydispersity in motility parameters) and confinement. A mixture of run-and-tumble particles with a distribution of tumbling rates (denoted generically by α) is considered on a 1D lattice. Particles facing each other cross at constant rate, rendering the lattice quasi-1D. To isolate the role of diversity, the global average α stays fixed. For a binary mixture with no particle crossing, the average cluster size (Lc) increases with the diversity as lower-α particles trap higher-α ones for longer. At finite crossing rate, particles escape from the clusters sooner, making Lc smaller and the diversity less important, even though crossing can enhance demixing of particle types between the cluster and gas phases. If the crossing rate is increased further, the clusters become controlled by particle crossing. We also consider an experiment-based continuous distribution of tumbling rates, revealing similar physics. Using parameters fitted from experiments with Escherichia coli bacteria, we predict that the error in estimating Lc without accounting for polydispersity is around 60%. We discuss how to find a binary system with the same CSD as the fully polydisperse mixture. An effective theory is developed and shown to give accurate expressions for the CSD, the effective α, and the average fraction of mobile particles. We give reasons why our qualitative results are expected to be valid for other active matter models and discuss the changes that would result from polydispersity in the active speed rather than in the tumbling rate.

Diversity of self-propulsion speeds reduces motility-induced clustering in confined active matter.

Self-propelled swimmers such as bacteria agglomerate into clusters as a result of their persistent motion. In 1D, those clusters do not coalesce macroscopically and the stationary cluster size distribution (CSD) takes an exponential form. We develop a minimal lattice model for active particles in narrow channels to study how clustering is affected by the interplay between self-propulsion speed diversity and confinement. A mixture of run-and-tumble particles with a distribution of self-propulsion speeds is simulated in 1D. Particles can swap positions at rates proportional to their relative self-propulsion speed. Without swapping, we find that the average cluster size Lc decreases with diversity and follows a non-arithmetic power mean of the single-component Lc's, unlike the case of tumbling-rate diversity previously studied. Effectively, the mixture is thus equivalent to a system of identical particles whose self-propulsion speed is the harmonic mean self-propulsion speed of the mixture. With swapping, particles escape more quickly from clusters. As a consequence, Lc decreases with swapping rates and depends less strongly on diversity. We derive a dynamical equilibrium theory for the CSDs of binary and fully polydisperse systems. Similarly to the clustering behaviour of one-component models, our qualitative results for mixtures are expected to be universal across active matter. Using literature experimental values for the self-propulsion speed diversity of unicellular swimmers known as choanoflagellates, which naturally differentiate into slower and faster cells, we predict that the error in estimating their Lcvia one-component models which use the conventional arithmetic mean self-propulsion speed is around 30%.

From conventional two-stage hepatectomy to ALPPS: Fifteen years of experience in a hepatobiliary surgery unit.

BACKGROUND: Hepatectomy in patients with large tumor load may result in postoperative liver failure and associated complications due to excessive liver parenchyma removal. Conventional two-stage hepatectomy (TSH) and associating liver partition and portal vein ligation for staged hepatectomy (ALPPS) technique are possible solutions to this problem. Colorectal liver metastases (CRLM) is the most frequent indication, and there is a need to assess outcomes for both techniques to improve surgical and long-term oncological outcomes in these patients. METHODS: A single-center retrospective study was designed to compare TSH with ALPPS in patients with initially unresectable bilateral liver tumors between January 2005 and January 2020. ALPPS was performed from January 2012 onwards as the technique of choice. Long-term overall survival (OS) and disease-free survival (DFS) were evaluated as primary outcome in CRLM patients. Postoperative morbidity, mortality and liver growth in all patients were also evaluated. RESULTS: A total of 38 staged hepatectomies were performed: 17 TSH and 21 ALPPS. Complete resection rate was 76.5% (n = 13) in the TSH group and 85.7% (n = 18) in the ALPPS group (P = 0.426). Overall major morbidity (Clavien-Dindo ≥ 3a) (stage 1 + stage 2) was 41.2% (n = 7) in TSH and 33.3% (n = 7) in ALPPS patients (P = 0.389), and perioperative 90-day mortalities were 11.8% (n = 2) vs. 19.0% (n = 4) in each group, respectively (P = 0.654). Intention-to-treat OS rates at 1 and 5 years in CRLM patients for TSH (n = 15) were 80% and 33%, and for ALPPS (n = 17) 76% and 35%, respectively. DFS rates at 1 and 5 years were 36% and 27% in the TSH group vs. 33% and 27% in the ALPPS group, respectively. CONCLUSIONS: ALPPS is an effective alternative to TSH in bilateral affecting liver tumors, allowing higher resection rate, but patients must be carefully selected. In CRLM patients similar long-term OS and DFS can be achieved with both techniques.

Caenorhabditis elegans Multi-Tracker Based on a Modified Skeleton Algorithm.

Automatic tracking of Caenorhabditis elegans (C. egans) in standard Petri dishes is challenging due to high-resolution image requirements when fully monitoring a Petri dish, but mainly due to potential losses of individual worm identity caused by aggregation of worms, overlaps and body contact. To date, trackers only automate tests for individual worm behaviors, canceling data when body contact occurs. However, essays automating contact behaviors still require solutions to this problem. In this work, we propose a solution to this difficulty using computer vision techniques. On the one hand, a skeletonization method is applied to extract skeletons in overlap and contact situations. On the other hand, new optimization methods are proposed to solve the identity problem during these situations. Experiments were performed with 70 tracks and 3779 poses (skeletons) of C. elegans. Several cost functions with different criteria have been evaluated, and the best results gave an accuracy of 99.42% in overlapping with other worms and noise on the plate using the modified skeleton algorithm and 98.73% precision using the classical skeleton algorithm.

Towards Lifespan Automation for Caenorhabditis elegans Based on Deep Learning: Analysing Convolutional and Recurrent Neural Networks for Dead or Live Classification.

The automation of lifespan assays with C. elegans in standard Petri dishes is a challenging problem because there are several problems hindering detection such as occlusions at the plate edges, dirt accumulation, and worm aggregations. Moreover, determining whether a worm is alive or dead can be complex as they barely move during the last few days of their lives. This paper proposes a method combining traditional computer vision techniques with a live/dead C. elegans classifier based on convolutional and recurrent neural networks from low-resolution image sequences. In addition to proposing a new method to automate lifespan, the use of data augmentation techniques is proposed to train the network in the absence of large numbers of samples. The proposed method achieved small error rates (3.54% ± 1.30% per plate) with respect to the manual curve, demonstrating its feasibility.

Rho-kinase pathway activation and apoptosis in circulating leucocytes in patients with heart failure with reduced ejection fraction.

BACKGROUND: Increased Rho-kinase activity in circulating leucocytes is observed in heart failure with reduced ejection fraction (HFrEF). However, there is little information in HFrEF regarding other Rho-kinase pathway components an on the relationship between Rho-kinase and apoptosis. Here, Rho-kinase activation levels and phosphorylation of major downstream molecules and apoptosis levels were measured for the first time both in HFrEF patients and healthy individuals. METHODS: Cross-sectional study comparing HFrEF patients (n = 20) and healthy controls (n = 19). Rho-kinase activity in circulating leucocytes (peripheral blood mononuclear cells, PBMCs) was determined by myosin light chain phosphatase 1 (MYPT1) and ezrin-radixin-moesin (ERM) phosphorylation. Rho-kinase cascade proteins phosphorylation p38-MAPK, myosin light chain-2, JAK and JNK were also analysed along with apoptosis. RESULTS: MYPT1 and ERM phosphorylation were significantly elevated in HFrEF patients, (3.9- and 4.8-fold higher than in controls, respectively). JAK phosphorylation was significantly increased by 300% over controls. Phosphorylation of downstream molecules p38-MAPK and myosin light chain-2 was significantly higher by 360% and 490%, respectively, while JNK phosphorylation was reduced by 60%. Catecholamine and angiotensin II levels were significantly higher in HFrEF patients, while angiotensin-(1-9) levels were lower. Apoptosis in circulating leucocytes was significantly increased in HFrEF patients by 2.8-fold compared with controls and significantly correlated with Rho-kinase activation. CONCLUSION: Rho-kinase pathway is activated in PMBCs from HFrEF patients despite optimal treatment, and it is closely associated with neurohormonal activation and with apoptosis. ROCK cascade inhibition might induce clinical benefits in HFrEF patients, and its assessment in PMBCs could be useful to evaluate reverse remodelling and disease regression.

Promoting the ethical use of safe and effective cell-based products: the Andalusian plan on regenerative medicine.

With regard to regenerative medicine, the expectations generated over the last two decades and the time involved in developing this type of therapies, together with the availability of devices that allow point-of-care treatments through the rapid isolation of cellular or plasma products from patients in the operating theater, represent the perfect breeding ground for the offering of unproven or unregulated therapies on a global scale. A multidisciplinary approach-one based on the collaboration of institutions that, from the perspective of their area of competence, can contribute to reversing this worrying situation-to this problem is essential. It is a priority for local health authorities to take measures that are adapted to the particular situation and regulatory framework of their respective territory. In this article, the authors present the regenerative medicine action plan promoted by the Andalusian Transplant Coordination (i.e., the action plan for the largest region in Spain), highlighting the aspects the authors believe are fundamental to its success. The authors describe, in summary form, the methodology, phases of the plan, actions designed, key collaborators, important milestones achieved and main lessons they have drawn from their experience so that this can serve as an example for other institutions interested in promoting the ethical use of this type of therapy.

Early left atrial dysfunction is associated with suboptimal cardiovascular health.

AIMS: Two-dimensional speckle-tracking echocardiography can assess left atrial (LA) function by measuring atrial volumes and deformation parameters (strain, strain rate). This cross-sectional analysis explores the association between ideal CV health (CVH), LA function, and systemic biomarkers in healthy individuals from the Chilean MAUCO Cohort. METHODS: We enrolled 95 MAUCO participants with different levels of CVH (mean age: 51 ± 8 years). We categorized participants into low or high CVH groups: A: 0-2, or B: 3-6 CVH risk factors. 2D echocardiography, glucose, insulin, total cholesterol, triglycerides, proBNP, hsCRP, insulin resistance index (HOMA), and right and left atrial strain (RASs and LASs, respectively) were determined. RESULTS: LASs was lower in Group A, while systolic and diastolic blood pressure (BP), body mass index (BMI), insulin, HOMA, total cholesterol, triglycerides, and LV and RV end-diastolic volume were significantly higher in Group A than Group B (P < .01). Change in LASs was inversely correlated with insulin (P = .040), HOMA (P = .013), total cholesterol (P = .039), glycemia (P = .018), and BMI (P = .0.037). CONCLUSION: LASs during the reservoir phase was diminished in subjects with a lower level of CVH. Higher insulin, HOMA, total cholesterol, glycemia, and BMI values were associated with decreased LA deformation during the reservoir phase. Morphofunctional alterations of the LA were also identified in the group with suboptimal CVH, as well as BP values in the range of hypertension. LA dysfunction in an asymptomatic population, along with metabolic syndrome, could be an early event in the continuum of CV damage.

Infrared Temperature Measurement Sensors of Overhead Power Conductors.

Efficiency in power lines operation is becoming more crucial as the electrification increases and more renewable energies are connected into the grid. New methods and sensors are being added to create smart grids to face these challenges and conductor temperature sensors are one of them. Contact temperature sensors have several problems regarding safety and electronic damage due to the electromagnetic fields induced on the conductors. The goal of this paper is to describe an infrared temperature measurement sensor and to compare contact and non-contact temperature measurements to estimate the temperature of power lines. Measurements were done for almost a year, storing around 150,000 measures of contact and infrared thermometers for many different weather and load conditions. The results conclude that the infrared system can be successfully used to control the temperature of the overhead conductor within a range of less than 4 ∘C difference with respect to contact temperature methods for the 88% of the samples and less than 6 ∘C for the 99%.

Reducing Results Variance in Lifespan Machines: An Analysis of the Influence of Vibrotaxis on Wild-Type Caenorhabditis elegans for the Death Criterion.

Nowadays, various artificial vision-based machines automate the lifespan assays of C. elegans. These automated machines present wider variability in results than manual assays because in the latter worms can be poked one by one to determine whether they are alive or not. Lifespan machines normally use a "dead or alive criterion" based on nematode position or pose changes, without poking worms. However, worms barely move on their last days of life, even though they are still alive. Therefore, a long monitoring period is necessary to observe motility in order to guarantee worms are actually dead, or a stimulus to prompt worm movement is required to reduce the lifespan variability measure. Here, a new automated vibrotaxis-based method for lifespan machines is proposed as a solution to prompt a motion response in all worms cultured on standard Petri plates in order to better distinguish between live and dead individuals. This simple automated method allows the stimulation of all animals through the whole plate at the same time and intensity, increasing the experiment throughput. The experimental results exhibited improved live-worm detection using this method, and most live nematodes (>93%) reacted to the vibration stimulus. This method increased machine sensitivity by decreasing results variance by approximately one half (from ±1 individual error per plate to ±0.6) and error in lifespan curve was reduced as well (from 2.6% to 1.2%).

Comparing perturbation models for evaluating stability of neuroimaging pipelines.

With an increase in awareness regarding a troubling lack of reproducibility in analytical software tools, the degree of validity in scientific derivatives and their downstream results has become unclear. The nature of reproducibility issues may vary across domains, tools, data sets, and computational infrastructures, but numerical instabilities are thought to be a core contributor. In neuroimaging, unexpected deviations have been observed when varying operating systems, software implementations, or adding negligible quantities of noise. In the field of numerical analysis, these issues have recently been explored through Monte Carlo Arithmetic, a method involving the instrumentation of floating-point operations with probabilistic noise injections at a target precision. Exploring multiple simulations in this context allows the characterization of the result space for a given tool or operation. In this article, we compare various perturbation models to introduce instabilities within a typical neuroimaging pipeline, including (i) targeted noise, (ii) Monte Carlo Arithmetic, and (iii) operating system variation, to identify the significance and quality of their impact on the resulting derivatives. We demonstrate that even low-order models in neuroimaging such as the structural connectome estimation pipeline evaluated here are sensitive to numerical instabilities, suggesting that stability is a relevant axis upon which tools are compared, alongside more traditional criteria such as biological feasibility, computational efficiency, or, when possible, accuracy. Heterogeneity was observed across participants which clearly illustrates a strong interaction between the tool and data set being processed, requiring that the stability of a given tool be evaluated with respect to a given cohort. We identify use cases for each perturbation method tested, including quality assurance, pipeline error detection, and local sensitivity analysis, and make recommendations for the evaluation of stability in a practical and analytically focused setting. Identifying how these relationships and recommendations scale to higher order computational tools, distinct data sets, and their implication on biological feasibility remain exciting avenues for future work.