Clay mineral continental amplifier for marine carbon sequestration in a greenhouse ocean

PubMed Central

Kennedy, Martin J.; Wagner, Thomas

2011-01-01

The majority of carbon sequestration at the Earth’s surface occurs in marine continental margin settings within fine-grained sediments whose mineral properties are a function of continental climatic conditions. We report very high mineral surface area (MSA) values of 300 and 570 m2 g in Late Cretaceous black shales from Ocean Drilling Program site 959 of the Deep Ivorian Basin that vary on subcentennial time scales corresponding with abrupt increases from approximately 3 to approximately 18% total organic carbon (TOC). The observed MSA changes with TOC across multiple scales of variability and on a sample-by-sample basis (centimeter scale), provides a rigorous test of a hypothesized influence on organic carbon burial by detrital clay mineral controlled MSA. Changes in TOC also correspond with geochemical and sedimentological evidence for water column anoxia. Bioturbated intervals show a lower organic carbon loading on mineral surface area of 0.1 mg-OC m-2 when compared to 0.4 mg-OC m-2 for laminated and sulfidic sediments. Although either anoxia or mineral surface protection may be capable of producing TOC of < 5%, when brought together they produced the very high TOC (10–18%) apparent in these sediments. This nonlinear response in carbon burial resulted from minor precession-driven changes of continental climate influencing clay mineral properties and runoff from the African continent. This study identifies a previously unrecognized land–sea connection among continental weathering, clay mineral production, and anoxia and a nonlinear effect on marine carbon sequestration during the Coniacian-Santonian Oceanic Anoxic Event 3 in the tropical eastern Atlantic. PMID:21576498

Synthesis of wavelet envelope in 2-D random media having power-law spectra: comparison with FD simulations

NASA Astrophysics Data System (ADS)

Sato, Haruo; Fehler, Michael C.

2016-10-01

The envelope broadening and the peak delay of the S-wavelet of a small earthquake with increasing travel distance are results of scattering by random velocity inhomogeneities in the earth medium. As a simple mathematical model, Sato proposed a new stochastic synthesis of the scalar wavelet envelope in 3-D von Kármán type random media when the centre wavenumber of the wavelet is in the power-law spectral range of the random velocity fluctuation. The essential idea is to split the random medium spectrum into two components using the centre wavenumber as a reference: the long-scale (low-wavenumber spectral) component produces the peak delay and the envelope broadening by multiple scattering around the forward direction; the short-scale (high-wavenumber spectral) component attenuates wave amplitude by wide angle scattering. The former is calculated by the Markov approximation based on the parabolic approximation and the latter is calculated by the Born approximation. Here, we extend the theory for the envelope synthesis of a wavelet in 2-D random media, which makes it easy to compare with finite difference (FD) simulation results. The synthetic wavelet envelope is analytically written by using the random medium parameters in the angular frequency domain. For the case that the power spectral density function of the random velocity fluctuation has a steep roll-off at large wavenumbers, the envelope broadening is small and frequency independent, and scattering attenuation is weak. For the case of a small roll-off, however, the envelope broadening is large and increases with frequency, and the scattering attenuation is strong and increases with frequency. As a preliminary study, we compare synthetic wavelet envelopes with the average of FD simulation wavelet envelopes in 50 synthesized random media, which are characterized by the RMS fractional velocity fluctuation ε = 0.05, correlation scale a = 5 km and the background wave velocity V0 = 4 km s-1. We use the radiation of a 2 Hz Ricker wavelet from a point source. For all the cases of von Kármán order κ = 0.1, 0.5 and 1, we find the synthetic wavelet envelopes are a good match to the characteristics of FD simulation wavelet envelopes in a time window starting from the onset through the maximum peak to the time when the amplitude decreases to half the peak amplitude.

Learning to Question: The Roles of Multiple Hypotheses, Successive Approximations, Balloons and Toilet Paper in University Science Programs of Southwestern Amazonia

ERIC Educational Resources Information Center

Brown, I. Foster

2008-01-01

Learning to question is essential for determining pathways of conservation and development in southwestern Amazonia during a time of rapid global environmental change. Teaching such an approach in graduate science programs in regional universities can be done using play-acting and simulation exercises. Multiple working hypotheses help students…

Evolution of a hybrid micro-macro entangled state of the qubit-oscillator system via the generalized rotating wave approximation

NASA Astrophysics Data System (ADS)

Chakrabarti, R.; Yogesh, V.

2016-04-01

We study the evolution of the hybrid entangled states in a bipartite (ultra) strongly coupled qubit-oscillator system. Using the generalized rotating wave approximation the reduced density matrices of the qubit and the oscillator are obtained. The reduced density matrix of the oscillator yields the phase space quasi probability distributions such as the diagonal P-representation, the Wigner W-distribution and the Husimi Q-function. In the strong coupling regime the Q-function evolves to uniformly separated macroscopically distinct Gaussian peaks representing ‘kitten’ states at certain specified times that depend on multiple time scales present in the interacting system. The ultrastrong coupling strength of the interaction triggers appearance of a large number of modes that quickly develop a randomization of their phase relationships. A stochastic averaging of the dynamical quantities sets in, and leads to the decoherence of the system. The delocalization in the phase space of the oscillator is studied by using the Wehrl entropy. The negativity of the W-distribution reflects the departure of the oscillator from the classical states, and allows us to study the underlying differences between various information-theoretic measures such as the Wehrl entropy and the Wigner entropy. Other features of nonclassicality such as the existence of the squeezed states and appearance of negative values of the Mandel parameter are realized during the course of evolution of the bipartite system. In the parametric regime studied here these properties do not survive in the time-averaged limit.

Multiscale skeletal representation of images via Voronoi diagrams

NASA Astrophysics Data System (ADS)

Marston, R. E.; Shih, Jian C.

1995-08-01

Polygonal approximations to skeletal or stroke-based representations of 2D objects may consume less storage and be sufficient to describe their shape for many applications. Multi- scale descriptions of object outlines are well established but corresponding methods for skeletal descriptions have been slower to develop. In this paper we offer a method of generating scale-based skeletal representation via the Voronoi diagram. The method has the advantages of less time complexity, a closer relationship between the skeletons at each scale and better control over simplification of the skeleton at lower scales. This is because the algorithm starts by generating the skeleton at the coarsest scale first, then it produces each finer scale, in an iterative manner, directly from the level below. The skeletal approximations produced by the algorithm also benefit from a strong relationship with the object outline, due to the structure of the Voronoi diagram.

Entanglement renormalization, quantum error correction, and bulk causality

NASA Astrophysics Data System (ADS)

Kim, Isaac H.; Kastoryano, Michael J.

2017-04-01

Entanglement renormalization can be viewed as an encoding circuit for a family of approximate quantum error correcting codes. The logical information becomes progres-sively more well-protected against erasure errors at larger length scales. In particular, an approximate variant of holographic quantum error correcting code emerges at low energy for critical systems. This implies that two operators that are largely separated in scales behave as if they are spatially separated operators, in the sense that they obey a Lieb-Robinson type locality bound under a time evolution generated by a local Hamiltonian.

Fast discovery and visualization of conserved regions in DNA sequences using quasi-alignment

PubMed Central

2013-01-01

Background Next Generation Sequencing techniques are producing enormous amounts of biological sequence data and analysis becomes a major computational problem. Currently, most analysis, especially the identification of conserved regions, relies heavily on Multiple Sequence Alignment and its various heuristics such as progressive alignment, whose run time grows with the square of the number and the length of the aligned sequences and requires significant computational resources. In this work, we present a method to efficiently discover regions of high similarity across multiple sequences without performing expensive sequence alignment. The method is based on approximating edit distance between segments of sequences using p-mer frequency counts. Then, efficient high-throughput data stream clustering is used to group highly similar segments into so called quasi-alignments. Quasi-alignments have numerous applications such as identifying species and their taxonomic class from sequences, comparing sequences for similarities, and, as in this paper, discovering conserved regions across related sequences. Results In this paper, we show that quasi-alignments can be used to discover highly similar segments across multiple sequences from related or different genomes efficiently and accurately. Experiments on a large number of unaligned 16S rRNA sequences obtained from the Greengenes database show that the method is able to identify conserved regions which agree with known hypervariable regions in 16S rRNA. Furthermore, the experiments show that the proposed method scales well for large data sets with a run time that grows only linearly with the number and length of sequences, whereas for existing multiple sequence alignment heuristics the run time grows super-linearly. Conclusion Quasi-alignment-based algorithms can detect highly similar regions and conserved areas across multiple sequences. Since the run time is linear and the sequences are converted into a compact clustering model, we are able to identify conserved regions fast or even interactively using a standard PC. Our method has many potential applications such as finding characteristic signature sequences for families of organisms and studying conserved and variable regions in, for example, 16S rRNA. PMID:24564200

Fast discovery and visualization of conserved regions in DNA sequences using quasi-alignment.

PubMed

Nagar, Anurag; Hahsler, Michael

2013-01-01

Next Generation Sequencing techniques are producing enormous amounts of biological sequence data and analysis becomes a major computational problem. Currently, most analysis, especially the identification of conserved regions, relies heavily on Multiple Sequence Alignment and its various heuristics such as progressive alignment, whose run time grows with the square of the number and the length of the aligned sequences and requires significant computational resources. In this work, we present a method to efficiently discover regions of high similarity across multiple sequences without performing expensive sequence alignment. The method is based on approximating edit distance between segments of sequences using p-mer frequency counts. Then, efficient high-throughput data stream clustering is used to group highly similar segments into so called quasi-alignments. Quasi-alignments have numerous applications such as identifying species and their taxonomic class from sequences, comparing sequences for similarities, and, as in this paper, discovering conserved regions across related sequences. In this paper, we show that quasi-alignments can be used to discover highly similar segments across multiple sequences from related or different genomes efficiently and accurately. Experiments on a large number of unaligned 16S rRNA sequences obtained from the Greengenes database show that the method is able to identify conserved regions which agree with known hypervariable regions in 16S rRNA. Furthermore, the experiments show that the proposed method scales well for large data sets with a run time that grows only linearly with the number and length of sequences, whereas for existing multiple sequence alignment heuristics the run time grows super-linearly. Quasi-alignment-based algorithms can detect highly similar regions and conserved areas across multiple sequences. Since the run time is linear and the sequences are converted into a compact clustering model, we are able to identify conserved regions fast or even interactively using a standard PC. Our method has many potential applications such as finding characteristic signature sequences for families of organisms and studying conserved and variable regions in, for example, 16S rRNA.

Wechsler Adult Intelligence Scale-Fourth Edition performance in relapsing-remitting multiple sclerosis.

PubMed

Ryan, Joseph J; Gontkovsky, Samuel T; Kreiner, David S; Tree, Heather A

2012-01-01

Forty patients with relapsing-remitting multiple sclerosis (MS) completed the 10 core Wechsler Adult Intelligence Scale-Fourth Edition (WAIS-IV) subtests. Means for age and education were 42.05 years (SD = 9.94) and 14.33 years (SD = 2.40). For all participants, the native language was English. The mean duration of MS diagnosis was 8.17 years (SD = 7.75), and the mean Expanded Disability Status Scale (EDSS; Kurtzke, 1983 ) score was 3.73 (SD = 1.41) with a range from 2.0 to 6.5. A control group of healthy individuals with similar demographic characteristics also completed the WAIS-IV and were provided by the test publisher. Compared to controls, patients with MS earned significantly lower subtest and composite scores. The patients' mean scores were consistently in the low-average to average range, and the patterns of performance across groups did not differ significantly, although there was a trend towards higher scores on the Verbal Comprehension Index (VCI) and lower scores on the Processing Speed Index (PSI). Approximately 78% of patients had actual Full Scale IQs that were significantly lower than preillness, demographically based IQ estimates.

Hierarchical Address Event Routing for Reconfigurable Large-Scale Neuromorphic Systems.

PubMed

Park, Jongkil; Yu, Theodore; Joshi, Siddharth; Maier, Christoph; Cauwenberghs, Gert

2017-10-01

We present a hierarchical address-event routing (HiAER) architecture for scalable communication of neural and synaptic spike events between neuromorphic processors, implemented with five Xilinx Spartan-6 field-programmable gate arrays and four custom analog neuromophic integrated circuits serving 262k neurons and 262M synapses. The architecture extends the single-bus address-event representation protocol to a hierarchy of multiple nested buses, routing events across increasing scales of spatial distance. The HiAER protocol provides individually programmable axonal delay in addition to strength for each synapse, lending itself toward biologically plausible neural network architectures, and scales across a range of hierarchies suitable for multichip and multiboard systems in reconfigurable large-scale neuromorphic systems. We show approximately linear scaling of net global synaptic event throughput with number of routing nodes in the network, at 3.6×10 7 synaptic events per second per 16k-neuron node in the hierarchy.

Correspondence: Reply to ‘Phantom phonon localization in relaxors’

DOE PAGES

Manley, Michael E.; Abernathy, Douglas L.; Budai, John D.

2017-12-05

The Correspondence by Gehring et al. mistakes Anderson phonon localization for the concept of an atomic-scale local mode. An atomic-scale local mode refers to a single atom vibrating on its own within a crystal. Such a local mode will have an almost flat intensity profile, but this is not the same as phonon localization. Anderson localization is a wave interference effect in a disordered system that results in waves becoming spatially localized. The length scale of the localized waves is set by the wavelength, which is approximately 2 nm in this case. This larger length scale in real space meansmore » narrower intensity profiles in reciprocal space. Here, we conclude that the claims in the Correspondence by Gehring et al. are incorrect because they mistakenly assume that the length scale for Anderson localization is atomic, and because the experimental observations rule out multiple scattering as the origin.« less

Correspondence: Reply to ‘Phantom phonon localization in relaxors’

DOE Office of Scientific and Technical Information (OSTI.GOV)

Manley, Michael E.; Abernathy, Douglas L.; Budai, John D.

The Correspondence by Gehring et al. mistakes Anderson phonon localization for the concept of an atomic-scale local mode. An atomic-scale local mode refers to a single atom vibrating on its own within a crystal. Such a local mode will have an almost flat intensity profile, but this is not the same as phonon localization. Anderson localization is a wave interference effect in a disordered system that results in waves becoming spatially localized. The length scale of the localized waves is set by the wavelength, which is approximately 2 nm in this case. This larger length scale in real space meansmore » narrower intensity profiles in reciprocal space. Here, we conclude that the claims in the Correspondence by Gehring et al. are incorrect because they mistakenly assume that the length scale for Anderson localization is atomic, and because the experimental observations rule out multiple scattering as the origin.« less

Effects of simplifying fracture network representation on inert chemical migration in fracture-controlled aquifers

USGS Publications Warehouse

Wellman, Tristan; Shapiro, Allen M.; Hill, Mary C.

2009-01-01

While it is widely recognized that highly permeable 'large-scale' fractures dominate chemical migration in many fractured aquifers, recent studies suggest that the pervasive 'small-scale' fracturing once considered of less significance can be equally important for characterizing the spatial extent and residence time associated with transport processes. A detailed examination of chemical migration through fracture-controlled aquifers is used to advance this conceptual understanding. The influence of fracture structure is evaluated by quantifying the effects to transport caused by a systematic removal of fractures from three-dimensional discrete fracture models whose attributes are derived from geologic and hydrologic conditions at multiple field sites. Results indicate that the effects to transport caused by network simplification are sensitive to the fracture network characteristics, degree of network simplification, and plume travel distance, but primarily in an indirect sense since correlation to individual attributes is limited. Transport processes can be 'enhanced' or 'restricted' from network simplification meaning that the elimination of fractures may increase or decrease mass migration, mean travel time, dispersion, and tailing of the concentration plume. The results demonstrate why, for instance, chemical migration may not follow the classic advection-dispersion equation where dispersion approximates the effect of the ignored geologic structure as a strictly additive process to the mean flow. The analyses further reveal that the prediction error caused by fracture network simplification is reduced by at least 50% using the median estimate from an ensemble of simplified fracture network models, and that the error from network simplification is at least 70% less than the stochastic variability from multiple realizations. Copyright 2009 by the American Geophysical Union.

The Influence of temporal sampling regime on the WFD classification of catchments within the Eden Demonstration Test Catchment Project

NASA Astrophysics Data System (ADS)

Jonczyk, Jennine; Haygarth, Phil; Quinn, Paul; Reaney, Sim

2014-05-01

A high temporal resolution data set from the Eden Demonstration Test Catchment (DTC) project is used to investigate the processes causing pollution and the influence of temporal sampling regime on the WFD classification of three catchments. This data highlights WFD standards may not be fit for purpose. The Eden DTC project is part of a UK government-funded project designed to provide robust evidence regarding how diffuse pollution can be cost-effectively controlled to improve and maintain water quality in rural river catchments. The impact of multiple water quality parameters on ecosystems and sustainable food production are being studied at the catchment scale. Three focus catchments approximately 10 km2 each, have been selected to represent the different farming practices and geophysical characteristics across the Eden catchment, Northern England. A field experimental programme has been designed to monitor the dynamics of agricultural diffuse pollution at multiple scales using state of the art sensors providing continuous real time data. The data set, which includes Total Phosphorus and Total Reactive Phosphorus, Nitrate, Ammonium, pH, Conductivity, Turbidity and Chlorophyll a reveals the frequency and duration of nutrient concentration target exceedance which arises from the prevalence of storm events of increasing magnitude. This data set is sub-sampled at different time intervals to explore how different sampling regimes affects our understanding of nutrient dynamics and the ramification of the different regimes to WFD chemical status. This presentation seeks to identify an optimum temporal resolution of data for effective catchment management and to question the usefulness of the WFD status metric for determining health of a system. Criteria based on high frequency short duration events needs to be accounted for.

Heat, temperature and Clausius inequality in a model for active Brownian particles

PubMed Central

Marconi, Umberto Marini Bettolo; Puglisi, Andrea; Maggi, Claudio

2017-01-01

Methods of stochastic thermodynamics and hydrodynamics are applied to a recently introduced model of active particles. The model consists of an overdamped particle subject to Gaussian coloured noise. Inspired by stochastic thermodynamics, we derive from the system’s Fokker-Planck equation the average exchanges of heat and work with the active bath and the associated entropy production. We show that a Clausius inequality holds, with the local (non-uniform) temperature of the active bath replacing the uniform temperature usually encountered in equilibrium systems. Furthermore, by restricting the dynamical space to the first velocity moments of the local distribution function we derive a hydrodynamic description where local pressure, kinetic temperature and internal heat fluxes appear and are consistent with the previous thermodynamic analysis. The procedure also shows under which conditions one obtains the unified coloured noise approximation (UCNA): such an approximation neglects the fast relaxation to the active bath and therefore yields detailed balance and zero entropy production. In the last part, by using multiple time-scale analysis, we provide a constructive method (alternative to UCNA) to determine the solution of the Kramers equation and go beyond the detailed balance condition determining negative entropy production. PMID:28429787

Heat, temperature and Clausius inequality in a model for active Brownian particles.

PubMed

Marconi, Umberto Marini Bettolo; Puglisi, Andrea; Maggi, Claudio

2017-04-21

Methods of stochastic thermodynamics and hydrodynamics are applied to a recently introduced model of active particles. The model consists of an overdamped particle subject to Gaussian coloured noise. Inspired by stochastic thermodynamics, we derive from the system's Fokker-Planck equation the average exchanges of heat and work with the active bath and the associated entropy production. We show that a Clausius inequality holds, with the local (non-uniform) temperature of the active bath replacing the uniform temperature usually encountered in equilibrium systems. Furthermore, by restricting the dynamical space to the first velocity moments of the local distribution function we derive a hydrodynamic description where local pressure, kinetic temperature and internal heat fluxes appear and are consistent with the previous thermodynamic analysis. The procedure also shows under which conditions one obtains the unified coloured noise approximation (UCNA): such an approximation neglects the fast relaxation to the active bath and therefore yields detailed balance and zero entropy production. In the last part, by using multiple time-scale analysis, we provide a constructive method (alternative to UCNA) to determine the solution of the Kramers equation and go beyond the detailed balance condition determining negative entropy production.

Electron diffusion region during magnetopause reconnection with an intermediate guide field: Magnetospheric multiscale observations

NASA Astrophysics Data System (ADS)

Chen, L.-J.; Hesse, M.; Wang, S.; Gershman, D.; Ergun, R. E.; Burch, J.; Bessho, N.; Torbert, R. B.; Giles, B.; Webster, J.; Pollock, C.; Dorelli, J.; Moore, T.; Paterson, W.; Lavraud, B.; Strangeway, R.; Russell, C.; Khotyaintsev, Y.; Lindqvist, P.-A.; Avanov, L.

2017-05-01

An electron diffusion region (EDR) in magnetic reconnection with a guide magnetic field approximately 0.2 times the reconnecting component is encountered by the four Magnetospheric Multiscale spacecraft at the Earth's magnetopause. The distinct substructures in the EDR on both sides of the reconnecting current sheet are visualized with electron distribution functions that are 2 orders of magnitude higher cadence than ever achieved to enable the following new findings: (1) Motion of the demagnetized electrons plays an important role to sustain the reconnection current and contributes to the dissipation due to the nonideal electric field, (2) the finite guide field dominates over the Hall magnetic field in an electron-scale region in the exhaust and modifies the electron flow dynamics in the EDR, (3) the reconnection current is in part carried by inflowing field-aligned electrons in the magnetosphere part of the EDR, and (4) the reconnection electric field measured by multiple spacecraft is uniform over at least eight electron skin depths and corresponds to a reconnection rate of approximately 0.1. The observations establish the first look at the structure of the EDR under a weak but not negligible guide field.

Detection of F-region electron density irregularities using incoherent-scatter radar

NASA Astrophysics Data System (ADS)

Gudivada, Krishna Prasad

Incoherent-scatter radar data from Poker Flat, Alaska has been used to determine size distributions of electron density structures in the evening time sector of the auroral zone. At high latitudes ionospheric plasma typically moves east-west with speeds of several hundred meters per second. Density irregularities that rapidly move through the radar beam are therefore observed as time-varying power fluctuations. The new phased array radar used for this study has been operated with several antenna directions with successive pulses transmitted in each direction. It is therefore possible to observe plasma Doppler velocities in multiple directions and determine the vector direction of the plasma motion. This near-simultaneous observation of the plasma velocity in conjunction with the electron density height profile data enable a new technique to determine the scale sizes of electron density fluctuations that move horizontally through the radar beam. The study focuses on the collision-less F-region ionosphere where the plasma drift is approximately constant with altitude. The experimental technique limits the range of scale sizes that may be studied to relatively large-scale sizes (i.e. greater than few tens of km). Results show that during magnetically disturbed conditions (Kp ≥ 4) when westward plasma velocities are relatively high (500-1000 m/s) the scale sizes of irregularities (often called plasma blobs) are in the range of 100-300 km and predominantly originate from the polar cap and are transported over long distances (˜1000 km) due to the long chemical recombination times (30-90 minutes). Some irregularities are caused by local auroral particle precipitation and have been identified with associated electron temperature enhancements. For cases of low magnetic activity (Kp ≤ 1), when the radar is located in a region of low plasma velocities (100-500 m/s) well south of the auroral oval (essentially a mid-latitude type ionosphere), the density distribution is always biased strongly toward small-scale sizes (less than 50 km).

Moiré patterns in doubly differential electron-momentum distributions in atomic ionization by mid-infrared lasers

NASA Astrophysics Data System (ADS)

Dran, Martín; Arbó, Diego G.

2018-05-01

We analyze the doubly differential electron momentum distribution in above-threshold ionization of atomic hydrogen by a linearly polarized mid-infrared laser pulse. We reproduce side rings in the momentum distribution with forward-backward symmetry previously observed by Lemell et al. [Phys. Rev. A 87, 013421 (2013), 10.1103/PhysRevA.87.013421], whose origin, as far as we know, has not been explained so far. By developing a Fourier theory of moiré patterns, we demonstrate that such structures stem from the interplay between intra- and intercycle interference patterns which work as two separate grids in the two-dimensional momentum domain. We use a three-dimensional (3D) description based on the saddle-point approximation (SPA) to unravel the nature of these structures. When the periods of the two grids (intra- and intercycle) are similar, principal moiré patterns arise symmetrically as concentric rings in the forward and backward directions at high electron kinetic energy. Higher order moiré patterns are observed and characterized when the period of one grid is multiple of the other. We find a scale law for the position (in momentum space) of the center of the moiré rings in the tunneling regime. We verify the SPA predictions by comparison with time-dependent distorted-wave strong-field approximation calculations and the solutions of the full 3D time-dependent Schrödinger equation.

Hierarchical Modeling and Robust Synthesis for the Preliminary Design of Large Scale Complex Systems

NASA Technical Reports Server (NTRS)

Koch, Patrick N.

1997-01-01

Large-scale complex systems are characterized by multiple interacting subsystems and the analysis of multiple disciplines. The design and development of such systems inevitably requires the resolution of multiple conflicting objectives. The size of complex systems, however, prohibits the development of comprehensive system models, and thus these systems must be partitioned into their constituent parts. Because simultaneous solution of individual subsystem models is often not manageable iteration is inevitable and often excessive. In this dissertation these issues are addressed through the development of a method for hierarchical robust preliminary design exploration to facilitate concurrent system and subsystem design exploration, for the concurrent generation of robust system and subsystem specifications for the preliminary design of multi-level, multi-objective, large-scale complex systems. This method is developed through the integration and expansion of current design techniques: Hierarchical partitioning and modeling techniques for partitioning large-scale complex systems into more tractable parts, and allowing integration of subproblems for system synthesis; Statistical experimentation and approximation techniques for increasing both the efficiency and the comprehensiveness of preliminary design exploration; and Noise modeling techniques for implementing robust preliminary design when approximate models are employed. Hierarchical partitioning and modeling techniques including intermediate responses, linking variables, and compatibility constraints are incorporated within a hierarchical compromise decision support problem formulation for synthesizing subproblem solutions for a partitioned system. Experimentation and approximation techniques are employed for concurrent investigations and modeling of partitioned subproblems. A modified composite experiment is introduced for fitting better predictive models across the ranges of the factors, and an approach for constructing partitioned response surfaces is developed to reduce the computational expense of experimentation for fitting models in a large number of factors. Noise modeling techniques are compared and recommendations are offered for the implementation of robust design when approximate models are sought. These techniques, approaches, and recommendations are incorporated within the method developed for hierarchical robust preliminary design exploration. This method as well as the associated approaches are illustrated through their application to the preliminary design of a commercial turbofan turbine propulsion system. The case study is developed in collaboration with Allison Engine Company, Rolls Royce Aerospace, and is based on the Allison AE3007 existing engine designed for midsize commercial, regional business jets. For this case study, the turbofan system-level problem is partitioned into engine cycle design and configuration design and a compressor modules integrated for more detailed subsystem-level design exploration, improving system evaluation. The fan and low pressure turbine subsystems are also modeled, but in less detail. Given the defined partitioning, these subproblems are investigated independently and concurrently, and response surface models are constructed to approximate the responses of each. These response models are then incorporated within a commercial turbofan hierarchical compromise decision support problem formulation. Five design scenarios are investigated, and robust solutions are identified. The method and solutions identified are verified by comparison with the AE3007 engine. The solutions obtained are similar to the AE3007 cycle and configuration, but are better with respect to many of the requirements.

Spin-locking of half-integer quadrupolar nuclei in nuclear magnetic resonance of solids: creation and evolution of coherences.

PubMed

Ashbrook, Sharon E; Wimperis, Stephen

2004-02-08

