Monte-Carlo simulation of Callisto's exosphere

We model Callisto's exosphere based on its ice as well as non-ice surface via the use of a Monte-Carlo exosphere model. For the ice component we implement two putative compositions that have been computed from two possible extreme formation scenarios of the satellite. One composition represents the oxidizing state and is based on the assumption that the building blocks of Callisto were formed in the protosolar nebula and the other represents the reducing state of the gas, based on the assumption that the satellite accreted from solids condensed in the jovian sub-nebula. For the non-ice component we implemented the compositions of typical CI as well as L type chondrites. Both chondrite types have been suggested to represent Callisto's non-ice composition best. As release processes we consider surface sublimation, ion sputtering and photon-stimulated desorption. Particles are followed on their individual trajectories until they either escape Callisto's gravitational attraction, return to the surface, are ionized, or are fragmented. Our density profiles show that whereas the sublimated species dominate close to the surface on the sun-lit side, their density profiles (with the exception of H and H-2) decrease much more rapidly than the sputtered particles. The Neutral gas and Ion Mass (NIM) spectrometer, which is part of the Particle Environment Package (PEP), will investigate Callisto's exosphere during the JUICE mission. Our simulations show that NIM will be able to detect sublimated and sputtered particles from both the ice and non-ice surface. NIM's measured chemical composition will allow us to distinguish between different formation scenarios. (C) 2015 Elsevier Inc. All rights reserved.

Cerebellar mechanisms for motor learning: Testing predictions from a large-scale computer simulation

The cerebellum is the major brain structure that contributes to our ability to improve movements through learning and experience. We have combined computer simulations with behavioral and lesion studies to investigate how modification of synaptic strength at two different sites within the cerebellum contributes to a simple form of motor learning—Pavlovian conditioning of the eyelid response. These studies are based on the wealth of knowledge about the intrinsic circuitry and physiology of the cerebellum and the straightforward manner in which this circuitry is engaged during eyelid conditioning. Thus, our simulations are constrained by the well-characterized synaptic organization of the cerebellum and further, the activity of cerebellar inputs during simulated eyelid conditioning is based on existing recording data. These simulations have allowed us to make two important predictions regarding the mechanisms underlying cerebellar function, which we have tested and confirmed with behavioral studies. The first prediction describes the mechanisms by which one of the sites of synaptic modification, the granule to Purkinje cell synapses (gr → Pkj) of the cerebellar cortex, could generate two time-dependent properties of eyelid conditioning—response timing and the ISI function. An empirical test of this prediction using small, electrolytic lesions of the cerebellar cortex revealed the pattern of results predicted by the simulations. The second prediction made by the simulations is that modification of synaptic strength at the other site of plasticity, the mossy fiber to deep nuclei synapses (mf → nuc), is under the control of Purkinje cell activity. The analysis predicts that this property should confer mf → nuc synapses with resistance to extinction. Thus, while extinction processes erase plasticity at the first site, residual plasticity at mf → nuc synapses remains. The residual plasticity at the mf → nuc site confers the cerebellum with the capability for rapid relearning long after the learned behavior has been extinguished. We confirmed this prediction using a lesion technique that reversibly disconnected the cerebellar cortex at various stages during extinction and reacquisition of eyelid responses. The results of these studies represent significant progress toward a complete understanding of how the cerebellum contributes to motor learning. ^

Visualizing the Simulation of 3-D Underwater Sensor Networks

The majority of sensor network research deals with land-based networks, which are essentially two-dimensional, and thus the majority of simulation and animation tools also only handle such networks. Underwater sensor networks on the other hand, are essentially 3D networks because the depth at which a sensor node is located needs to be considered as well. Due to that additional dimension, specialized tools need to be used when conducting simulations for experimentation. The School of Engineering’s Underwater Sensor Network (UWSN) lab is conducting research on underwater sensor networks and requires simulation tools for 3D networks. The lab has extended NS-2, a widely used network simulator, so that it can simulate three-dimensional networks. However, NAM, a widely used network animator, currently only supports two-dimensional networks and no extensions have been implemented to give it three-dimensional capabilities. In this project, we develop a network visualization tool that functions similarly to NAM but is able to render network environments in full 3-D. It is able to take as input a NS-2 trace file (the same file taken as input by NAM), create the environment, position the sensor nodes, and animate the events of the simulation. Further, the visualization tool is easy to use, especially friendly to NAM users, as it is designed to follow the interfaces and functions similar to NAM. So far, the development has fulfilled the basic functionality. Future work includes fully functional capabilities for visualization and much improved user interfaces.

Functional architecture of mammalian horizontal cells

In the rabbit retina, there are two kinds of horizontal cells (HCs). The A-type HC is a large axonless cell which contacts cones exclusively. The B-type HC is an axon bearing cell. While the somatic dendrites of B-type HCs also contact cones, the axon expands into an elaborately branched structure, the axon terminal (AT), which contacts a large number of rods. It is difficult to label the different HCs selectively by immunochemical methods. Therefore, we developed dye injection methods to label each type of HC. Then it was possible, (1) to describe the detailed structure of the AT (2) to identify the glutamate receptors mediating cone input to A and B-type HCs and rod input to ATs and (3) to test the hypothesis that the B-type HCs are coupled via Cx57 gap junctions. ^ To obtain well filled examples of single HCs, it was necessary to block gap junction coupling to stop the spread of Neurobiotin through the network. We used dye coupling in A-type HCs to screen a series of potential gap junction antagonists. One of these compounds, meclofenamic acid (MFA), was potent, water soluble and easily reversible. This compound may be a useful tool to manipulate gap junction coupling. ^ In the presence of MFA, Neurobiotin passed down the axon of B-type HCs to reveal the detailed structure of the AT. We observed that only one AT ending entered each rod spherule invagination. This observation was confirmed by calculation and two dye injections. ^ Glutamate is the neurotransmitter used by both rods and cones. AMPA receptors were colocalized with the dendrites of A and B-type HCs at each cone pedicle. In addition, AMPA receptors were located on the AT ending at each rod spherule. Thus rod and cone input to HCs is mediated by AMPA receptors. ^ A-type and B-type HCs may express different connexins because they have different dye-coupling properties. Recently, we found that connexin50 (Cx50) is expressed by A-type HCs. B-type HCs and B-type ATs are also independently coupled. Cx57 was expressed in the OPL and double label studies showed that Cx 57 was colocalized with the AT matrix but not with the somatic dendrites of B-type HCs. ^ In summary, we have identified a useful gap junction antagonist, MFA. There is one AT ending at each rod spherule, rods inputs to ATs is mediated by AMPA receptors and coupling in the AT matrix is mediated by Cx57. This confirms that HCs with different properties use distinct connexins. The properties of ATs described in this research are consistent. The connections and properties reported here suggest that ATs functions as rod HCs and provide a negative feedback signal to rods. ^

Outcomes of tuberculosis positive non-injection drug users in Harris County from October 1995 through September 2004

Background. The population-based Houston Tuberculosis Initiative (HTI) study has enrolled and gathered demographic, social, behavioral, and disease related data on more than 80% of all reported Mycobacterium Tuberculosis (MTB) cases and 90% of all culture positive patients in Houston/Harris County over a 9 year period (from October 1995-September 2004). During this time period 33% (n=1210) of HTI MTB cases have reported a history of drug use. Of those MTB cases reporting a history of drug use, a majority of them (73.6%), are non-injection drug users (NIDUs). ^ Other than HIV, drug use is the single most important risk factor for progression from latent to infectious tuberculosis (TB). In addition, drug use is associated with increased transmission of active TB, as seen by the increased number of clonally related strains or clusters (see definition on page 30) found in this population. The deregulatory effects of drug use on immune function are well documented. Associations between drug use and increased morbidity have been reported since the late 1970's. However, limited research focused on the immunological consequence of non-injection drug use and its relation to tuberculosis infection among TB patients is available. ^ Methods. TB transmission patterns, symptoms, and prevalence of co-morbidities were a focus of this project. Smoking is known to suppress Nitric Oxide (NO) production and interfere with immune function. In order to limit any possible confounding due to smoking two separate analyses were done. Non-injection drug user smokers (NIDU-S) were compared to non-drug user smokers (NDU-S) and non-injection drug user non-smokers (NIDU-NS) were compared to non-drug user non-smokers (NDU-NS) individually. Specifically proportions, chi-square p-values, and (where appropriate) odds ratios with 95% confidence intervals were calculated to assess characteristics and potential associations of co-morbidities and symptoms of TB among NIDUs HTI TB cases. ^ Results. Significant differences in demographic characteristics and risk factors were found. In addition drug users were found to have a decreased risk for cancer, diabetes mellitus, and chronic pulmonary disease. They were at increased risk of having HIV/AIDS diagnosis, liver disease, and trauma related morbidities. Drug users were more likely to have pulmonary TB disease, and a significantly increased amount of clonally related strains of TB or "clusters" were seen in both smokers and non-smoker drug users when compared to their non-drug user counterparts. Drug users are more likely to belong to print groups (clonally related TB strains with matching spoligotypes) including print one and print three and the Beijing family group, s1. Drug users were found to be no more likely to experience drug resistance to TB therapy and were likely to be cured of disease upon completion of therapy. ^ Conclusion. Drug users demographic and behavioral risk factors put them at an increased risk contracting and spreading TB disease throughout the community. Their increased levels of clustering are evidence of recent transmission and the significance of certain print groups among this population indicate the transmission is from within the social family. For these reasons a focus on this "at risk population" is critical to the success of future public health interventions. Successful completion of directly observed therapy (DOT), the tracking of TB outbreaks and incidence through molecular characterization, and increased diagnostic strategies have led to the stabilization of TB incidence in Houston, Harris County over the past 9 years and proven that the Houston Tuberculosis Initiative has played a critical role in the control and prevention of TB transmission. ^