Spin-locking of half-integer quadrupolar nuclei, such as 23Na (I=3/2) and 27Al (I=5/2), is of renewed interest owing to the development of variants of the multiple-quantum and satellite-transition magic angle spinning (MAS) nuclear magnetic resonance experiments that either utilize spin-locking directly or offer the possibility that spin-locked states may arise. However, the large magnitude and, under MAS, the time dependence of the quadrupolar interaction often result in complex spin-locking phenomena that are not widely understood. Here we show that, following the application of a spin-locking pulse, a variety of coherence transfer processes occur on a time scale of approximately 1/omegaQ before the spin system settles down into a spin-locked state which may itself be time dependent if MAS is performed. We show theoretically for both spin I=3/2 and 5/2 nuclei that the spin-locked state created by this initial rapid dephasing typically consists of a variety of single- and multiple-quantum coherences and nonequilibrium population states and we discuss the subsequent evolution of these under MAS. In contrast to previous work, we consider spin-locking using a wide range of radio frequency field strengths, i.e., a range that covers both the "strong-field" (omega1 > omegaQPAS and "weak-field" (omega1 < omegaQPAS limits. Single- and multiple-quantum filtered spin-locking experiments on NaNO2, NaNO3, and Al(acac)3, under both static and MAS conditions, are used to illustrate and confirm the results of the theoretical discussion.

Danube Delta Coastline Dynamics in the Last 160 Years

NASA Astrophysics Data System (ADS)

Tătui, Florin; Vespremeanu-Stroe, Alfred; Constantinescu, Ştefan; Zăinescu, Florin

2017-04-01

Wave-dominated deltaic coasts depend on the balance between wave climate and sediment supply, which controls the medium and long-term shoreline evolution. Interestingly, the common plan shapes of the wave-dominated lobes impose different wave exposures and longshore sediment transport magnitudes on the lobe flanks, characterized by ever changing aspects which make these sandy coasts some of the most mobile world coastlines. The Danube Delta coast consists of approximately 220 km (both Romanian and Ukrainian sectors) of tideless, medium-energy low-lying sandy beaches interrupted by multiple river mouths and, sometimes, by engineering structures (Sulina jetties and Midia harbour). The objective of this study is to examine and explain the factors which have driven the Danube Delta coastline dynamics at multi-annual to multi-decadal and centennial time-scales. Our analysis is based on multiple shorelines extracted from historical and modern maps (since mid-19th century), recent medium to high resolution satellite images (since 1984), aerial photos (since 1969), GPS surveys (available after 1990) and LIDAR data (2011), which were comparatively analysed by means of GIS techniques. Nowadays, more than half ( 55%) of the Romanian Danube Delta shoreline (disposed in five littoral cells) is affected by erosion. The present coastline configuration is the result of the long-term evolution of this deltaic coast. Depending on the temporal and spatial scales taken into consideration, different driving forces changed the leading role in the dynamics of Danube Delta shoreline in the last 160 years. At centennial time-scale, the threefold decrease of Danube sediment discharge in the last century (especially after 1950, as a result of dams` construction in the Danube watershed) explains the significantly higher shoreline migration rates and area changes between 1856 and 1961/1979 in comparison with the subsequent period, especially along the accumulative sectors. For the Chilia prograding lobe, this resulted in the decrease with more than 75% of the progradation rates and with approximately 90% of the corresponding area change rates, marking its transition, since mid-20th century, from fluvial-dominated morphology to wave-influenced aspect and behaviour. Also, since the beginning of the 20th century, the asymmetric Sf. Gheorghe lobe (the other active lobe of the Danube), experienced dramatic changes of its millennial prograding pattern expressed by the complete cessation of the updrift coastal progradation and the prevalence of erosion in front of the river mouth, whose sediments are feeding far-positioned downdrift depocentres. These changes are reflected by the recent (1930s-present) river mouth dynamics, characterized by cessation of its long-term seaward expansion in favour of downdrift migration, indicating the transition of the Sf. Gheorghe mouth from an asymmetric to a deflected wave-influenced delta morphology. At multi-decadal scale, different modes of climate variability (e.g. North Atlantic Oscillation) control the storminess variations along the Danube Delta coast. Hence, active storminess during 1961-1979 time interval determined very high shoreline dynamics, with two-three times higher shoreline migration rates than afterwards, when a decrease in storminess favoured less dynamic coastlines (on both prograding and erosive sectors). At inter-annual scale, waterline mobility is influenced by storm regime and river floods. Our findings should support the sustainable coastal management and planning, providing a better understanding of past and present coastal processes along the Danube Delta coast.

Leo Szilard Lectureship Award Talk - Universal Scaling Laws from Cells to Cities; A Physicist's Search for Quantitative, Unified Theories of Biological and Social Structure and Dynamics

NASA Astrophysics Data System (ADS)

West, Geoffrey

2013-04-01

Many of the most challenging, exciting and profound questions facing science and society, from the origins of life to global sustainability, fall under the banner of ``complex adaptive systems.'' This talk explores how scaling can be used to begin to develop physics-inspired quantitative, predictive, coarse-grained theories for understanding their structure, dynamics and organization based on underlying mathematisable principles. Remarkably, most physiological, organisational and life history phenomena in biology and socio-economic systems scale in a simple and ``universal'' fashion: metabolic rate scales approximately as the 3/4-power of mass over 27 orders of magnitude from complex molecules to the largest organisms. Time-scales (such as lifespans and growth-rates) and sizes (such as genome lengths and RNA densities) scale with exponents which are typically simple multiples of 1/4, suggesting that fundamental constraints underlie much of the generic structure and dynamics of living systems. These scaling laws follow from dynamical and geometrical properties of space-filling, fractal-like, hierarchical branching networks, presumed optimised by natural selection. This leads to a general framework that potentially captures essential features of diverse systems including vasculature, ontogenetic growth, cancer, aging and mortality, sleep, cell size, and DNA nucleotide substitution rates. Cities and companies also scale: wages, profits, patents, crime, disease, pollution, road lengths scale similarly across the globe, reflecting underlying universal social network dynamics which point to general principles of organization transcending their individuality. These have dramatic implications for global sustainability: innovation and wealth creation that fuel social systems, left unchecked, potentially sow the seeds for their inevitable collapse.

Optimizing hydraulic retention times in denitrifying woodchip bioreactors treating recirculating aquaculture system wastewater

USDA-ARS?s Scientific Manuscript database

The performance of wood-based denitrifying bioreactors to treat high-nitrate wastewaters from aquaculture systems has not previously been demonstrated. Four pilot-scale woodchip bioreactors (approximately 1:10 scale) were constructed and operated for 268 d to determine the optimal range of design hy...

Validity and reliability of the multidimensional assessment of fatigue scale in Iranian patients with relapsing-remitting subtype of multiple sclerosis.

PubMed

Behrangrad, Shabnam; Kordi Yoosefinejad, Amin

2018-03-01

The purpose of this study is to investigate the validity and reliability of the Persian version of the Multidimensional Assessment of Fatigue Scale (MAFS) in an Iranian population with multiple sclerosis. A self-reported survey on fatigue including the MAFS, Fatigue Impact Scale and demographic measures was completed by 130 patients with multiple sclerosis and 60 healthy persons sampled with a convenience method. Test-retest reliability and validity were evaluated 3 days apart. Construct validity of the MAFS was assessed with the Fatigue Impact Scale. The MAFS had high internal consistency (Cronbach's alpha >0.9) and 3-d test-retest reliability (intraclass correlation coefficient = 0.99). Correlation between the Fatigue Impact Scale and MAFS was high (r = 0.99). Correlation between MAFS scores and the Expanded Disability Status Scale was also strong (r = 0.85). Questionnaire items showed acceptable item-scale correlation (0.968-0.993). The Persian version of the MAFS appears to be a valid and reliable questionnaire. It is an appropriate short multidimensional instrument to assess fatigue in patients with multiple sclerosis in clinical practice and research. Implications for Rehabilitation The Persian version of Multidimensional Assessment of Fatigue is a valid and reliable instrument for the assessment and monitoring the fatigue in Persian-language patients with multiple sclerosis. It is very easy to administer and a time efficient scale in comparison to other instruments evaluating fatigue in patients with multiple sclerosis.

Nonlinear Dynamics and Nucleation Kinetics in Near-Critical Liquids

NASA Technical Reports Server (NTRS)

Patashinski, Alexander Z.; Ratner, Mark A.; Pines, Vladimir

1996-01-01

The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.

Low rank approximation methods for MR fingerprinting with large scale dictionaries.

PubMed

Yang, Mingrui; Ma, Dan; Jiang, Yun; Hamilton, Jesse; Seiberlich, Nicole; Griswold, Mark A; McGivney, Debra

2018-04-01

This work proposes new low rank approximation approaches with significant memory savings for large scale MR fingerprinting (MRF) problems. We introduce a compressed MRF with randomized singular value decomposition method to significantly reduce the memory requirement for calculating a low rank approximation of large sized MRF dictionaries. We further relax this requirement by exploiting the structures of MRF dictionaries in the randomized singular value decomposition space and fitting them to low-degree polynomials to generate high resolution MRF parameter maps. In vivo 1.5T and 3T brain scan data are used to validate the approaches. T 1 , T 2 , and off-resonance maps are in good agreement with that of the standard MRF approach. Moreover, the memory savings is up to 1000 times for the MRF-fast imaging with steady-state precession sequence and more than 15 times for the MRF-balanced, steady-state free precession sequence. The proposed compressed MRF with randomized singular value decomposition and dictionary fitting methods are memory efficient low rank approximation methods, which can benefit the usage of MRF in clinical settings. They also have great potentials in large scale MRF problems, such as problems considering multi-component MRF parameters or high resolution in the parameter space. Magn Reson Med 79:2392-2400, 2018. © 2017 International Society for Magnetic Resonance in Medicine. © 2017 International Society for Magnetic Resonance in Medicine.

A rotor optimization using regression analysis

NASA Technical Reports Server (NTRS)

Giansante, N.

1984-01-01

The design and development of helicopter rotors is subject to the many design variables and their interactions that effect rotor operation. Until recently, selection of rotor design variables to achieve specified rotor operational qualities has been a costly, time consuming, repetitive task. For the past several years, Kaman Aerospace Corporation has successfully applied multiple linear regression analysis, coupled with optimization and sensitivity procedures, in the analytical design of rotor systems. It is concluded that approximating equations can be developed rapidly for a multiplicity of objective and constraint functions and optimizations can be performed in a rapid and cost effective manner; the number and/or range of design variables can be increased by expanding the data base and developing approximating functions to reflect the expanded design space; the order of the approximating equations can be expanded easily to improve correlation between analyzer results and the approximating equations; gradients of the approximating equations can be calculated easily and these gradients are smooth functions reducing the risk of numerical problems in the optimization; the use of approximating functions allows the problem to be started easily and rapidly from various initial designs to enhance the probability of finding a global optimum; and the approximating equations are independent of the analysis or optimization codes used.

Multiple scales in metapopulations of public goods producers

NASA Astrophysics Data System (ADS)

Bauer, Marianne; Frey, Erwin

2018-04-01

Multiple scales in metapopulations can give rise to paradoxical behavior: in a conceptual model for a public goods game, the species associated with a fitness cost due to the public good production can be stabilized in the well-mixed limit due to the mere existence of these scales. The scales in this model involve a length scale corresponding to separate patches, coupled by mobility, and separate time scales for reproduction and interaction with a local environment. Contrary to the well-mixed high mobility limit, we find that for low mobilities, the interaction rate progressively stabilizes this species due to stochastic effects, and that the formation of spatial patterns is not crucial for this stabilization.

Structural design using equilibrium programming formulations

NASA Technical Reports Server (NTRS)

Scotti, Stephen J.

1995-01-01

Solutions to increasingly larger structural optimization problems are desired. However, computational resources are strained to meet this need. New methods will be required to solve increasingly larger problems. The present approaches to solving large-scale problems involve approximations for the constraints of structural optimization problems and/or decomposition of the problem into multiple subproblems that can be solved in parallel. An area of game theory, equilibrium programming (also known as noncooperative game theory), can be used to unify these existing approaches from a theoretical point of view (considering the existence and optimality of solutions), and be used as a framework for the development of new methods for solving large-scale optimization problems. Equilibrium programming theory is described, and existing design techniques such as fully stressed design and constraint approximations are shown to fit within its framework. Two new structural design formulations are also derived. The first new formulation is another approximation technique which is a general updating scheme for the sensitivity derivatives of design constraints. The second new formulation uses a substructure-based decomposition of the structure for analysis and sensitivity calculations. Significant computational benefits of the new formulations compared with a conventional method are demonstrated.

RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.

PubMed

Stamatakis, Alexandros

2006-11-01

RAxML-VI-HPC (randomized axelerated maximum likelihood for high performance computing) is a sequential and parallel program for inference of large phylogenies with maximum likelihood (ML). Low-level technical optimizations, a modification of the search algorithm, and the use of the GTR+CAT approximation as replacement for GTR+Gamma yield a program that is between 2.7 and 52 times faster than the previous version of RAxML. A large-scale performance comparison with GARLI, PHYML, IQPNNI and MrBayes on real data containing 1000 up to 6722 taxa shows that RAxML requires at least 5.6 times less main memory and yields better trees in similar times than the best competing program (GARLI) on datasets up to 2500 taxa. On datasets > or =4000 taxa it also runs 2-3 times faster than GARLI. RAxML has been parallelized with MPI to conduct parallel multiple bootstraps and inferences on distinct starting trees. The program has been used to compute ML trees on two of the largest alignments to date containing 25,057 (1463 bp) and 2182 (51,089 bp) taxa, respectively. icwww.epfl.ch/~stamatak

the observation, simulation and evaluation of lake-air interaction process over a high altitude small lake on the Tibetan Plateau

NASA Astrophysics Data System (ADS)

Wang, Binbin; Ma, Yaoming; Ma, Weiqiang; Su, Bob

2017-04-01

Lakes are an important part of the landscape on the Tibetan Plateau. The area that contains most of the plateau lakes has been expanding in recent years, but the impact of lakes on lake-atmosphere energy and water interactions is poorly understood because of a lack of observational data and adequate modeling systems. Furthermore, Precise measurements of evaporation and understanding of the physical controls on turbulent heat flux over lakes at different time scales have fundamental significance for catchment-scale water balance analysis and local-scale climate modeling. To test the performance of lake-air turbulent exchange models over high-altitude lakes and to understanding the driving forces for turbulent heat flux and obtain the actual evaporation over the small high-altitude lakes, an eddy covariance observational system was built above the water surface of the small Nam Co Lake (with an altitude of 4715 m and an area of approximately 1 km2) in April 2012. Firstly, we proposed the proper Charnock coefficient (0.031) and the roughness Reynolds number (0.54) for simulation using turbulent data in 2012, and validated the results using data in 2013 independently; secondly, wind speed shows significance at half-hourly time scales, whereas water vapor and temperature gradients have higher correlations over daily and monthly time scales in lake-air turbulent heat exchange; thirdly, the total evaporation in this small lake (812 mm) is approximately 200 mm larger than that from adjacent Nam Co (approximately 627 mm) during their ice-free seasons. Moreover, the energy stored during April to June is mainly released during September to November, suggesting an energy balance closure value of 0.97 over the entire ice-free season; lastly, 10 evaporation estimation methods are evaluated with the prepared datasets.

Prediction of Autism at 3 Years from Behavioural and Developmental Measures in High-Risk Infants: A Longitudinal Cross-Domain Classifier Analysis

ERIC Educational Resources Information Center

Bussu, G.; Jones, E. J. H.; Charman, T.; Johnson, M. H.; Buitelaar, J. K.; Baron-Cohen, S.; Bedford, R.; Bolton, P.; Blasi, A.; Chandler, S.; Cheung, C.; Davies, K.; Elsabbagh, M.; Fernandes, J.; Gammer, I.; Garwood, H.; Gliga, T.; Guiraud, J.; Hudry, K.; Liew, M.; Lloyd-Fox, S.; Maris, H.; O'Hara, L.; Pasco, G.; Pickles, A.; Ribeiro, H.; Salomone, E.; Tucker, L.; Volein, A.

2018-01-01

We integrated multiple behavioural and developmental measures from multiple time-points using machine learning to improve early prediction of individual Autism Spectrum Disorder (ASD) outcome. We examined Mullen Scales of Early Learning, Vineland Adaptive Behavior Scales, and early ASD symptoms between 8 and 36 months in high-risk siblings (HR; n…

[Multiple time scales analysis of spatial differentiation characteristics of non-point source nitrogen loss within watershed].

PubMed

Liu, Mei-bing; Chen, Xing-wei; Chen, Ying

2015-07-01

Identification of the critical source areas of non-point source pollution is an important means to control the non-point source pollution within the watershed. In order to further reveal the impact of multiple time scales on the spatial differentiation characteristics of non-point source nitrogen loss, a SWAT model of Shanmei Reservoir watershed was developed. Based on the simulation of total nitrogen (TN) loss intensity of all 38 subbasins, spatial distribution characteristics of nitrogen loss and critical source areas were analyzed at three time scales of yearly average, monthly average and rainstorms flood process, respectively. Furthermore, multiple linear correlation analysis was conducted to analyze the contribution of natural environment and anthropogenic disturbance on nitrogen loss. The results showed that there were significant spatial differences of TN loss in Shanmei Reservoir watershed at different time scales, and the spatial differentiation degree of nitrogen loss was in the order of monthly average > yearly average > rainstorms flood process. TN loss load mainly came from upland Taoxi subbasin, which was identified as the critical source area. At different time scales, land use types (such as farmland and forest) were always the dominant factor affecting the spatial distribution of nitrogen loss, while the effect of precipitation and runoff on the nitrogen loss was only taken in no fertilization month and several processes of storm flood at no fertilization date. This was mainly due to the significant spatial variation of land use and fertilization, as well as the low spatial variability of precipitation and runoff.

Orbital and millennial-scale features of atmospheric CH4 over the past 800,000 years.

PubMed

Loulergue, Laetitia; Schilt, Adrian; Spahni, Renato; Masson-Delmotte, Valérie; Blunier, Thomas; Lemieux, Bénédicte; Barnola, Jean-Marc; Raynaud, Dominique; Stocker, Thomas F; Chappellaz, Jérôme

2008-05-15

Atmospheric methane is an important greenhouse gas and a sensitive indicator of climate change and millennial-scale temperature variability. Its concentrations over the past 650,000 years have varied between approximately 350 and approximately 800 parts per 10(9) by volume (p.p.b.v.) during glacial and interglacial periods, respectively. In comparison, present-day methane levels of approximately 1,770 p.p.b.v. have been reported. Insights into the external forcing factors and internal feedbacks controlling atmospheric methane are essential for predicting the methane budget in a warmer world. Here we present a detailed atmospheric methane record from the EPICA Dome C ice core that extends the history of this greenhouse gas to 800,000 yr before present. The average time resolution of the new data is approximately 380 yr and permits the identification of orbital and millennial-scale features. Spectral analyses indicate that the long-term variability in atmospheric methane levels is dominated by approximately 100,000 yr glacial-interglacial cycles up to approximately 400,000 yr ago with an increasing contribution of the precessional component during the four more recent climatic cycles. We suggest that changes in the strength of tropical methane sources and sinks (wetlands, atmospheric oxidation), possibly influenced by changes in monsoon systems and the position of the intertropical convergence zone, controlled the atmospheric methane budget, with an additional source input during major terminations as the retreat of the northern ice sheet allowed higher methane emissions from extending periglacial wetlands. Millennial-scale changes in methane levels identified in our record as being associated with Antarctic isotope maxima events are indicative of ubiquitous millennial-scale temperature variability during the past eight glacial cycles.

Muscles innervated by a single motor neuron exhibit divergent synaptic properties on multiple time scales.

PubMed

Blitz, Dawn M; Pritchard, Amy E; Latimer, John K; Wakefield, Andrew T

2017-04-01

Adaptive changes in the output of neural circuits underlying rhythmic behaviors are relayed to muscles via motor neuron activity. Presynaptic and postsynaptic properties of neuromuscular junctions can impact the transformation from motor neuron activity to muscle response. Further, synaptic plasticity occurring on the time scale of inter-spike intervals can differ between multiple muscles innervated by the same motor neuron. In rhythmic behaviors, motor neuron bursts can elicit additional synaptic plasticity. However, it is unknown whether plasticity regulated by the longer time scale of inter-burst intervals also differs between synapses from the same neuron, and whether any such distinctions occur across a physiological activity range. To address these issues, we measured electrical responses in muscles innervated by a chewing circuit neuron, the lateral gastric (LG) motor neuron, in a well-characterized small motor system, the stomatogastric nervous system (STNS) of the Jonah crab, Cancer borealis In vitro and in vivo , sensory, hormonal and modulatory inputs elicit LG bursting consisting of inter-spike intervals of 50-250 ms and inter-burst intervals of 2-24 s. Muscles expressed similar facilitation measured with paired stimuli except at the shortest inter-spike interval. However, distinct decay time constants resulted in differences in temporal summation. In response to bursting activity, augmentation occurred to different extents and saturated at different inter-burst intervals. Further, augmentation interacted with facilitation, resulting in distinct intra-burst facilitation between muscles. Thus, responses of multiple target muscles diverge across a physiological activity range as a result of distinct synaptic properties sensitive to multiple time scales. © 2017. Published by The Company of Biologists Ltd.

High-speed femtosecond pump-probe spectroscopy with a smart pixel detector array.

PubMed

Bourquin, S; Prasankumar, R P; Kärtner, F X; Fujimoto, J G; Lasser, T; Salathé, R P

2003-09-01

A new femtosecond pump-probe spectroscopy technique is demonstrated that permits the high-speed, parallel acquisition of pump-probe measurements at multiple wavelengths. This is made possible by use of a novel, two-dimensional smart pixel detector array that performs amplitude demodulation in real time on each pixel. This detector array can not only achieve sensitivities comparable with lock-in amplification but also simultaneously performs demodulation of probe transmission signals at multiple wavelengths, thus permitting rapid time- and wavelength-resolved femtosecond pump-probe spectroscopy. Measurements on a thin sample of bulk GaAs are performed across 58 simultaneous wavelengths. Differential probe transmission changes as small as approximately 2 x 10(-4) can be measured over a 5-ps delay scan in only approximately 3 min. This technology can be applied to a wide range of pump-probe measurements in condensed matter, chemistry, and biology.

A method to approximate a closest loadability limit using multiple load flow solutions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yorino, Naoto; Harada, Shigemi; Cheng, Haozhong

A new method is proposed to approximate a closest loadability limit (CLL), or closest saddle node bifurcation point, using a pair of multiple load flow solutions. More strictly, the obtainable points by the method are the stationary points including not only CLL but also farthest and saddle points. An operating solution and a low voltage load flow solution are used to efficiently estimate the node injections at a CLL as well as the left and right eigenvectors corresponding to the zero eigenvalue of the load flow Jacobian. They can be used in monitoring loadability margin, in identification of weak spotsmore » in a power system and in the examination of an optimal control against voltage collapse. Most of the computation time of the proposed method is taken in calculating the load flow solution pair. The remaining computation time is less than that of an ordinary load flow.« less

Generalization of Wilemski-Fixman-Weiss decoupling approximation to the case involving multiple sinks of different sizes, shapes, and reactivities.

PubMed

Uhm, Jesik; Lee, Jinuk; Eun, Changsun; Lee, Sangyoub

2006-08-07

We generalize the Wilemski-Fixman-Weiss decoupling approximation to calculate the transient rate of absorption of point particles into multiple sinks of different sizes, shapes, and reactivities. As an application we consider the case involving two spherical sinks. We obtain a Laplace-transform expression for the transient rate that is in excellent agreement with computer simulations. The long-time steady-state rate has a relatively simple expression, which clearly shows the dependence on the diffusion constant of the particles and on the sizes and reactivities of sinks, and its numerical result is in good agreement with the known exact result that is given in terms of recursion relations.

The quantum-field renormalization group in the problem of a growing phase boundary

DOE Office of Scientific and Technical Information (OSTI.GOV)

Antonov, N.V.; Vasil`ev, A.N.

1995-09-01

Within the quantum-field renormalization-group approach we examine the stochastic equation discussed by S.I. Pavlik in describing a randomly growing phase boundary. We show that, in contrast to Pavlik`s assertion, the model is not multiplicatively renormalizable and that its consistent renormalization-group analysis requires introducing an infinite number of counterterms and the respective coupling constants ({open_quotes}charge{close_quotes}). An explicit calculation in the one-loop approximation shows that a two-dimensional surface of renormalization-group points exits in the infinite-dimensional charge space. If the surface contains an infrared stability region, the problem allows for scaling with the nonuniversal critical dimensionalities of the height of the phase boundarymore » and time, {delta}{sub h} and {delta}{sub t}, which satisfy the exact relationship 2 {delta}{sub h}= {delta}{sub t} + d, where d is the dimensionality of the phase boundary. 23 refs., 1 tab.« less

Online Visualization and Analysis of Global Half-Hourly Infrared Satellite Data

NASA Technical Reports Server (NTRS)

Liu, Zhong; Ostrenga, Dana; Leptoukh, Gregory

2011-01-01

nfrared (IR) images (approximately 11-micron channel) recorded by satellite sensors have been widely used in weather forecasting, research, and classroom education since the Nimbus program. Unlike visible images, IR imagery can reveal cloud features without sunlight illumination; therefore, they can be used to monitor weather phenomena day and night. With geostationary satellites deployed around the globe, it is possible to monitor weather events 24/7 at a temporal resolution that polar-orbiting satellites cannot achieve at the present time. When IR data from multiple geostationary satellites are merged to form a single product--also known as a merged product--it allows for observing weather on a global scale. Its high temporal resolution (e.g., every half hour) also makes it an ideal ancillary dataset for supporting other satellite missions, such as the Tropical Rainfall Measuring Mission (TRMM), etc., by providing additional background information about weather system evolution.

Production of dense plasmas in a hypocycloidal pinch apparatus

NASA Technical Reports Server (NTRS)

Lee, J. H.; Mcfarland, D. R.; Hohl, F.

1977-01-01

A high-power pinch apparatus consisting of disk electrodes was developed, and diagnostic measurements to study its mechanism of dense plasma production have been made. The collapse fronts of the current sheets are well organized, and dense plasma foci are produced on the axis with radial stability in excess of 5 microsec. A plasma density greater than 10 to the 18th power per cu cm is determined with Stark broadening and CO2 laser absorption. Essentially complete absorption of a high-energy CO2 laser beam has been observed. A plasma temperature of approximately 1 keV is measured with differential transmission of soft X-rays through thin foils. The advantages of this apparatus over the coaxial plasma focus are improvements in (1) plasma volume, (2) stability, (3) containment time, (4) access to additional heating by laser or electron beams, and (5) the possibility of scaling up to a multiple array for high-power operation.

Dense plasma focus production in a hypocycloidal pinch

NASA Technical Reports Server (NTRS)

Lee, J. H.; Mcfarland, D. R.; Hohl, F.

1975-01-01

A type of high-power pinch apparatus consisting of disk electrodes was developed, and diagnostic measurements to study its mechanism of dense plasma production were made. The collapse fronts of the current sheets are well organized, and dense plasma focuses are produced on the axis with radial stability in excess of 5 microns. A plasma density greater than 10 to the 18th power/cubic cm was determined with Stark broadening and CO2 laser absorption. A plasma temperature of approximately 1 keV was measured with differential transmission of soft X-rays through thin foils. Essentially complete absorption of a high-energy CO2 laser beam was observed. The advantages of this apparatus over the coaxial plasma focus are in (1) the plasma volume, (2) the stability, (3) the containment time, (4) the easy access to additional heating by laser or electron beams, and (5) the possibility of scaling up to a multiple array for high-power operation.

Bacterial chromosome organization and segregation

PubMed Central

Badrinarayanan, Anjana; Le, Tung BK; Laub, Michael T

2016-01-01

If fully stretched out, a typical bacterial chromosome would be nearly one millimeter long, or approximately 1000 times the length of a cell. Not only must cells massively compact their genetic material, but they must also organize their DNA in a manner that is compatible with a range of cellular processes, including DNA replication, DNA repair, homologous recombination, and horizontal gene transfer. Recent work, driven in part by technological advances, has begun to reveal the general principles of chromosome organization in bacteria. Here, drawing on studies of many different organisms, we review the emerging picture of how bacterial chromosomes are structured at multiple length-scales, highlighting the functions of various DNA-binding proteins and impact of physical forces. Additionally, we discuss the spatial dynamics of chromosomes, particularly during their segregation to daughter cells. Although there has been tremendous progress, we also highlight gaps that remain in understanding chromosome organization and segregation. PMID:26566111

Multiple Streaming and the Probability Distribution of Density in Redshift Space

NASA Astrophysics Data System (ADS)

Hui, Lam; Kofman, Lev; Shandarin, Sergei F.

2000-07-01

We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple streaming using the Zeldovich approximation (ZA), and compute the average number of streams in both real and redshift space. We find that multiple streaming can be significant in redshift space but negligible in real space, even at moderate values of the linear fluctuation amplitude (σl<~1). Moreover, unlike their real-space counterparts, redshift-space multiple streams can flow past each other with minimal interactions. Such nonlinear redshift-space effects, which are physically distinct from the fingers-of-God due to small-scale virialized motions, might in part explain the well-known departure of redshift distortions from the classic linear prediction by Kaiser, even at relatively large scales where the corresponding density field in real space is well described by linear perturbation theory. We also compute, using the ZA, the probability distribution function (PDF) of the density, as well as S3, in real and redshift space, and compare it with the PDF measured from N-body simulations. The role of caustics in defining the character of the high-density tail is examined. We find that (non-Lagrangian) smoothing, due to both finite resolution or discreteness and small-scale velocity dispersions, is very effective in erasing caustic structures, unless the initial power spectrum is sufficiently truncated.

Laboratory studies of scales for measuring helicopter noise

NASA Technical Reports Server (NTRS)

Ollerhead, J. B.

1982-01-01

The adequacy of the effective perceived noise level (EPNL) procedure for rating helicopter noise annoyance was investigated. Recordings of 89 helicopters and 30 fixed wing aircraft (CTOL) flyover sounds were rated with respect to annoyance by groups of approximately 40 subjects. The average annoyance scores were transformed to annoyance levels defined as the equally annoying sound levels of a fixed reference sound. The sound levels of the test sounds were measured on various scales, with and without corrections for duration, tones, and impulsiveness. On average, the helicopter sounds were judged equally annoying to CTOL sounds when their duration corrected levels are approximately 2 dB higher. Multiple regression analysis indicated that, provided the helicopter/CTOL difference of about 2 dB is taken into account, the particular linear combination of level, duration, and tone corrections inherent in EPNL is close to optimum. The results reveal no general requirement for special EPNL correction terms to penalize helicopter sounds which are particularly impulsive; impulsiveness causes spectral and temporal changes which themselves adequately amplify conventionally measured sound levels.

First passage times for a tracer particle in single file diffusion and fractional Brownian motion.

PubMed

Sanders, Lloyd P; Ambjörnsson, Tobias

2012-05-07

We investigate the full functional form of the first passage time density (FPTD) of a tracer particle in a single-file diffusion (SFD) system whose population is: (i) homogeneous, i.e., all particles having the same diffusion constant and (ii) heterogeneous, with diffusion constants drawn from a heavy-tailed power-law distribution. In parallel, the full FPTD for fractional Brownian motion [fBm-defined by the Hurst parameter, H ∈ (0, 1)] is studied, of interest here as fBm and SFD systems belong to the same universality class. Extensive stochastic (non-Markovian) SFD and fBm simulations are performed and compared to two analytical Markovian techniques: the method of images approximation (MIA) and the Willemski-Fixman approximation (WFA). We find that the MIA cannot approximate well any temporal scale of the SFD FPTD. Our exact inversion of the Willemski-Fixman integral equation captures the long-time power-law exponent, when H ≥ 1/3, as predicted by Molchan [Commun. Math. Phys. 205, 97 (1999)] for fBm. When H < 1/3, which includes homogeneous SFD (H = 1/4), and heterogeneous SFD (H < 1/4), the WFA fails to agree with any temporal scale of the simulations and Molchan's long-time result. SFD systems are compared to their fBm counter parts; and in the homogeneous system both scaled FPTDs agree on all temporal scales including also, the result by Molchan, thus affirming that SFD and fBm dynamics belong to the same universality class. In the heterogeneous case SFD and fBm results for heterogeneity-averaged FPTDs agree in the asymptotic time limit. The non-averaged heterogeneous SFD systems display a lack of self-averaging. An exponential with a power-law argument, multiplied by a power-law pre-factor is shown to describe well the FPTD for all times for homogeneous SFD and sub-diffusive fBm systems.

Scenario-based modeling for multiple allocation hub location problem under disruption risk: multiple cuts Benders decomposition approach

NASA Astrophysics Data System (ADS)

Yahyaei, Mohsen; Bashiri, Mahdi

2017-12-01

The hub location problem arises in a variety of domains such as transportation and telecommunication systems. In many real-world situations, hub facilities are subject to disruption. This paper deals with the multiple allocation hub location problem in the presence of facilities failure. To model the problem, a two-stage stochastic formulation is developed. In the proposed model, the number of scenarios grows exponentially with the number of facilities. To alleviate this issue, two approaches are applied simultaneously. The first approach is to apply sample average approximation to approximate the two stochastic problem via sampling. Then, by applying the multiple cuts Benders decomposition approach, computational performance is enhanced. Numerical studies show the effective performance of the SAA in terms of optimality gap for small problem instances with numerous scenarios. Moreover, performance of multi-cut Benders decomposition is assessed through comparison with the classic version and the computational results reveal the superiority of the multi-cut approach regarding the computational time and number of iterations.

Sound velocity and absorption in a coarsening foam.

PubMed

Mujica, Nicolás; Fauve, Stéphan

2002-08-01

We present experimental measurements of sound velocity and absorption in a commercial shaving foam. We observe that both quantities evolve with time as the foam coarsens increasing its mean bubble radius . By varying the acoustic frequency we probe the foam from the large wavelength regime, lambda approximately 1500, down to the scale lambda approximately 20. Sound absorption alpha varies significantly with both the foam age and the excitation frequency. After an initial transition time of 20 min, the attenuation per wavelength, alphalambda, varies linearly with the foam age. In addition, for evolution times smaller than approximately 90 min, we observe that alphalambda scales linearly with both foam age and frequency. From these scalings we show that the thermal dissipation mechanism is the dominant one. Sound velocity c is initially frequency independent but the medium becomes slightly dispersive as the foam coarsens. We observe that sound velocity depends on the evolution of the structure of the foam, even in the large wavelength regime. After 2 h of foam coarsening, c decreases at least by a factor of 20%, due to the softening of the foam. These facts are explained by considering the liquid matrix elasticity, due to the presence of surfactant molecules. A simple model of foam structure, combined with results of Biot's theory for porous media, gives both good qualitative and quantitative agreement with our experimental results in the low frequency regime.

Chaotic dynamics of Comet 1P/Halley: Lyapunov exponent and survival time expectancy

NASA Astrophysics Data System (ADS)

Muñoz-Gutiérrez, M. A.; Reyes-Ruiz, M.; Pichardo, B.

2015-03-01

The orbital elements of Comet Halley are known to a very high precision, suggesting that the calculation of its future dynamical evolution is straightforward. In this paper we seek to characterize the chaotic nature of the present day orbit of Comet Halley and to quantify the time-scale over which its motion can be predicted confidently. In addition, we attempt to determine the time-scale over which its present day orbit will remain stable. Numerical simulations of the dynamics of test particles in orbits similar to that of Comet Halley are carried out with the MERCURY 6.2 code. On the basis of these we construct survival time maps to assess the absolute stability of Halley's orbit, frequency analysis maps to study the variability of the orbit, and we calculate the Lyapunov exponent for the orbit for variations in initial conditions at the level of the present day uncertainties in our knowledge of its orbital parameters. On the basis of our calculations of the Lyapunov exponent for Comet Halley, the chaotic nature of its motion is demonstrated. The e-folding time-scale for the divergence of initially very similar orbits is approximately 70 yr. The sensitivity of the dynamics on initial conditions is also evident in the self-similarity character of the survival time and frequency analysis maps in the vicinity of Halley's orbit, which indicates that, on average, it is unstable on a time-scale of hundreds of thousands of years. The chaotic nature of Halley's present day orbit implies that a precise determination of its motion, at the level of the present-day observational uncertainty, is difficult to predict on a time-scale of approximately 100 yr. Furthermore, we also find that the ejection of Halley from the Solar system or its collision with another body could occur on a time-scale as short as 10 000 yr.

Models of inertial range spectra of interplanetary magnetohydrodynamic turbulence

NASA Technical Reports Server (NTRS)

Zhou, YE; Matthaeus, William H.

1990-01-01

A framework based on turbulence theory is presented to develop approximations for the local turbulence effects that are required in transport models. An approach based on Kolmogoroff-style dimensional analysis is presented as well as one based on a wave-number diffusion picture. Particular attention is given to the case of MHD turbulence with arbitrary cross helicity and with arbitrary ratios of the Alfven time scale and the nonlinear time scale.

Measurement of Thunderstorm Cloud-Top Parameters Using High-Frequency Satellite Imagery

DTIC Science & Technology

1978-01-01

short wave was present well to the south of this system approximately 2000 ka west of Baja California. Two distinct flow patterns were present, one...view can be observed in near real time whereas radar observations, although excellent for local purposes, involve substantial errors when composited...on a large scale. The time delay in such large scale compositing is critical when attempting to monitor convective cloud systems for a potential

The relationship between challenging behaviours, mood and interest/pleasure in adults with severe and profound intellectual disabilities.

PubMed

Bernstein, A M; Visconti, K J; Csorba, J; Radvanyi, K; Rojahn, J

2015-11-01

We investigated whether current mood and interest/pleasure ratings in adults with moderate to profound intellectual disabilities were predictive of challenging behaviour [self-injurious behaviour (SIB), aggressive/destructive behaviour and stereotypic behaviour] and vice versa. In this combined cross-sectional and longitudinal study, staff members of a Hungarian residential facility completed translated versions of the Behaviour Problems Inventory-Short Form (BPI-S), the Challenging Behaviour Interview (CBI) and the Mood, Interest and Pleasure Questionnaire-Short Form (MIPQ-S) for 50 participants at two time points, approximately 4 to 5 months apart. Bivariate correlations from data concurrently assessed at Time-1 showed significant linear relationships between the SIB (both frequency and severity scores) and Interest/Pleasure sub-scales, and the Aggressive/Destructive Behaviour (severity scores) and the MIPQ-S Mood sub-scales (unadjusted for multiple correlations). All of these effects were found with the BPI-S data, but not with the CBI. Multiple regression analyses revealed that (1) low interest/pleasure assessed at Time-1 predicted high SIB (frequency and severity) at Time-2. (2) Interest/pleasure was not predictive of aggressive or stereotypic behaviour. (3) Mood at Time-1 did not predict any of the three types of behaviour problems at Time-2. (4) In reverse, high SIB (frequency and severity) at Time-1 predicted low interest/pleasure ratings at Time-2. (5) Surprisingly, frequent aggressive/destructive behaviour predicted high interest/pleasure. (6) Stereotypic behaviour scores at Time-1 did not predict interest/pleasure ratings at Time-2. Again, all of these effects were only found with the BPI-S data, but not with the CBI. Internal consistency, test-retest reliability and concurrent validity of the Hungarian versions of all three questionnaires had generally satisfactory outcomes. The fact that increasingly frequent and severe SIB was predicted by declining measures of interest/pleasure is consistent with previous studies. Contrary to those earlier studies, however, we found that SIB was not predicted by mood and that aggressive/destructive behaviour actually predicted future elevated mood. Implications for future research regarding the directional relationship between affective states such as mood and interest and pleasure, on the one hand, and challenging behaviour, on the other, were discussed. © 2015 MENCAP and International Association of the Scientific Study of Intellectual and Developmental Disabilities and John Wiley & Sons Ltd.

Centrality and energy dependence of charged-particle multiplicities in heavy ion collisions in the context of elementary reactions

NASA Astrophysics Data System (ADS)

Back, B. B.; Baker, M. D.; Barton, D. S.; Betts, R. R.; Ballintijn, M.; Bickley, A. A.; Bindel, R.; Budzanowski, A.; Busza, W.; Carroll, A.; Decowski, M. P.; García, E.; George, N.; Gulbrandsen, K.; Gushue, S.; Halliwell, C.; Hamblen, J.; Heintzelman, G. A.; Henderson, C.; Hofman, D. J.; Hollis, R. S.; Hołyński, R.; Holzman, B.; Iordanova, A.; Johnson, E.; Kane, J. L.; Katzy, J.; Khan, N.; Kucewicz, W.; Kulinich, P.; Kuo, C. M.; Lin, W. T.; Manly, S.; McLeod, D.; Michałowski, J.; Mignerey, A. C.; Nouicer, R.; Olszewski, A.; Pak, R.; Park, I. C.; Pernegger, H.; Reed, C.; Remsberg, L. P.; Reuter, M.; Roland, C.; Roland, G.; Rosenberg, L.; Sagerer, J.; Sarin, P.; Sawicki, P.; Skulski, W.; Steadman, S. G.; Steinberg, P.; Stephans, G. S. F.; Stodulski, M.; Sukhanov, A.; Tang, J.-L.; Teng, R.; Trzupek, A.; Vale, C.; Nieuwenhuizen, G. J. Van; Verdier, R.; Veres, G. I.; Wadsworth, B.; Wolfs, F. L. H.; Wosiek, B.; Woźniak, K.; Wuosmaa, A. H.; Wysłouch, B.

2006-08-01

The PHOBOS experiment at the BNL Relativistic Heavy Ion Collider has measured the total multiplicity of primary charged particles as a function of collision centrality in Au+Au collisions at sNN= 19.6, 130, and 200 GeV. An approximate independence of / on the number of participating nucleons is observed, reminiscent of “wounded nucleon” scaling (Nch∝Npart) observed in proton-nucleus collisions. Unlike p+A, the constant of proportionality does not seem to be set by the pp/p¯p data at the same energy. Rather, there seems to be a surprising correspondence with the total multiplicity measured in e+e- annihilations, as well as the rapidity shape measured over a large range. The energy dependence of the integrated multiplicity per participant pair shows that e+e- and A+A data agree over a large range of center-of-mass energies (s>20 GeV), and pp/p¯p data can be brought to agree approximately with the e+e- data by correcting for the typical energy taken away by leading particles. This is suggestive of a mechanism for soft particle production that depends mainly on the amount of available energy. It is conjectured that the dominant distinction between A+A and p+p collisions is the multiple collisions per participant, which appears to be sufficient to substantially reduce the energy taken away by leading particles.

Fractal structure of the interplanetary magnetic field

NASA Technical Reports Server (NTRS)

Burlaga, L. F.; Klein, L. W.

1985-01-01

Under some conditions, time series of the interplanetary magnetic field strength and components have the properties of fractal curves. Magnetic field measurements made near 8.5 AU by Voyager 2 from June 5 to August 24, 1981 were self-similar over time scales from approximately 20 sec to approximately 3 x 100,000 sec, and the fractal dimension of the time series of the strength and components of the magnetic field was D = 5/3, corresponding to a power spectrum P(f) approximately f sup -5/3. Since the Kolmogorov spectrum for homogeneous, isotropic, stationary turbulence is also f sup -5/3, the Voyager 2 measurements are consistent with the observation of an inertial range of turbulence extending over approximately four decades in frequency. Interaction regions probably contributed most of the power in this interval. As an example, one interaction region is discussed in which the magnetic field had a fractal dimension D = 5/3.

Single-grain 40Ar-39Ar ages of glauconies: implications for the geologic time scale and global sea level variations

PubMed

Smith; Evensen; York; Odin

1998-03-06

The mineral series glaucony supplies 40% of the absolute-age database for the geologic time scale of the last 250 million years. However, glauconies have long been suspected of giving young potassium-argon ages on bulk samples. Laser-probe argon-argon dating shows that glaucony populations comprise grains with a wide range of ages, suggesting a period of genesis several times longer ( approximately 5 million years) than previously thought. An estimate of the age of their enclosing sediments (and therefore of time scale boundaries) is given by the oldest nonrelict grains in the glaucony populations, whereas the formation times of the younger grains appear to be modulated by global sea level.

Methods and optical fibers that decrease pulse degradation resulting from random chromatic dispersion

DOEpatents

Chertkov, Michael; Gabitov, Ildar

2004-03-02

The present invention provides methods and optical fibers for periodically pinning an actual (random) accumulated chromatic dispersion of an optical fiber to a predicted accumulated dispersion of the fiber through relatively simple modifications of fiber-optic manufacturing methods or retrofitting of existing fibers. If the pinning occurs with sufficient frequency (at a distance less than or are equal to a correlation scale), pulse degradation resulting from random chromatic dispersion is minimized. Alternatively, pinning may occur quasi-periodically, i.e., the pinning distance is distributed between approximately zero and approximately two to three times the correlation scale.

GAPPARD: a computationally efficient method of approximating gap-scale disturbance in vegetation models

NASA Astrophysics Data System (ADS)

Scherstjanoi, M.; Kaplan, J. O.; Thürig, E.; Lischke, H.

2013-02-01

Models of vegetation dynamics that are designed for application at spatial scales larger than individual forest gaps suffer from several limitations. Typically, either a population average approximation is used that results in unrealistic tree allometry and forest stand structure, or models have a high computational demand because they need to simulate both a series of age-based cohorts and a number of replicate patches to account for stochastic gap-scale disturbances. The detail required by the latter method increases the number of calculations by two to three orders of magnitude compared to the less realistic population average approach. In an effort to increase the efficiency of dynamic vegetation models without sacrificing realism, and to explore patterns of spatial scaling in forests, we developed a new method for simulating stand-replacing disturbances that is both accurate and 10-50x faster than approaches that use replicate patches. The GAPPARD (approximating GAP model results with a Probabilistic Approach to account for stand Replacing Disturbances) method works by postprocessing the output of deterministic, undisturbed simulations of a cohort-based vegetation model by deriving the distribution of patch ages at any point in time on the basis of a disturbance probability. With this distribution, the expected value of any output variable can be calculated from the output values of the deterministic undisturbed run at the time corresponding to the patch age. To account for temporal changes in model forcing, e.g., as a result of climate change, GAPPARD performs a series of deterministic simulations and interpolates between the results in the postprocessing step. We integrated the GAPPARD method in the forest models LPJ-GUESS and TreeM-LPJ, and evaluated these in a series of simulations along an altitudinal transect of an inner-alpine valley. With GAPPARD applied to LPJ-GUESS results were insignificantly different from the output of the original model LPJ-GUESS using 100 replicate patches, but simulation time was reduced by approximately the factor 10. Our new method is therefore highly suited rapidly approximating LPJ-GUESS results, and provides the opportunity for future studies over large spatial domains, allows easier parameterization of tree species, faster identification of areas of interesting simulation results, and comparisons with large-scale datasets and forest models.

A neuromorphic implementation of multiple spike-timing synaptic plasticity rules for large-scale neural networks

PubMed Central

Wang, Runchun M.; Hamilton, Tara J.; Tapson, Jonathan C.; van Schaik, André

2015-01-01

We present a neuromorphic implementation of multiple synaptic plasticity learning rules, which include both Spike Timing Dependent Plasticity (STDP) and Spike Timing Dependent Delay Plasticity (STDDP). We present a fully digital implementation as well as a mixed-signal implementation, both of which use a novel dynamic-assignment time-multiplexing approach and support up to 226 (64M) synaptic plasticity elements. Rather than implementing dedicated synapses for particular types of synaptic plasticity, we implemented a more generic synaptic plasticity adaptor array that is separate from the neurons in the neural network. Each adaptor performs synaptic plasticity according to the arrival times of the pre- and post-synaptic spikes assigned to it, and sends out a weighted or delayed pre-synaptic spike to the post-synaptic neuron in the neural network. This strategy provides great flexibility for building complex large-scale neural networks, as a neural network can be configured for multiple synaptic plasticity rules without changing its structure. We validate the proposed neuromorphic implementations with measurement results and illustrate that the circuits are capable of performing both STDP and STDDP. We argue that it is practical to scale the work presented here up to 236 (64G) synaptic adaptors on a current high-end FPGA platform. PMID:26041985

GAPPARD: a computationally efficient method of approximating gap-scale disturbance in vegetation models

NASA Astrophysics Data System (ADS)

Scherstjanoi, M.; Kaplan, J. O.; Thürig, E.; Lischke, H.

2013-09-01

Models of vegetation dynamics that are designed for application at spatial scales larger than individual forest gaps suffer from several limitations. Typically, either a population average approximation is used that results in unrealistic tree allometry and forest stand structure, or models have a high computational demand because they need to simulate both a series of age-based cohorts and a number of replicate patches to account for stochastic gap-scale disturbances. The detail required by the latter method increases the number of calculations by two to three orders of magnitude compared to the less realistic population average approach. In an effort to increase the efficiency of dynamic vegetation models without sacrificing realism, we developed a new method for simulating stand-replacing disturbances that is both accurate and faster than approaches that use replicate patches. The GAPPARD (approximating GAP model results with a Probabilistic Approach to account for stand Replacing Disturbances) method works by postprocessing the output of deterministic, undisturbed simulations of a cohort-based vegetation model by deriving the distribution of patch ages at any point in time on the basis of a disturbance probability. With this distribution, the expected value of any output variable can be calculated from the output values of the deterministic undisturbed run at the time corresponding to the patch age. To account for temporal changes in model forcing (e.g., as a result of climate change), GAPPARD performs a series of deterministic simulations and interpolates between the results in the postprocessing step. We integrated the GAPPARD method in the vegetation model LPJ-GUESS, and evaluated it in a series of simulations along an altitudinal transect of an inner-Alpine valley. We obtained results very similar to the output of the original LPJ-GUESS model that uses 100 replicate patches, but simulation time was reduced by approximately the factor 10. Our new method is therefore highly suited for rapidly approximating LPJ-GUESS results, and provides the opportunity for future studies over large spatial domains, allows easier parameterization of tree species, faster identification of areas of interesting simulation results, and comparisons with large-scale datasets and results of other forest models.

Large-scale phenomena, chapter 3, part D

NASA Technical Reports Server (NTRS)

1975-01-01

Oceanic phenomena with horizontal scales from approximately 100 km up to the widths of the oceans themselves are examined. Data include: shape of geoid, quasi-stationary anomalies due to spatial variations in sea density and steady current systems, and the time dependent variations due to tidal and meteorological forces and to varying currents.

Accuracy of the adiabatic-impulse approximation for closed and open quantum systems

NASA Astrophysics Data System (ADS)

Tomka, Michael; Campos Venuti, Lorenzo; Zanardi, Paolo

2018-03-01

We study the adiabatic-impulse approximation (AIA) as a tool to approximate the time evolution of quantum states when driven through a region of small gap. Such small-gap regions are a common situation in adiabatic quantum computing and having reliable approximations is important in this context. The AIA originates from the Kibble-Zurek theory applied to continuous quantum phase transitions. The Kibble-Zurek mechanism was developed to predict the power-law scaling of the defect density across a continuous quantum phase transition. Instead, here we quantify the accuracy of the AIA via the trace norm distance with respect to the exact evolved state. As expected, we find that for short times or fast protocols, the AIA outperforms the simple adiabatic approximation. However, for large times or slow protocols, the situation is actually reversed and the AIA provides a worse approximation. Nevertheless, we found a variation of the AIA that can perform better than the adiabatic one. This counterintuitive modification consists in crossing the region of small gap twice. Our findings are illustrated by several examples of driven closed and open quantum systems.

Black holes from large N singlet models

NASA Astrophysics Data System (ADS)

Amado, Irene; Sundborg, Bo; Thorlacius, Larus; Wintergerst, Nico

2018-03-01

The emergent nature of spacetime geometry and black holes can be directly probed in simple holographic duals of higher spin gravity and tensionless string theory. To this end, we study time dependent thermal correlation functions of gauge invariant observables in suitably chosen free large N gauge theories. At low temperature and on short time scales the correlation functions encode propagation through an approximate AdS spacetime while interesting departures emerge at high temperature and on longer time scales. This includes the existence of evanescent modes and the exponential decay of time dependent boundary correlations, both of which are well known indicators of bulk black holes in AdS/CFT. In addition, a new time scale emerges after which the correlation functions return to a bulk thermal AdS form up to an overall temperature dependent normalization. A corresponding length scale was seen in equal time correlation functions in the same models in our earlier work.

1SXPS: A Deep Swift X-Ray Telescope Point Source Catalog with Light Curves and Spectra

NASA Technical Reports Server (NTRS)

Evans, P. A.; Osborne, J. P.; Beardmore, A. P.; Page, K. L.; Willingale, R.; Mountford, C. J.; Pagani, C.; Burrows, D. N.; Kennea, J. A.; Perri, M.;

The Long-term Light Curves of X-ray Binaries Contain Simultaneous Periodic and Random Components

NASA Technical Reports Server (NTRS)

White, Nicholas E. (Technical Monitor); Boyd, Patricia T.; Smale, Alan P.

2002-01-01

LMC X-3 and Cyg X-2 show large amplitude X-ray fluctuations that have been attributed to a warped accretion disk. Cyg X-3 displays high amplitude, apparently non-periodic oscillations. We reanalyze these systems using RXTE ASM data and time-frequency decomposition techniques. We find that the long-term variations in Cyg X-2 can be completely characterized by excursions whose durations are integer multiples of the orbital period, including one essentially identical to the reported "period" of 78 days. Cyg X-3 can be characterized in terms of integer multiples of a 71-day fundamental period unrelated to the 4.8 day orbital period, but suggestively close to the approximately equal to greater than 60 day reported precession period of the relativistic jet inferred from recent radio observations. The long-term excursions of LMC X-3 are related to each other by rational fractions, suggesting the characteristic time scale is 10.594 days, shorter than any observed excursion to date. We explore the phase space evolution of the light curves using a natural embedding and find that all three systems possess two rotation centers that organize the phase space trajectories, one of low luminosity and the other of high luminosity. The implications of this repeatable behavior on generic models of accretion disk dynamics and mass transfer variability are explored.

Dynamical tuning for MPC using population games: A water supply network application.

PubMed

Barreiro-Gomez, Julian; Ocampo-Martinez, Carlos; Quijano, Nicanor

2017-07-01

Model predictive control (MPC) is a suitable strategy for the control of large-scale systems that have multiple design requirements, e.g., multiple physical and operational constraints. Besides, an MPC controller is able to deal with multiple control objectives considering them within the cost function, which implies to determine a proper prioritization for each of the objectives. Furthermore, when the system has time-varying parameters and/or disturbances, the appropriate prioritization might vary along the time as well. This situation leads to the need of a dynamical tuning methodology. This paper addresses the dynamical tuning issue by using evolutionary game theory. The advantages of the proposed method are highlighted and tested over a large-scale water supply network with periodic time-varying disturbances. Finally, results are analyzed with respect to a multi-objective MPC controller that uses static tuning. Copyright © 2017 ISA. Published by Elsevier Ltd. All rights reserved.

The effects of vent location, event scale and time forecasts on pyroclastic density current hazard maps at Campi Flegrei caldera (Italy)

NASA Astrophysics Data System (ADS)

Bevilacqua, Andrea; Neri, Augusto; Bisson, Marina; Esposti Ongaro, Tomaso; Flandoli, Franco; Isaia, Roberto; Rosi, Mauro; Vitale, Stefano

2017-09-01

This study presents a new method for producing long-term hazard maps for pyroclastic density currents (PDC) originating at Campi Flegrei caldera. Such method is based on a doubly stochastic approach and is able to combine the uncertainty assessments on the spatial location of the volcanic vent, the size of the flow and the expected time of such an event. The results are obtained by using a Monte Carlo approach and adopting a simplified invasion model based on the box model integral approximation. Temporal assessments are modelled through a Cox-type process including self-excitement effects, based on the eruptive record of the last 15 kyr. Mean and percentile maps of PDC invasion probability are produced, exploring their sensitivity to some sources of uncertainty and to the effects of the dependence between PDC scales and the caldera sector where they originated. Conditional maps representative of PDC originating inside limited zones of the caldera, or of PDC with a limited range of scales are also produced. Finally, the effect of assuming different time windows for the hazard estimates is explored, also including the potential occurrence of a sequence of multiple events. Assuming that the last eruption of Monte Nuovo (A.D. 1538) marked the beginning of a new epoch of activity similar to the previous ones, results of the statistical analysis indicate a mean probability of PDC invasion above 5% in the next 50 years on almost the entire caldera (with a probability peak of 25% in the central part of the caldera). In contrast, probability values reduce by a factor of about 3 if the entire eruptive record is considered over the last 15 kyr, i.e. including both eruptive epochs and quiescent periods.

Pattern analysis of total item score and item response of the Kessler Screening Scale for Psychological Distress (K6) in a nationally representative sample of US adults

PubMed Central

Kawasaki, Yohei; Ide, Kazuki; Akutagawa, Maiko; Yamada, Hiroshi; Yutaka, Ono; Furukawa, Toshiaki A.

2017-01-01

Background Several recent studies have shown that total scores on depressive symptom measures in a general population approximate an exponential pattern except for the lower end of the distribution. Furthermore, we confirmed that the exponential pattern is present for the individual item responses on the Center for Epidemiologic Studies Depression Scale (CES-D). To confirm the reproducibility of such findings, we investigated the total score distribution and item responses of the Kessler Screening Scale for Psychological Distress (K6) in a nationally representative study. Methods Data were drawn from the National Survey of Midlife Development in the United States (MIDUS), which comprises four subsamples: (1) a national random digit dialing (RDD) sample, (2) oversamples from five metropolitan areas, (3) siblings of individuals from the RDD sample, and (4) a national RDD sample of twin pairs. K6 items are scored using a 5-point scale: “none of the time,” “a little of the time,” “some of the time,” “most of the time,” and “all of the time.” The pattern of total score distribution and item responses were analyzed using graphical analysis and exponential regression model. Results The total score distributions of the four subsamples exhibited an exponential pattern with similar rate parameters. The item responses of the K6 approximated a linear pattern from “a little of the time” to “all of the time” on log-normal scales, while “none of the time” response was not related to this exponential pattern. Discussion The total score distribution and item responses of the K6 showed exponential patterns, consistent with other depressive symptom scales. PMID:28289560

A simple approximation for larval retention around reefs

NASA Astrophysics Data System (ADS)

Cetina-Heredia, Paulina; Connolly, Sean R.

2011-09-01

Estimating larval retention at individual reefs by local scale three-dimensional flows is a significant problem for understanding, and predicting, larval dispersal. Determining larval dispersal commonly involves the use of computationally demanding and expensively calibrated/validated hydrodynamic models that resolve reef wake eddies. This study models variation in larval retention times for a range of reef shapes and circulation regimes, using a reef-scale three-dimensional hydrodynamic model. It also explores how well larval retention time can be estimated based on the "Island Wake Parameter", a measure of the degree of flow turbulence in the wake of reefs that is a simple function of flow speed, reef dimension, and vertical diffusion. The mean residence times found in the present study (0.48-5.64 days) indicate substantial potential for self-recruitment of species whose larvae are passive, or weak swimmers, for the first several days after release. Results also reveal strong and significant relationships between the Island Wake Parameter and mean residence time, explaining 81-92% of the variability in retention among reefs across a range of unidirectional flow speeds and tidal regimes. These findings suggest that good estimates of larval retention may be obtained from relatively coarse-scale characteristics of the flow, and basic features of reef geomorphology. Such approximations may be a valuable tool for modeling connectivity and meta-population dynamics over large spatial scales, where explicitly characterizing fine-scale flows around reef requires a prohibitive amount of computation and extensive model calibration.

Local Approximation and Hierarchical Methods for Stochastic Optimization

NASA Astrophysics Data System (ADS)

Cheng, Bolong

In this thesis, we present local and hierarchical approximation methods for two classes of stochastic optimization problems: optimal learning and Markov decision processes. For the optimal learning problem class, we introduce a locally linear model with radial basis function for estimating the posterior mean of the unknown objective function. The method uses a compact representation of the function which avoids storing the entire history, as is typically required by nonparametric methods. We derive a knowledge gradient policy with the locally parametric model, which maximizes the expected value of information. We show the policy is asymptotically optimal in theory, and experimental works suggests that the method can reliably find the optimal solution on a range of test functions. For the Markov decision processes problem class, we are motivated by an application where we want to co-optimize a battery for multiple revenue, in particular energy arbitrage and frequency regulation. The nature of this problem requires the battery to make charging and discharging decisions at different time scales while accounting for the stochastic information such as load demand, electricity prices, and regulation signals. Computing the exact optimal policy becomes intractable due to the large state space and the number of time steps. We propose two methods to circumvent the computation bottleneck. First, we propose a nested MDP model that structure the co-optimization problem into smaller sub-problems with reduced state space. This new model allows us to understand how the battery behaves down to the two-second dynamics (that of the frequency regulation market). Second, we introduce a low-rank value function approximation for backward dynamic programming. This new method only requires computing the exact value function for a small subset of the state space and approximate the entire value function via low-rank matrix completion. We test these methods on historical price data from the PJM Interconnect and show that it outperforms the baseline approach used in the industry.

On the Role of Multi-Scale Processes in CO2 Storage Security and Integrity

NASA Astrophysics Data System (ADS)

Pruess, K.; Kneafsey, T. J.

2008-12-01

Consideration of multiple scales in subsurface processes is usually referred to the spatial domain, where we may attempt to relate process descriptions and parameters from pore and bench (Darcy) scale to much larger field and regional scales. However, multiple scales occur also in the time domain, and processes extending over a broad range of time scales may be very relevant to CO2 storage and containment. In some cases, such as in the convective instability induced by CO2 dissolution in saline waters, space and time scales are coupled in the sense that perturbations induced by CO2 injection will grow concurrently over many orders of magnitude in both space and time. In other cases, CO2 injection may induce processes that occur on short time scales, yet may affect large regions. Possible examples include seismicity that may be triggered by CO2 injection, or hypothetical release events such as "pneumatic eruptions" that may discharge substantial amounts of CO2 over a short time period. This paper will present recent advances in our experimental and modeling studies of multi-scale processes. Specific examples that will be discussed include (1) the process of CO2 dissolution-diffusion-convection (DDC), that can greatly accelerate the rate at which free-phase CO2 is stored as aqueous solute; (2) self- enhancing and self-limiting processes during CO2 leakage through faults, fractures, or improperly abandoned wells; and (3) porosity and permeability reduction from salt precipitation near CO2 injection wells, and mitigation of corresponding injectivity loss. This work was supported by the Office of Basic Energy Sciences and by the Zero Emission Research and Technology project (ZERT) under Contract No. DE-AC02-05CH11231 with the U.S. Department of Energy.

Gradient design for liquid chromatography using multi-scale optimization.

PubMed

López-Ureña, S; Torres-Lapasió, J R; Donat, R; García-Alvarez-Coque, M C

2018-01-26

In reversed phase-liquid chromatography, the usual solution to the "general elution problem" is the application of gradient elution with programmed changes of organic solvent (or other properties). A correct quantification of chromatographic peaks in liquid chromatography requires well resolved signals in a proper analysis time. When the complexity of the sample is high, the gradient program should be accommodated to the local resolution needs of each analyte. This makes the optimization of such situations rather troublesome, since enhancing the resolution for a given analyte may imply a collateral worsening of the resolution of other analytes. The aim of this work is to design multi-linear gradients that maximize the resolution, while fulfilling some restrictions: all peaks should be eluted before a given maximal time, the gradient should be flat or increasing, and sudden changes close to eluting peaks are penalized. Consequently, an equilibrated baseline resolution for all compounds is sought. This goal is achieved by splitting the optimization problem in a multi-scale framework. In each scale κ, an optimization problem is solved with N κ  ≈ 2 κ variables that are used to build the gradients. The N κ variables define cubic splines written in terms of a B-spline basis. This allows expressing gradients as polygonals of M points approximating the splines. The cubic splines are built using subdivision schemes, a technique of fast generation of smooth curves, compatible with the multi-scale framework. Owing to the nature of the problem and the presence of multiple local maxima, the algorithm used in the optimization problem of each scale κ should be "global", such as the pattern-search algorithm. The multi-scale optimization approach is successfully applied to find the best multi-linear gradient for resolving a mixture of amino acid derivatives. Copyright © 2017 Elsevier B.V. All rights reserved.

Statistical significance approximation in local trend analysis of high-throughput time-series data using the theory of Markov chains.

PubMed

Xia, Li C; Ai, Dongmei; Cram, Jacob A; Liang, Xiaoyi; Fuhrman, Jed A; Sun, Fengzhu

2015-09-21

Local trend (i.e. shape) analysis of time series data reveals co-changing patterns in dynamics of biological systems. However, slow permutation procedures to evaluate the statistical significance of local trend scores have limited its applications to high-throughput time series data analysis, e.g., data from the next generation sequencing technology based studies. By extending the theories for the tail probability of the range of sum of Markovian random variables, we propose formulae for approximating the statistical significance of local trend scores. Using simulations and real data, we show that the approximate p-value is close to that obtained using a large number of permutations (starting at time points >20 with no delay and >30 with delay of at most three time steps) in that the non-zero decimals of the p-values obtained by the approximation and the permutations are mostly the same when the approximate p-value is less than 0.05. In addition, the approximate p-value is slightly larger than that based on permutations making hypothesis testing based on the approximate p-value conservative. The approximation enables efficient calculation of p-values for pairwise local trend analysis, making large scale all-versus-all comparisons possible. We also propose a hybrid approach by integrating the approximation and permutations to obtain accurate p-values for significantly associated pairs. We further demonstrate its use with the analysis of the Polymouth Marine Laboratory (PML) microbial community time series from high-throughput sequencing data and found interesting organism co-occurrence dynamic patterns. The software tool is integrated into the eLSA software package that now provides accelerated local trend and similarity analysis pipelines for time series data. The package is freely available from the eLSA website: http://bitbucket.org/charade/elsa.

On the use of variability time-scales as an early classifier of radio transients and variables

NASA Astrophysics Data System (ADS)

Pietka, M.; Staley, T. D.; Pretorius, M. L.; Fender, R. P.

2017-11-01

We have shown previously that a broad correlation between the peak radio luminosity and the variability time-scales, approximately L ∝ τ5, exists for variable synchrotron emitting sources and that different classes of astrophysical sources occupy different regions of luminosity and time-scale space. Based on those results, we investigate whether the most basic information available for a newly discovered radio variable or transient - their rise and/or decline rate - can be used to set initial constraints on the class of events from which they originate. We have analysed a sample of ≈800 synchrotron flares, selected from light curves of ≈90 sources observed at 5-8 GHz, representing a wide range of astrophysical phenomena, from flare stars to supermassive black holes. Selection of outbursts from the noisy radio light curves has been done automatically in order to ensure reproducibility of results. The distribution of rise/decline rates for the selected flares is modelled as a Gaussian probability distribution for each class of object, and further convolved with estimated areal density of that class in order to correct for the strong bias in our sample. We show in this way that comparing the measured variability time-scale of a radio transient/variable of unknown origin can provide an early, albeit approximate, classification of the object, and could form part of a suite of measurements used to provide early categorization of such events. Finally, we also discuss the effect scintillating sources will have on our ability to classify events based on their variability time-scales.

Optical properties of electrohydrodynamic convection patterns: rigorous and approximate methods.

PubMed

Bohley, Christian; Heuer, Jana; Stannarius, Ralf

2005-12-01

We analyze the optical behavior of two-dimensionally periodic structures that occur in electrohydrodynamic convection (EHC) patterns in nematic sandwich cells. These structures are anisotropic, locally uniaxial, and periodic on the scale of micrometers. For the first time, the optics of these structures is investigated with a rigorous method. The method used for the description of the electromagnetic waves interacting with EHC director patterns is a numerical approach that discretizes directly the Maxwell equations. It works as a space-grid-time-domain method and computes electric and magnetic fields in time steps. This so-called finite-difference-time-domain (FDTD) method is able to generate the fields with arbitrary accuracy. We compare this rigorous method with earlier attempts based on ray-tracing and analytical approximations. Results of optical studies of EHC structures made earlier based on ray-tracing methods are confirmed for thin cells, when the spatial periods of the pattern are sufficiently large. For the treatment of small-scale convection structures, the FDTD method is without alternatives.

Multiple continental radiations and correlates of diversification in Lupinus (Leguminosae): testing for key innovation with incomplete taxon sampling.

PubMed

Drummond, Christopher S; Eastwood, Ruth J; Miotto, Silvia T S; Hughes, Colin E

2012-05-01

Replicate radiations provide powerful comparative systems to address questions about the interplay between opportunity and innovation in driving episodes of diversification and the factors limiting their subsequent progression. However, such systems have been rarely documented at intercontinental scales. Here, we evaluate the hypothesis of multiple radiations in the genus Lupinus (Leguminosae), which exhibits some of the highest known rates of net diversification in plants. Given that incomplete taxon sampling, background extinction, and lineage-specific variation in diversification rates can confound macroevolutionary inferences regarding the timing and mechanisms of cladogenesis, we used Bayesian relaxed clock phylogenetic analyses as well as MEDUSA and BiSSE birth-death likelihood models of diversification, to evaluate the evolutionary patterns of lineage accumulation in Lupinus. We identified 3 significant shifts to increased rates of net diversification (r) relative to background levels in the genus (r = 0.18-0.48 lineages/myr). The primary shift occurred approximately 4.6 Ma (r = 0.48-1.76) in the montane regions of western North America, followed by a secondary shift approximately 2.7 Ma (r = 0.89-3.33) associated with range expansion and diversification of allopatrically distributed sister clades in the Mexican highlands and Andes. We also recovered evidence for a third independent shift approximately 6.5 Ma at the base of a lower elevation eastern South American grassland and campo rupestre clade (r = 0.36-1.33). Bayesian ancestral state reconstructions and BiSSE likelihood analyses of correlated diversification indicated that increased rates of speciation are strongly associated with the derived evolution of perennial life history and invasion of montane ecosystems. Although we currently lack hard evidence for "replicate adaptive radiations" in the sense of convergent morphological and ecological trajectories among species in different clades, these results are consistent with the hypothesis that iteroparity functioned as an adaptive key innovation, providing a mechanism for range expansion and rapid divergence in upper elevation regions across much of the New World.

Framework to trade optimality for local processing in large-scale wavefront reconstruction problems.

PubMed

Haber, Aleksandar; Verhaegen, Michel

2016-11-15

We show that the minimum variance wavefront estimation problems permit localized approximate solutions, in the sense that the wavefront value at a point (excluding unobservable modes, such as the piston mode) can be approximated by a linear combination of the wavefront slope measurements in the point's neighborhood. This enables us to efficiently compute a wavefront estimate by performing a single sparse matrix-vector multiplication. Moreover, our results open the possibility for the development of wavefront estimators that can be easily implemented in a decentralized/distributed manner, and in which the estimate optimality can be easily traded for computational efficiency. We numerically validate our approach on Hudgin wavefront sensor geometries, and the results can be easily generalized to Fried geometries.

[Analyses of pesticides in drinking water from small-scale water supplies in Schleswig-Holstein, Germany].

PubMed

Hippelein, M; Matthiessen, A; Kolychalow, O; Ostendorp, G

2012-12-01

In rural areas of Schleswig-Holstein, Germany, drinking water for about 37 000 people is provided by approximately 10 000 small-scale water supplies. For those wells data on pesticides in the drinking water are rare. In this study 100 small-scale water supplies, mainly situated in areas with intensive agriculture, fruit-growing or tree-nursery, were selected and the drinking water was analysed for pesticides. In 68 samples at least one agent or metabolite was detectable, 38 samples showed multiple contaminations. The metabolites dimethylsulfamide and chloridazone-desphenyl were found in nearly 40% of the wells in concentrations up to 42 µg/L. Bentazone was the most frequently detected biocidal agent. These data show that pesticides in drinking water from small-scale supplies are a notable issue in preventive public health. © Georg Thieme Verlag KG Stuttgart · New York.

A walk through the approximations of ab initio multiple spawning

NASA Astrophysics Data System (ADS)

Mignolet, Benoit; Curchod, Basile F. E.

2018-04-01

Full multiple spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the ab initio multiple spawning method which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the ab initio multiple spawning approximations on different observables and (ii) the convergence of the ab initio multiple spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying ab initio multiple spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.

A walk through the approximations of ab initio multiple spawning.

PubMed

Mignolet, Benoit; Curchod, Basile F E

2018-04-07

Full multiple spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the ab initio multiple spawning method which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the ab initio multiple spawning approximations on different observables and (ii) the convergence of the ab initio multiple spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying ab initio multiple spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.

A result about scale transformation families in approximation

NASA Astrophysics Data System (ADS)

Apprato, Dominique; Gout, Christian

2000-06-01

Scale transformations are common in approximation. In surface approximation from rapidly varying data, one wants to suppress, or at least dampen the oscillations of the approximation near steep gradients implied by the data. In that case, scale transformations can be used to give some control over overshoot when the surface has large variations of its gradient. Conversely, in image analysis, scale transformations are used in preprocessing to enhance some features present on the image or to increase jumps of grey levels before segmentation of the image. In this paper, we establish the convergence of an approximation method which allows some control over the behavior of the approximation. More precisely, we study the convergence of an approximation from a data set of , while using scale transformations on the values before and after classical approximation. In addition, the construction of scale transformations is also given. The algorithm is presented with some numerical examples.

Microfilament-Eruption Mechanism for Solar Spicules

NASA Technical Reports Server (NTRS)

Sterling, Alphonse C.; Moore, Ronald L.

2017-01-01

Recent studies indicate that solar coronal jets result from eruption of small-scale filaments, or "minifilaments" (Sterling et al. 2015, Nature, 523, 437; Panesar et al. ApJL, 832L, 7). In many aspects, these coronal jets appear to be small-scale versions of long-recognized large-scale solar eruptions that are often accompanied by eruption of a large-scale filament and that produce solar flares and coronal mass ejections (CMEs). In coronal jets, a jet-base bright point (JBP) that is often observed to accompany the jet and that sits on the magnetic neutral line from which the minifilament erupts, corresponds to the solar flare of larger-scale eruptions that occurs at the neutral line from which the large-scale filament erupts. Large-scale eruptions are relatively uncommon (approximately 1 per day) and occur with relatively large-scale erupting filaments (approximately 10 (sup 5) kilometers long). Coronal jets are more common (approximately 100s per day), but occur from erupting minifilaments of smaller size (approximately 10 (sup 4) kilometers long). It is known that solar spicules are much more frequent (many millions per day) than coronal jets. Just as coronal jets are small-scale versions of large-scale eruptions, here we suggest that solar spicules might in turn be small-scale versions of coronal jets; we postulate that the spicules are produced by eruptions of "microfilaments" of length comparable to the width of observed spicules (approximately 300 kilometers). A plot of the estimated number of the three respective phenomena (flares/CMEs, coronal jets, and spicules) occurring on the Sun at a given time, against the average sizes of erupting filaments, minifilaments, and the putative microfilaments, results in a size distribution that can be fitted with a power-law within the estimated uncertainties. The counterparts of the flares of large-scale eruptions and the JBPs of jets might be weak, pervasive, transient brightenings observed in Hinode/CaII images, and the production of spicules by microfilament eruptions might explain why spicules spin, as do coronal jets. The expected small-scale neutral lines from which the microfilaments would be expected to erupt would be difficult to detect reliably with current instrumentation, but might be apparent with instrumentation of the near future. A full report on this work appears in Sterling and Moore 2016, ApJL, 829, L9.

Web-based physiotherapy for people moderately affected with Multiple Sclerosis; quantitative and qualitative data from a randomized, controlled pilot study.

PubMed

Paul, Lorna; Coulter, Elaine H; Miller, Linda; McFadyen, Angus; Dorfman, Joe; Mattison, Paul George G

2014-09-01

To explore the effectiveness and participant experience of web-based physiotherapy for people moderately affected with Multiple Sclerosis (MS) and to provide data to establish the sample size required for a fully powered, definitive randomized controlled study. A randomized controlled pilot study. Rehabilitation centre and participants' homes. Thirty community dwelling adults moderately affected by MS (Expanded Disability Status Scale 5-6.5). Twelve weeks of individualised web-based physiotherapy completed twice per week or usual care (control). Online exercise diaries were monitored; participants were telephoned weekly by the physiotherapist and exercise programmes altered remotely by the physiotherapist as required. The following outcomes were completed at baseline and after 12 weeks; 25 Foot Walk, Berg Balance Scale, Timed Up and Go, Multiple Sclerosis Impact Scale, Leeds MS Quality of Life Scale, MS-Related Symptom Checklist and Hospital Anxiety and Depression Scale. The intervention group also completed a website evaluation questionnaire and interviews. Participants reported that website was easy to use, convenient, and motivating and would be happy to use in the future. There was no statistically significant difference in the primary outcome measure, the timed 25ft walk in the intervention group (P=0.170), or other secondary outcome measures, except the Multiple Sclerosis Impact Scale (P=0.048). Effect sizes were generally small to moderate. People with MS were very positive about web-based physiotherapy. The results suggested that 80 participants, 40 in each group, would be sufficient for a fully powered, definitive randomized controlled trial. © The Author(s) 2014.

Viscous-resistive layer in Rayleigh-Taylor instability

NASA Astrophysics Data System (ADS)

Silveira, F. E. M.; Orlandi, H. I.

2017-03-01

In this work, new scaling laws of the time growth rate γ of the Rayleigh-Taylor instability with the plasma resistivity η, kinematic viscosity ν, and electron number density ne are derived. A viscosity scale is defined in terms of the time decay of the perturbative fluid flow perpendicular to the equilibrium magnetic field, at the quasi-static approximation. Such a scale provides the identification of a viscous layer that can be combined with the resistive layer to produce a viscous-resistive layer. The latter, in turn, is found to satisfy an algebraic biquadratic equation. When viscous effects are negligible, it is shown that the viscous-resistive layer is given by the resistive layer. Somewhat surprisingly, when viscous effects cannot be neglected, it is shown that the viscous-resistive layer is given by the geometric mean of the resistive and viscous layers. A dispersion relation for the time growth rate is derived in terms of the viscous-resistive layer. When viscous effects cannot be neglected, two new scaling laws are found. At the quasi-static approximation, it is shown that γ ˜ (ην)1/4. However, on account of a finite electron mass, it is shown that γ˜(ν/ne ) 1 /3 . Further developments of our formulation are addressed in connection with a finite compressibility in the perturbative flow.

GCM simulations of volcanic aerosol forcing. I - Climate changes induced by steady-state perturbations

NASA Technical Reports Server (NTRS)

Pollack, James B.; Rind, David; Lacis, Andrew; Hansen, James E.; Sato, Makiko; Ruedy, Reto

1993-01-01

The response of the climate system to a temporally and spatially constant amount of volcanic particles is simulated using a general circulation model (GCM). The optical depth of the aerosols is chosen so as to produce approximately the same amount of forcing as results from doubling the present CO2 content of the atmosphere and from the boundary conditions associated with the peak of the last ice age. The climate changes produced by long-term volcanic aerosol forcing are obtained by differencing this simulation and one made for the present climate with no volcanic aerosol forcing. The simulations indicate that a significant cooling of the troposphere and surface can occur at times of closely spaced multiple sulfur-rich volcanic explosions that span time scales of decades to centuries. The steady-state climate response to volcanic forcing includes a large expansion of sea ice, especially in the Southern Hemisphere; a resultant large increase in surface and planetary albedo at high latitudes; and sizable changes in the annually and zonally averaged air temperature.

Fast computation of an optimal controller for large-scale adaptive optics.

PubMed

Massioni, Paolo; Kulcsár, Caroline; Raynaud, Henri-François; Conan, Jean-Marc

2011-11-01

The linear quadratic Gaussian regulator provides the minimum-variance control solution for a linear time-invariant system. For adaptive optics (AO) applications, under the hypothesis of a deformable mirror with instantaneous response, such a controller boils down to a minimum-variance phase estimator (a Kalman filter) and a projection onto the mirror space. The Kalman filter gain can be computed by solving an algebraic Riccati matrix equation, whose computational complexity grows very quickly with the size of the telescope aperture. This "curse of dimensionality" makes the standard solvers for Riccati equations very slow in the case of extremely large telescopes. In this article, we propose a way of computing the Kalman gain for AO systems by means of an approximation that considers the turbulence phase screen as the cropped version of an infinite-size screen. We demonstrate the advantages of the methods for both off- and on-line computational time, and we evaluate its performance for classical AO as well as for wide-field tomographic AO with multiple natural guide stars. Simulation results are reported.

Universal Shapes formed by Interacting Cracks

NASA Astrophysics Data System (ADS)

Fender, Melissa; Lechenault, Frederic; Daniels, Karen

2011-03-01

Brittle failure through multiple cracks occurs in a wide variety of contexts, from microscopic failures in dental enamel and cleaved silicon to geological faults and planetary ice crusts. In each of these situations, with complicated curvature and stress geometries, pairwise interactions between approaching cracks nonetheless produce characteristically curved fracture paths known in the geologic literature as en passant cracks. While the fragmentation of solids via many interacting cracks has seen wide investigation, less attention has been paid to the details of individual crack-crack interactions. We investigate the origins of this widely observed crack pattern using a rectangular elastic plate which is notched on each long side and then subjected to quasistatic uniaxial strain from the short side. The two cracks propagate along approximately straight paths until the pass each other, after which they curve and release a lenticular fragment. We find that, for materials with diverse mechanical properties, the shape of this fragment has an aspect ratio of 2:1, with the length scale set by the initial cracks offset s and the time scale set by the ratio of s to the pulling velocity. The cracks have a universal square root shape, which we understand by using a simple geometric model and the crack-crack interaction.

Hopelessness as a predictor of attempted suicide among first admission patients with psychosis: A 10-year cohort study

PubMed Central

Klonsky, E. David; Kotov, Roman; Bakst, Shelly; Rabinowitz, Jonathan; Bromet, Evelyn J.

2011-01-01

Little is known about the longitudinal relationship of hopelessness to attempted suicide in psychotic disorders. The present study addresses this gap by assessing hopelessness and attempted suicide at multiple time-points over 10-years in a first-admission cohort with psychosis (n=414). Approximately 1 in 5 participants attempted suicide during the 10-year follow-up, and those who attempted suicide scored significantly higher at baseline on the Beck Hopelessness Scale. In general, a given assessment of hopelessness (i.e., baseline, 6-months, 24-months, and 48-months) reliably predicted attempted suicide up to 4–6 years later, but not beyond. Structural equation modeling indicated that hopelessness prospectively predicted attempted suicide even when controlling for previous attempts. Notably, a cut-point of 3 or greater on the Beck Hopelessness Scale yielded sensitivity and specificity values similar to those found in non-psychotic populations using a cut-point of 9. Results suggest that hopelessness in individuals with psychotic disorders confers information about suicide risk above and beyond history of attempted suicide. Moreover, in comparison to non-psychotic populations, even relatively modest levels of hopelessness appear to confer risk for suicide in psychotic disorders. PMID:22320192

A parallel implementation of an off-lattice individual-based model of multicellular populations

NASA Astrophysics Data System (ADS)

Harvey, Daniel G.; Fletcher, Alexander G.; Osborne, James M.; Pitt-Francis, Joe

2015-07-01

As computational models of multicellular populations include ever more detailed descriptions of biophysical and biochemical processes, the computational cost of simulating such models limits their ability to generate novel scientific hypotheses and testable predictions. While developments in microchip technology continue to increase the power of individual processors, parallel computing offers an immediate increase in available processing power. To make full use of parallel computing technology, it is necessary to develop specialised algorithms. To this end, we present a parallel algorithm for a class of off-lattice individual-based models of multicellular populations. The algorithm divides the spatial domain between computing processes and comprises communication routines that ensure the model is correctly simulated on multiple processors. The parallel algorithm is shown to accurately reproduce the results of a deterministic simulation performed using a pre-existing serial implementation. We test the scaling of computation time, memory use and load balancing as more processes are used to simulate a cell population of fixed size. We find approximate linear scaling of both speed-up and memory consumption on up to 32 processor cores. Dynamic load balancing is shown to provide speed-up for non-regular spatial distributions of cells in the case of a growing population.

Magnetospheric MultiScale (MMS) System Manager

NASA Technical Reports Server (NTRS)

Schiff, Conrad; Maher, Francis Alfred; Henely, Sean Philip; Rand, David

2014-01-01

The Magnetospheric MultiScale (MMS) mission is an ambitious NASA space science mission in which 4 spacecraft are flown in tight formation about a highly elliptical orbit. Each spacecraft has multiple instruments that measure particle and field compositions in the Earths magnetosphere. By controlling the members relative motion, MMS can distinguish temporal and spatial fluctuations in a way that a single spacecraft cannot.To achieve this control, 2 sets of four maneuvers, distributed evenly across the spacecraft must be performed approximately every 14 days. Performing a single maneuver on an individual spacecraft is usually labor intensive and the complexity becomes clearly increases with four. As a result, the MMS flight dynamics team turned to the System Manager to put the routine or error-prone under machine control freeing the analysts for activities that require human judgment.The System Manager is an expert system that is capable of handling operations activities associated with performing MMS maneuvers. As an expert system, it can work off a known schedule, launching jobs based on a one-time occurrence or on a set reoccurring schedule. It is also able to detect situational changes and use event-driven programming to change schedules, adapt activities, or call for help.

A hybrid gyrokinetic ion and isothermal electron fluid code for astrophysical plasma

NASA Astrophysics Data System (ADS)

Kawazura, Y.; Barnes, M.

2018-05-01

This paper describes a new code for simulating astrophysical plasmas that solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal electron fluid (ITEF) Schekochihin et al. (2009) [9]. This model captures ion kinetic effects that are important near the ion gyro-radius scale while electron kinetic effects are ordered out by an electron-ion mass ratio expansion. The code is developed by incorporating the ITEF approximation into AstroGK, an Eulerian δf gyrokinetics code specialized to a slab geometry Numata et al. (2010) [41]. The new code treats the linear terms in the ITEF equations implicitly while the nonlinear terms are treated explicitly. We show linear and nonlinear benchmark tests to prove the validity and applicability of the simulation code. Since the fast electron timescale is eliminated by the mass ratio expansion, the Courant-Friedrichs-Lewy condition is much less restrictive than in full gyrokinetic codes; the present hybrid code runs ∼ 2√{mi /me } ∼ 100 times faster than AstroGK with a single ion species and kinetic electrons where mi /me is the ion-electron mass ratio. The improvement of the computational time makes it feasible to execute ion scale gyrokinetic simulations with a high velocity space resolution and to run multiple simulations to determine the dependence of turbulent dynamics on parameters such as electron-ion temperature ratio and plasma beta.

Body size distributions signal a regime shift in a lake ...

EPA Pesticide Factsheets

Communities of organisms, from mammals to microorganisms, have discontinuous distributions of body size. This pattern of size structuring is a conservative trait of community organization and is a product of processes that occur at multiple spatial and temporal scales. In this study, we assessed whether body size patterns serve as an indicator of a threshold between alternative regimes. Over the past 7000 years, the biological communities of Foy Lake (Montana,USA) have undergone a major regime shift owing to climate change. We used a palaeoecological record of diatom communities to estimate diatom sizes, and then analysed the discontinuous distribution of organism sizes over time. We used Bayesian classification and regression tree models to determine that all time intervals exhibited aggregations of sizes separated by gaps in the distribution and found a significant change in diatom body size distributions approximately 150 years before the identified ecosystem regime shift. We suggest that discontinuity analysis is a useful addition to the suite of tools for the detection of early warning signals of regime shifts. Communities of organisms from mammals to microorganisms have discontinuous distributions of body size. This pattern of size structuring is a conservative trait of community organization and is a product of processes that occur at discrete spatial and temporal scales within ecosystems. Here, a paleoecological record of diatom community change is use

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thomas, D; O’Connell, D; Lamb, J

Purpose: To demonstrate real-time dose calculation of free-breathing MRI guided Co−60 treatments, using a motion model and Monte-Carlo dose calculation to accurately account for the interplay between irregular breathing motion and an IMRT delivery. Methods: ViewRay Co-60 dose distributions were optimized on ITVs contoured from free-breathing CT images of lung cancer patients. Each treatment plan was separated into 0.25s segments, accounting for the MLC positions and beam angles at each time point. A voxel-specific motion model derived from multiple fast-helical free-breathing CTs and deformable registration was calculated for each patient. 3D images for every 0.25s of a simulated treatment weremore » generated in real time, here using a bellows signal as a surrogate to accurately account for breathing irregularities. Monte-Carlo dose calculation was performed every 0.25s of the treatment, with the number of histories in each calculation scaled to give an overall 1% statistical uncertainty. Each dose calculation was deformed back to the reference image using the motion model and accumulated. The static and real-time dose calculations were compared. Results: Image generation was performed in real time at 4 frames per second (GPU). Monte-Carlo dose calculation was performed at approximately 1frame per second (CPU), giving a total calculation time of approximately 30 minutes per treatment. Results show both cold- and hot-spots in and around the ITV, and increased dose to contralateral lung as the tumor moves in and out of the beam during treatment. Conclusion: An accurate motion model combined with a fast Monte-Carlo dose calculation allows almost real-time dose calculation of a free-breathing treatment. When combined with sagittal 2D-cine-mode MRI during treatment to update the motion model in real time, this will allow the true delivered dose of a treatment to be calculated, providing a useful tool for adaptive planning and assessing the effectiveness of gated treatments.« less

Enhanced sampling by multiple molecular dynamics trajectories: carbonmonoxy myoglobin 10 micros A0-->A(1-3) transition from ten 400 picosecond simulations.

PubMed

Loccisano, Anne E; Acevedo, Orlando; DeChancie, Jason; Schulze, Brita G; Evanseck, Jeffrey D

2004-05-01

The utility of multiple trajectories to extend the time scale of molecular dynamics simulations is reported for the spectroscopic A-states of carbonmonoxy myoglobin (MbCO). Experimentally, the A0-->A(1-3) transition has been observed to be 10 micros at 300 K, which is beyond the time scale of standard molecular dynamics simulations. To simulate this transition, 10 short (400 ps) and two longer time (1.2 ns) molecular dynamics trajectories, starting from five different crystallographic and solution phase structures with random initial velocities centered in a 37 A radius sphere of water, have been used to sample the native-fold of MbCO. Analysis of the ensemble of structures gathered over the cumulative 5.6 ns reveals two biomolecular motions involving the side chains of His64 and Arg45 to explain the spectroscopic states of MbCO. The 10 micros A0-->A(1-3) transition involves the motion of His64, where distance between His64 and CO is found to vary up to 8.8 +/- 1.0 A during the transition of His64 from the ligand (A(1-3)) to bulk solvent (A0). The His64 motion occurs within a single trajectory only once, however the multiple trajectories populate the spectroscopic A-states fully. Consequently, multiple independent molecular dynamics simulations have been found to extend biomolecular motion from 5 ns of total simulation to experimental phenomena on the microsecond time scale.

Measuring the impact of multiple sclerosis: Enhancing the measurement performance of the Multiple Sclerosis Impact Scale (MSIS-29) using Rasch Measurement Theory (RMT)

PubMed Central

Cleanthous, Sophie; Kinter, Elizabeth; Marquis, Patrick; Petrillo, Jennifer; You, Xiaojun; Wakeford, Craig; Sabatella, Guido

2017-01-01

Background Study objectives were to evaluate the Multiple Sclerosis Impact Scale (MSIS-29) and explore an optimized scoring structure based on empirical post-hoc analyses of data from the Phase III ADVANCE clinical trial. Methods ADVANCE MSIS-29 data from six time-points were analyzed in a sample of patients with relapsing–remitting multiple sclerosis (RRMS). Rasch Measurement Theory (RMT) analysis was undertaken to examine three broad areas: sample-to-scale targeting, measurement scale properties, and sample measurement validity. Interpretation of results led to an alternative MSIS-29 scoring structure, further evaluated alongside responsiveness of the original and revised scales at Week 48. Results RMT analysis provided mixed evidence for Physical and Psychological Impact scales that were sub-optimally targeted at the lower functioning end of the scales. Their conceptual basis could also stand to improve based on item fit results. The revised MSIS-29 rescored scales improved but did not resolve the measurement scale properties and targeting of the MSIS-29. In two out of three revised scales, responsiveness analysis indicated strengthened ability to detect change. Conclusion The revised MSIS-29 provides an initial evidence-based improved patient-reported outcome (PRO) instrument for evaluating the impact of MS. Revised scoring improves conceptual clarity and interpretation of scores by refining scale structure to include Symptoms, Psychological Impact, and General Limitations. Clinical trial ADVANCE (ClinicalTrials.gov identifier NCT00906399). PMID:29104758

A Multi-Scale Distribution Model for Non-Equilibrium Populations Suggests Resource Limitation in an Endangered Rodent

PubMed Central

Bean, William T.; Stafford, Robert; Butterfield, H. Scott; Brashares, Justin S.

2014-01-01

Species distributions are known to be limited by biotic and abiotic factors at multiple temporal and spatial scales. Species distribution models, however, frequently assume a population at equilibrium in both time and space. Studies of habitat selection have repeatedly shown the difficulty of estimating resource selection if the scale or extent of analysis is incorrect. Here, we present a multi-step approach to estimate the realized and potential distribution of the endangered giant kangaroo rat. First, we estimate the potential distribution by modeling suitability at a range-wide scale using static bioclimatic variables. We then examine annual changes in extent at a population-level. We define “available” habitat based on the total suitable potential distribution at the range-wide scale. Then, within the available habitat, model changes in population extent driven by multiple measures of resource availability. By modeling distributions for a population with robust estimates of population extent through time, and ecologically relevant predictor variables, we improved the predictive ability of SDMs, as well as revealed an unanticipated relationship between population extent and precipitation at multiple scales. At a range-wide scale, the best model indicated the giant kangaroo rat was limited to areas that received little to no precipitation in the summer months. In contrast, the best model for shorter time scales showed a positive relation with resource abundance, driven by precipitation, in the current and previous year. These results suggest that the distribution of the giant kangaroo rat was limited to the wettest parts of the drier areas within the study region. This multi-step approach reinforces the differing relationship species may have with environmental variables at different scales, provides a novel method for defining “available” habitat in habitat selection studies, and suggests a way to create distribution models at spatial and temporal scales relevant to theoretical and applied ecologists. PMID:25237807

Scaled CMOS Technology Reliability Users Guide

NASA Technical Reports Server (NTRS)

White, Mark

2010-01-01

The desire to assess the reliability of emerging scaled microelectronics technologies through faster reliability trials and more accurate acceleration models is the precursor for further research and experimentation in this relevant field. The effect of semiconductor scaling on microelectronics product reliability is an important aspect to the high reliability application user. From the perspective of a customer or user, who in many cases must deal with very limited, if any, manufacturer's reliability data to assess the product for a highly-reliable application, product-level testing is critical in the characterization and reliability assessment of advanced nanometer semiconductor scaling effects on microelectronics reliability. A methodology on how to accomplish this and techniques for deriving the expected product-level reliability on commercial memory products are provided.Competing mechanism theory and the multiple failure mechanism model are applied to the experimental results of scaled SDRAM products. Accelerated stress testing at multiple conditions is applied at the product level of several scaled memory products to assess the performance degradation and product reliability. Acceleration models are derived for each case. For several scaled SDRAM products, retention time degradation is studied and two distinct soft error populations are observed with each technology generation: early breakdown, characterized by randomly distributed weak bits with Weibull slope (beta)=1, and a main population breakdown with an increasing failure rate. Retention time soft error rates are calculated and a multiple failure mechanism acceleration model with parameters is derived for each technology. Defect densities are calculated and reflect a decreasing trend in the percentage of random defective bits for each successive product generation. A normalized soft error failure rate of the memory data retention time in FIT/Gb and FIT/cm2 for several scaled SDRAM generations is presented revealing a power relationship. General models describing the soft error rates across scaled product generations are presented. The analysis methodology may be applied to other scaled microelectronic products and their key parameters.

Variability of Disk Emission in Pre-Main Sequence and Related Stars. II. Variability in the Gas and Dust Emission of the Herbig Fe Star SAO 206462

NASA Technical Reports Server (NTRS)

Sitko, Michael L.; Day, Amanda N.; Kimes, Robin L.; Beerman, Lori C.; Martus, Cameron; Lynch, David K.; Russell, Ray W.; Grady, Carol A.; Schneider, Glenn; Lisse, Carey M.;

al3c: high-performance software for parameter inference using Approximate Bayesian Computation.

PubMed

Stram, Alexander H; Marjoram, Paul; Chen, Gary K

2015-11-01

The development of Approximate Bayesian Computation (ABC) algorithms for parameter inference which are both computationally efficient and scalable in parallel computing environments is an important area of research. Monte Carlo rejection sampling, a fundamental component of ABC algorithms, is trivial to distribute over multiple processors but is inherently inefficient. While development of algorithms such as ABC Sequential Monte Carlo (ABC-SMC) help address the inherent inefficiencies of rejection sampling, such approaches are not as easily scaled on multiple processors. As a result, current Bayesian inference software offerings that use ABC-SMC lack the ability to scale in parallel computing environments. We present al3c, a C++ framework for implementing ABC-SMC in parallel. By requiring only that users define essential functions such as the simulation model and prior distribution function, al3c abstracts the user from both the complexities of parallel programming and the details of the ABC-SMC algorithm. By using the al3c framework, the user is able to scale the ABC-SMC algorithm in parallel computing environments for his or her specific application, with minimal programming overhead. al3c is offered as a static binary for Linux and OS-X computing environments. The user completes an XML configuration file and C++ plug-in template for the specific application, which are used by al3c to obtain the desired results. Users can download the static binaries, source code, reference documentation and examples (including those in this article) by visiting https://github.com/ahstram/al3c. astram@usc.edu Supplementary data are available at Bioinformatics online. © The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.

The scaling of model test results to predict intake hot gas reingestion for STOVL aircraft with augmented vectored thrust engines

NASA Technical Reports Server (NTRS)

Penrose, C. J.

1987-01-01

The difficulties of modeling the complex recirculating flow fields produced by multiple jet STOVL aircraft close to the ground have led to extensive use of experimental model tests to predict intake Hot Gas Reingestion (HGR). Model test results reliability is dependent on a satisfactory set of scaling rules which must be validated by fully comparable full scale tests. Scaling rules devised in the U.K. in the mid 60's gave good model/full scale agreement for the BAe P1127 aircraft. Until recently no opportunity has occurred to check the applicability of the rules to the high energy exhaust of current ASTOVL aircraft projects. Such an opportunity has arisen following tests on a Tethered Harrier. Comparison of this full scale data and results from tests on a model configuration approximating to the full scale aircraft geometry has shown discrepancies between HGR levels. These discrepancies although probably due to geometry and other model/scale differences indicate some reexamination of the scaling rules is needed. Therefore the scaling rules are reviewed, further scaling studies planned are described and potential areas for further work are suggested.

Multiple Streaming and the Probability Distribution of Density in Redshift Space

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hui, Lam; Kofman, Lev; Shandarin, Sergei F.

2000-07-01

We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple streaming using the Zeldovich approximation (ZA), and compute the average number of streams in both real and redshift space. We find that multiple streaming can be significant in redshift space but negligible in real space, even at moderate values of the linear fluctuation amplitude ({sigma}{sub l}(less-or-similar sign)1). Moreover, unlike their real-space counterparts, redshift-space multiple streams can flow past each other with minimal interactions. Such nonlinear redshift-space effects, which aremore » physically distinct from the fingers-of-God due to small-scale virialized motions, might in part explain the well-known departure of redshift distortions from the classic linear prediction by Kaiser, even at relatively large scales where the corresponding density field in real space is well described by linear perturbation theory. We also compute, using the ZA, the probability distribution function (PDF) of the density, as well as S{sub 3}, in real and redshift space, and compare it with the PDF measured from N-body simulations. The role of caustics in defining the character of the high-density tail is examined. We find that (non-Lagrangian) smoothing, due to both finite resolution or discreteness and small-scale velocity dispersions, is very effective in erasing caustic structures, unless the initial power spectrum is sufficiently truncated. (c) 2000 The American Astronomical Society.« less

Multi-scale enhancement of climate prediction over land by improving the model sensitivity to vegetation variability

NASA Astrophysics Data System (ADS)

Alessandri, A.; Catalano, F.; De Felice, M.; Hurk, B. V. D.; Doblas-Reyes, F. J.; Boussetta, S.; Balsamo, G.; Miller, P. A.

2017-12-01

Here we demonstrate, for the first time, that the implementation of a realistic representation of vegetation in Earth System Models (ESMs) can significantly improve climate simulation and prediction across multiple time-scales. The effective sub-grid vegetation fractional coverage vary seasonally and at interannual time-scales in response to leaf-canopy growth, phenology and senescence. Therefore it affects biophysical parameters such as the surface resistance to evapotranspiration, albedo, roughness lenght, and soil field capacity. To adequately represent this effect in the EC-Earth ESM, we included an exponential dependence of the vegetation cover on the Leaf Area Index.By comparing two sets of simulations performed with and without the new variable fractional-coverage parameterization, spanning from centennial (20th Century) simulations and retrospective predictions to the decadal (5-years), seasonal (2-4 months) and weather (4 days) time-scales, we show for the first time a significant multi-scale enhancement of vegetation impacts in climate simulation and prediction over land. Particularly large effects at multiple time scales are shown over boreal winter middle-to-high latitudes over Canada, West US, Eastern Europe, Russia and eastern Siberia due to the implemented time-varying shadowing effect by tree-vegetation on snow surfaces. Over Northern Hemisphere boreal forest regions the improved representation of vegetation-cover consistently correct the winter warm biases, improves the climate change sensitivity, the decadal potential predictability as well as the skill of forecasts at seasonal and weather time-scales. Significant improvements of the prediction of 2m temperature and rainfall are also shown over transitional land surface hot spots. Both the potential predictability at decadal time-scale and seasonal-forecasts skill are enhanced over Sahel, North American Great Plains, Nordeste Brazil and South East Asia, mainly related to improved performance in the surface evapotranspiration.Above results are discussed in a peer-review paper just being accepted for publication on Climate Dynamics (Alessandri et al., 2017; doi:10.1007/s00382-017-3766-y).

X-ray observations of a flare in NGC 4151 from OSO 8

NASA Technical Reports Server (NTRS)

Mushotzky, R. F.; Holt, S. S.; Serlemitsos, P. J.

1978-01-01

The 2-60-keV flux from NGC 4151 has been observed to change by a factor of 2 on a time scale of 1.5 days. No fluctuations in excess of a factor of 3 are detected on time scales less than 4 hours. During a total observation of approximately 11 days there were no statistically significant changes in spectral shape. The spectrum can be fitted by a power law with photon index of about 1.42 + or - 0.06 and a hydrogen column density of approximately 7.5 + or - 0.5 x 10 to the 22nd power atoms/sq cm. A 2-sigma residual to this fit implies fluorescent Fe line emission with an equivalent width of about 240 eV. Both synchrotron self-Compton and thermal Compton models are consistent with the X-ray data.

Large-Scale Cubic-Scaling Random Phase Approximation Correlation Energy Calculations Using a Gaussian Basis.

PubMed

Wilhelm, Jan; Seewald, Patrick; Del Ben, Mauro; Hutter, Jürg

2016-12-13

We present an algorithm for computing the correlation energy in the random phase approximation (RPA) in a Gaussian basis requiring [Formula: see text] operations and [Formula: see text] memory. The method is based on the resolution of the identity (RI) with the overlap metric, a reformulation of RI-RPA in the Gaussian basis, imaginary time, and imaginary frequency integration techniques, and the use of sparse linear algebra. Additional memory reduction without extra computations can be achieved by an iterative scheme that overcomes the memory bottleneck of canonical RPA implementations. We report a massively parallel implementation that is the key for the application to large systems. Finally, cubic-scaling RPA is applied to a thousand water molecules using a correlation-consistent triple-ζ quality basis.

A multiple-time-scale turbulence model based on variable partitioning of turbulent kinetic energy spectrum

NASA Technical Reports Server (NTRS)

Kim, S.-W.; Chen, C.-P.

1987-01-01

A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

A multiple-time-scale turbulence model based on variable partitioning of the turbulent kinetic energy spectrum

NASA Technical Reports Server (NTRS)

Kim, S.-W.; Chen, C.-P.

1989-01-01

A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

Estimating False Discovery Proportion Under Arbitrary Covariance Dependence*

PubMed Central

Fan, Jianqing; Han, Xu; Gu, Weijie

2012-01-01

Multiple hypothesis testing is a fundamental problem in high dimensional inference, with wide applications in many scientific fields. In genome-wide association studies, tens of thousands of tests are performed simultaneously to find if any SNPs are associated with some traits and those tests are correlated. When test statistics are correlated, false discovery control becomes very challenging under arbitrary dependence. In the current paper, we propose a novel method based on principal factor approximation, which successfully subtracts the common dependence and weakens significantly the correlation structure, to deal with an arbitrary dependence structure. We derive an approximate expression for false discovery proportion (FDP) in large scale multiple testing when a common threshold is used and provide a consistent estimate of realized FDP. This result has important applications in controlling FDR and FDP. Our estimate of realized FDP compares favorably with Efron (2007)’s approach, as demonstrated in the simulated examples. Our approach is further illustrated by some real data applications. We also propose a dependence-adjusted procedure, which is more powerful than the fixed threshold procedure. PMID:24729644

Influence of orbital precession on the polar methane accumulation on Titan

NASA Astrophysics Data System (ADS)

Liu, J.; Schneider, T.

2014-12-01

Data collected by Cassini Spacecraft indicate that lakes on Titan are primarily found in the polar regions, preferentially in the north. It has been suggested that the hemispherical asymmetry in lake distribution is related to Saturn's orbital precession, which changes the seasonal distribution of solar radiation on Titan, but not the annual mean (Aharonson et al., 2009; Schneider et al., 2012). Saturn's current longitude of perihelion is near northern winter solstice. Hence, the northern summer on Titan is longer and less intense than the southern summer. The longer northern summer leads to greater net precipitation in the annual mean and the methane accumulation over the northern polar region (Schneider et al. 2012). Saturn's perihelion precesses over an approximately 45-kyr period, so the solar radiation at the top of Titan's atmosphere varies on this time scale. Here we investigate how the orbital precession influences the polar methane accumulation with a three-dimensional atmospheric model coupled to a dynamic surface reservoir of methane (Schneider et al. 2012). We find that methane accumulation is closely tied to Saturn's orbital precession. At the time when Saturn's longitude of perihelion is 180 degree away from the present day value, methane is mainly accumulated in the southern polar region due to the stronger annual-mean precipitation there induced by the longer southern summer. The annual-mean evaporation is largely unchanged with orbital precession, since it scales with the annual-mean insolation, which does not change under orbital precession. When Saturn's longitude of perihelion is close to equinox, methane is approximately evenly distributed in the northern and southern polar regions, and the lake dichotomy disappears. The timescale of methane redistribution from one pole to the other is short compared with the timescale of orbital precession, so the surface methane distribution can be viewed as being approximately in equilibrium with the solar forcing at any given historic time. These results indicate that the methane lake distribution on Titan likely varies over an approximately 45 Kyr time scale.

Exploration of EEG features of Alzheimer's disease using continuous wavelet transform.

PubMed

Ghorbanian, Parham; Devilbiss, David M; Hess, Terry; Bernstein, Allan; Simon, Adam J; Ashrafiuon, Hashem

2015-09-01

We have developed a novel approach to elucidate several discriminating EEG features of Alzheimer's disease. The approach is based on the use of a variety of continuous wavelet transforms, pairwise statistical tests with multiple comparison correction, and several decision tree algorithms, in order to choose the most prominent EEG features from a single sensor. A pilot study was conducted to record EEG signals from Alzheimer's disease (AD) patients and healthy age-matched control (CTL) subjects using a single dry electrode device during several eyes-closed (EC) and eyes-open (EO) resting conditions. We computed the power spectrum distribution properties and wavelet and sample entropy of the wavelet coefficients time series at scale ranges approximately corresponding to the major brain frequency bands. A predictive index was developed using the results from statistical tests and decision tree algorithms to identify the most reliable significant features of the AD patients when compared to healthy controls. The three most dominant features were identified as larger absolute mean power and larger standard deviation of the wavelet scales corresponding to 4-8 Hz (θ) during EO and lower wavelet entropy of the wavelet scales corresponding to 8-12 Hz (α) during EC, respectively. The fourth reliable set of distinguishing features of AD patients was lower relative power of the wavelet scales corresponding to 12-30 Hz (β) followed by lower skewness of the wavelet scales corresponding to 2-4 Hz (upper δ), both during EO. In general, the results indicate slowing and lower complexity of EEG signal in AD patients using a very easy-to-use and convenient single dry electrode device.

Global hydrodynamic modelling of flood inundation in continental rivers: How can we achieve it?

NASA Astrophysics Data System (ADS)

Yamazaki, D.

2016-12-01

Global-scale modelling of river hydrodynamics is essential for understanding global hydrological cycle, and is also required in interdisciplinary research fields . Global river models have been developed continuously for more than two decades, but modelling river flow at a global scale is still a challenging topic because surface water movement in continental rivers is a multi-spatial-scale phenomena. We have to consider the basin-wide water balance (>1000km scale), while hydrodynamics in river channels and floodplains is regulated by much smaller-scale topography (<100m scale). For example, heavy precipitation in upstream regions may later cause flooding in farthest downstream reaches. In order to realistically simulate the timing and amplitude of flood wave propagation for a long distance, consideration of detailed local topography is unavoidable. I have developed the global hydrodynamic model CaMa-Flood to overcome this scale-discrepancy of continental river flow. The CaMa-Flood divides river basins into multiple "unit-catchments", and assumes the water level is uniform within each unit-catchment. One unit-catchment is assigned to each grid-box defined at the typical spatial resolution of global climate models (10 100 km scale). Adopting a uniform water level in a >10km river segment seems to be a big assumption, but it is actually a good approximation for hydrodynamic modelling of continental rivers. The number of grid points required for global hydrodynamic simulations is largely reduced by this "unit-catchment assumption". Alternative to calculating 2-dimensional floodplain flows as in regional flood models, the CaMa-Flood treats floodplain inundation in a unit-catchment as a sub-grid physics. The water level and inundated area in each unit-catchment are diagnosed from water volume using topography parameters derived from high-resolution digital elevation models. Thus, the CaMa-Flood is at least 1000 times computationally more efficient compared to regional flood inundation models while the reality of simulated flood dynamics is kept. I will explain in detail how the CaMa-Flood model has been constructed from high-resolution topography datasets, and how the model can be used for various interdisciplinary applications.

The effect of rhythmic-cued motor imagery on walking, fatigue and quality of life in people with multiple sclerosis: A randomised controlled trial.

PubMed

Seebacher, Barbara; Kuisma, Raija; Glynn, Angela; Berger, Thomas

2017-02-01

Motor imagery and rhythmic auditory stimulation are physiotherapy strategies for walking rehabilitation. To investigate the effect of motor imagery combined with rhythmic cueing on walking, fatigue and quality of life (QoL) in people with multiple sclerosis (MS). Individuals with MS and Expanded Disability Status Scale scores of 1.5-4.5 were randomised into one of three groups: 17 minutes of motor imagery, six times per week, for 4 weeks, with music (A) or metronome cues (B), both with verbal cueing, and (C) controls. Primary outcomes were walking speed (Timed 25-Foot Walk) and distance (6-Minute Walk Test). Secondary outcomes were walking perception (Multiple Sclerosis Walking Scale-12), fatigue (Modified Fatigue Impact Scale) and QoL (Short Form-36 Health Survey, Multiple Sclerosis Impact Scale-29, Euroquol-5D-3L Questionnaire). Of the 112 participants randomised, 101 completed the study. Compared to controls, both interventions significantly improved walking speed, distance and perception. Significant improvements in cognitive but not psychosocial fatigue were seen in the intervention groups, and physical fatigue improved only in the music-based group. Both interventions improved QoL; however, music-cued motor imagery was superior at improving health-related QoL. Rhythmic-cued motor imagery improves walking, fatigue and QoL in people with MS, with music-cued motor imagery being more effective.

Dynamical evolution and spatial mixing of multiple population globular clusters

NASA Astrophysics Data System (ADS)

Vesperini, Enrico; McMillan, Stephen L. W.; D'Antona, Francesca; D'Ercole, Annibale

2013-03-01

Numerous spectroscopic and photometric observational studies have provided strong evidence for the widespread presence of multiple stellar populations in globular clusters. In this paper, we study the long-term dynamical evolution of multiple population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars. In previous studies, we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we explore the time-scales and the dynamics of the spatial mixing of the FG and the SG populations and their dependence on the SG initial concentration. Our simulations show that, as the evolution proceeds, the radial profile of the SG/FG number ratio, NSG/NFG, is characterized by three regions: (1) a flat inner part; (2) a declining part in which FG stars are increasingly dominant and (3) an outer region where the NSG/NFG profile flattens again (the NSG/NFG profile may rise slightly again in the outermost cluster regions). Until mixing is complete and the NSG/NFG profile is flat over the entire cluster, the radial variation of NSG/NFG implies that the fraction of SG stars determined by observations covering a limited range of radial distances is not, in general, equal to the SG global fraction, (NSG/NFG)glob. The distance at which NSG/NFG equals (NSG/NFG)glob is approximately between 1 and 2 cluster half-mass radii. The time-scale for complete mixing depends on the SG initial concentration, but in all cases complete mixing is expected only for clusters in advanced evolutionary phases, having lost at least 60-70 per cent of their mass due to two-body relaxation (in addition to the early FG loss due to the cluster expansion triggered by SNII ejecta and gas expulsion).The results of our simulations suggest that in many Galactic globular clusters the SG should still be more spatially concentrated than the FG.

Sediment dynamics over multiple time scales in Dyke Marsh Preserve (Potomac River, VA)

NASA Astrophysics Data System (ADS)

Palinkas, C. M.; Walters, D.

2010-12-01

Tidal freshwater marshes are critical components of fluvial and estuarine ecosystems, yet sediment dynamics within them have not received as much attention as their saltwater counterparts. This study examines sedimentation in Dyke Marsh Preserve, located on the Potomac River (VA), focusing on understanding the spatial variability present over multiple time scales. Bimonthly sediment data were collected using ceramic tiles, and seasonal- and decadal-scale sedimentation was determined via 7Be (half-life 53.3 days) and 210Pb (half-life 22.3 years), respectively. Results were also compared to SET data collected by the National Park Service since 2006. Preliminary data indicate that sites at lower elevations have higher sedimentation rates, likely related to their close proximity to the sediment source. Mass accumulation rates generally decreased with increasing time scale, such that the seasonal rates were greater than the SET-derived accretion rates, which were in turn greater than the decadal-scale rates. However, the bimonthly rates were the lowest observed, probably because the sampling period (May-October 2010) did not include the main depositional period of the year, which would be integrated by the other techniques.

Changes in photochemically significant solar UV spectral irradiance as estimated by the composite Mg II index and scale factors

NASA Technical Reports Server (NTRS)

Deland, Matthew T.; Cebula, Richard P.

1994-01-01

Quantitative assessment of the impact of solar ultraviolet irradiance variations on stratospheric ozone abundances currently requires the use of proxy indicators. The Mg II core-to-wing index has been developed as an indicator of solar UV activity between 175-400 nm that is independent of most instrument artifacts, and measures solar variability on both rotational and solar cycle time scales. Linear regression fits have been used to merge the individual Mg II index data sets from the Nimbus-7, NOAA-9, and NOAA-11 instruments onto a single reference scale. The change in 27-dayrunning average of the composite Mg II index from solar maximum to solar minimum is approximately 8 percent for solar cycle 21, and approximately 9 percent for solar cycle 22 through January 1992. Scaling factors based on the short-term variations in the Mg II index and solar irradiance data sets have been developed to estimate solar variability at mid-UV and near-UV wavelengths. Near 205 nm, where solar irradiance variations are important for stratospheric photo-chemistry and dynamics, the estimated change in irradiance during solar cycle 22 is approximately 10 percent using the composite Mg II index and scale factors.

Dissipation of Turbulence in the Solar Wind as Measured by Cluster

NASA Technical Reports Server (NTRS)

Goldstein, Melvyn

2012-01-01

Turbulence in fluids and plasmas is a scale-dependent process that generates fluctuations towards ever-smaller scales until dissipation occurs. Recent Cluster observations in the solar wind demonstrate the existence of a cascade of magnetic energy from the scale of the proton Larmor radius, where kinetic properties of ions invalidate fluid approximations, down to the electron Larmor radius, where electrons become demagnetized. The cascade is quasi-two-dimensional and has been interpreted as consisting of highly oblique kinetic Alfvenic fluctuations that dissipate near at the electron gyroradius scale via proton and electron Landau damping. Here we investigate for the first time the spatial properties of the turbulence at these scales. We report the presence of thin current sheets and discontinuities with spatial sizes greater than or approximately equal to the proton Larmor radius. These isolated structures may be manifestations of intermittency, and such would localize sites of turbulent dissipation. Studying the relationship between turbulent dissipation, reconnection and intermittency is crucial for understanding the dynamics of laboratory and astrophysical plasmas.

Detecting Genomic Clustering of Risk Variants from Sequence Data: Cases vs. Controls

PubMed Central

Schaid, Daniel J.; Sinnwell, Jason P.; McDonnell, Shannon K.; Thibodeau, Stephen N.

2013-01-01

As the ability to measure dense genetic markers approaches the limit of the DNA sequence itself, taking advantage of possible clustering of genetic variants in, and around, a gene would benefit genetic association analyses, and likely provide biological insights. The greatest benefit might be realized when multiple rare variants cluster in a functional region. Several statistical tests have been developed, one of which is based on the popular Kulldorff scan statistic for spatial clustering of disease. We extended another popular spatial clustering method – Tango’s statistic – to genomic sequence data. An advantage of Tango’s method is that it is rapid to compute, and when single test statistic is computed, its distribution is well approximated by a scaled chi-square distribution, making computation of p-values very rapid. We compared the Type-I error rates and power of several clustering statistics, as well as the omnibus sequence kernel association test (SKAT). Although our version of Tango’s statistic, which we call “Kernel Distance” statistic, took approximately half the time to compute than the Kulldorff scan statistic, it had slightly less power than the scan statistic. Our results showed that the Ionita-Laza version of Kulldorff’s scan statistic had the greatest power over a range of clustering scenarios. PMID:23842950

Initial eccentricity and constituent quark number scaling of elliptic flow in ideal and viscous dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chaudhuri, A. K.

2010-04-15

In the Israel-Stewart theory of dissipative hydrodynamics, the scaling properties of elliptic flow in Au+Au collisions are studied. The initial energy density of the fluid was fixed to reproduce STAR data on phi-meson multiplicity in 0-5% Au+Au collisions such that, irrespective of fluid viscosity, entropy at the freeze-out is similar in ideal or in viscous evolution. The initial eccentricity or constituent quark number scaling is only approximate in ideal or minimally viscous (eta/s=1/4pi) fluid. Eccentricity scaling becomes nearly exact in more viscous fluid (eta/s>=0.12). However, in more viscous fluid, constituent quark number scaled elliptic flow for mesons and baryons splitsmore » into separate scaling functions. Simulated flows also do not exhibit 'universal scaling'; that is, elliptic flow scaled by the constituent quark number and charged particles v{sub 2} is not a single function of transverse kinetic energy scaled by the quark number. From a study of the violation of universal scaling, we obtain an estimate of quark-gluon plasma viscosity, eta/s=0.12+-0.03. The error is statistical only. The systematic error in eta/s could be as large.« less

Summary of types of radiation belt electron precipitation observed by BARREL

NASA Astrophysics Data System (ADS)

Halford, Alexa

2016-07-01

The Balloon Array for Relativistic Radiation belt Electron Loss (BARREL) was able to infer precipitation of radiation belt electrons on multiple time scales and due to multiple loss mechanisms. One storm will be specifically highlighted which occurred on 26 January 2013 when a solar wind shock hit the Earth. Although MeV electrons were observed to be lost due to an EMIC wave event [Zhang et al in prep], and multiple periods of electron loss during substorms were observed [Rae et al submitted JGR, Mann et al in prep], we will consider an event period where loss associated with multiple time scales, and thus possibly different loss mechanisms was observed from 1000 - 1200 UT on 26 January 2013. At about 1005 UT on 26 January 2013 an injection of radiation belt electrons followed by drift echoes for energies of ˜80 - 400 keV. BARREL observed X-rays with energies less than 180 keV associated with multiple temporal structures during the drift echo event period. The Van Allen Probes were at similar L-values but upwards of 2 hours away in MLT. Upper band chorus and ULF waves were observed during the event period. Throughout the beginning of the event period, microbursts were clearly observed. During this time lower band chorus waves as well as time domain structures were observed at Van Allen Probe A located upwards of 2 hours away in MLT. This large difference in MLT meant that neither potential loss mechanism was able to be clearly associated with the microbursts. As the lower band chorus and time domain structures were observed to recede, the microbursts were also observed to subside. ULF time scale modulation of the X-rays was also observed throughout most of the event period. We will examine if the ULF waves are the cause of the precipitation themselves, or are modulating the loss of particles from a secondary loss mechanism [Brito et al 2015 JGR, Rae et al Submitted JGR]. Although the 100s ms and ULF time scales are clearly observed, there is an ˜20 minute overarching structure observed in the X-rays at BARREL. This longer time scale appears to match the drift period of the ˜300 keV electrons observed by the Van Allen probes. However the inferred energy of the precipitating electrons is ˜150 keV. It is unclear what may be causing the ˜20 minute structure in the X-rays. At the time of writing this abstract, it is unclear if the drifting of the 300 keV electrons is related to the precipitation of the lower energy electrons (< 180 keV) or if it is just coincidence that they have the same temporal structure.

Hamiltonian chaos acts like a finite energy reservoir: accuracy of the Fokker-Planck approximation.

PubMed

Riegert, Anja; Baba, Nilüfer; Gelfert, Katrin; Just, Wolfram; Kantz, Holger

2005-02-11

The Hamiltonian dynamics of slow variables coupled to fast degrees of freedom is modeled by an effective stochastic differential equation. Formal perturbation expansions, involving a Markov approximation, yield a Fokker-Planck equation in the slow subspace which respects conservation of energy. A detailed numerical and analytical analysis of suitable model systems demonstrates the feasibility of obtaining the system specific drift and diffusion terms and the accuracy of the stochastic approximation on all time scales. Non-Markovian and non-Gaussian features of the fast variables are negligible.

Multiple time scale analysis of sediment and runoff changes in the Lower Yellow River

NASA Astrophysics Data System (ADS)

Chi, Kaige; Gang, Zhao; Pang, Bo; Huang, Ziqian

2018-06-01

Sediment and runoff changes of seven hydrological stations along the Lower Yellow River (LYR) (Huayuankou Station, Jiahetan Station, Gaocun Station, Sunkou Station, Ai Shan Station, Qikou Station and Lijin Station) from 1980 to 2003 were alanyzed at multiple time scale. The maximum value of monthly, daily and hourly sediment load and runoff conservations were also analyzed with the annually mean value. Mann-Kendall non-parametric mathematics correlation test and Hurst coefficient method were adopted in the study. Research results indicate that (1) the runoff of seven hydrological stations was significantly reduced in the study period at different time scales. However, the trends of sediment load in these stations were not obvious. The sediment load of Huayuankou, Jiahetan and Aishan stations even slightly increased with the runoff decrease. (2) The trends of the sediment load with different time scale showed differences at Luokou and Lijin stations. Although the annually and monthly sediment load were broadly flat, the maximum hourly sediment load showed decrease trend. (3) According to the Hurst coefficients, the trend of sediment and runoff will be continue without taking measures, which proved the necessary of runoff-sediment regulation scheme.

Pilates exercise training vs. physical therapy for improving walking and balance in people with multiple sclerosis: a randomized controlled trial.

PubMed

Kalron, Alon; Rosenblum, Uri; Frid, Lior; Achiron, Anat

2017-03-01

Evaluate the effects of a Pilates exercise programme on walking and balance in people with multiple sclerosis and compare this exercise approach to conventional physical therapy sessions. Randomized controlled trial. Multiple Sclerosis Center, Sheba Medical Center, Tel-Hashomer, Israel. Forty-five people with multiple sclerosis, 29 females, mean age (SD) was 43.2 (11.6) years; mean Expanded Disability Status Scale (S.D) was 4.3 (1.3). Participants received 12 weekly training sessions of either Pilates ( n=22) or standardized physical therapy ( n=23) in an outpatient basis. Spatio-temporal parameters of walking and posturography parameters during static stance. Functional tests included the Time Up and Go Test, 2 and 6-minute walk test, Functional Reach Test, Berg Balance Scale and the Four Square Step Test. In addition, the following self-report forms included the Multiple Sclerosis Walking Scale and Modified Fatigue Impact Scale. At the termination, both groups had significantly increased their walking speed ( P=0.021) and mean step length ( P=0.023). According to the 2-minute and 6-minute walking tests, both groups at the end of the intervention program had increased their walking speed. Mean (SD) increase in the Pilates and physical therapy groups were 39.1 (78.3) and 25.3 (67.2) meters, respectively. There was no effect of group X time in all instrumented and clinical balance and gait measures. Pilates is a possible treatment option for people with multiple sclerosis in order to improve their walking and balance capabilities. However, this approach does not have any significant advantage over standardized physical therapy.

Spatial correlation and irradiance statistics in a multiple-beam terrestrial free-space optical communication link.

PubMed

Anguita, Jaime A; Neifeld, Mark A; Vasic, Bane V

2007-09-10

By means of numerical simulations we analyze the statistical properties of the power fluctuations induced by the incoherent superposition of multiple transmitted laser beams in a terrestrial free-space optical communication link. The measured signals arising from different transmitted optical beams are found to be statistically correlated. This channel correlation increases with receiver aperture and propagation distance. We find a simple scaling rule for the spatial correlation coefficient in terms of the propagation distance and we are able to predict the scintillation reduction in previously reported experiments with good accuracy. We propose an approximation to the probability density function of the received power of a spatially correlated multiple-beam system in terms of the parameters of the single-channel gamma-gamma function. A bit-error-rate evaluation is also presented to demonstrate the improvement of a multibeam system over its single-beam counterpart.

In situ redox manipulation of subsurface sediments from Fort Lewis, Washington: Iron reduction and TCE dechlorination mechanisms

DOE Office of Scientific and Technical Information (OSTI.GOV)

JE Szecsody; JS Fruchter; DS Sklarew

2000-03-21

Pacific Northwest National Laboratory (PNNL) conducted a bench-scale study to determine how effective chemically treated Ft. Lewis sediments can degrade trichloroethylene (TCE). The objectives of this experimental study were to quantify: (1) sediment reduction and oxidation reactions, (2) TCE degradation reactions, and (3) other significant geochemical changes that occurred. Sediment reduction and oxidation were investigated to determine the mass of reducible iron in the Ft. Lewis sediments and the rate of this reduction and subsequent oxidation at different temperatures. The temperature dependence was needed to be able to predict field-scale reduction in the relatively cold ({approximately}11 C) Ft. Lewis aquifer.more » Results of these experiments were used in conjunction with other geochemical and hydraulic characterization to design the field-scale injection experiment and predict barrier longevity. For example, the sediment reduction rate controls the amount of time required for the dithionite solution to fully react with sediments. Sediment oxidation experiments were additionally conducted to determine the oxidation rate and provide a separate measure of the mass of reduced iron. Laboratory experiments that were used to meet these objectives included: (1) sediment reduction in batch (static) systems, (2) sediment reduction in 1-D columns, and (3) sediment oxidation in 1-D columns. Multiple reaction modeling was conducted to quantify the reactant masses and reaction rates.« less

Low-Cost Ultra-High Spatial and Temporal Resolution Mapping of Intertidal Rock Platforms

NASA Astrophysics Data System (ADS)

Bryson, M.; Johnson-Roberson, M.; Murphy, R.

2012-07-01

Intertidal ecosystems have primarily been studied using field-based sampling; remote sensing offers the ability to collect data over large areas in a snapshot of time which could compliment field-based sampling methods by extrapolating them into the wider spatial and temporal context. Conventional remote sensing tools (such as satellite and aircraft imaging) provide data at relatively course, sub-meter resolutions or with limited temporal resolutions and relatively high costs for small-scale environmental science and ecology studies. In this paper, we describe a low-cost, kite-based imaging system and photogrammetric pipeline that was developed for constructing highresolution, 3D, photo-realistic terrain models of intertidal rocky shores. The processing pipeline uses automatic image feature detection and matching, structure-from-motion and photo-textured terrain surface reconstruction algorithms that require minimal human input and only a small number of ground control points and allow the use of cheap, consumer-grade digital cameras. The resulting maps combine colour and topographic information at sub-centimeter resolutions over an area of approximately 100m, thus enabling spatial properties of the intertidal environment to be determined across a hierarchy of spatial scales. Results of the system are presented for an intertidal rock platform at Cape Banks, Sydney, Australia. Potential uses of this technique include mapping of plant (micro- and macro-algae) and animal (e.g. gastropods) assemblages at multiple spatial and temporal scales.

Influence of configuration effects on multiple burst simulation testing

NASA Technical Reports Server (NTRS)

Emanuely, J. L.; Cantaloube, M.

1991-01-01

During the initial phase of a lightning strike attachment on an aircraft, fast current pulses (rise time approximately 100 ns, I(sub max) approximately few kA) were measured, which can create equipment upsets or disturbances. This threat, made of repetitive pulses and usually called 'multiple bursts', can be reproduced at the equipment interfaces assuming that the transfer function of the structure was determined. The normalized waveform H (10 kA - 100 ns rise time) is the reference for one of these pulses. The importance of the coaxial return path termination for the injection of the wave H is emphasized. According to the constitutive materials of the test bed, and the adaptation of the line, the natural oscillations of the structure and the internal coupling mechanisms can be modified. As a conclusion, various test configurations in relation with the nature of the test bed and the characteristics of the generator are detailed, for a more accurate ground simulation of the attachment phase.

Accelerating Electrostatic Surface Potential Calculation with Multiscale Approximation on Graphics Processing Units

PubMed Central

Anandakrishnan, Ramu; Scogland, Tom R. W.; Fenley, Andrew T.; Gordon, John C.; Feng, Wu-chun; Onufriev, Alexey V.

2010-01-01

Tools that compute and visualize biomolecular electrostatic surface potential have been used extensively for studying biomolecular function. However, determining the surface potential for large biomolecules on a typical desktop computer can take days or longer using currently available tools and methods. Two commonly used techniques to speed up these types of electrostatic computations are approximations based on multi-scale coarse-graining and parallelization across multiple processors. This paper demonstrates that for the computation of electrostatic surface potential, these two techniques can be combined to deliver significantly greater speed-up than either one separately, something that is in general not always possible. Specifically, the electrostatic potential computation, using an analytical linearized Poisson Boltzmann (ALPB) method, is approximated using the hierarchical charge partitioning (HCP) multiscale method, and parallelized on an ATI Radeon 4870 graphical processing unit (GPU). The implementation delivers a combined 934-fold speed-up for a 476,040 atom viral capsid, compared to an equivalent non-parallel implementation on an Intel E6550 CPU without the approximation. This speed-up is significantly greater than the 42-fold speed-up for the HCP approximation alone or the 182-fold speed-up for the GPU alone. PMID:20452792

Canopy wake measurements using multiple scanning wind LiDARs

NASA Astrophysics Data System (ADS)

Markfort, C. D.; Carbajo Fuertes, F.; Iungo, V.; Stefan, H. G.; Porte-Agel, F.

2014-12-01

Canopy wakes have been shown, in controlled wind tunnel experiments, to significantly affect the fluxes of momentum, heat and other scalars at the land and water surface over distances of ˜O(1 km), see Markfort et al. (EFM, 2013). However, there are currently no measurements of the velocity field downwind of a full-scale forest canopy. Point-based anemometer measurements of wake turbulence provide limited insight into the extent and details of the wake structure, whereas scanning Doppler wind LiDARs can provide information on how the wake evolves in space and varies over time. For the first time, we present measurements of the velocity field in the wake of a tall patch of forest canopy. The patch consists of two uniform rows of 40-meter tall deciduous, plane trees, which border either side of the Allée de Dorigny, near the EPFL campus. The canopy is approximately 250 m long, and it is approximately 40 m wide, along the direction of the wind. A challenge faced while making field measurements is that the wind rarely intersects a canopy normal to the edge. The resulting wake flow may be deflected relative to the mean inflow. Using multiple LiDARs, we measure the evolution of the wake due to an oblique wind blowing over the canopy. One LiDAR is positioned directly downwind of the canopy to measure the flow along the mean wind direction and the other is positioned near the canopy to evaluate the transversal component of the wind and how it varies with downwind distance from the canopy. Preliminary results show that the open trunk space near the base of the canopy results in a surface jet that can be detected just downwind of the canopy and farther downwind dissipates as it mixes with the wake flow above. A time-varying recirculation zone can be detected by the periodic reversal of the velocity near the surface, downwind of the canopy. The implications of canopy wakes for measurement and modeling of surface fluxes will be discussed.

BME Estimation of Residential Exposure to Ambient PM10 and Ozone at Multiple Time Scales

PubMed Central

Yu, Hwa-Lung; Chen, Jiu-Chiuan; Christakos, George; Jerrett, Michael

2009-01-01

Background Long-term human exposure to ambient pollutants can be an important contributing or etiologic factor of many chronic diseases. Spatiotemporal estimation (mapping) of long-term exposure at residential areas based on field observations recorded in the U.S. Environmental Protection Agency’s Air Quality System often suffer from missing data issues due to the scarce monitoring network across space and the inconsistent recording periods at different monitors. Objective We developed and compared two upscaling methods: UM1 (data aggregation followed by exposure estimation) and UM2 (exposure estimation followed by data aggregation) for the long-term PM10 (particulate matter with aerodynamic diameter ≤ 10 μm) and ozone exposure estimations and applied them in multiple time scales to estimate PM and ozone exposures for the residential areas of the Health Effects of Air Pollution on Lupus (HEAPL) study. Method We used Bayesian maximum entropy (BME) analysis for the two upscaling methods. We performed spatiotemporal cross-validations at multiple time scales by UM1 and UM2 to assess the estimation accuracy across space and time. Results Compared with the kriging method, the integration of soft information by the BME method can effectively increase the estimation accuracy for both pollutants. The spatiotemporal distributions of estimation errors from UM1 and UM2 were similar. The cross-validation results indicated that UM2 is generally better than UM1 in exposure estimations at multiple time scales in terms of predictive accuracy and lack of bias. For yearly PM10 estimations, both approaches have comparable performance, but the implementation of UM1 is associated with much lower computation burden. Conclusion BME-based upscaling methods UM1 and UM2 can assimilate core and site-specific knowledge bases of different formats for long-term exposure estimation. This study shows that UM1 can perform reasonably well when the aggregation process does not alter the spatiotemporal structure of the original data set; otherwise, UM2 is preferable. PMID:19440491

Deviations from uniform power law scaling in nonstationary time series

NASA Technical Reports Server (NTRS)

Viswanathan, G. M.; Peng, C. K.; Stanley, H. E.; Goldberger, A. L.

1997-01-01

A classic problem in physics is the analysis of highly nonstationary time series that typically exhibit long-range correlations. Here we test the hypothesis that the scaling properties of the dynamics of healthy physiological systems are more stable than those of pathological systems by studying beat-to-beat fluctuations in the human heart rate. We develop techniques based on the Fano factor and Allan factor functions, as well as on detrended fluctuation analysis, for quantifying deviations from uniform power-law scaling in nonstationary time series. By analyzing extremely long data sets of up to N = 10(5) beats for 11 healthy subjects, we find that the fluctuations in the heart rate scale approximately uniformly over several temporal orders of magnitude. By contrast, we find that in data sets of comparable length for 14 subjects with heart disease, the fluctuations grow erratically, indicating a loss of scaling stability.

Toward automatic time-series forecasting using neural networks.

PubMed

Yan, Weizhong

2012-07-01

Over the past few decades, application of artificial neural networks (ANN) to time-series forecasting (TSF) has been growing rapidly due to several unique features of ANN models. However, to date, a consistent ANN performance over different studies has not been achieved. Many factors contribute to the inconsistency in the performance of neural network models. One such factor is that ANN modeling involves determining a large number of design parameters, and the current design practice is essentially heuristic and ad hoc, this does not exploit the full potential of neural networks. Systematic ANN modeling processes and strategies for TSF are, therefore, greatly needed. Motivated by this need, this paper attempts to develop an automatic ANN modeling scheme. It is based on the generalized regression neural network (GRNN), a special type of neural network. By taking advantage of several GRNN properties (i.e., a single design parameter and fast learning) and by incorporating several design strategies (e.g., fusing multiple GRNNs), we have been able to make the proposed modeling scheme to be effective for modeling large-scale business time series. The initial model was entered into the NN3 time-series competition. It was awarded the best prediction on the reduced dataset among approximately 60 different models submitted by scholars worldwide.

Calibrating a method for simulated long-term ageing of biochar

NASA Astrophysics Data System (ADS)

Sohi, Saran; Cross, Andrew

2013-04-01

We recently established a procedure that imposes oxidatiave ageing to biochar and charcoal samples over a short time-frame, that provided carbon mass loss in the range projected for wild-fire charcoal in soil over a period of approximately 100 years. The stability of biochar samples in soil (relative to charcoal) range from 45-98% could be determined repeatably with high precision. Initial tests to understand the kinetics of the accelerated ageing method showed progressive increase in surface O concentration when examined by X-ray photoelectron spectroscopy (XPS) that slowly reached equilibrium. These trends resembled patterns observed in climate-for-time studies elsewhere, on centennial time-frame. We have extended this work to a preliminary direct calibration by matching progressive oxidation achieved in the laboratory to the surface composition of charcoal fragments recovered from the environment after periods of hundred to thousands of years. We have also applied artificial ageing to the same sets of naturally pre-aged charcoal fragments, and to recreated fresh charcoal. In this presentation of the first approach to quantifiably relate a laboratory test for biochar carbon stability to field data covering multiple time scales, we report on both the process and the implications for the stability of carbon stored in biochar under different climates and diverse agro-ecosystems.

An Incompressible, Depth-Averaged Lattice Boltzmann Method for Liquid Flow in Microfluidic Devices with Variable Aperture

DOE PAGES

Laleian, Artin; Valocchi, Albert J.; Werth, Charles J.

2015-11-24

Two-dimensional (2D) pore-scale models have successfully simulated microfluidic experiments of aqueous-phase flow with mixing-controlled reactions in devices with small aperture. A standard 2D model is not generally appropriate when the presence of mineral precipitate or biomass creates complex and irregular three-dimensional (3D) pore geometries. We modify the 2D lattice Boltzmann method (LBM) to incorporate viscous drag from the top and bottom microfluidic device (micromodel) surfaces, typically excluded in a 2D model. Viscous drag from these surfaces can be approximated by uniformly scaling a steady-state 2D velocity field at low Reynolds number. We demonstrate increased accuracy by approximating the viscous dragmore » with an analytically-derived body force which assumes a local parabolic velocity profile across the micromodel depth. Accuracy of the generated 2D velocity field and simulation permeability have not been evaluated in geometries with variable aperture. We obtain permeabilities within approximately 10% error and accurate streamlines from the proposed 2D method relative to results obtained from 3D simulations. Additionally, the proposed method requires a CPU run time approximately 40 times less than a standard 3D method, representing a significant computational benefit for permeability calculations.« less

Development of multiple chemical sensitivities in laborers after acute gasoline fume exposure in an underground tunneling operation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Davidoff, A.L.; Keyl, P.M.; Meggs, W.

1998-05-01

In this article, investigators report on the presence and nature of chemical sensitivities and other indices of illness in a cohort of workers excavating a new subway tunnel located under a former gasoline station. The workers were exposed to gasoline fumes for up to approximately 2 mo when they inadvertently dug into soil contaminated by gasoline. The cohort was unique in several ways: (a) contact with gasoline was made by the workers at a time when no one had complained of multiple chemical sensitivities syndrome; (b) all were males of low socioeconomic status; (c) the exposure was well documented; (d)more » the cohort could be considered naive because, at the time of the study, the men were not members of support groups and were not being seen by clinical ecologists, and they were not labeled, either by self or others, as having multiple chemical sensitivities syndrome or any related diagnosis; and (e) at the time of interview, all workers contacted appeared to be either gainfully employed or laid off temporarily and seeking gainful employment. The authors explored the health status of the workers at two different times: (1) soon after the tunnel was closed as a result of high, measured benzene-exposure levels and (2) 10--13 mo after the tunnel was closed. The workers were chronically overexposed to gasoline fumes, after which approximately one-fourth of their random sample of relatively naive, low-socio-economic-status male laborers--although neither disabled nor generally litigious--reported the new onset of chemical hypersensitivities and other characteristics that fit conservative criteria for multiple chemical sensitivities syndrome.« less

The multiple time scales of sleep dynamics as a challenge for modelling the sleeping brain.

PubMed

Olbrich, Eckehard; Claussen, Jens Christian; Achermann, Peter

2011-10-13

A particular property of the sleeping brain is that it exhibits dynamics on very different time scales ranging from the typical sleep oscillations such as sleep spindles and slow waves that can be observed in electroencephalogram (EEG) segments of several seconds duration over the transitions between the different sleep stages on a time scale of minutes to the dynamical processes involved in sleep regulation with typical time constants in the range of hours. There is an increasing body of work on mathematical and computational models addressing these different dynamics, however, usually considering only processes on a single time scale. In this paper, we review and present a new analysis of the dynamics of human sleep EEG at the different time scales and relate the findings to recent modelling efforts pointing out both the achievements and remaining challenges.

Large-scale disturbance legacies and the climate sensitivity of primary Picea abies forests.

PubMed

Schurman, Jonathan S; Trotsiuk, Volodymyr; Bače, Radek; Čada, Vojtěch; Fraver, Shawn; Janda, Pavel; Kulakowski, Dominik; Labusova, Jana; Mikoláš, Martin; Nagel, Thomas A; Seidl, Rupert; Synek, Michal; Svobodová, Kristýna; Chaskovskyy, Oleh; Teodosiu, Marius; Svoboda, Miroslav

2018-05-01

Determining the drivers of shifting forest disturbance rates remains a pressing global change issue. Large-scale forest dynamics are commonly assumed to be climate driven, but appropriately scaled disturbance histories are rarely available to assess how disturbance legacies alter subsequent disturbance rates and the climate sensitivity of disturbance. We compiled multiple tree ring-based disturbance histories from primary Picea abies forest fragments distributed throughout five European landscapes spanning the Bohemian Forest and the Carpathian Mountains. The regional chronology includes 11,595 tree cores, with ring dates spanning the years 1750-2000, collected from 560 inventory plots in 37 stands distributed across a 1,000 km geographic gradient, amounting to the largest disturbance chronology yet constructed in Europe. Decadal disturbance rates varied significantly through time and declined after 1920, resulting in widespread increases in canopy tree age. Approximately 75% of current canopy area recruited prior to 1900. Long-term disturbance patterns were compared to an historical drought reconstruction, and further linked to spatial variation in stand structure and contemporary disturbance patterns derived from LANDSAT imagery. Historically, decadal Palmer drought severity index minima corresponded to higher rates of canopy removal. The severity of contemporary disturbances increased with each stand's estimated time since last major disturbance, increased with mean diameter, and declined with increasing within-stand structural variability. Reconstructed spatial patterns suggest that high small-scale structural variability has historically acted to reduce large-scale susceptibility and climate sensitivity of disturbance. Reduced disturbance rates since 1920, a potential legacy of high 19th century disturbance rates, have contributed to a recent region-wide increase in disturbance susceptibility. Increasingly common high-severity disturbances throughout primary Picea forests of Central Europe should be reinterpreted in light of both legacy effects (resulting in increased susceptibility) and climate change (resulting in increased exposure to extreme events). © 2018 John Wiley & Sons Ltd.

Intrinsic vs. spurious long-range memory in high-frequency records of environmental radioactivity. Critical re-assessment and application to indoor 222Rn concentrations from Coimbra, Portugal

NASA Astrophysics Data System (ADS)

Donner, R. V.; Potirakis, S. M.; Barbosa, S. M.; Matos, J. A. O.; Pereira, A. J. S. C.; Neves, L. J. P. F.

2015-05-01

The presence or absence of long-range correlations in the environmental radioactivity fluctuations has recently attracted considerable interest. Among a multiplicity of practically relevant applications, identifying and disentangling the environmental factors controlling the variable concentrations of the radioactive noble gas radon is important for estimating its effect on human health and the efficiency of possible measures for reducing the corresponding exposition. In this work, we present a critical re-assessment of a multiplicity of complementary methods that have been previously applied for evaluating the presence of long-range correlations and fractal scaling in environmental radon variations with a particular focus on the specific properties of the underlying time series. As an illustrative case study, we subsequently re-analyze two high-frequency records of indoor radon concentrations from Coimbra, Portugal, each of which spans several weeks of continuous measurements at a high temporal resolution of five minutes.Our results reveal that at the study site, radon concentrations exhibit complex multi-scale dynamics with qualitatively different properties at different time-scales: (i) essentially white noise in the high-frequency part (up to time-scales of about one hour), (ii) spurious indications of a non-stationary, apparently long-range correlated process (at time scales between some hours and one day) arising from marked periodic components, and (iii) low-frequency variability indicating a true long-range dependent process. In the presence of such multi-scale variability, common estimators of long-range memory in time series are prone to fail if applied to the raw data without previous separation of time-scales with qualitatively different dynamics.

Atomoxetine pharmacokinetics in healthy Chinese subjects and effect of the CYP2D6*10 allele.

PubMed

Cui, Yi M; Teng, Choo H; Pan, Alan X; Yuen, Eunice; Yeo, Kwee P; Zhou, Ying; Zhao, Xia; Long, Amanda J; Bangs, Mark E; Wise, Stephen D

2007-10-01

To characterize atomoxetine pharmacokinetics, explore the effect of the homozygous CYP2D6*10 genotype on atomoxetine pharmacokinetics and evaluate the tolerability of atomoxetine, in healthy Chinese subjects. Twenty-four subjects, all CYP2D6 extensive metabolizers (EM), were randomized to receive atomoxetine (40 mg qd for 3 days, then 80 mg qd for 7 days) or matching placebo (2 : 1 ratio) in a double-blind fashion. Atomoxetine serum concentrations were measured following single (40 mg) and multiple (80 mg) doses. Adverse events, clinical safety laboratory data and vital signs were assessed during the study. Atomoxetine was rapidly absorbed with median time to maximum serum concentrations of approximately 1.5 h after single and multiple doses. Atomoxetine concentrations appeared to decrease monoexponentially with a mean apparent terminal half-life (t(1/2)) of approximately 4 h. The apparent clearance, apparent volume of distribution and t(1/2) following single and multiple doses were similar, suggesting linear pharmacokinetics with respect to time. Homozygous CYP2D6*10 subjects had 50% lower clearances compared with other EM subjects, resulting in twofold higher mean exposures. No clinically significant changes or abnormalities were noted in laboratory data and vital signs. The pharmacokinetics of atomoxetine in healthy Chinese subjects appears comparable to other ethnic populations. Multiple dosing of 80 mg qd atomoxetine was well tolerated in this study.

Slip-rate measurements on the Karakorum Fault may imply secular variations in fault motion.

PubMed

Chevalier, M-L; Ryerson, F J; Tapponnier, P; Finkel, R C; Van Der Woerd, J; Haibing, Li; Qing, Liu

2005-01-21

Beryllium-10 surface exposure dating of offset moraines on one branch of the Karakorum Fault west of the Gar basin yields a long-term (140- to 20-thousand-year) right-lateral slip rate of approximately 10.7 +/- 0.7 millimeters per year. This rate is 10 times larger than that inferred from recent InSAR analyses ( approximately 1 +/- 3 millimeters per year) that span approximately 8 years and sample all branches of the fault. The difference in slip-rate determinations suggests that large rate fluctuations may exist over centennial or millennial time scales. Such fluctuations would be consistent with mechanical coupling between the seismogenic, brittle-creep, and ductile shear sections of faults that reach deep into the crust.

Quantitative model of diffuse speckle contrast analysis for flow measurement.

PubMed

Liu, Jialin; Zhang, Hongchao; Lu, Jian; Ni, Xiaowu; Shen, Zhonghua

2017-07-01

Diffuse speckle contrast analysis (DSCA) is a noninvasive optical technique capable of monitoring deep tissue blood flow. However, a detailed study of the speckle contrast model for DSCA has yet to be presented. We deduced the theoretical relationship between speckle contrast and exposure time and further simplified it to a linear approximation model. The feasibility of this linear model was validated by the liquid phantoms which demonstrated that the slope of this linear approximation was able to rapidly determine the Brownian diffusion coefficient of the turbid media at multiple distances using multiexposure speckle imaging. Furthermore, we have theoretically quantified the influence of optical property on the measurements of the Brownian diffusion coefficient which was a consequence of the fact that the slope of this linear approximation was demonstrated to be equal to the inverse of correlation time of the speckle.

Subsampled Hessian Newton Methods for Supervised Learning.

PubMed

Wang, Chien-Chih; Huang, Chun-Heng; Lin, Chih-Jen

2015-08-01

Newton methods can be applied in many supervised learning approaches. However, for large-scale data, the use of the whole Hessian matrix can be time-consuming. Recently, subsampled Newton methods have been proposed to reduce the computational time by using only a subset of data for calculating an approximation of the Hessian matrix. Unfortunately, we find that in some situations, the running speed is worse than the standard Newton method because cheaper but less accurate search directions are used. In this work, we propose some novel techniques to improve the existing subsampled Hessian Newton method. The main idea is to solve a two-dimensional subproblem per iteration to adjust the search direction to better minimize the second-order approximation of the function value. We prove the theoretical convergence of the proposed method. Experiments on logistic regression, linear SVM, maximum entropy, and deep networks indicate that our techniques significantly reduce the running time of the subsampled Hessian Newton method. The resulting algorithm becomes a compelling alternative to the standard Newton method for large-scale data classification.

Behavioral responses of wolves to roads: scale-dependent ambivalence

PubMed Central

Nelson, Lindsey; Wabakken, Petter; Sand, Håkan; Liberg, Olof

2014-01-01

Throughout their recent recovery in several industrialized countries, large carnivores have had to cope with a changed landscape dominated by human infrastructure. Population growth depends on the ability of individuals to adapt to these changes by making use of new habitat features and at the same time to avoid increased risks of mortality associated with human infrastructure. We analyzed the summer movements of 19 GPS-collared resident wolves (Canis lupus L.) from 14 territories in Scandinavia in relation to roads. We used resource and step selection functions, including >12000 field-checked GPS-positions and 315 kill sites. Wolves displayed ambivalent responses to roads depending on the spatial scale, road type, time of day, behavioral state, and reproductive status. At the site scale (approximately 0.1 km2), they selected for roads when traveling, nearly doubling their travel speed. Breeding wolves moved the fastest. At the patch scale (10 km2), house density rather than road density was a significant negative predictor of wolf patch selection. At the home range scale (approximately 1000 km2), breeding wolves increased gravel road use with increasing road availability, although at a lower rate than expected. Wolves have adapted to use roads for ease of travel, but at the same time developed a cryptic behavior to avoid human encounters. This behavioral plasticity may have been important in allowing the successful recovery of wolf populations in industrialized countries. However, we emphasize the role of roads as a potential cause of increased human-caused mortality. PMID:25419085

Behavioral responses of wolves to roads: scale-dependent ambivalence.

PubMed

Zimmermann, Barbara; Nelson, Lindsey; Wabakken, Petter; Sand, Håkan; Liberg, Olof

2014-11-01

Throughout their recent recovery in several industrialized countries, large carnivores have had to cope with a changed landscape dominated by human infrastructure. Population growth depends on the ability of individuals to adapt to these changes by making use of new habitat features and at the same time to avoid increased risks of mortality associated with human infrastructure. We analyzed the summer movements of 19 GPS-collared resident wolves ( Canis lupus L.) from 14 territories in Scandinavia in relation to roads. We used resource and step selection functions, including >12000 field-checked GPS-positions and 315 kill sites. Wolves displayed ambivalent responses to roads depending on the spatial scale, road type, time of day, behavioral state, and reproductive status. At the site scale (approximately 0.1 km 2 ), they selected for roads when traveling, nearly doubling their travel speed. Breeding wolves moved the fastest. At the patch scale (10 km 2 ), house density rather than road density was a significant negative predictor of wolf patch selection. At the home range scale (approximately 1000 km 2 ), breeding wolves increased gravel road use with increasing road availability, although at a lower rate than expected. Wolves have adapted to use roads for ease of travel, but at the same time developed a cryptic behavior to avoid human encounters. This behavioral plasticity may have been important in allowing the successful recovery of wolf populations in industrialized countries. However, we emphasize the role of roads as a potential cause of increased human-caused mortality.

Memory: Ironing Out a Wrinkle in Time.

PubMed

Miller, Adam M P; Frankland, Paul W; Josselyn, Sheena A

2018-05-21

Individual hippocampal neurons encode time over seconds, whereas large-scale changes in population activity of hippocampal neurons encode time over minutes and days. New research shows how the hippocampus represents these multiple timescales simultaneously. Copyright © 2018 Elsevier Ltd. All rights reserved.

Efficiency and flexibility using implicit methods within atmosphere dycores

NASA Astrophysics Data System (ADS)

Evans, K. J.; Archibald, R.; Norman, M. R.; Gardner, D. J.; Woodward, C. S.; Worley, P.; Taylor, M.

2016-12-01

A suite of explicit and implicit methods are evaluated for a range of configurations of the shallow water dynamical core within the spectral-element Community Atmosphere Model (CAM-SE) to explore their relative computational performance. The configurations are designed to explore the attributes of each method under different but relevant model usage scenarios including varied spectral order within an element, static regional refinement, and scaling to large problem sizes. The limitations and benefits of using explicit versus implicit, with different discretizations and parameters, are discussed in light of trade-offs such as MPI communication, memory, and inherent efficiency bottlenecks. For the regionally refined shallow water configurations, the implicit BDF2 method is about the same efficiency as an explicit Runge-Kutta method, without including a preconditioner. Performance of the implicit methods with the residual function executed on a GPU is also presented; there is speed up for the residual relative to a CPU, but overwhelming transfer costs motivate moving more of the solver to the device. Given the performance behavior of implicit methods within the shallow water dynamical core, the recommendation for future work using implicit solvers is conditional based on scale separation and the stiffness of the problem. The strong growth of linear iterations with increasing resolution or time step size is the main bottleneck to computational efficiency. Within the hydrostatic dynamical core, of CAM-SE, we present results utilizing approximate block factorization preconditioners implemented using the Trilinos library of solvers. They reduce the cost of linear system solves and improve parallel scalability. We provide a summary of the remaining efficiency considerations within the preconditioner and utilization of the GPU, as well as a discussion about the benefits of a time stepping method that provides converged and stable solutions for a much wider range of time step sizes. As more complex model components, for example new physics and aerosols, are connected in the model, having flexibility in the time stepping will enable more options for combining and resolving multiple scales of behavior.

A time series of urban extent in China using DSMP/OLS nighttime light data

PubMed Central

Chen, Dongsheng; Chen, Le; Wang, Huan; Guan, Qingfeng

2018-01-01

Urban extent data play an important role in urban management and urban studies, such as monitoring the process of urbanization and changes in the spatial configuration of urban areas. Traditional methods of extracting urban-extent information are primarily based on manual investigations and classifications using remote sensing images, and these methods have such problems as large costs in labor and time and low precision. This study proposes an improved, simplified and flexible method for extracting urban extents over multiple scales and the construction of spatiotemporal models using DMSP/OLS nighttime light (NTL) for practical situations. This method eliminates the regional temporal and spatial inconsistency of thresholding NTL in large-scale and multi-temporal scenes. Using this method, we have extracted the urban extents and calculated the corresponding areas on the county, municipal and provincial scales in China from 2000 to 2012. In addition, validation with the data of reference data shows that the overall accuracy (OA), Kappa and F1 Scores were 0.996, 0.793, and 0.782, respectively. We increased the spatial resolution of the urban extent to 500 m (approximately four times finer than the results of previous studies). Based on the urban extent dataset proposed above, we analyzed changes in urban extents over time and observed that urban sprawl has grown in all of the counties of China. We also identified three patterns of urban sprawl: Early Urban Growth, Constant Urban Growth and Recent Urban Growth. In addition, these trends of urban sprawl are consistent with the western, eastern and central cities of China, respectively, in terms of their spatial distribution, socioeconomic characteristics and historical background. Additionally, the urban extents display the spatial configurations of urban areas intuitively. The proposed urban extent dataset is available for download and can provide reference data and support for future studies of urbanization and urban planning. PMID:29795685

Solution of the Average-Passage Equations for the Incompressible Flow through Multiple-Blade-Row Turbomachinery

DTIC Science & Technology

1994-02-01

numerical treatment. An explicit numerical procedure based on Runqe-Kutta time stepping for cell-centered, hexahedral finite volumes is...An explicit numerical procedure based on Runge-Kutta time stepping for cell-centered, hexahedral finite volumes is outlined for the approximate...Discretization 16 3.1 Cell-Centered Finite -Volume Discretization in Space 16 3.2 Artificial Dissipation 17 3.3 Time Integration 21 3.4 Convergence

Chorus Whistler Wave Source Scales As Determined From Multipoint Van Allen Probe Measurements

NASA Technical Reports Server (NTRS)

Agapitov, O.; Blum, L. W.; Mozer, F. S.; Bonnell, J. W.; Wygant, J.

2017-01-01

Whistler mode chorus waves are particularly important in outer radiation belt dynamics due to their key role in controlling the acceleration and scattering of electrons over a very wide energy range. The key parameters for both nonlinear and quasi-linear treatment of wave-particle interactions are the temporal and spatial scales of the wave source region and coherence of the wave field perturbations. Neither the source scale nor the coherence scale is well established experimentally, mostly because of a lack of multipoint VLF waveform measurements. We present an unprecedentedly long interval of coordinated VLF waveform measurements (sampled at 16384 s(exp -1)) aboard the two Van Allen Probes spacecraft-9 h (0800-1200 UT and 1700-2200 UT) during two consecutive apogees on 15 July 2014. The spacecraft separations varied from about 100 to 5000 km (mostly radially); measurements covered an L shell range from 3 to 6; magnetic local time 0430-0900, and magnetic latitudes were approximately 15 and approximately 5 deg during the two orbits. Using time-domain correlation techniques, the single chorus source spatial extent transverse to the background magnetic field has been determined to be about 550-650 km for upper band chorus waves with amplitudes less than 100 pT and up to 800 km for larger amplitude, lower band chorus waves. The ratio between wave amplitudes measured on the two spacecraft is also examined to reveal that the wave amplitude distribution within a single chorus element generation area can be well approximated by a Gaussian exp(-0.5 x r (exp 2)/r(sub 0)(exp 2)), with the characteristic scale r(sub 0) around 300 km. Waves detected by the two spacecraft were found to be coherent in phase at distances up to 400 km.

Investigating lithium-ion battery materials during overcharge-induced thermal runaway: an operando and multi-scale X-ray CT study.

PubMed

Finegan, Donal P; Scheel, Mario; Robinson, James B; Tjaden, Bernhard; Di Michiel, Marco; Hinds, Gareth; Brett, Dan J L; Shearing, Paul R

2016-11-16

Catastrophic failure of lithium-ion batteries occurs across multiple length scales and over very short time periods. A combination of high-speed operando tomography, thermal imaging and electrochemical measurements is used to probe the degradation mechanisms leading up to overcharge-induced thermal runaway of a LiCoO 2 pouch cell, through its interrelated dynamic structural, thermal and electrical responses. Failure mechanisms across multiple length scales are explored using a post-mortem multi-scale tomography approach, revealing significant morphological and phase changes in the LiCoO 2 electrode microstructure and location dependent degradation. This combined operando and multi-scale X-ray computed tomography (CT) technique is demonstrated as a comprehensive approach to understanding battery degradation and failure.

Linear Approximation SAR Azimuth Processing Study

NASA Technical Reports Server (NTRS)

Lindquist, R. B.; Masnaghetti, R. K.; Belland, E.; Hance, H. V.; Weis, W. G.

1979-01-01

A segmented linear approximation of the quadratic phase function that is used to focus the synthetic antenna of a SAR was studied. Ideal focusing, using a quadratic varying phase focusing function during the time radar target histories are gathered, requires a large number of complex multiplications. These can be largely eliminated by using linear approximation techniques. The result is a reduced processor size and chip count relative to ideally focussed processing and a correspondingly increased feasibility for spaceworthy implementation. A preliminary design and sizing for a spaceworthy linear approximation SAR azimuth processor meeting requirements similar to those of the SEASAT-A SAR was developed. The study resulted in a design with approximately 1500 IC's, 1.2 cubic feet of volume, and 350 watts of power for a single look, 4000 range cell azimuth processor with 25 meters resolution.

Computing the universe: how large-scale simulations illuminate galaxies and dark energy

NASA Astrophysics Data System (ADS)

O'Shea, Brian

2015-04-01

High-performance and large-scale computing is absolutely to understanding astronomical objects such as stars, galaxies, and the cosmic web. This is because these are structures that operate on physical, temporal, and energy scales that cannot be reasonably approximated in the laboratory, and whose complexity and nonlinearity often defies analytic modeling. In this talk, I show how the growth of computing platforms over time has facilitated our understanding of astrophysical and cosmological phenomena, focusing primarily on galaxies and large-scale structure in the Universe.

Time-localized wavelet multiple regression and correlation

NASA Astrophysics Data System (ADS)

Fernández-Macho, Javier

2018-02-01

This paper extends wavelet methodology to handle comovement dynamics of multivariate time series via moving weighted regression on wavelet coefficients. The concept of wavelet local multiple correlation is used to produce one single set of multiscale correlations along time, in contrast with the large number of wavelet correlation maps that need to be compared when using standard pairwise wavelet correlations with rolling windows. Also, the spectral properties of weight functions are investigated and it is argued that some common time windows, such as the usual rectangular rolling window, are not satisfactory on these grounds. The method is illustrated with a multiscale analysis of the comovements of Eurozone stock markets during this century. It is shown how the evolution of the correlation structure in these markets has been far from homogeneous both along time and across timescales featuring an acute divide across timescales at about the quarterly scale. At longer scales, evidence from the long-term correlation structure can be interpreted as stable perfect integration among Euro stock markets. On the other hand, at intramonth and intraweek scales, the short-term correlation structure has been clearly evolving along time, experiencing a sharp increase during financial crises which may be interpreted as evidence of financial 'contagion'.

On the use of finite difference matrix-vector products in Newton-Krylov solvers for implicit climate dynamics with spectral elements

DOE PAGES

Woodward, Carol S.; Gardner, David J.; Evans, Katherine J.

2015-01-01

Efficient solutions of global climate models require effectively handling disparate length and time scales. Implicit solution approaches allow time integration of the physical system with a step size governed by accuracy of the processes of interest rather than by stability of the fastest time scales present. Implicit approaches, however, require the solution of nonlinear systems within each time step. Usually, a Newton's method is applied to solve these systems. Each iteration of the Newton's method, in turn, requires the solution of a linear model of the nonlinear system. This model employs the Jacobian of the problem-defining nonlinear residual, but thismore » Jacobian can be costly to form. If a Krylov linear solver is used for the solution of the linear system, the action of the Jacobian matrix on a given vector is required. In the case of spectral element methods, the Jacobian is not calculated but only implemented through matrix-vector products. The matrix-vector multiply can also be approximated by a finite difference approximation which may introduce inaccuracy in the overall nonlinear solver. In this paper, we review the advantages and disadvantages of finite difference approximations of these matrix-vector products for climate dynamics within the spectral element shallow water dynamical core of the Community Atmosphere Model.« less

Reusable Launch Vehicle Control in Multiple Time Scale Sliding Modes

NASA Technical Reports Server (NTRS)

Shtessel, Yuri

1999-01-01

A reusable launch vehicle control problem during ascent is addressed via multiple-time scaled continuous sliding mode control. The proposed sliding mode controller utilizes a two-loop structure and provides robust, de-coupled tracking of both orientation angle command profiles and angular rate command profiles in the presence of bounded external disturbances and plant uncertainties. Sliding mode control causes the angular rate and orientation angle tracking error dynamics to be constrained to linear, de-coupled, homogeneous, and vector valued differential equations with desired eigenvalues placement. The dual-time scale sliding mode controller was designed for the X-33 technology demonstration sub-orbital launch vehicle in the launch mode. 6DOF simulation results show that the designed controller provides robust, accurate, de-coupled tracking of the orientation angle command profiles in presence of external disturbances and vehicle inertia uncertainties. It creates possibility to operate the X-33 vehicle in an aircraft-like mode with reduced pre-launch adjustment of the control system.

Linear regulator design for stochastic systems by a multiple time scales method

NASA Technical Reports Server (NTRS)

Teneketzis, D.; Sandell, N. R., Jr.

1976-01-01

A hierarchically-structured, suboptimal controller for a linear stochastic system composed of fast and slow subsystems is considered. The controller is optimal in the limit as the separation of time scales of the subsystems becomes infinite. The methodology is illustrated by design of a controller to suppress the phugoid and short period modes of the longitudinal dynamics of the F-8 aircraft.

Combined action of time-delay and colored cross-associated multiplicative and additive noises on stability and stochastic resonance for a stochastic metapopulation system

NASA Astrophysics Data System (ADS)

Wang, Kang-Kang; Zong, De-Cai; Wang, Ya-Jun; Li, Sheng-Hong

2016-05-01

In this paper, the transition between the stable state of a big density and the extinction state and stochastic resonance (SR) for a time-delayed metapopulation system disturbed by colored cross-correlated noises are investigated. By applying the fast descent method, the small time-delay approximation and McNamara and Wiesenfeld's SR theory, we investigate the impacts of time-delay, the multiplicative, additive noises and colored cross-correlated noise on the SNR and the shift between the two states of the system. Numerical results show that the multiplicative, additive noises and time-delay can all speed up the transition from the stable state to the extinction state, while the correlation noise and its correlation time can slow down the extinction process of the population system. With respect to SNR, the multiplicative noise always weakens the SR effect, while noise correlation time plays a dual role in motivating the SR phenomenon. Meanwhile, time-delay mainly plays a negative role in stimulating the SR phenomenon. Conversely, it could motivate the SR effect to increase the strength of the cross-correlation noise in the SNR-β plot, while the increase of additive noise intensity will firstly excite SR, and then suppress the SR effect.

Doubly stochastic Poisson process models for precipitation at fine time-scales

NASA Astrophysics Data System (ADS)

Ramesh, Nadarajah I.; Onof, Christian; Xie, Dichao

2012-09-01

This paper considers a class of stochastic point process models, based on doubly stochastic Poisson processes, in the modelling of rainfall. We examine the application of this class of models, a neglected alternative to the widely-known Poisson cluster models, in the analysis of fine time-scale rainfall intensity. These models are mainly used to analyse tipping-bucket raingauge data from a single site but an extension to multiple sites is illustrated which reveals the potential of this class of models to study the temporal and spatial variability of precipitation at fine time-scales.

Large-scale structure in brane-induced gravity. I. Perturbation theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Scoccimarro, Roman

2009-11-15

We study the growth of subhorizon perturbations in brane-induced gravity using perturbation theory. We solve for the linear evolution of perturbations taking advantage of the symmetry under gauge transformations along the extra-dimension to decouple the bulk equations in the quasistatic approximation, which we argue may be a better approximation at large scales than thought before. We then study the nonlinearities in the bulk and brane equations, concentrating on the workings of the Vainshtein mechanism by which the theory becomes general relativity (GR) at small scales. We show that at the level of the power spectrum, to a good approximation, themore » effect of nonlinearities in the modified gravity sector may be absorbed into a renormalization of the gravitational constant. Since the relation between the lensing potential and density perturbations is entirely unaffected by the extra physics in these theories, the modified gravity can be described in this approximation by a single function, an effective gravitational constant for nonrelativistic motion that depends on space and time. We develop a resummation scheme to calculate it, and provide predictions for the nonlinear power spectrum. At the level of the large-scale bispectrum, the leading order corrections are obtained by standard perturbation theory techniques, and show that the suppression of the brane-bending mode leads to characteristic signatures in the non-Gaussianity generated by gravity, generic to models that become GR at small scales through second-derivative interactions. We compare the predictions in this work to numerical simulations in a companion paper.« less

Accuracy & Computational Considerations for Wide--Angle One--way Seismic Propagators and Multiple Scattering by Invariant Embedding

NASA Astrophysics Data System (ADS)

Thomson, C. J.

2004-12-01

Pseudodifferential operators (PSDOs) yield in principle exact one--way seismic wave equations, which are attractive both conceptually and for their promise of computational efficiency. The one--way operators can be extended to include multiple--scattering effects, again in principle exactly. In practice approximations must be made and, as an example, the variable--wavespeed Helmholtz equation for scalar waves in two space dimensions is here factorized to give the one--way wave equation. This simple case permits clear identification of a sequence of physically reasonable approximations to be used when the mathematically exact PSDO one--way equation is implemented on a computer. As intuition suggests, these approximations hinge on the medium gradients in the direction transverse to the main propagation direction. A key point is that narrow--angle approximations are to be avoided in the interests of accuracy. Another key consideration stems from the fact that the so--called ``standard--ordering'' PSDO indicates how lateral interpolation of the velocity structure can significantly reduce computational costs associated with the Fourier or plane--wave synthesis lying at the heart of the calculations. The decision on whether a slow or a fast Fourier transform code should be used rests upon how many lateral model parameters are truly distinct. A third important point is that the PSDO theory shows what approximations are necessary in order to generate an exponential one--way propagator for the laterally varying case, representing the intuitive extension of classical integral--transform solutions for a laterally homogeneous medium. This exponential propagator suggests the use of larger discrete step sizes, and it can also be used to approach phase--screen like approximations (though the latter are not the main interest here). Numerical comparisons with finite--difference solutions will be presented in order to assess the approximations being made and to gain an understanding of computation time differences. The ideas described extend to the three--dimensional, generally anisotropic case and to multiple scattering by invariant embedding.

Greater sage-grouse population trends across Wyoming

USGS Publications Warehouse

Edmunds, David; Aldridge, Cameron L.; O'Donnell, Michael; Monroe, Adrian

2018-01-01

The scale at which analyses are performed can have an effect on model results and often one scale does not accurately describe the ecological phenomena of interest (e.g., population trends) for wide-ranging species: yet, most ecological studies are performed at a single, arbitrary scale. To best determine local and regional trends for greater sage-grouse (Centrocercus urophasianus) in Wyoming, USA, we modeled density-independent and -dependent population growth across multiple spatial scales relevant to management and conservation (Core Areas [habitat encompassing approximately 83% of the sage-grouse population on ∼24% of surface area in Wyoming], local Working Groups [7 regional areas for which groups of local experts are tasked with implementing Wyoming's statewide sage-grouse conservation plan at the local level], Core Area status (Core Area vs. Non-Core Area) by Working Groups, and Core Areas by Working Groups). Our goal was to determine the influence of fine-scale population trends (Core Areas) on larger-scale populations (Working Group Areas). We modeled the natural log of change in population size ( peak M lek counts) by time to calculate the finite rate of population growth (λ) for each population of interest from 1993 to 2015. We found that in general when Core Area status (Core Area vs. Non-Core Area) was investigated by Working Group Area, the 2 populations trended similarly and agreed with the overall trend of the Working Group Area. However, at the finer scale where Core Areas were analyzed separately, Core Areas within the same Working Group Area often trended differently and a few large Core Areas could influence the overall Working Group Area trend and mask trends occurring in smaller Core Areas. Relatively close fine-scale populations of sage-grouse can trend differently, indicating that large-scale trends may not accurately depict what is occurring across the landscape (e.g., local effects of gas and oil fields may be masked by increasing larger populations). 

Perturbation solutions of combustion instability problems

NASA Technical Reports Server (NTRS)

Googerdy, A.; Peddieson, J., Jr.; Ventrice, M.

1979-01-01

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Protein Conformational Dynamics Probed by Single-Molecule Electron Transfer

NASA Astrophysics Data System (ADS)

Yang, Haw; Luo, Guobin; Karnchanaphanurach, Pallop; Louie, Tai-Man; Rech, Ivan; Cova, Sergio; Xun, Luying; Xie, X. Sunney

2003-10-01

Electron transfer is used as a probe for angstrom-scale structural changes in single protein molecules. In a flavin reductase, the fluorescence of flavin is quenched by a nearby tyrosine residue by means of photo-induced electron transfer. By probing the fluorescence lifetime of the single flavin on a photon-by-photon basis, we were able to observe the variation of flavin-tyrosine distance over time. We could then determine the potential of mean force between the flavin and the tyrosine, and a correlation analysis revealed conformational fluctuation at multiple time scales spanning from hundreds of microseconds to seconds. This phenomenon suggests the existence of multiple interconverting conformers related to the fluctuating catalytic reactivity.

Multiscale Modeling in the Clinic: Drug Design and Development

DOE Office of Scientific and Technical Information (OSTI.GOV)

Clancy, Colleen E.; An, Gary; Cannon, William R.

A wide range of length and time scales are relevant to pharmacology, especially in drug development, drug design and drug delivery. Therefore, multi-scale computational modeling and simulation methods and paradigms that advance the linkage of phenomena occurring at these multiple scales have become increasingly important. Multi-scale approaches present in silico opportunities to advance laboratory research to bedside clinical applications in pharmaceuticals research. This is achievable through the capability of modeling to reveal phenomena occurring across multiple spatial and temporal scales, which are not otherwise readily accessible to experimentation. The resultant models, when validated, are capable of making testable predictions tomore » guide drug design and delivery. In this review we describe the goals, methods, and opportunities of multi-scale modeling in drug design and development. We demonstrate the impact of multiple scales of modeling in this field. We indicate the common mathematical techniques employed for multi-scale modeling approaches used in pharmacology and present several examples illustrating the current state-of-the-art regarding drug development for: Excitable Systems (Heart); Cancer (Metastasis and Differentiation); Cancer (Angiogenesis and Drug Targeting); Metabolic Disorders; and Inflammation and Sepsis. We conclude with a focus on barriers to successful clinical translation of drug development, drug design and drug delivery multi-scale models.« less

A multiple scales approach to sound generation by vibrating bodies

NASA Technical Reports Server (NTRS)

Geer, James F.; Pope, Dennis S.

1992-01-01

The problem of determining the acoustic field in an inviscid, isentropic fluid generated by a solid body whose surface executes prescribed vibrations is formulated and solved as a multiple scales perturbation problem, using the Mach number M based on the maximum surface velocity as the perturbation parameter. Following the idea of multiple scales, new 'slow' spacial scales are introduced, which are defined as the usual physical spacial scale multiplied by powers of M. The governing nonlinear differential equations lead to a sequence of linear problems for the perturbation coefficient functions. However, it is shown that the higher order perturbation functions obtained in this manner will dominate the lower order solutions unless their dependence on the slow spacial scales is chosen in a certain manner. In particular, it is shown that the perturbation functions must satisfy an equation similar to Burgers' equation, with a slow spacial scale playing the role of the time-like variable. The method is illustrated by a simple one-dimenstional example, as well as by three different cases of a vibrating sphere. The results are compared with solutions obtained by purely numerical methods and some insights provided by the perturbation approach are discussed.

Concurrent temporal channels for auditory processing: Oscillatory neural entrainment reveals segregation of function at different scales

PubMed Central

Tian, Xing; Rowland, Jess; Poeppel, David

2017-01-01

Natural sounds convey perceptually relevant information over multiple timescales, and the necessary extraction of multi-timescale information requires the auditory system to work over distinct ranges. The simplest hypothesis suggests that temporal modulations are encoded in an equivalent manner within a reasonable intermediate range. We show that the human auditory system selectively and preferentially tracks acoustic dynamics concurrently at 2 timescales corresponding to the neurophysiological theta band (4–7 Hz) and gamma band ranges (31–45 Hz) but, contrary to expectation, not at the timescale corresponding to alpha (8–12 Hz), which has also been found to be related to auditory perception. Listeners heard synthetic acoustic stimuli with temporally modulated structures at 3 timescales (approximately 190-, approximately 100-, and approximately 30-ms modulation periods) and identified the stimuli while undergoing magnetoencephalography recording. There was strong intertrial phase coherence in the theta band for stimuli of all modulation rates and in the gamma band for stimuli with corresponding modulation rates. The alpha band did not respond in a similar manner. Classification analyses also revealed that oscillatory phase reliably tracked temporal dynamics but not equivalently across rates. Finally, mutual information analyses quantifying the relation between phase and cochlear-scaled correlations also showed preferential processing in 2 distinct regimes, with the alpha range again yielding different patterns. The results support the hypothesis that the human auditory system employs (at least) a 2-timescale processing mode, in which lower and higher perceptual sampling scales are segregated by an intermediate temporal regime in the alpha band that likely reflects different underlying computations. PMID:29095816

Simulated fault injection - A methodology to evaluate fault tolerant microprocessor architectures

NASA Technical Reports Server (NTRS)

Choi, Gwan S.; Iyer, Ravishankar K.; Carreno, Victor A.

1990-01-01

A simulation-based fault-injection method for validating fault-tolerant microprocessor architectures is described. The approach uses mixed-mode simulation (electrical/logic analysis), and injects transient errors in run-time to assess the resulting fault impact. As an example, a fault-tolerant architecture which models the digital aspects of a dual-channel real-time jet-engine controller is used. The level of effectiveness of the dual configuration with respect to single and multiple transients is measured. The results indicate 100 percent coverage of single transients. Approximately 12 percent of the multiple transients affect both channels; none result in controller failure since two additional levels of redundancy exist.

Spatiotemporal analysis of land use and land cover change in the Brazilian Amazon

PubMed Central

Li, Guiying; Moran, Emilio; Hetrick, Scott

2013-01-01

This paper provides a comparative analysis of land use and land cover (LULC) changes among three study areas with different biophysical environments in the Brazilian Amazon at multiple scales, from per-pixel, polygon, census sector, to study area. Landsat images acquired in the years of 1990/1991, 1999/2000, and 2008/2010 were used to examine LULC change trajectories with the post-classification comparison approach. A classification system composed of six classes – forest, savanna, other-vegetation (secondary succession and plantations), agro-pasture, impervious surface, and water, was designed for this study. A hierarchical-based classification method was used to classify Landsat images into thematic maps. This research shows different spatiotemporal change patterns, composition and rates among the three study areas and indicates the importance of analyzing LULC change at multiple scales. The LULC change analysis over time for entire study areas provides an overall picture of change trends, but detailed change trajectories and their spatial distributions can be better examined at a per-pixel scale. The LULC change at the polygon scale provides the information of the changes in patch sizes over time, while the LULC change at census sector scale gives new insights on how human-induced activities (e.g., urban expansion, roads, and land use history) affect LULC change patterns and rates. This research indicates the necessity to implement change detection at multiple scales for better understanding the mechanisms of LULC change patterns and rates. PMID:24127130

Hybrid Grid and Basis Set Approach to Quantum Chemistry DMRG

NASA Astrophysics Data System (ADS)

Stoudenmire, Edwin Miles; White, Steven

We present a new approach for using DMRG for quantum chemistry that combines the advantages of a basis set with that of a grid approximation. Because DMRG scales linearly for quasi-one-dimensional systems, it is feasible to approximate the continuum with a fine grid in one direction while using a standard basis set approach for the transverse directions. Compared to standard basis set methods, we reach larger systems and achieve better scaling when approaching the basis set limit. The flexibility and reduced costs of our approach even make it feasible to incoporate advanced DMRG techniques such as simulating real-time dynamics. Supported by the Simons Collaboration on the Many-Electron Problem.

Subparsec-scale structure and evolution of Centaurus A (NGC5128).

PubMed Central

Jauncey, D L; Tingay, S J; Preston, R A; Reynolds, J E; Lovell, J E; McCulloch, P M; Tzioumis, A K; Costa, M E; Murphy, D W; Meier, D L; Jones, D L; Amy, S W; Biggs, J D; Blair, D G; Clay, R W; Edwards, P G; Ellingsen, S P; Ferris, R H; Gough, R G; Harbison, P; Jones, P A; King, E A; Kemball, A J; Migenes, V; Nicolson, G D; Sinclair, M W; Van Ommen, T; Wark, R M; White, G L

1995-01-01

We present a series of 8.4-GHz very-long-baseline radio interferometry images of the nucleus of Centaurus A (NGC5128) made with a Southern Hemisphere array, representing a 3.3-year monitoring effort. The nuclear radio jet is approximately 50 milliarcseconds in extent, or at the 3.5-megaparsec distance of NGC5128, approximately 1 parsec in length. Subluminal motion is seen and structural changes are observed on time scales shorter than 4 months. High-resolution observations at 4.8 and 8.4 GHz made in November 1992 reveal a complex morphology and allow us to unambiguously identify the self-absorbed core located at the southwestern end of the jet. PMID:11607599

Subparsec-scale structure and evolution of Centaurus A (NGC5128).

PubMed

Jauncey, D L; Tingay, S J; Preston, R A; Reynolds, J E; Lovell, J E; McCulloch, P M; Tzioumis, A K; Costa, M E; Murphy, D W; Meier, D L; Jones, D L; Amy, S W; Biggs, J D; Blair, D G; Clay, R W; Edwards, P G; Ellingsen, S P; Ferris, R H; Gough, R G; Harbison, P; Jones, P A; King, E A; Kemball, A J; Migenes, V; Nicolson, G D; Sinclair, M W; Van Ommen, T; Wark, R M; White, G L

1995-12-05

We present a series of 8.4-GHz very-long-baseline radio interferometry images of the nucleus of Centaurus A (NGC5128) made with a Southern Hemisphere array, representing a 3.3-year monitoring effort. The nuclear radio jet is approximately 50 milliarcseconds in extent, or at the 3.5-megaparsec distance of NGC5128, approximately 1 parsec in length. Subluminal motion is seen and structural changes are observed on time scales shorter than 4 months. High-resolution observations at 4.8 and 8.4 GHz made in November 1992 reveal a complex morphology and allow us to unambiguously identify the self-absorbed core located at the southwestern end of the jet.

Multi-scale occupancy estimation and modelling using multiple detection methods

USGS Publications Warehouse

Nichols, James D.; Bailey, Larissa L.; O'Connell, Allan F.; Talancy, Neil W.; Grant, Evan H. Campbell; Gilbert, Andrew T.; Annand, Elizabeth M.; Husband, Thomas P.; Hines, James E.

2008-01-01

Occupancy estimation and modelling based on detection–nondetection data provide an effective way of exploring change in a species’ distribution across time and space in cases where the species is not always detected with certainty. Today, many monitoring programmes target multiple species, or life stages within a species, requiring the use of multiple detection methods. When multiple methods or devices are used at the same sample sites, animals can be detected by more than one method.We develop occupancy models for multiple detection methods that permit simultaneous use of data from all methods for inference about method-specific detection probabilities. Moreover, the approach permits estimation of occupancy at two spatial scales: the larger scale corresponds to species’ use of a sample unit, whereas the smaller scale corresponds to presence of the species at the local sample station or site.We apply the models to data collected on two different vertebrate species: striped skunks Mephitis mephitis and red salamanders Pseudotriton ruber. For striped skunks, large-scale occupancy estimates were consistent between two sampling seasons. Small-scale occupancy probabilities were slightly lower in the late winter/spring when skunks tend to conserve energy, and movements are limited to males in search of females for breeding. There was strong evidence of method-specific detection probabilities for skunks. As anticipated, large- and small-scale occupancy areas completely overlapped for red salamanders. The analyses provided weak evidence of method-specific detection probabilities for this species.Synthesis and applications. Increasingly, many studies are utilizing multiple detection methods at sampling locations. The modelling approach presented here makes efficient use of detections from multiple methods to estimate occupancy probabilities at two spatial scales and to compare detection probabilities associated with different detection methods. The models can be viewed as another variation of Pollock's robust design and may be applicable to a wide variety of scenarios where species occur in an area but are not always near the sampled locations. The estimation approach is likely to be especially useful in multispecies conservation programmes by providing efficient estimates using multiple detection devices and by providing device-specific detection probability estimates for use in survey design.

Balance as a measurement of fatigue in postcall residents.

PubMed

Cuthbertson, David W; Bershad, Eric M; Sangi-Haghpeykar, Haleh; Cohen, Helen S

2015-02-01

To determine the degree to which balance is impaired by fatigue in postcall medical residents and to assess the validity of posturography in quantifying that impairment. Prospective study. Medical residents underwent the following tests before and immediately after overnight call: 1) computerized dynamic posturography (CDP) using the Equitest force platform (Natus Medical Inc., Clackamas, OR); 2) Clinical Test of Sensory Interaction and Balance, in which the residents stood on a compliant foam mat with closed eyes and flexed and extended their necks; 3) a computerized reaction time test (available at www.humanbenchmark.com). All tests were performed before and after call at approximately the same time of day to exclude circadian rhythm influence on fatigue. The residents also completed questionnaires, including the Stanford Sleepiness Scale, describing their level of fatigue before and after call. A statistically significant difference was found between a resident's pre- and postcall test with regard to reaction time (P = .01), Stanford Sleepiness Scale score (P = .001), and preference for the visual system as measured by CDP (P = .05). Linear regression analysis showed a statistically significant relationship between sleep deprivation and multiple performance variables measured by CDP. With the growing attention paid to resident duty hours, there is an increasing need for research involving fatigue and practical ways to measure it. This study shows that residents who are measurably fatigued (both objectively and subjectively) may have difficulty utilizing vestibular input during quiet standing but can compensate by means of somatosensory and visual input. NA. © 2014 The American Laryngological, Rhinological and Otological Society, Inc.

Long-term urban carbon dioxide observations reveal spatial and temporal dynamics related to urban characteristics and growth

DOE PAGES

Mitchell, Logan E.; Lin, John C.; Bowling, David R.; ...

2018-03-05

Cities are concentrated areas of CO 2 emissions and have become the foci of policies for mitigation actions. However, atmospheric measurement networks suitable for evaluating urban emissions over time are scarce. Here we present a unique long-term (decadal) record of CO 2 mole fractions from five sites across Utah’s metropolitan Salt Lake Valley. We examine “excess” CO 2 above background conditions resulting from local emissions and meteorological conditions. We ascribe CO 2 trends to changes in emissions, since we did not find longterm trends in atmospheric mixing proxies. Three contrasting CO 2 trends emerged across urban types: negative trends atmore » a residentialindustrial site, positive trends at a site surrounded by rapid suburban growth, and relatively constant CO 2 over time at multiple sites in the established, residential, and commercial urban core. Analysis of populationwithin the atmospheric footprints of the different sites reveals approximately equal increases in population influencing the observed CO 2, implying a nonlinear relationshipwith CO 2 emissions: Population growth in rural areas that experienced suburban development was associated with increasing emissions while population growth in the developed urban core was associated with stable emissions. Four state-of-the-art global-scale emission inventories also have a nonlinear relationship with population density across the city; however, in contrast to our observations, they all have nearly constant emissions over time. Our results indicate that decadal scale changes in urban CO 2 emissions are detectable through monitoring networks and constitute a valuable approach to evaluate emission inventories and studies of urban carbon cycles.« less

Long-term urban carbon dioxide observations reveal spatial and temporal dynamics related to urban characteristics and growth

NASA Astrophysics Data System (ADS)

Mitchell, Logan E.; Lin, John C.; Bowling, David R.; Pataki, Diane E.; Strong, Courtenay; Schauer, Andrew J.; Bares, Ryan; Bush, Susan E.; Stephens, Britton B.; Mendoza, Daniel; Mallia, Derek; Holland, Lacey; Gurney, Kevin R.; Ehleringer, James R.

2018-03-01

Cities are concentrated areas of CO2 emissions and have become the foci of policies for mitigation actions. However, atmospheric measurement networks suitable for evaluating urban emissions over time are scarce. Here we present a unique long-term (decadal) record of CO2 mole fractions from five sites across Utah’s metropolitan Salt Lake Valley. We examine “excess” CO2 above background conditions resulting from local emissions and meteorological conditions. We ascribe CO2 trends to changes in emissions, since we did not find long-term trends in atmospheric mixing proxies. Three contrasting CO2 trends emerged across urban types: negative trends at a residential-industrial site, positive trends at a site surrounded by rapid suburban growth, and relatively constant CO2 over time at multiple sites in the established, residential, and commercial urban core. Analysis of population within the atmospheric footprints of the different sites reveals approximately equal increases in population influencing the observed CO2, implying a nonlinear relationship with CO2 emissions: Population growth in rural areas that experienced suburban development was associated with increasing emissions while population growth in the developed urban core was associated with stable emissions. Four state-of-the-art global-scale emission inventories also have a nonlinear relationship with population density across the city; however, in contrast to our observations, they all have nearly constant emissions over time. Our results indicate that decadal scale changes in urban CO2 emissions are detectable through monitoring networks and constitute a valuable approach to evaluate emission inventories and studies of urban carbon cycles.

Long-term urban carbon dioxide observations reveal spatial and temporal dynamics related to urban characteristics and growth

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mitchell, Logan E.; Lin, John C.; Bowling, David R.

Cities are concentrated areas of CO 2 emissions and have become the foci of policies for mitigation actions. However, atmospheric measurement networks suitable for evaluating urban emissions over time are scarce. Here we present a unique long-term (decadal) record of CO 2 mole fractions from five sites across Utah’s metropolitan Salt Lake Valley. We examine “excess” CO 2 above background conditions resulting from local emissions and meteorological conditions. We ascribe CO 2 trends to changes in emissions, since we did not find longterm trends in atmospheric mixing proxies. Three contrasting CO 2 trends emerged across urban types: negative trends atmore » a residentialindustrial site, positive trends at a site surrounded by rapid suburban growth, and relatively constant CO 2 over time at multiple sites in the established, residential, and commercial urban core. Analysis of populationwithin the atmospheric footprints of the different sites reveals approximately equal increases in population influencing the observed CO 2, implying a nonlinear relationshipwith CO 2 emissions: Population growth in rural areas that experienced suburban development was associated with increasing emissions while population growth in the developed urban core was associated with stable emissions. Four state-of-the-art global-scale emission inventories also have a nonlinear relationship with population density across the city; however, in contrast to our observations, they all have nearly constant emissions over time. Our results indicate that decadal scale changes in urban CO 2 emissions are detectable through monitoring networks and constitute a valuable approach to evaluate emission inventories and studies of urban carbon cycles.« less

Direct measurement of local material properties within living embryonic tissues

NASA Astrophysics Data System (ADS)

Serwane, Friedhelm; Mongera, Alessandro; Rowghanian, Payam; Kealhofer, David; Lucio, Adam; Hockenbery, Zachary; Campàs, Otger

The shaping of biological matter requires the control of its mechanical properties across multiple scales, ranging from single molecules to cells and tissues. Despite their relevance, measurements of the mechanical properties of sub-cellular, cellular and supra-cellular structures within living embryos pose severe challenges to existing techniques. We have developed a technique that uses magnetic droplets to measure the mechanical properties of complex fluids, including in situ and in vivo measurements within living embryos ,across multiple length and time scales. By actuating the droplets with magnetic fields and recording their deformation we probe the local mechanical properties, at any length scale we choose by varying the droplets' diameter. We use the technique to determine the subcellular mechanics of individual blastomeres of zebrafish embryos, and bridge the gap to the tissue scale by measuring the local viscosity and elasticity of zebrafish embryonic tissues. Using this technique, we show that embryonic zebrafish tissues are viscoelastic with a fluid-like behavior at long time scales. This technique will enable mechanobiology and mechano-transduction studies in vivo, including the study of diseases correlated with tissue stiffness, such as cancer.

Multiprocessing on supercomputers for computational aerodynamics

NASA Technical Reports Server (NTRS)

Yarrow, Maurice; Mehta, Unmeel B.

1990-01-01

Very little use is made of multiple processors available on current supercomputers (computers with a theoretical peak performance capability equal to 100 MFLOPs or more) in computational aerodynamics to significantly improve turnaround time. The productivity of a computer user is directly related to this turnaround time. In a time-sharing environment, the improvement in this speed is achieved when multiple processors are used efficiently to execute an algorithm. The concept of multiple instructions and multiple data (MIMD) through multi-tasking is applied via a strategy which requires relatively minor modifications to an existing code for a single processor. Essentially, this approach maps the available memory to multiple processors, exploiting the C-FORTRAN-Unix interface. The existing single processor code is mapped without the need for developing a new algorithm. The procedure for building a code utilizing this approach is automated with the Unix stream editor. As a demonstration of this approach, a Multiple Processor Multiple Grid (MPMG) code is developed. It is capable of using nine processors, and can be easily extended to a larger number of processors. This code solves the three-dimensional, Reynolds averaged, thin-layer and slender-layer Navier-Stokes equations with an implicit, approximately factored and diagonalized method. The solver is applied to generic oblique-wing aircraft problem on a four processor Cray-2 computer. A tricubic interpolation scheme is developed to increase the accuracy of coupling of overlapped grids. For the oblique-wing aircraft problem, a speedup of two in elapsed (turnaround) time is observed in a saturated time-sharing environment.

Intrinsic vs. spurious long-range memory in high-frequency records of environmental radioactivity - Critical re-assessment and application to indoor 222Rn concentrations from Coimbra, Portugal

NASA Astrophysics Data System (ADS)

Donner, Reik V.; Potirakis, Stelios M.; Barbosa, Susana M.; Matos, Jose A. O.

2015-04-01

The presence or absence of long-range correlations in environmental radioactivity fluctuations has recently attracted considerable interest. Among a multiplicity of practically relevant applications, identifying and disentangling the environmental factors controlling the variable concentrations of the radioactive noble gas Radon is important for estimating its effect on human health and the efficiency of possible measures for reducing the corresponding exposition. In this work, we present a critical re-assessment of a multiplicity of complementary methods that have been previously applied for evaluating the presence of long-range correlations and fractal scaling in environmental Radon variations with a particular focus on the specific properties of the underlying time series. As an illustrative case study, we subsequently re-analyze two high-frequency records of indoor Radon concentrations from Coimbra, Portugal, each of which spans several months of continuous measurements at a high temporal resolution of five minutes. Our results reveal that at the study site, Radon concentrations exhibit complex multi-scale dynamics with qualitatively different properties at different time-scales: (i) essentially white noise in the high-frequency part (up to time-scales of about one hour), (ii) spurious indications of a non-stationary, apparently long-range correlated process (at time scales between hours and one day) arising from marked periodic components probably related to tidal frequencies, and (iii) low-frequency variability indicating a true long-range dependent process, which might be dominated by a response to meteorological drivers. In the presence of such multi-scale variability, common estimators of long-range memory in time series are necessarily prone to fail if applied to the raw data without previous separation of time-scales with qualitatively different dynamics. We emphasize that similar properties can be found in other types of geophysical time series (for example, tide gauge records), calling for a careful application of time series analysis tools when studying such data.

Computer-assisted 3D kinematic analysis of all leg joints in walking insects.

PubMed

Bender, John A; Simpson, Elaine M; Ritzmann, Roy E

2010-10-26

High-speed video can provide fine-scaled analysis of animal behavior. However, extracting behavioral data from video sequences is a time-consuming, tedious, subjective task. These issues are exacerbated where accurate behavioral descriptions require analysis of multiple points in three dimensions. We describe a new computer program written to assist a user in simultaneously extracting three-dimensional kinematics of multiple points on each of an insect's six legs. Digital video of a walking cockroach was collected in grayscale at 500 fps from two synchronized, calibrated cameras. We improved the legs' visibility by painting white dots on the joints, similar to techniques used for digitizing human motion. Compared to manual digitization of 26 points on the legs over a single, 8-second bout of walking (or 106,496 individual 3D points), our software achieved approximately 90% of the accuracy with 10% of the labor. Our experimental design reduced the complexity of the tracking problem by tethering the insect and allowing it to walk in place on a lightly oiled glass surface, but in principle, the algorithms implemented are extensible to free walking. Our software is free and open-source, written in the free language Python and including a graphical user interface for configuration and control. We encourage collaborative enhancements to make this tool both better and widely utilized.

Evaluating the status of individuals and populations: Advantages of multiple approaches and time scales: Chapter 6

USGS Publications Warehouse

Monson, Daniel H.; Bowen, Lizabeth

2015-01-01

Overall, a variety of indices used to measure population status throughout the sea otter’s range have provided insights for understanding the mechanisms driving the trajectory of various sea otter populations, which a single index could not, and we suggest using multiple methods to measure a population’s status at multiple spatial and temporal scales. The work described here also illustrates the usefulness of long-term data sets and/or approaches that can be used to assess population status retrospectively, providing information otherwise not available. While not all systems will be as amenable to using all the approaches presented here, we expect innovative researchers could adapt analogous multi-scale methods to a broad range of habitats and species including apex predators occupying the top trophic levels, which are often of conservation concern.

Approximation of reliabilities for multiple-trait model with maternal effects.

PubMed

Strabel, T; Misztal, I; Bertrand, J K

2001-04-01

Reliabilities for a multiple-trait maternal model were obtained by combining reliabilities obtained from single-trait models. Single-trait reliabilities were obtained using an approximation that supported models with additive and permanent environmental effects. For the direct effect, the maternal and permanent environmental variances were assigned to the residual. For the maternal effect, variance of the direct effect was assigned to the residual. Data included 10,550 birth weight, 11,819 weaning weight, and 3,617 postweaning gain records of Senepol cattle. Reliabilities were obtained by generalized inversion and by using single-trait and multiple-trait approximation methods. Some reliabilities obtained by inversion were negative because inbreeding was ignored in calculating the inverse of the relationship matrix. The multiple-trait approximation method reduced the bias of approximation when compared with the single-trait method. The correlations between reliabilities obtained by inversion and by multiple-trait procedures for the direct effect were 0.85 for birth weight, 0.94 for weaning weight, and 0.96 for postweaning gain. Correlations for maternal effects for birth weight and weaning weight were 0.96 to 0.98 for both approximations. Further improvements can be achieved by refining the single-trait procedures.

Naming Game with Multiple Hearers

NASA Astrophysics Data System (ADS)

Li, Bing; Chen, Guanrong; Chow, Tommy W. S.

2013-05-01

A new model called Naming Game with Multiple Hearers (NGMH) is proposed in this paper. A naming game over a population of individuals aims to reach consensus on the name of an object through pair-wise local interactions among all the individuals. The proposed NGMH model describes the learning process of a new word, in a population with one speaker and multiple hearers, at each interaction towards convergence. The characteristics of NGMH are examined on three types of network topologies, namely ER random-graph network, WS small-world network, and BA scale-free network. Comparative analysis on the convergence time is performed, revealing that the topology with a larger average (node) degree can reach consensus faster than the others over the same population. It is found that, for a homogeneous network, the average degree is the limiting value of the number of hearers, which reduces the individual ability of learning new words, consequently decreasing the convergence time; for a scale-free network, this limiting value is the deviation of the average degree. It is also found that a network with a larger clustering coefficient takes longer time to converge; especially a small-word network with smallest rewiring possibility takes longest time to reach convergence. As more new nodes are being added to scale-free networks with different degree distributions, their convergence time appears to be robust against the network-size variation. Most new findings reported in this paper are different from that of the single-speaker/single-hearer naming games documented in the literature.

Short-term cognitive improvement in schizophrenics treated with typical and atypical neuroleptics.

PubMed

Rollnik, Jens D; Borsutzky, Marthias; Huber, Thomas J; Mogk, Hannu; Seifert, Jürgen; Emrich, Hinderk M; Schneider, Udo

2002-01-01

Atypical neuroleptics seem to be more beneficial than typical ones with respect to long-term neuropsychological functioning. Thus, most studies focus on the long-term effects of neuroleptics. We were interested in whether atypical neuroleptic treatment is also superior to typical drugs over relatively short periods of time. We studied 20 schizophrenic patients [10 males, mean age 35.5 years, mean Brief Psychiatric Rating Scale (BPRS) score at entry 58.9] admitted to our hospital with acute psychotic exacerbation. Nine of them were treated with typical and 11 with atypical neuroleptics. In addition, 14 healthy drug-free subjects (6 males, mean age 31.2 years) were enrolled in the study and compared to the patients. As neuropsychological tools, a divided attention test, the Vienna reaction time test, the Benton visual retention test, digit span and a Multiple Choice Word Fluency Test (MWT-B) were used during the first week after admission, within the third week and before discharge (approximately 3 months). Patients scored significantly worse than healthy controls on nearly all tests (except Vienna reaction time). Clinical ratings [BPRS and Positive and Negative Symptom Scale for Schizophrenia (PANSS)] improved markedly (p < 0.01), without a significant difference between typical and atypical medication. Clinical improvement (PANSS total score) correlated with less mistakes on the Benton test (r = 0.762, p = 0.017) and an improvement on the divided attention task (r = 0.705, p = 0.034). Neuropsychological functioning (explicit memory, p < 0.01; divided attention, p < 0.05) moderately improved for both groups under treatment but without a significant difference between atypical and typical antipsychotic drugs. Over short periods of time (3 months), neuropsychological disturbances in schizophrenia seem to be moderately responsive to both typical and atypical neuroleptics. Copyright 2002 S. Karger AG, Basel

Thermal architecture for the SPIDER flight cryostat

NASA Astrophysics Data System (ADS)

Gudmundsson, J. E.; Ade, P. A. R.; Amiri, M.; Benton, S. J.; Bihary, R.; Bock, J. J.; Bond, J. R.; Bonetti, J. A.; Bryan, S. A.; Burger, B.; Chiang, H. C.; Contaldi, C. R.; Crill, B. P.; Doré, O.; Farhang, M.; Filippini, J.; Fissel, L. M.; Gandilo, N. N.; Golwala, S. R.; Halpern, M.; Hasselfield, M.; Hilton, G.; Holmes, W.; Hristov, V. V.; Irwin, K. D.; Jones, W. C.; Kuo, C. L.; MacTavish, C. J.; Mason, P. V.; Montroy, T. E.; Morford, T. A.; Netterfield, C. B.; O'Dea, D. T.; Rahlin, A. S.; Reintsema, C. D.; Ruhl, J. E.; Runyan, M. C.; Schenker, M. A.; Shariff, J. A.; Soler, J. D.; Trangsrud, A.; Tucker, C.; Tucker, R. S.; Turner, A. D.

2010-07-01

We describe the cryogenic system for SPIDER, a balloon-borne microwave polarimeter that will map 8% of the sky with degree-scale angular resolution. The system consists of a 1284 L liquid helium cryostat and a 16 L capillary-filled superfluid helium tank, which provide base operating temperatures of 4 K and 1.5 K, respectively. Closed-cycle 3He adsorption refrigerators supply sub-Kelvin cooling power to multiple focal planes, which are housed in monochromatic telescope inserts. The main helium tank is suspended inside the vacuum vessel with thermally insulating fiberglass flexures, and shielded from thermal radiation by a combination of two vapor cooled shields and multi-layer insulation. This system allows for an extremely low instrumental background and a hold time in excess of 25 days. The total mass of the cryogenic system, including cryogens, is approximately 1000 kg. This enables conventional long duration balloon flights. We will discuss the design, thermal analysis, and qualification of the cryogenic system.

A compact holographic optical tweezers instrument

NASA Astrophysics Data System (ADS)

Gibson, G. M.; Bowman, R. W.; Linnenberger, A.; Dienerowitz, M.; Phillips, D. B.; Carberry, D. M.; Miles, M. J.; Padgett, M. J.

2012-11-01

Holographic optical tweezers have found many applications including the construction of complex micron-scale 3D structures and the control of tools and probes for position, force, and viscosity measurement. We have developed a compact, stable, holographic optical tweezers instrument which can be easily transported and is compatible with a wide range of microscopy techniques, making it a valuable tool for collaborative research. The instrument measures approximately 30×30×35 cm and is designed around a custom inverted microscope, incorporating a fibre laser operating at 1070 nm. We designed the control software to be easily accessible for the non-specialist, and have further improved its ease of use with a multi-touch iPad interface. A high-speed camera allows multiple trapped objects to be tracked simultaneously. We demonstrate that the compact instrument is stable to 0.5 nm for a 10 s measurement time by plotting the Allan variance of the measured position of a trapped 2 μm silica bead. We also present a range of objects that have been successfully manipulated.

Relationship between long working hours and depression: a 3-year longitudinal study of clerical workers.

PubMed

Amagasa, Takashi; Nakayama, Takeo

2013-08-01

To clarify how long working hours affect the likelihood of current and future depression. Using data from four repeated measurements collected from 218 clerical workers, four models associating work-related factors to the depressive mood scale were established. The final model was constructed after comparing and testing the goodness-of-fit index using structural equation modeling. Multiple logistic regression analysis was also performed. The final model showed the best fit (normed fit index = 0.908; goodness-of-fit index = 0.936; root-mean-square error of approximation = 0.018). Its standardized total effect indicated that long working hours affected depression at the time of evaluation and 1 to 3 years later. The odds ratio for depression risk was 14.7 in employees who were not long-hours overworked according to the initial survey but who were long-hours overworked according to the second survey. Long working hours increase current and future risks of depression.

DEVELOPMENT AND PEER REVIEW OF TIME-TO-EFFECT MODELS FOR THE ANALYSIS OF NEUROTOXICITY AND OTHER TIME DEPENDENT DATA

EPA Science Inventory

Neurobehavioral studies pose unique challenges for dose-response modeling, including small sample size and relatively large intra-subject variation, repeated measurements over time, multiple endpoints with both continuous and ordinal scales, and time dependence of risk characteri...

A theory of cerebellar cortex and adaptive motor control based on two types of universal function approximation capability.

PubMed

Fujita, Masahiko

2016-03-01

Lesions of the cerebellum result in large errors in movements. The cerebellum adaptively controls the strength and timing of motor command signals depending on the internal and external environments of movements. The present theory describes how the cerebellar cortex can control signals for accurate and timed movements. A model network of the cerebellar Golgi and granule cells is shown to be equivalent to a multiple-input (from mossy fibers) hierarchical neural network with a single hidden layer of threshold units (granule cells) that receive a common recurrent inhibition (from a Golgi cell). The weighted sum of the hidden unit signals (Purkinje cell output) is theoretically analyzed regarding the capability of the network to perform two types of universal function approximation. The hidden units begin firing as the excitatory inputs exceed the recurrent inhibition. This simple threshold feature leads to the first approximation theory, and the network final output can be any continuous function of the multiple inputs. When the input is constant, this output becomes stationary. However, when the recurrent unit activity is triggered to decrease or the recurrent inhibition is triggered to increase through a certain mechanism (metabotropic modulation or extrasynaptic spillover), the network can generate any continuous signals for a prolonged period of change in the activity of recurrent signals, as the second approximation theory shows. By incorporating the cerebellar capability of two such types of approximations to a motor system, in which learning proceeds through repeated movement trials with accompanying corrections, accurate and timed responses for reaching the target can be adaptively acquired. Simple models of motor control can solve the motor error vs. sensory error problem, as well as the structural aspects of credit (or error) assignment problem. Two physiological experiments are proposed for examining the delay and trace conditioning of eyelid responses, as well as saccade adaptation, to investigate this novel idea of cerebellar processing. Copyright © 2015 Elsevier Ltd. All rights reserved.

Network Modeling and Energy-Efficiency Optimization for Advanced Machine-to-Machine Sensor Networks

PubMed Central

Jung, Sungmo; Kim, Jong Hyun; Kim, Seoksoo

2012-01-01

Wireless machine-to-machine sensor networks with multiple radio interfaces are expected to have several advantages, including high spatial scalability, low event detection latency, and low energy consumption. Here, we propose a network model design method involving network approximation and an optimized multi-tiered clustering algorithm that maximizes node lifespan by minimizing energy consumption in a non-uniformly distributed network. Simulation results show that the cluster scales and network parameters determined with the proposed method facilitate a more efficient performance compared to existing methods. PMID:23202190

Metamodeling and the Critic-based approach to multi-level optimization.

PubMed

Werbos, Ludmilla; Kozma, Robert; Silva-Lugo, Rodrigo; Pazienza, Giovanni E; Werbos, Paul J

2012-08-01

Large-scale networks with hundreds of thousands of variables and constraints are becoming more and more common in logistics, communications, and distribution domains. Traditionally, the utility functions defined on such networks are optimized using some variation of Linear Programming, such as Mixed Integer Programming (MIP). Despite enormous progress both in hardware (multiprocessor systems and specialized processors) and software (Gurobi) we are reaching the limits of what these tools can handle in real time. Modern logistic problems, for example, call for expanding the problem both vertically (from one day up to several days) and horizontally (combining separate solution stages into an integrated model). The complexity of such integrated models calls for alternative methods of solution, such as Approximate Dynamic Programming (ADP), which provide a further increase in the performance necessary for the daily operation. In this paper, we present the theoretical basis and related experiments for solving the multistage decision problems based on the results obtained for shorter periods, as building blocks for the models and the solution, via Critic-Model-Action cycles, where various types of neural networks are combined with traditional MIP models in a unified optimization system. In this system architecture, fast and simple feed-forward networks are trained to reasonably initialize more complicated recurrent networks, which serve as approximators of the value function (Critic). The combination of interrelated neural networks and optimization modules allows for multiple queries for the same system, providing flexibility and optimizing performance for large-scale real-life problems. A MATLAB implementation of our solution procedure for a realistic set of data and constraints shows promising results, compared to the iterative MIP approach. Copyright © 2012 Elsevier Ltd. All rights reserved.

Stacking Multiple Ion Captures in The High Performance Antiproton Trap (HiPAT)

NASA Technical Reports Server (NTRS)

Martin, James J.; Lewis, Raymond A.; Sims, William H.; Chakrabarti, Suman; Pearson, Boise; Fant, Wallace E.

2004-01-01

The High performance Antiproton Trap (HiPAT) research project was initiated by the Marshall Space Flight Center's propulsion Research Center to examining the fundamental behavior of low energy antiprotons. Stored antiproton would ultimately be used for experimental demonstration of basic propulsive concepts. Matter-antimatter annihilation produces approximately 10(exp 8) MJ/g nearly 10 orders of magnitude more energy per unit mass than chemical based combustion, hence NASA's interest. To achieve containment, HiPAT utilizes a type of electromagnetic bottle know as a Penning trap positioned within an ultrahigh vacuum test section. Recently, the HiPAT hardware configuration has been enhanced to facilitate the capture of multiple normal matter ion burst. This endeavor is often referred to as "stacking" and used to increasing the number of captured particles. A prior normal matter experimental effort, successfully demonstrated the effectiveness of single burst capture. The stacking process is accomplished by manipulating the electric field generated by the confinement electrodes i.e. adjusting the well potential depth. These potential well values are initially configured to maximize the quantity of captured ions per burst; shallow wells with a depth of 100 volt or less (referenced to the incoming ion beam energy) are typically selected. Once captured, a cooling interval is required to reduce the energy of trapped particles below the lower extent of the "trap door" (or leading electrode) ion emitting potential. This is necessary such that a new burst of hot ions can be introduced while preventing those already inside from escaping. The cooling time is driven by a combination of mechanisms such as synchrotron radiation, background gas scattering, and resistive damping in a time scale on the order of minutes. A potential for reducing this hold period is to actively manipulate the electric field shape, using the power supply control system, to produce a deeper potential well. The trapped ions ride down the morphing well to a lower potential energy, leaving margin to once again cycle the "trap door" capturing a new burst while limiting the number of escaping ions. By adjusting the magnitude and timing of the potential well shaping, subsequent shots from the ion beam can be captured on a time scale shorter then that dictated by the standard inject, capture and cool technique. Currently, experimental tests are under way examining the functionally of this system for stacking multiple ions bursts within the HiPAT system.

Large-Scale Flows and Magnetic Fields Produced by Rotating Convection in a Quasi-Geostrophic Model of Planetary Cores

NASA Astrophysics Data System (ADS)

Guervilly, C.; Cardin, P.

2017-12-01

Convection is the main heat transport process in the liquid cores of planets. The convective flows are thought to be turbulent and constrained by rotation (corresponding to high Reynolds numbers Re and low Rossby numbers Ro). Under these conditions, and in the absence of magnetic fields, the convective flows can produce coherent Reynolds stresses that drive persistent large-scale zonal flows. The formation of large-scale flows has crucial implications for the thermal evolution of planets and the generation of large-scale magnetic fields. In this work, we explore this problem with numerical simulations using a quasi-geostrophic approximation to model convective and zonal flows at Re 104 and Ro 10-4 for Prandtl numbers relevant for liquid metals (Pr 0.1). The formation of intense multiple zonal jets strongly affects the convective heat transport, leading to the formation of a mean temperature staircase. We also study the generation of magnetic fields by the quasi-geostrophic flows at low magnetic Prandtl numbers.

MMS Observations of Electron-Scale Filamentary Currents in the Reconnection Exhaust and Near the X Line

NASA Technical Reports Server (NTRS)

Phan, T. D.; Eastwood, J. P.; Cassak, P. A.; Oieroset, M.; Gosling, J. T.; Gershman, D. J.; Mozer, F. S.; Shay, M. A.; Fujimoto, M.; Daughton, W.;

Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging

DOEpatents

Anderson, R.N.; Boulanger, A.; Bagdonas, E.P.; Xu, L.; He, W.

1996-12-17

The invention utilizes 3-D and 4-D seismic surveys as a means of deriving information useful in petroleum exploration and reservoir management. The methods use both single seismic surveys (3-D) and multiple seismic surveys separated in time (4-D) of a region of interest to determine large scale migration pathways within sedimentary basins, and fine scale drainage structure and oil-water-gas regions within individual petroleum producing reservoirs. Such structure is identified using pattern recognition tools which define the regions of interest. The 4-D seismic data sets may be used for data completion for large scale structure where time intervals between surveys do not allow for dynamic evolution. The 4-D seismic data sets also may be used to find variations over time of small scale structure within individual reservoirs which may be used to identify petroleum drainage pathways, oil-water-gas regions and, hence, attractive drilling targets. After spatial orientation, and amplitude and frequency matching of the multiple seismic data sets, High Amplitude Event (HAE) regions consistent with the presence of petroleum are identified using seismic attribute analysis. High Amplitude Regions are grown and interconnected to establish plumbing networks on the large scale and reservoir structure on the small scale. Small scale variations over time between seismic surveys within individual reservoirs are identified and used to identify drainage patterns and bypassed petroleum to be recovered. The location of such drainage patterns and bypassed petroleum may be used to site wells. 22 figs.

Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging

DOEpatents

Anderson, Roger N.; Boulanger, Albert; Bagdonas, Edward P.; Xu, Liqing; He, Wei

1996-01-01

The invention utilizes 3-D and 4-D seismic surveys as a means of deriving information useful in petroleum exploration and reservoir management. The methods use both single seismic surveys (3-D) and multiple seismic surveys separated in time (4-D) of a region of interest to determine large scale migration pathways within sedimentary basins, and fine scale drainage structure and oil-water-gas regions within individual petroleum producing reservoirs. Such structure is identified using pattern recognition tools which define the regions of interest. The 4-D seismic data sets may be used for data completion for large scale structure where time intervals between surveys do not allow for dynamic evolution. The 4-D seismic data sets also may be used to find variations over time of small scale structure within individual reservoirs which may be used to identify petroleum drainage pathways, oil-water-gas regions and, hence, attractive drilling targets. After spatial orientation, and amplitude and frequency matching of the multiple seismic data sets, High Amplitude Event (HAE) regions consistent with the presence of petroleum are identified using seismic attribute analysis. High Amplitude Regions are grown and interconnected to establish plumbing networks on the large scale and reservoir structure on the small scale. Small scale variations over time between seismic surveys within individual reservoirs are identified and used to identify drainage patterns and bypassed petroleum to be recovered. The location of such drainage patterns and bypassed petroleum may be used to site wells.

An Approximate Markov Model for the Wright-Fisher Diffusion and Its Application to Time Series Data.

PubMed

Ferrer-Admetlla, Anna; Leuenberger, Christoph; Jensen, Jeffrey D; Wegmann, Daniel

2016-06-01

The joint and accurate inference of selection and demography from genetic data is considered a particularly challenging question in population genetics, since both process may lead to very similar patterns of genetic diversity. However, additional information for disentangling these effects may be obtained by observing changes in allele frequencies over multiple time points. Such data are common in experimental evolution studies, as well as in the comparison of ancient and contemporary samples. Leveraging this information, however, has been computationally challenging, particularly when considering multilocus data sets. To overcome these issues, we introduce a novel, discrete approximation for diffusion processes, termed mean transition time approximation, which preserves the long-term behavior of the underlying continuous diffusion process. We then derive this approximation for the particular case of inferring selection and demography from time series data under the classic Wright-Fisher model and demonstrate that our approximation is well suited to describe allele trajectories through time, even when only a few states are used. We then develop a Bayesian inference approach to jointly infer the population size and locus-specific selection coefficients with high accuracy and further extend this model to also infer the rates of sequencing errors and mutations. We finally apply our approach to recent experimental data on the evolution of drug resistance in influenza virus, identifying likely targets of selection and finding evidence for much larger viral population sizes than previously reported. Copyright © 2016 by the Genetics Society of America.

Review of subjective measures of human response to aircraft noise

NASA Technical Reports Server (NTRS)

Cawthorn, J. M.; Mayes, W. H.

1976-01-01

The development of aircraft noise rating scales and indexes is reviewed up to the present time. Single event scales, multiple event indexes, and their interrelation with each other, are considered. Research requirements for further refinement and development of aircraft noise rating quantification factors are discussed.

The evolution of structure in the universe from axions

NASA Technical Reports Server (NTRS)

Stecker, F. W.; Shafi, Q.

1982-01-01

A scenario where axions provide the dark matter in the universe is considered. Fluctuations in the axion field density produced by domain walls and strings cause the appearance of axion clumps of masses of order 10 to the 6th power solar mass which most likely collapse to black holes by or at the time that the universe becomes axion dominated at T is approximately 10 eV. These objects form the building blocks for a clustering hierarchy theory of galaxy and supercluster formation on scales up to approximately 10 Mpc and approximately 10 to the 15th power solar mass.

Time-Frequency Analysis of Beach Bacteria Variations and its Implication for Recreational Water Quality Modeling

EPA Science Inventory

This paper explores the potential of time-frequency wavelet analysis in resolving beach bacteria concentration and possible explanatory variables across multiple time scales with temporal information still preserved. The wavelet scalograms of E. coli concentrations and the explan...

Cross-Scale Modelling of Subduction from Minute to Million of Years Time Scale

NASA Astrophysics Data System (ADS)

Sobolev, S. V.; Muldashev, I. A.

2015-12-01

Subduction is an essentially multi-scale process with time-scales spanning from geological to earthquake scale with the seismic cycle in-between. Modelling of such process constitutes one of the largest challenges in geodynamic modelling today.Here we present a cross-scale thermomechanical model capable of simulating the entire subduction process from rupture (1 min) to geological time (millions of years) that employs elasticity, mineral-physics-constrained non-linear transient viscous rheology and rate-and-state friction plasticity. The model generates spontaneous earthquake sequences. The adaptive time-step algorithm recognizes moment of instability and drops the integration time step to its minimum value of 40 sec during the earthquake. The time step is then gradually increased to its maximal value of 5 yr, following decreasing displacement rates during the postseismic relaxation. Efficient implementation of numerical techniques allows long-term simulations with total time of millions of years. This technique allows to follow in details deformation process during the entire seismic cycle and multiple seismic cycles. We observe various deformation patterns during modelled seismic cycle that are consistent with surface GPS observations and demonstrate that, contrary to the conventional ideas, the postseismic deformation may be controlled by viscoelastic relaxation in the mantle wedge, starting within only a few hours after the great (M>9) earthquakes. Interestingly, in our model an average slip velocity at the fault closely follows hyperbolic decay law. In natural observations, such deformation is interpreted as an afterslip, while in our model it is caused by the viscoelastic relaxation of mantle wedge with viscosity strongly varying with time. We demonstrate that our results are consistent with the postseismic surface displacement after the Great Tohoku Earthquake for the day-to-year time range. We will also present results of the modeling of deformation of the upper plate during multiple earthquake cycles at times of hundred thousand and million years and discuss effect of great earthquakes in changing long-term stress field in the upper plate.

Fine-scale detection of pollutants by a benthic marine jellyfish.

PubMed

Epstein, Hannah E; Templeman, Michelle A; Kingsford, Michael J

2016-06-15

Local sources of pollution can vary immensely on small geographic scales and short time frames due to differences in runoff and adjacent land use. This study examined the rate of uptake and retention of trace metals in Cassiopea maremetens, a benthic marine jellyfish, over a short time frame and in the presence of multiple pollutants. This study also validated the ability of C. maremetens to uptake metals in the field. Experimental manipulation demonstrated that metal accumulation in jellyfish tissue began within 24h of exposure to treated water and trended for higher accumulation in the presence of multiple pollutants. C. maremetens was found to uptake trace metals in the field and provide unique signatures among locations. This fine-scale detection and rapid accumulation of metals in jellyfish tissue can have major implications for both biomonitoring and the trophic transfer of pollutants through local ecosystems. Copyright © 2016 Elsevier Ltd. All rights reserved.

Falcon: Visual analysis of large, irregularly sampled, and multivariate time series data in additive manufacturing

DOE PAGES

Steed, Chad A.; Halsey, William; Dehoff, Ryan; ...

2017-02-16

Flexible visual analysis of long, high-resolution, and irregularly sampled time series data from multiple sensor streams is a challenge in several domains. In the field of additive manufacturing, this capability is critical for realizing the full potential of large-scale 3D printers. Here, we propose a visual analytics approach that helps additive manufacturing researchers acquire a deep understanding of patterns in log and imagery data collected by 3D printers. Our specific goals include discovering patterns related to defects and system performance issues, optimizing build configurations to avoid defects, and increasing production efficiency. We introduce Falcon, a new visual analytics system thatmore » allows users to interactively explore large, time-oriented data sets from multiple linked perspectives. Falcon provides overviews, detailed views, and unique segmented time series visualizations, all with adjustable scale options. To illustrate the effectiveness of Falcon at providing thorough and efficient knowledge discovery, we present a practical case study involving experts in additive manufacturing and data from a large-scale 3D printer. The techniques described are applicable to the analysis of any quantitative time series, though the focus of this paper is on additive manufacturing.« less

Falcon: Visual analysis of large, irregularly sampled, and multivariate time series data in additive manufacturing

DOE Office of Scientific and Technical Information (OSTI.GOV)

Steed, Chad A.; Halsey, William; Dehoff, Ryan

Flexible visual analysis of long, high-resolution, and irregularly sampled time series data from multiple sensor streams is a challenge in several domains. In the field of additive manufacturing, this capability is critical for realizing the full potential of large-scale 3D printers. Here, we propose a visual analytics approach that helps additive manufacturing researchers acquire a deep understanding of patterns in log and imagery data collected by 3D printers. Our specific goals include discovering patterns related to defects and system performance issues, optimizing build configurations to avoid defects, and increasing production efficiency. We introduce Falcon, a new visual analytics system thatmore » allows users to interactively explore large, time-oriented data sets from multiple linked perspectives. Falcon provides overviews, detailed views, and unique segmented time series visualizations, all with adjustable scale options. To illustrate the effectiveness of Falcon at providing thorough and efficient knowledge discovery, we present a practical case study involving experts in additive manufacturing and data from a large-scale 3D printer. The techniques described are applicable to the analysis of any quantitative time series, though the focus of this paper is on additive manufacturing.« less

The freshwater transport and dynamics of the western Maine coastal current

USGS Publications Warehouse

Geyer, W.R.; Signell, R.P.; Fong, D.A.; Wang, Jingyuan; Anderson, D.M.; Keafer, B.A.

2004-01-01

Observations in the Gulf of Maine, USA, were used to characterize the freshwater transport, temporal variability and dynamics of the western Maine coastal current. These observations included moored measurements, multiple hydrographic surveys, and drifter releases during April–July of 1993 and 1994. There is a strong seasonal signal in salinity and along-shore velocity of the coastal current, caused by the freshwater inputs of the rivers entering the western Gulf. Surface salinity within the coastal current during the spring freshet is typically 2 psu below ambient, and along-shore currents in the surface layer are directed southwestward at speeds of 0.10–0.20 m s−1, occasionally reaching 0.50 m s−1. The plume thickness is typically 10–20 m in water depths of 50–100 m, thus it is well isolated from the bottom over most of its areal extent. The along-coast freshwater transport within the plume varies considerably due to variations in wind stress, but on time scales of weeks to months it follows the variations of riverine input, with a time lag consistent with the advective velocity. Less than half of the transport of the coastal current is explained by the baroclinic gradient; the barotropic forcing associated with the larger-scale dynamics of the Gulf of Maine accounts for about 60% of the transport. The volume of freshwater transport in the coastal current exceeds the local riverine input of fresh water by 30%, suggesting a significant contribution of freshwater transport from the St. John River, 500 km northeastward. The measurements within the western Maine coastal current, however, indicate a significant decrease in the baroclinic transport of fresh water along the coast, with an e-folding scale of approximately 200 km.

Long-term changes in river-floodplain dynamics: implications for salmonid habitat in the Interior Columbia Basin, USA.

PubMed

Tomlinson, Matthew J; Gergel, Sarah E; Beechie, Timothy J; McClure, Michelle M

2011-07-01

Rivers and their associated floodplains are among the world's most highly altered ecosystems, resulting in billions of dollars in restoration expenditures. Successful restoration of these systems requires information at multiple spatial scales (from localized reaches to broader-scale watersheds), as well as information spanning long time frames. Here, we develop a suite of historical landscape indicators of riverine status, primarily from the perspective of salmonid management, using a case study in the Interior Columbia Basin, Washington, USA. We use a combination of historical and modern aerial photography to quantify changes in land cover and reach type, as well as potential fish habitat within channel and off-channel floodplain areas. As of 1949, 55% of the Wenatchee River floodplain had been converted to agriculture. By 2006, 62% had been modified by anthropogenic development, of which 20% was due to urban expansion. The historical percentage of agricultural land in the watershed and the contemporary percentage of urban area surpass thresholds in land cover associated with deleterious impacts on river systems. In addition, the abundance of reach types associated with the highest quality salmonid habitat (island braided and meandering reaches) has declined due to conversion to straight reach types. The area occupied by fish habitats associated with channel migration (slow/stagnant channels and dry channels) has declined approximately 25-30%. Along highly modified rivers, these habitats have also become increasingly fragmented. Caveats related to visual quality and seasonal timing of historical photographs were important considerations in the interpretation of changes witnessed for headwater island braided systems, as well as for floodplain ponds. Development of rigorous, long-term, multi-scale monitoring techniques is necessary to guide the management and restoration of river-floodplain systems for the diversity of ecosystem services they provide.

Improved scaling of temperature-accelerated dynamics using localization

NASA Astrophysics Data System (ADS)

Shim, Yunsic; Amar, Jacques G.

2016-07-01

While temperature-accelerated dynamics (TAD) is a powerful method for carrying out non-equilibrium simulations of systems over extended time scales, the computational cost of serial TAD increases approximately as N3 where N is the number of atoms. In addition, although a parallel TAD method based on domain decomposition [Y. Shim et al., Phys. Rev. B 76, 205439 (2007)] has been shown to provide significantly improved scaling, the dynamics in such an approach is only approximate while the size of activated events is limited by the spatial decomposition size. Accordingly, it is of interest to develop methods to improve the scaling of serial TAD. As a first step in understanding the factors which determine the scaling behavior, we first present results for the overall scaling of serial TAD and its components, which were obtained from simulations of Ag/Ag(100) growth and Ag/Ag(100) annealing, and compare with theoretical predictions. We then discuss two methods based on localization which may be used to address two of the primary "bottlenecks" to the scaling of serial TAD with system size. By implementing both of these methods, we find that for intermediate system-sizes, the scaling is improved by almost a factor of N1/2. Some additional possible methods to improve the scaling of TAD are also discussed.

AntiClustal: Multiple Sequence Alignment by antipole clustering and linear approximate 1-median computation.

PubMed

Di Pietro, C; Di Pietro, V; Emmanuele, G; Ferro, A; Maugeri, T; Modica, E; Pigola, G; Pulvirenti, A; Purrello, M; Ragusa, M; Scalia, M; Shasha, D; Travali, S; Zimmitti, V

2003-01-01

In this paper we present a new Multiple Sequence Alignment (MSA) algorithm called AntiClusAl. The method makes use of the commonly use idea of aligning homologous sequences belonging to classes generated by some clustering algorithm, and then continue the alignment process ina bottom-up way along a suitable tree structure. The final result is then read at the root of the tree. Multiple sequence alignment in each cluster makes use of the progressive alignment with the 1-median (center) of the cluster. The 1-median of set S of sequences is the element of S which minimizes the average distance from any other sequence in S. Its exact computation requires quadratic time. The basic idea of our proposed algorithm is to make use of a simple and natural algorithmic technique based on randomized tournaments which has been successfully applied to large size search problems in general metric spaces. In particular a clustering algorithm called Antipole tree and an approximate linear 1-median computation are used. Our algorithm compared with Clustal W, a widely used tool to MSA, shows a better running time results with fully comparable alignment quality. A successful biological application showing high aminoacid conservation during evolution of Xenopus laevis SOD2 is also cited.

Unified double- and single-sided homogeneous Green’s function representations

PubMed Central

van der Neut, Joost; Slob, Evert

2016-01-01

In wave theory, the homogeneous Green’s function consists of the impulse response to a point source, minus its time-reversal. It can be represented by a closed boundary integral. In many practical situations, the closed boundary integral needs to be approximated by an open boundary integral because the medium of interest is often accessible from one side only. The inherent approximations are acceptable as long as the effects of multiple scattering are negligible. However, in case of strongly inhomogeneous media, the effects of multiple scattering can be severe. We derive double- and single-sided homogeneous Green’s function representations. The single-sided representation applies to situations where the medium can be accessed from one side only. It correctly handles multiple scattering. It employs a focusing function instead of the backward propagating Green’s function in the classical (double-sided) representation. When reflection measurements are available at the accessible boundary of the medium, the focusing function can be retrieved from these measurements. Throughout the paper, we use a unified notation which applies to acoustic, quantum-mechanical, electromagnetic and elastodynamic waves. We foresee many interesting applications of the unified single-sided homogeneous Green’s function representation in holographic imaging and inverse scattering, time-reversed wave field propagation and interferometric Green’s function retrieval. PMID:27436983

Unified double- and single-sided homogeneous Green's function representations

NASA Astrophysics Data System (ADS)

Wapenaar, Kees; van der Neut, Joost; Slob, Evert

2016-06-01

In wave theory, the homogeneous Green's function consists of the impulse response to a point source, minus its time-reversal. It can be represented by a closed boundary integral. In many practical situations, the closed boundary integral needs to be approximated by an open boundary integral because the medium of interest is often accessible from one side only. The inherent approximations are acceptable as long as the effects of multiple scattering are negligible. However, in case of strongly inhomogeneous media, the effects of multiple scattering can be severe. We derive double- and single-sided homogeneous Green's function representations. The single-sided representation applies to situations where the medium can be accessed from one side only. It correctly handles multiple scattering. It employs a focusing function instead of the backward propagating Green's function in the classical (double-sided) representation. When reflection measurements are available at the accessible boundary of the medium, the focusing function can be retrieved from these measurements. Throughout the paper, we use a unified notation which applies to acoustic, quantum-mechanical, electromagnetic and elastodynamic waves. We foresee many interesting applications of the unified single-sided homogeneous Green's function representation in holographic imaging and inverse scattering, time-reversed wave field propagation and interferometric Green's function retrieval.

Evaluation of scaling invariance embedded in short time series.

PubMed

Pan, Xue; Hou, Lei; Stephen, Mutua; Yang, Huijie; Zhu, Chenping

2014-01-01

Scaling invariance of time series has been making great contributions in diverse research fields. But how to evaluate scaling exponent from a real-world series is still an open problem. Finite length of time series may induce unacceptable fluctuation and bias to statistical quantities and consequent invalidation of currently used standard methods. In this paper a new concept called correlation-dependent balanced estimation of diffusion entropy is developed to evaluate scale-invariance in very short time series with length ~10(2). Calculations with specified Hurst exponent values of 0.2,0.3,...,0.9 show that by using the standard central moving average de-trending procedure this method can evaluate the scaling exponents for short time series with ignorable bias (≤0.03) and sharp confidential interval (standard deviation ≤0.05). Considering the stride series from ten volunteers along an approximate oval path of a specified length, we observe that though the averages and deviations of scaling exponents are close, their evolutionary behaviors display rich patterns. It has potential use in analyzing physiological signals, detecting early warning signals, and so on. As an emphasis, the our core contribution is that by means of the proposed method one can estimate precisely shannon entropy from limited records.

Evaluation of Scaling Invariance Embedded in Short Time Series

PubMed Central

Pan, Xue; Hou, Lei; Stephen, Mutua; Yang, Huijie; Zhu, Chenping

2014-01-01

Scaling invariance of time series has been making great contributions in diverse research fields. But how to evaluate scaling exponent from a real-world series is still an open problem. Finite length of time series may induce unacceptable fluctuation and bias to statistical quantities and consequent invalidation of currently used standard methods. In this paper a new concept called correlation-dependent balanced estimation of diffusion entropy is developed to evaluate scale-invariance in very short time series with length . Calculations with specified Hurst exponent values of show that by using the standard central moving average de-trending procedure this method can evaluate the scaling exponents for short time series with ignorable bias () and sharp confidential interval (standard deviation ). Considering the stride series from ten volunteers along an approximate oval path of a specified length, we observe that though the averages and deviations of scaling exponents are close, their evolutionary behaviors display rich patterns. It has potential use in analyzing physiological signals, detecting early warning signals, and so on. As an emphasis, the our core contribution is that by means of the proposed method one can estimate precisely shannon entropy from limited records. PMID:25549356

3D-structured illumination microscopy reveals clustered DNA double-strand break formation in widespread γH2AX foci after high LET heavy-ion particle radiation.

PubMed

Hagiwara, Yoshihiko; Niimi, Atsuko; Isono, Mayu; Yamauchi, Motohiro; Yasuhara, Takaaki; Limsirichaikul, Siripan; Oike, Takahiro; Sato, Hiro; Held, Kathryn D; Nakano, Takashi; Shibata, Atsushi

2017-12-12

DNA double-strand breaks (DSBs) induced by ionising radiation are considered the major cause of genotoxic mutations and cell death. While DSBs are dispersed throughout chromatin after X-rays or γ-irradiation, multiple types of DNA damage including DSBs, single-strand breaks and base damage can be generated within 1-2 helical DNA turns, defined as a complex DNA lesion, after high Linear Energy Transfer (LET) particle irradiation. In addition to the formation of complex DNA lesions, recent evidence suggests that multiple DSBs can be closely generated along the tracks of high LET particle irradiation. Herein, by using three dimensional (3D)-structured illumination microscopy, we identified the formation of 3D widespread γH2AX foci after high LET carbon-ion irradiation. The large γH2AX foci in G 2 -phase cells encompassed multiple foci of replication protein A (RPA), a marker of DSBs undergoing resection during homologous recombination. Furthermore, we demonstrated by 3D analysis that the distance between two individual RPA foci within γH2AX foci was approximately 700 nm. Together, our findings suggest that high LET heavy-ion particles induce clustered DSB formation on a scale of approximately 1 μm 3 . These closely localised DSBs are considered to be a risk for the formation of chromosomal rearrangement after heavy-ion irradiation.

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

NASA Technical Reports Server (NTRS)

Buchert, T.; Melott, A. L.; Weiss, A. G.

1993-01-01

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Multiscale Granger causality

NASA Astrophysics Data System (ADS)

Faes, Luca; Nollo, Giandomenico; Stramaglia, Sebastiano; Marinazzo, Daniele

2017-10-01

In the study of complex physical and biological systems represented by multivariate stochastic processes, an issue of great relevance is the description of the system dynamics spanning multiple temporal scales. While methods to assess the dynamic complexity of individual processes at different time scales are well established, multiscale analysis of directed interactions has never been formalized theoretically, and empirical evaluations are complicated by practical issues such as filtering and downsampling. Here we extend the very popular measure of Granger causality (GC), a prominent tool for assessing directed lagged interactions between joint processes, to quantify information transfer across multiple time scales. We show that the multiscale processing of a vector autoregressive (AR) process introduces a moving average (MA) component, and describe how to represent the resulting ARMA process using state space (SS) models and to combine the SS model parameters for computing exact GC values at arbitrarily large time scales. We exploit the theoretical formulation to identify peculiar features of multiscale GC in basic AR processes, and demonstrate with numerical simulations the much larger estimation accuracy of the SS approach compared to pure AR modeling of filtered and downsampled data. The improved computational reliability is exploited to disclose meaningful multiscale patterns of information transfer between global temperature and carbon dioxide concentration time series, both in paleoclimate and in recent years.

Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO 2 in fractured rocks

DOE PAGES

Zhou, Quanlin; Oldenburg, Curtis M.; Spangler, Lee H.; ...

2017-01-05

Analytical solutions with infinite exponential series are available to calculate the rate of diffusive transfer between low-permeability blocks and high-permeability zones in the subsurface. Truncation of these series is often employed by neglecting the early-time regime. Here in this paper, we present unified-form approximate solutions in which the early-time and the late-time solutions are continuous at a switchover time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the first coefficient dependent only on the dimensionless area-to-volume ratio. The last two coefficients are either determined analytically for isotropic blocks (e.g., spheresmore » and slabs) or obtained by fitting the exact solutions, and they solely depend on the aspect ratios for rectangular columns and parallelepipeds. For the late-time solutions, only the leading exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic rectangular blocks. The optimal switchover time is between 0.157 and 0.229, with highest relative approximation error less than 0.2%. The solutions are used to demonstrate the storage of dissolved CO 2 in fractured reservoirs with low-permeability matrix blocks of single and multiple shapes and sizes. These approximate solutions are building blocks for development of analytical and numerical tools for hydraulic, solute, and thermal diffusion processes in low-permeability matrix blocks.« less

Robust decentralized hybrid adaptive output feedback fuzzy control for a class of large-scale MIMO nonlinear systems and its application to AHS.

PubMed

Huang, Yi-Shao; Liu, Wel-Ping; Wu, Min; Wang, Zheng-Wu

2014-09-01

This paper presents a novel observer-based decentralized hybrid adaptive fuzzy control scheme for a class of large-scale continuous-time multiple-input multiple-output (MIMO) uncertain nonlinear systems whose state variables are unmeasurable. The scheme integrates fuzzy logic systems, state observers, and strictly positive real conditions to deal with three issues in the control of a large-scale MIMO uncertain nonlinear system: algorithm design, controller singularity, and transient response. Then, the design of the hybrid adaptive fuzzy controller is extended to address a general large-scale uncertain nonlinear system. It is shown that the resultant closed-loop large-scale system keeps asymptotically stable and the tracking error converges to zero. The better characteristics of our scheme are demonstrated by simulations. Copyright © 2014. Published by Elsevier Ltd.

Small scale noise and wind tunnel tests of upper surface blowing nozzle flap concepts. Volume 1. Aerodynamic test results

NASA Technical Reports Server (NTRS)

Renselaer, D. J.; Nishida, R. S.; Wilkin, C. A.

1975-01-01

The results and analyses of aerodynamic and acoustic studies conducted on the small scale noise and wind tunnel tests of upper surface blowing nozzle flap concepts are presented. Various types of nozzle flap concepts were tested. These are an upper surface blowing concept with a multiple slot arrangement with seven slots (seven slotted nozzle), an upper surface blowing type with a large nozzle exit at approximately mid-chord location in conjunction with a powered trailing edge flap with multiple slots (split flow or partially slotted nozzle). In addition, aerodynamic tests were continued on a similar multi-slotted nozzle flap, but with 14 slots. All three types of nozzle flap concepts tested appear to be about equal in overall aerodynamic performance but with the split flow nozzle somewhat better than the other two nozzle flaps in the landing approach mode. All nozzle flaps can be deflected to a large angle to increase drag without significant loss in lift. The nozzle flap concepts appear to be viable aerodynamic drag modulation devices for landing.

Analysis of remote sensing data for evaluation of vegetation resources

NASA Technical Reports Server (NTRS)

1970-01-01

Research has centered around: (1) completion of a study on the use of remote sensing techniques as an aid to multiple use management; (2) determination of the information transfer at various image resolution levels for wildland areas; and (3) determination of the value of small scale multiband, multidate photography for the analysis of vegetation resources. In addition, a substantial effort was made to upgrade the automatic image classification and spectral signature acquisition capabilities of the laboratory. It was found that: (1) Remote sensing techniques should be useful in multiple use management to provide a first-cut analysis of an area. (2) Imagery with 400-500 feet ground resolvable distance (GRD), such as that expected from ERTS-1, should allow discriminations to be made between woody vegetation, grassland, and water bodies with approximately 80% accuracy. (3) Barley and wheat acreages in Maricopa County, Arizona could be estimated with acceptable accuracies using small scale multiband, multidate photography. Sampling errors for acreages of wheat, barley, small grains (wheat and barley combined), and all cropland were 13%, 11%, 8% and 3% respectively.

Dietary human exposure to mercury in two artisanal small-scale gold mining communities of northwestern Colombia.

PubMed

Salazar-Camacho, Carlos; Salas-Moreno, Manuel; Marrugo-Madrid, Siday; Marrugo-Negrete, José; Díez, Sergi

2017-10-01

Artisanal and small-scale gold mining (ASGM) is the largest anthropogenic source of mercury pollution worldwide, posing a grave threat to human health. The present study identifies current levels of mercury in the human population from mining areas of the Chocó Department, Colombia, through total mercury (THg) and methylmercury (MeHg) measurements in human hair. Mercury exposure of the local population was assessed in two towns affected by ASGM and was related to different variables of interest. Concentrations of THg in human hair ranged from 0.06 to 17.54ppm and the mean value for the subjects under study was 2.48ppm. Men had significantly higher levels than women in both towns (3.29ppm vs. 0.77ppm). Fish consumption was related to a marked increase of THg in hair, with mean values close to five times higher in frequent fish consumers (5-7 times/week) than in non-fish consumers (4.80ppm vs. 0.90ppm). A multiple linear regression model was fitted successfully (R=0.671) and reveals that gender, fish consumption and location of residence were significant indicators of Hg levels in hair, while no significant relationship was found for age. Approximately 60% of subjects tested had THg levels that exceeded the U.S. Environmental Protection Agency reference dose of 1.0ppm, while 25% surpassed that of the World Health Organization (2.2ppm). Copyright © 2017 Elsevier Ltd. All rights reserved.

Using analogy to learn about phenomena at scales outside human perception.

PubMed

Resnick, Ilyse; Davatzes, Alexandra; Newcombe, Nora S; Shipley, Thomas F

2017-01-01

Understanding and reasoning about phenomena at scales outside human perception (for example, geologic time) is critical across science, technology, engineering, and mathematics. Thus, devising strong methods to support acquisition of reasoning at such scales is an important goal in science, technology, engineering, and mathematics education. In two experiments, we examine the use of analogical principles in learning about geologic time. Across both experiments we find that using a spatial analogy (for example, a time line) to make multiple alignments, and keeping all unrelated components of the analogy held constant (for example, keep the time line the same length), leads to better understanding of the magnitude of geologic time. Effective approaches also include hierarchically and progressively aligning scale information (Experiment 1) and active prediction in making alignments paired with immediate feedback (Experiments 1 and 2).

Local conformational dynamics in alpha-helices measured by fast triplet transfer.

PubMed

Fierz, Beat; Reiner, Andreas; Kiefhaber, Thomas

2009-01-27

Coupling fast triplet-triplet energy transfer (TTET) between xanthone and naphthylalanine to the helix-coil equilibrium in alanine-based peptides allowed the observation of local equilibrium fluctuations in alpha-helices on the nanoseconds to microseconds time scale. The experiments revealed faster helix unfolding in the terminal regions compared with the central parts of the helix with time constants varying from 250 ns to 1.4 micros at 5 degrees C. Local helix formation occurs with a time constant of approximately 400 ns, independent of the position in the helix. Comparing the experimental data with simulations using a kinetic Ising model showed that the experimentally observed dynamics can be explained by a 1-dimensional boundary diffusion with position-independent elementary time constants of approximately 50 ns for the addition and of approximately 65 ns for the removal of an alpha-helical segment. The elementary time constant for helix growth agrees well with previously measured time constants for formation of short loops in unfolded polypeptide chains, suggesting that helix elongation is mainly limited by a conformational search.

Solitary waves in shallow water hydrodynamics and magnetohydrodynamics in rotating spherical coordinates

NASA Astrophysics Data System (ADS)

London, Steven D.

2018-01-01

In a recent paper (London, Geophys. Astrophys. Fluid Dyn. 2017, vol. 111, pp. 115-130, referred to as L1), we considered a perfect electrically conducting rotating fluid in the presence of an ambient toroidal magnetic field, governed by the shallow water magnetohydrodynamic (MHD) equations in a modified equatorial ?-plane approximation. In conjunction with a WKB type approximation, we used a multiple scale asymptotic scheme, previously developed by Boyd (J. Phys. Oceanogr. 1980, vol. 10, pp. 1699-1717) for equatorial solitary hydrodynamic waves, and found solitary MHD waves. In this paper, as in L1, we apply a WKB type approximation in order to extend the results of L1 from the modified ?-plane to the full spherical geometry. We have included differential rotation in the analysis in order to make the results more relevant to the solar case. In addition, we consider the case of hydrodynamic waves on the rotating sphere in the presence of a differential rotation intended to roughly model the varying large scale currents in the oceans and atmosphere. In the hydrodynamic case, we find the usual equatorial solitary waves as found by Boyd, as well as waves in bands away from the equator for sufficiently strong currents. In the MHD case, we find basically the same equatorial waves found in L1. L1 also found non-equatorial modes; no such modes are found in the full spherical geometry.

A scale-invariant internal representation of time.

PubMed

Shankar, Karthik H; Howard, Marc W

2012-01-01

We propose a principled way to construct an internal representation of the temporal stimulus history leading up to the present moment. A set of leaky integrators performs a Laplace transform on the stimulus function, and a linear operator approximates the inversion of the Laplace transform. The result is a representation of stimulus history that retains information about the temporal sequence of stimuli. This procedure naturally represents more recent stimuli more accurately than less recent stimuli; the decrement in accuracy is precisely scale invariant. This procedure also yields time cells that fire at specific latencies following the stimulus with a scale-invariant temporal spread. Combined with a simple associative memory, this representation gives rise to a moment-to-moment prediction that is also scale invariant in time. We propose that this scale-invariant representation of temporal stimulus history could serve as an underlying representation accessible to higher-level behavioral and cognitive mechanisms. In order to illustrate the potential utility of this scale-invariant representation in a variety of fields, we sketch applications using minimal performance functions to problems in classical conditioning, interval timing, scale-invariant learning in autoshaping, and the persistence of the recency effect in episodic memory across timescales.

The temporal representation of the delay of dynamic iterated rippled noise with positive and negative gain by single units in the ventral cochlear nucleus.

PubMed

Sayles, Mark; Winter, Ian Michael

2007-09-26

Spike trains were recorded from single units in the ventral cochlear nucleus of the anaesthetised guinea-pig in response to dynamic iterated rippled noise with positive and negative gain. The short-term running waveform autocorrelation functions of these stimuli show peaks at integer multiples of the time-varying delay when the gain is +1, and troughs at odd-integer multiples and peaks at even-integer multiples of the time-varying delay when the gain is -1. In contrast, the short-term autocorrelation of the Hilbert envelope shows peaks at integer multiples of the time-varying delay for both positive and negative gain stimuli. A running short-term all-order interspike interval analysis demonstrates the ability of single units to represent the modulated pitch contour in their short-term interval statistics. For units with low best frequency (approximate < or = 1.1 kHz) the temporal discharge pattern reflected the waveform fine structure regardless of unit classification (Primary-like, Chopper). For higher best frequency units the pattern of response varied according to unit type. Chopper units with best frequency approximate > or = 1.1 kHz responded to envelope modulation; showing no difference between their response to stimuli with positive and negative gain. Primary-like units with best frequencies in the range 1-3 kHz were still able to represent the difference in the temporal fine structure between dynamic rippled noise with positive and negative gain. No unit with a best frequency above 3 kHz showed a response to the temporal fine structure. Chopper units in this high frequency group showed significantly greater representation of envelope modulation relative to primary-like units with the same range of best frequencies. These results show that at the level of the cochlear nucleus there exists sufficient information in the time domain to represent the time-varying pitch associated with dynamic iterated rippled noise.

Trellis Tone Modulation Multiple-Access for Peer Discovery in D2D Networks

PubMed Central

Lim, Chiwoo; Kim, Sang-Hyo

2018-01-01

In this paper, a new non-orthogonal multiple-access scheme, trellis tone modulation multiple-access (TTMMA), is proposed for peer discovery of distributed device-to-device (D2D) communication. The range and capacity of discovery are important performance metrics in peer discovery. The proposed trellis tone modulation uses single-tone transmission and achieves a long discovery range due to its low Peak-to-Average Power Ratio (PAPR). The TTMMA also exploits non-orthogonal resource assignment to increase the discovery capacity. For the multi-user detection of superposed multiple-access signals, a message-passing algorithm with supplementary schemes are proposed. With TTMMA and its message-passing demodulation, approximately 1.5 times the number of devices are discovered compared to the conventional frequency division multiple-access (FDMA)-based discovery. PMID:29673167

Trellis Tone Modulation Multiple-Access for Peer Discovery in D2D Networks.

PubMed

Lim, Chiwoo; Jang, Min; Kim, Sang-Hyo

2018-04-17

In this paper, a new non-orthogonal multiple-access scheme, trellis tone modulation multiple-access (TTMMA), is proposed for peer discovery of distributed device-to-device (D2D) communication. The range and capacity of discovery are important performance metrics in peer discovery. The proposed trellis tone modulation uses single-tone transmission and achieves a long discovery range due to its low Peak-to-Average Power Ratio (PAPR). The TTMMA also exploits non-orthogonal resource assignment to increase the discovery capacity. For the multi-user detection of superposed multiple-access signals, a message-passing algorithm with supplementary schemes are proposed. With TTMMA and its message-passing demodulation, approximately 1.5 times the number of devices are discovered compared to the conventional frequency division multiple-access (FDMA)-based discovery.

MIMO equalization with adaptive step size for few-mode fiber transmission systems.

PubMed

van Uden, Roy G H; Okonkwo, Chigo M; Sleiffer, Vincent A J M; de Waardt, Hugo; Koonen, Antonius M J

2014-01-13

Optical multiple-input multiple-output (MIMO) transmission systems generally employ minimum mean squared error time or frequency domain equalizers. Using an experimental 3-mode dual polarization coherent transmission setup, we show that the convergence time of the MMSE time domain equalizer (TDE) and frequency domain equalizer (FDE) can be reduced by approximately 50% and 30%, respectively. The criterion used to estimate the system convergence time is the time it takes for the MIMO equalizer to reach an average output error which is within a margin of 5% of the average output error after 50,000 symbols. The convergence reduction difference between the TDE and FDE is attributed to the limited maximum step size for stable convergence of the frequency domain equalizer. The adaptive step size requires a small overhead in the form of a lookup table. It is highlighted that the convergence time reduction is achieved without sacrificing optical signal-to-noise ratio performance.

Do patient and proxy agree? Long-term changes in multiple sclerosis physical impact and walking ability on patient-reported outcome scales.

PubMed

Sonder, Judith M; Balk, Lisanne J; Bosma, Libertje V A E; Polman, Chris H; Uitdehaag, Bernard M J

2014-10-01

Patient-reported outcome scales (PROs) are useful in monitoring changes in multiple sclerosis (MS) over time. Although these scales are reliable and valid measures in longitudinal studies in MS patients, it is unknown what the impact is when obtaining longitudinal data from proxies. The objective of this paper is to compare longitudinal changes in patient and proxy responses on PROs assessing physical impact of MS and walking ability. In a prospective observational study, data on the Multiple Sclerosis Impact Scale (MSIS-29 physical) and Multiple Sclerosis Walking Scale (MSWS-12) were obtained from 137 patient-proxy couples at baseline and at two-year follow-up. Demographic and disease-related variables explaining agreement or disagreement between patients and proxies were investigated using linear regression analyses. Full agreement was found in 56% (MSIS) and 62% (MSWS) of the patient-proxy couples. Complete disagreement was very rare for both scales (2% MSIS, 5% MSWS). When patients were more positive than proxies, a higher age, longer disease duration, longer patient-proxy relationship and increased levels of depression, anxiety and caregiver burden in proxies were observed. In the majority of the patient-proxy couples there was agreement. Proxies can serve as a valuable source of information, but caution remains essential when using scores from proxies. © The Author(s), 2014.

Simultaneous small- and wide-angle scattering at high X-ray energies.

PubMed

Daniels, J E; Pontoni, D; Hoo, Rui Ping; Honkimäki, V

2010-07-01

Combined small- and wide-angle X-ray scattering (SAXS/WAXS) is a powerful technique for the study of materials at length scales ranging from atomic/molecular sizes (a few angstroms) to the mesoscopic regime ( approximately 1 nm to approximately 1 microm). A set-up to apply this technique at high X-ray energies (E > 50 keV) has been developed. Hard X-rays permit the execution of at least three classes of investigations that are significantly more difficult to perform at standard X-ray energies (8-20 keV): (i) in situ strain analysis revealing anisotropic strain behaviour both at the atomic (WAXS) as well as at the mesoscopic (SAXS) length scales, (ii) acquisition of WAXS patterns to very large q (>20 A(-1)) thus allowing atomic pair distribution function analysis (SAXS/PDF) of micro- and nano-structured materials, and (iii) utilization of complex sample environments involving thick X-ray windows and/or samples that can be penetrated only by high-energy X-rays. Using the reported set-up a time resolution of approximately two seconds was demonstrated. It is planned to further improve this time resolution in the near future.

keV-Scale sterile neutrino sensitivity estimation with time-of-flight spectroscopy in KATRIN using self-consistent approximate Monte Carlo

NASA Astrophysics Data System (ADS)

Steinbrink, Nicholas M. N.; Behrens, Jan D.; Mertens, Susanne; Ranitzsch, Philipp C.-O.; Weinheimer, Christian

2018-03-01

We investigate the sensitivity of the Karlsruhe Tritium Neutrino Experiment (KATRIN) to keV-scale sterile neutrinos, which are promising dark matter candidates. Since the active-sterile mixing would lead to a second component in the tritium β-spectrum with a weak relative intensity of order sin ^2θ ≲ 10^{-6}, additional experimental strategies are required to extract this small signature and to eliminate systematics. A possible strategy is to run the experiment in an alternative time-of-flight (TOF) mode, yielding differential TOF spectra in contrast to the integrating standard mode. In order to estimate the sensitivity from a reduced sample size, a new analysis method, called self-consistent approximate Monte Carlo (SCAMC), has been developed. The simulations show that an ideal TOF mode would be able to achieve a statistical sensitivity of sin ^2θ ˜ 5 × 10^{-9} at one σ , improving the standard mode by approximately a factor two. This relative benefit grows significantly if additional exemplary systematics are considered. A possible implementation of the TOF mode with existing hardware, called gated filtering, is investigated, which, however, comes at the price of a reduced average signal rate.

Scaling of electromagnetic transducers for shunt damping and energy harvesting

NASA Astrophysics Data System (ADS)

Elliott, Stephen J.; Zilletti, Michele

2014-04-01

In order for an electromagnetic transducer to operate well as either a mechanical shunt damper or as a vibration energy harvester, it must have good electromechanical coupling. A simple two-port analysis is used to derive a non-dimensional measure of electromechanical coupling, which must be large compared with unity for efficient operation in both of these applications. The two-port parameters for an inertial electromagnetic transducer are derived, from which this non-dimensional coupling parameter can be evaluated. The largest value that this parameter takes is approximately equal to the square of the magnetic flux density times the length of wire in the field, divided by the mechanical damping times the electrical resistance. This parameter is found to be only of the order of one for voice coil devices that weigh approximately 1 kg, and so such devices are generally not efficient, within the definition used here, in either of these applications. The non-dimensional coupling parameter is found to scale in approximate proportion to the device's characteristic length, however, and so although miniaturised devices are less efficient, greater efficiency can be obtained with large devices, such as those used to control civil engineering structures.

Changes in channel morphology over human time scales [Chapter 32

Treesearch

John M. Buffington

2012-01-01

Rivers are exposed to changing environmental conditions over multiple spatial and temporal scales, with the imposed environmental conditions and response potential of the river modulated to varying degrees by human activity and our exploitation of natural resources. Watershed features that control river morphology include topography (valley slope and channel...

IRIS Arrays: Observing Wavefields at Multiple Scales and Frequencies

NASA Astrophysics Data System (ADS)

Sumy, D. F.; Woodward, R.; Frassetto, A.

2014-12-01

The Incorporated Research Institutions for Seismology (IRIS) provides instruments for creating and operating seismic arrays at a wide range of scales. As an example, for over thirty years the IRIS PASSCAL program has provided instruments to individual Principal Investigators to deploy arrays of all shapes and sizes on every continent. These arrays have ranged from just a few sensors to hundreds or even thousands of sensors, covering areas with dimensions of meters to thousands of kilometers. IRIS also operates arrays directly, such as the USArray Transportable Array (TA) as part of the EarthScope program. Since 2004, the TA has rolled across North America, at any given time spanning a swath of approximately 800 km by 2,500 km, and thus far sampling 2% of the Earth's surface. This achievement includes all of the lower-48 U.S., southernmost Canada, and now parts of Alaska. IRIS has also facilitated specialized arrays in polar environments and on the seafloor. In all cases, the data from these arrays are freely available to the scientific community. As the community of scientists who use IRIS facilities and data look to the future they have identified a clear need for new array capabilities. In particular, as part of its Wavefields Initiative, IRIS is exploring new technologies that can enable large, dense array deployments to record unaliased wavefields at a wide range of frequencies. Large-scale arrays might utilize multiple sensor technologies to best achieve observing objectives and optimize equipment and logistical costs. Improvements in packaging and power systems can provide equipment with reduced size, weight, and power that will reduce logistical constraints for large experiments, and can make a critical difference for deployments in harsh environments or other situations where rapid deployment is required. We will review the range of existing IRIS array capabilities with an overview of previous and current deployments and examples of data and results. We will review existing IRIS projects that explore new array capabilities and highlight future directions for IRIS instrumentation facilities.

High-accuracy self-mixing interferometer based on multiple reflections using a simple external reflecting mirror

NASA Astrophysics Data System (ADS)

Wang, Xiu-lin; Wei, Zheng; Wang, Rui; Huang, Wen-cai

2018-05-01

A self-mixing interferometer (SMI) with resolution twenty times higher than that of a conventional interferometer is developed by multiple reflections. Only by employing a simple external reflecting mirror, the multiple-pass optical configuration can be constructed. The advantage of the configuration is simple and easy to make the light re-injected back into the laser cavity. Theoretical analysis shows that the resolution of measurement is scalable by adjusting the number of reflections. The experiment shows that the proposed method has the optical resolution of approximate λ/40. The influence of displacement sensitivity gain ( G) is further analyzed and discussed in practical experiments.

Phase Coupling Between Spectral Components of Collapsing Langmuir Solitons in Solar Type III Radio Bursts

NASA Technical Reports Server (NTRS)

Thejappa, G.; MacDowall, R. J.; Bergamo, M.

2012-01-01

We present the high time resolution observations of one of the Langmuir wave packets obtained in the source region of a solar type III radio burst. This wave packet satisfies the threshold condition of the supersonic modulational instability, as well as the criterion of a collapsing Langmuir soliton, i.e., the spatial scale derived from its peak intensity is less than that derived from its short time scale. The spectrum of t his wave packet contains an intense spectral peak at local electron plasma frequency, f(sub pe) and relatively weaker peaks at 2f(sub pe) and 3f(sub pe). We apply the wavelet based bispectral analysis technique on this wave packet and compute the bicoherence between its spectral components. It is found that the bicoherence exhibits two peaks at (approximately f(sub pe), approximately f(sub pe)) and (approximately f(sub pe) approximately 2f(sub pe)), which strongly suggest that the spectral peak at 2f(sub pe) probably corresponds to the second harmonic radio emission, generated as a result of the merging of antiparallel propagating Langmuir waves trapped in the collapsing Langmuir soliton, and, the spectral peak at 3f(sub pe) probably corresponds to the third harmonic radio emission, generated as a result of merging of a trapped Langmuir wave and a second harmonic electromagnetic wave.

Influence of anisotropy on anomalous scaling of a passive scalar advected by the Navier-Stokes velocity field.

PubMed

Jurcisinová, E; Jurcisin, M; Remecký, R

2009-10-01

The influence of weak uniaxial small-scale anisotropy on the stability of the scaling regime and on the anomalous scaling of the single-time structure functions of a passive scalar advected by the velocity field governed by the stochastic Navier-Stokes equation is investigated by the field theoretic renormalization group and operator-product expansion within one-loop approximation of a perturbation theory. The explicit analytical expressions for coordinates of the corresponding fixed point of the renormalization-group equations as functions of anisotropy parameters are found, the stability of the three-dimensional Kolmogorov-like scaling regime is demonstrated, and the dependence of the borderline dimension d(c) is an element of (2,3] between stable and unstable scaling regimes is found as a function of the anisotropy parameters. The dependence of the turbulent Prandtl number on the anisotropy parameters is also briefly discussed. The influence of weak small-scale anisotropy on the anomalous scaling of the structure functions of a passive scalar field is studied by the operator-product expansion and their explicit dependence on the anisotropy parameters is present. It is shown that the anomalous dimensions of the structure functions, which are the same (universal) for the Kraichnan model, for the model with finite time correlations of the velocity field, and for the model with the advection by the velocity field driven by the stochastic Navier-Stokes equation in the isotropic case, can be distinguished by the assumption of the presence of the small-scale anisotropy in the systems even within one-loop approximation. The corresponding comparison of the anisotropic anomalous dimensions for the present model with that obtained within the Kraichnan rapid-change model is done.

Early-Onset Multiple Sclerosis in Isfahan, Iran: Report of the Demographic and Clinical Features of 221 Patients.

PubMed

Etemadifar, Masoud; Nourian, Sayed-Mohammadamin; Nourian, Niloofaralsadat; Abtahi, Seyed-Hossein; Sayahi, Farnaz; Saraf, Zahra; Fereidan-Esfahani, Mahboobeh

2016-06-01

It is estimated that early-onset multiple sclerosis multiple sclerosis (early-onset multiple sclerosis) approximately incorporates 3-5% of the multiple sclerosis population. In this report on early-onset multiple sclerosis, the authors aimed to define demographic, clinical and imaging features in a case-series of true-childhood multiple sclerosis and to compare its characteristics with juvenile multiple sclerosis. The authors inspected the records of multiple sclerosis patients who were registered by Isfahan MS Society. Clinical and demographic data of children with less than 16 years of age were reviewed retrospectively. Out of 4536 multiple sclerosis patients referred to the authors' center, 221 patients (4.8%) had multiple sclerosis starting at the age of 16 or less (11 true-childhood multiple sclerosis vs 210 juvenile-onset multiple sclerosis); the female to male ratio was 4.81:1. In the mean follow-up period of 6.2 years, 22 patients (10.5%) had positive family history of multiple sclerosis, 196 (88.6%) patients were classified as relapsing-remitting multiple sclerosis, the mean (± SD Expanded Disability Status Scale) was 1.5 ± 1.1 at the last evaluation. The most common initial presentation was optic nerve involvement (36.1%) and cerebellar sign and symptoms (14.6%). In all, 13 patients (5.8%) had experienced seizure in the course of multiple sclerosis. This study indicated that early-onset multiple sclerosis is not rare condition and overwhelmingly affects girls even at prepubertal onset. Physicians should consider multiple sclerosis in suspicious pediatric cases. © The Author(s) 2016.

An accurate and efficient laser-envelope solver for the modeling of laser-plasma accelerators

DOE PAGES

Benedetti, C.; Schroeder, C. B.; Geddes, C. G. R.; ...

2017-10-17

Detailed and reliable numerical modeling of laser-plasma accelerators (LPAs), where a short and intense laser pulse interacts with an underdense plasma over distances of up to a meter, is a formidably challenging task. This is due to the great disparity among the length scales involved in the modeling, ranging from the micron scale of the laser wavelength to the meter scale of the total laser-plasma interaction length. The use of the time-averaged ponderomotive force approximation, where the laser pulse is described by means of its envelope, enables efficient modeling of LPAs by removing the need to model the details ofmore » electron motion at the laser wavelength scale. Furthermore, it allows simulations in cylindrical geometry which captures relevant 3D physics at 2D computational cost. A key element of any code based on the time-averaged ponderomotive force approximation is the laser envelope solver. In this paper we present the accurate and efficient envelope solver used in the code INF & RNO (INtegrated Fluid & paRticle simulatioN cOde). The features of the INF & RNO laser solver enable an accurate description of the laser pulse evolution deep into depletion even at a reasonably low resolution, resulting in significant computational speed-ups.« less

An accurate and efficient laser-envelope solver for the modeling of laser-plasma accelerators

DOE Office of Scientific and Technical Information (OSTI.GOV)

Benedetti, C.; Schroeder, C. B.; Geddes, C. G. R.

Detailed and reliable numerical modeling of laser-plasma accelerators (LPAs), where a short and intense laser pulse interacts with an underdense plasma over distances of up to a meter, is a formidably challenging task. This is due to the great disparity among the length scales involved in the modeling, ranging from the micron scale of the laser wavelength to the meter scale of the total laser-plasma interaction length. The use of the time-averaged ponderomotive force approximation, where the laser pulse is described by means of its envelope, enables efficient modeling of LPAs by removing the need to model the details ofmore » electron motion at the laser wavelength scale. Furthermore, it allows simulations in cylindrical geometry which captures relevant 3D physics at 2D computational cost. A key element of any code based on the time-averaged ponderomotive force approximation is the laser envelope solver. In this paper we present the accurate and efficient envelope solver used in the code INF & RNO (INtegrated Fluid & paRticle simulatioN cOde). The features of the INF & RNO laser solver enable an accurate description of the laser pulse evolution deep into depletion even at a reasonably low resolution, resulting in significant computational speed-ups.« less

An accurate and efficient laser-envelope solver for the modeling of laser-plasma accelerators

NASA Astrophysics Data System (ADS)

Benedetti, C.; Schroeder, C. B.; Geddes, C. G. R.; Esarey, E.; Leemans, W. P.

2018-01-01

Detailed and reliable numerical modeling of laser-plasma accelerators (LPAs), where a short and intense laser pulse interacts with an underdense plasma over distances of up to a meter, is a formidably challenging task. This is due to the great disparity among the length scales involved in the modeling, ranging from the micron scale of the laser wavelength to the meter scale of the total laser-plasma interaction length. The use of the time-averaged ponderomotive force approximation, where the laser pulse is described by means of its envelope, enables efficient modeling of LPAs by removing the need to model the details of electron motion at the laser wavelength scale. Furthermore, it allows simulations in cylindrical geometry which captures relevant 3D physics at 2D computational cost. A key element of any code based on the time-averaged ponderomotive force approximation is the laser envelope solver. In this paper we present the accurate and efficient envelope solver used in the code INF&RNO (INtegrated Fluid & paRticle simulatioN cOde). The features of the INF&RNO laser solver enable an accurate description of the laser pulse evolution deep into depletion even at a reasonably low resolution, resulting in significant computational speed-ups.

Dynamic simulation of concentrated macromolecular solutions with screened long-range hydrodynamic interactions: Algorithm and limitations

PubMed Central

Ando, Tadashi; Chow, Edmond; Skolnick, Jeffrey

2013-01-01

Hydrodynamic interactions exert a critical effect on the dynamics of macromolecules. As the concentration of macromolecules increases, by analogy to the behavior of semidilute polymer solutions or the flow in porous media, one might expect hydrodynamic screening to occur. Hydrodynamic screening would have implications both for the understanding of macromolecular dynamics as well as practical implications for the simulation of concentrated macromolecular solutions, e.g., in cells. Stokesian dynamics (SD) is one of the most accurate methods for simulating the motions of N particles suspended in a viscous fluid at low Reynolds number, in that it considers both far-field and near-field hydrodynamic interactions. This algorithm traditionally involves an O(N3) operation to compute Brownian forces at each time step, although asymptotically faster but more complex SD methods are now available. Motivated by the idea of hydrodynamic screening, the far-field part of the hydrodynamic matrix in SD may be approximated by a diagonal matrix, which is equivalent to assuming that long range hydrodynamic interactions are completely screened. This approximation allows sparse matrix methods to be used, which can reduce the apparent computational scaling to O(N). Previously there were several simulation studies using this approximation for monodisperse suspensions. Here, we employ newly designed preconditioned iterative methods for both the computation of Brownian forces and the solution of linear systems, and consider the validity of this approximation in polydisperse suspensions. We evaluate the accuracy of the diagonal approximation method using an intracellular-like suspension. The diffusivities of particles obtained with this approximation are close to those with the original method. However, this approximation underestimates intermolecular correlated motions, which is a trade-off between accuracy and computing efficiency. The new method makes it possible to perform large-scale and long-time simulation with an approximate accounting of hydrodynamic interactions. PMID:24089734