Ligand-stimulated internalization of a Dictyostelium discoideum G protein-coupled cAMP receptor

The social amoeba, Dictyostelium discoideum, undergoes a remarkable starvation-induced program of development that transforms a population of unicellular amoebae into a fruiting body composed of resistant spores suspended on a stalk. During this development, secreted cAMP drives chemotaxis of the amoebae, leading to their aggregation, and subsequent differentiation and morphogenesis. Four sequentially expressed G protein-coupled receptors (GPCRs) for cAMP play critical roles in this process. The first of these, cAR1, is essential for aggregation as it mediates chemotaxis as well as the propagation of secreted cAMP waves throughout aggregating populations. Ligand-induced internalization has been shown to regulate a variety of GPCRs. However, little was known at the outset of this study about the role of internalization in the regulation of cAR1 function or, for that matter, in developmental systems in general. For this study, cAMP-induced cAR1 internalization was assessed by measuring (1) the reduction of cell surface binding sites for [ 3H]cAMP and (2) the redistribution of YFP-tagged receptors to the cell's interior, cAMP was found to induce little or no loss of ligand binding (LLB) in vegetative cells. However, the ability to induce LLB increased progressively over the initial 6 hrs of development, reaching ∼70% in cells undergoing aggregation. Despite these reductions in surface binding, detectable cAR1-YFP redistribution could be induced by cAMP only after the cells reached the mound stage (10 hrs) and was found to occur naturally by the ensuing slug stage (18 hrs). Site-directed substitution of a cluster of 5 serines in the receptor's cytoplasmic tail that was previously shown to be the principal site of cAMP-induced cAR1 phosphorylation impaired both LLB and receptor redistribution and furthermore resulted in mound-stage developmental arrest, suggesting that phosphorylation of cAR1 is a prerequisite for its internalization and that cAR1 internalization is required for post-aggregative development. To assess the involvement of clathrin mediated endocytosis, Dictyostelium cells lacking the clathrin light chain gene (clc-) or either of two dynamin genes were examined and found to be defective in LLB and, in the case of clc- cells, also cAR1 redistribution and turnover. Furthermore, cAR1 overexpression in clc- cells (like the serine mutant in wild-type cells) promoted developmental arrest in mounds. The mound-arrest phenotype was also recapitulated in a wild-type background by the specific expression of cAR1 in prestalk cells (but not prespore cells), suggesting that development depends critically on internalization and clearance of cAR1 from these cells. Persistent cAR1 expression following aggregation was found to be associated with aberrant expression of prestalk and prespore genes, which may adversely affect development in the prestalk cell lineage. The PI3 kinase-TORC2 signal transduction pathway, known to be important for Dictyostelium chemotaxis and internalization of yeast pheromone receptors, was examined using chemical inhibitors and null cells and found to be necessary for cAR1 internalization. In conclusion, cAR1 was shown to be similar to other GPCRs in that its internalization depends on phosphorylation of cytoplasmic domain serines, utilizes clathrin and dynamin, and involves the TORC2 complex. In addition, the findings presented here that cAR1 internalization is both developmentally regulated and required for normal development represent a novel regulatory paradigm that might pertain to other GPCRs known to play important roles in the development of humans and other metazoans. ^

A simulation study of the standard design, the rolling six design, the CRM, and the modified CRM in Phase I clinical trials

The Phase I clinical trial is considered the "first in human" study in medical research to examine the toxicity of a new agent. It determines the maximum tolerable dose (MTD) of a new agent, i.e., the highest dose in which toxicity is still acceptable. Several phase I clinical trial designs have been proposed in the past 30 years. The well known standard method, so called the 3+3 design, is widely accepted by clinicians since it is the easiest to implement and it does not need a statistical calculation. Continual reassessment method (CRM), a design uses Bayesian method, has been rising in popularity in the last two decades. Several variants of the CRM design have also been suggested in numerous statistical literatures. Rolling six is a new method introduced in pediatric oncology in 2008, which claims to shorten the trial duration as compared to the 3+3 design. The goal of the present research was to simulate clinical trials and compare these phase I clinical trial designs. Patient population was created by discrete event simulation (DES) method. The characteristics of the patients were generated by several distributions with the parameters derived from a historical phase I clinical trial data review. Patients were then selected and enrolled in clinical trials, each of which uses the 3+3 design, the rolling six, or the CRM design. Five scenarios of dose-toxicity relationship were used to compare the performance of the phase I clinical trial designs. One thousand trials were simulated per phase I clinical trial design per dose-toxicity scenario. The results showed the rolling six design was not superior to the 3+3 design in terms of trial duration. The time to trial completion was comparable between the rolling six and the 3+3 design. However, they both shorten the duration as compared to the two CRM designs. Both CRMs were superior to the 3+3 design and the rolling six in accuracy of MTD estimation. The 3+3 design and rolling six tended to assign more patients to undesired lower dose levels. The toxicities were slightly greater in the CRMs.^

A simulation study of the effects of small center enrollment in multi -center clinical trials

Multi-center clinical trials are very common in the development of new drugs and devices. One concern in such trials, is the effect of individual investigational sites enrolling small numbers of patients on the overall result. Can the presence of small centers cause an ineffective treatment to appear effective when treatment-by-center interaction is not statistically significant?^ In this research, simulations are used to study the effect that centers enrolling few patients may have on the analysis of clinical trial data. A multi-center clinical trial with 20 sites is simulated to investigate the effect of a new treatment in comparison to a placebo treatment. Twelve of these 20 investigational sites are considered small, each enrolling less than four patients per treatment group. Three clinical trials are simulated with sample sizes of 100, 170 and 300. The simulated data is generated with various characteristics, one in which treatment should be considered effective and another where treatment is not effective. Qualitative interactions are also produced within the small sites to further investigate the effect of small centers under various conditions.^ Standard analysis of variance methods and the "sometimes-pool" testing procedure are applied to the simulated data. One model investigates treatment and center effect and treatment-by-center interaction. Another model investigates treatment effect alone. These analyses are used to determine the power to detect treatment-by-center interactions, and the probability of type I error.^ We find it is difficult to detect treatment-by-center interactions when only a few investigational sites enrolling a limited number of patients participate in the interaction. However, we find no increased risk of type I error in these situations. In a pooled analysis, when the treatment is not effective, the probability of finding a significant treatment effect in the absence of significant treatment-by-center interaction is well within standard limits of type I error. ^

Integrated Data Sets of the EU FP5 Research Project ORFOIS: Origin and fate of biogenic particle fluxes in the ocean and their interactions with atmospheric CO2 concentrations as well as the marine sediment

For a reliable simulation of the time and space dependent CO2 redistribution between ocean and atmosphere an appropriate time dependent simulation of particle dynamics processes is essential but has not been carried out so far. The major difficulties were the lack of suitable modules for particle dynamics and early diagenesis (in order to close the carbon and nutrient budget) in ocean general circulation models, and the lack of an understanding of biogeochemical processes, such as the partial dissolution of calcareous particles in oversaturated water. The main target of ORFOIS was to fill in this gap in our knowledge and prediction capability infrastructure. This goal has been achieved step by step. At first comprehensive data bases (already existing data) of observations of relevance for the three major types of biogenic particles, organic carbon (POC), calcium carbonate (CaCO3), and biogenic silica (BSi or opal), as well as for refractory particles of terrestrial origin were collated and made publicly available.

Tab. 1 Mole % CO2 and 13C values of CO2 at baseline and at the final monitoring event for the eight observation wells

During CO2 storage operations in mature oilfields or saline aquifers it is desirable to trace the movement of injected CO2 for verification and safety purposes. We demonstrate the successful use of carbon isotope abundance ratios for tracing the movement of CO2 injected at the Cardium CO2 Storage Monitoring project in Alberta between 2005 and 2007. Injected CO2 had a d13C value of -4.6±1.1 per mil that was more than 10 per mil higher than the carbon isotope ratios of casing gas CO2 prior to CO2 injection with average d13C values ranging from -15.9 to -23.5 per mil. After commencement of CO2 injection, d13C values of casing gas CO2 increased in all observation wells towards those of the injected CO2 consistent with a two-source end-member mixing model. At four wells located in a NE-SW trend with respect to the injection wells, breakthrough of injected CO2 was registered chemically (>50 mol % CO2) and isotopically 1-6 months after commencement of CO2 injection resulting in cumulative CO2 fluxes exceeding 100000 m**3 during the observation period. At four other wells, casing gas CO2 contents remained below 5 mol % resulting in low cumulative CO2 fluxes (<2000 m**3) throughout the entire observation period, but carbon isotope ratios indicated contributions between <30 and 80% of injected CO2. Therefore, we conclude that monitoring the movement of CO2 in the injection reservoir with geochemical and isotopic techniques is an effective approach to determine plume expansion and to identify potential preferential flow paths provided that the isotopic composition of injected CO2 is constant and distinct from that of baseline CO2.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Table 1. Mean geochemical values for Baseline (pre-injection) and Monitor-1 (post-injection) samples, and the percent change between the two weighted means

The Weyburn Oil Field, Saskatchewan is the site of a large (5000 tonnes/day of CO2) CO2-EOR injection project By EnCana Corporation. Pre- and post-injection samples (Baseline and Monitor-1, respectively) of produced fluids from approximately 45 vertical wells were taken and chemically analyzed to determine changes in the fluid chemistry and isotope composition between August 2000 and March 2001. After 6 months of CO2 injection, geochemical parameters including pH, [HCO3], [Ca], [Mg], and ?13CO2(g) point to areas in which injected CO2 dissolution and reservoir carbonate mineral dissolution have occurred. Pre-injection fluid compositions suggest that the reservoir brine in the injection area may be capable of storing as much as 100 million tonnes of dissolved CO2. Modeling of water-rock reactions show that clay minerals and feldspar, although volumetrically insignificant, may be capable of acting as pH buffers, allowing injected CO2 to be stored as bicarbonate in the formation water or as newly precipitated carbonate minerals, given favorable reaction kinetics.

A coupled model for carbon and radiocarbon evolution during the last deglaciation, link to model results

Changes in the ventilation of the Southern Ocean are thought to play an important role on deglacial carbon and radiocarbon evolution, but have not been tested within a coupled climate-carbon model. Here, we present such a simulation based on a simple scenario of transient deglacial sinking of brines - sea-ice salt rejections - around Antarctica, which modulates Southern Ocean ventilation. This experiment is able to reproduce deglacial atmospheric changes in carbon and radiocarbon but also ocean radiocarbon records measured in the Atlantic, Southern and Pacific Oceans. Simulated for the first time in a fully coupled climate-carbon model including radiocarbon, our modeling results suggest that the deglacial changes in atmospheric carbon dioxide and radiocarbon were achieved by means of a breakdown in the glacial brine-induced stratification of the Southern Ocean.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

Simulation Model of Benthic Antarctic Assemblages (SIMBAA), link to PC-executable

SIMBAA is a spatially explicit, individual-based simulation model. It was developed to analyse the response of populations of Antarctic benthic species and their diversity to iceberg scouring. This disturbance is causing a high local mortality providing potential space for new colonisation. Traits can be attributed to model species, e.g. in terms of reproduction, dispersal, and life span. Physical disturbances can be designed in space and time, e.g. in terms of size, shape, and frequency. Environmental heterogeneity can be considered by cell-specific capacities to host a certain number of individuals. When grid cells become empty (after a disturbance event or due to natural mortality of of an individual), a lottery decides which individual from which species stored in a pool of candidates (for this cell) will recruit in that cell. After a defined period the individuals become mature and their offspring are dispersed and stored in the pool of candidates. The biological parameters and disturbance regimes decide on how long an individual lives. Temporal development of single populations of species as well as Shannon diversity are depicted in the main window graphically and primary values are listed. Examples for simulations can be loaded and saved as sgf-files. The results are also shown in an additional window in a dimensionless area with 50 x 50 cells, which contain single individuals depicted as circles; their colour indicates the assignment to the self-designed model species and the size represents their age. Dominant species per cell and disturbed areas can also be depicted. Output of simulation runs can be saved as images, which can be assembled to video-clips by standard computer programs (see GIF-examples of which "Demo 1" represents the response of the Antarctic benthos to iceberg scouring and "Demo 2" represents a simulation of a deep-sea benthic habitat).