The inverse problem of the calculus of variations for discrete systems

NASA Astrophysics Data System (ADS)

Barbero-Liñán, María; Farré Puiggalí, Marta; Ferraro, Sebastián; Martín de Diego, David

2018-05-01

We develop a geometric version of the inverse problem of the calculus of variations for discrete mechanics and constrained discrete mechanics. The geometric approach consists of using suitable Lagrangian and isotropic submanifolds. We also provide a transition between the discrete and the continuous problems and propose variationality as an interesting geometric property to take into account in the design and computer simulation of numerical integrators for constrained systems. For instance, nonholonomic mechanics is generally non variational but some special cases admit an alternative variational description. We apply some standard nonholonomic integrators to such an example to study which ones conserve this property.

NO TIMING VARIATIONS OBSERVED IN THIRD TRANSIT OF SNOW-LINE EXOPLANET KEPLER-421b

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

Dalba, Paul A.; Muirhead, Philip S., E-mail: pdalba@bu.edu

2016-07-20

We observed Kepler-421 during the anticipated third transit of the snow-line exoplanet Kepler-421b in order to constrain the existence and extent of transit timing variations (TTVs). Previously, the Kepler spacecraft only observed two transits of Kepler-421b, leaving the planet’s transit ephemeris unconstrained. Our visible light, time-series observations from the 4.3 m Discovery Channel Telescope were designed to capture pre-transit baseline and the partial transit of Kepler-421b, barring significant TTVs. We use the light curves to assess the probabilities of various transit models using both the posterior odds ratio and the Bayesian Information Criterion, and find that a transit model withmore » no TTVs is favored to 3.6 σ confidence. These observations suggest that Kepler-421b is either alone in its system or is only experiencing minor dynamic interactions with an unseen companion. With the Kepler-421b ephemeris constrained, we calculate future transit times and discuss the opportunity to characterize the atmosphere of this cold, long-period exoplanet via transmission spectroscopy. Our investigation emphasizes the difficulties associated with observing long-period exoplanet transits and the consequences that arise from failing to refine transit ephemerides.« less

Synthetic Constraint of Ecosystem C Models Using Radiocarbon and Net Primary Production (NPP) in New Zealand Grazing Land

NASA Astrophysics Data System (ADS)

Baisden, W. T.

2011-12-01

Time-series radiocarbon measurements have substantial ability to constrain the size and residence time of the soil C pools commonly represented in ecosystem models. Radiocarbon remains unique in the ability to constrain the large stabilized C pool with decadal residence times. Radiocarbon also contributes usefully to constraining the size and turnover rate of the passive pool, but typically struggles to constrain pools with residence times less than a few years. Overall, the number of pools and associated turnover rates that can be constrained depends upon the number of time-series samples available, the appropriateness of chemical or physical fractions to isolate unequivocal pools, and the utility of additional C flux data to provide additional constraints. In New Zealand pasture soils, we demonstrate the ability to constrain decadal turnover times with in a few years for the stabilized pool and reasonably constrain the passive fraction. Good constraint is obtained with two time-series samples spaced 10 or more years apart after 1970. Three or more time-series samples further improve the level of constraint. Work within this context shows that a two-pool model does explain soil radiocarbon data for the most detailed profiles available (11 time-series samples), and identifies clear and consistent differences in rates of C turnover and passive fraction in Andisols vs Non-Andisols. Furthermore, samples from multiple horizons can commonly be combined, yielding consistent residence times and passive fraction estimates that are stable with, or increase with, depth in different sites. Radiocarbon generally fails to quantify rapid C turnover, however. Given that the strength of radiocarbon is estimating the size and turnover of the stabilized (decadal) and passive (millennial) pools, the magnitude of fast cycling pool(s) can be estimated by subtracting the radiocarbon-based estimates of turnover within stabilized and passive pools from total estimates of NPP. In grazing land, these estimates can be derived primarily from measured aboveground NPP and calculated belowground NPP. Results suggest that only 19-36% of heterotrophic soil respiration is derived from the soil C with rapid turnover times. A final logical step in synthesis is the analysis of temporal variation in NPP, primarily due to climate, as driver of changes in plant inputs and resulting in dynamic changes in rapid and decadal soil C pools. In sites with good time series samples from 1959-1975, we examine the apparent impacts of measured or modelled (Biome-BGC) NPP on soil Δ14C. Ultimately, these approaches have the ability to empirically constrain, and provide limited verification, of the soil C cycle as commonly depicted ecosystem biogeochemistry models.

Detecting aseismic strain transients from seismicity data

USGS Publications Warehouse

Llenos, A.L.; McGuire, J.J.

2011-01-01

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate- and state-dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space-time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M???1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor. Copyright 2011 by the American Geophysical Union.

X-ray wind tomography of IGR J17252-3616

NASA Astrophysics Data System (ADS)

Manousakis, Antonios; Walter, Roland

IGR J17252-3616 is an heavily absorbed and eclipsing High Mass X-ray Binary with an ab-sorbing hydrogen column density >1023 cm-2 . We have observed it with XMM-Newton to understand the geometry of the absorbing material. Observations were scheduled in order to cover as many orbital phases as possible. Timing analysis is constraining the orbital solution and the physical parameters of the system. Spectral analysis reveals remarkable variations of the absorbing column density and of the Iron Kα fluorescence line around the eclipse. These variations allow to map the geometry of the absorbing and reflection material. Very large accretion structures could be imaged for the first time.

Constraining a Possible Variation of G with Type Ia Supernovae

NASA Astrophysics Data System (ADS)

Mould, Jeremy; Uddin, Syed A.

2014-03-01

Astrophysical cosmology constrains the variation of Newton's Constant in a manner complementary to laboratory experiments, such as the celebrated lunar laser ranging campaign. Supernova cosmology is an example of the former and has attained campaign status, following planning by a Dark Energy Task Force in 2005. In this paper, we employ the full SNIa data set to the end of 2013 to set a limit on G variation. In our approach, we adopt the standard candle delineation of the redshift distance relation. We set an upper limit on its rate of change |dot{G}/G| of 0.1 parts per billion per year over 9 Gyrs. By contrast, lunar laser ranging tests variation of G over the last few decades. Conversely, one may adopt the laboratory result as a prior and constrain the effect of variable G in dark energy equation of state experiments to δw < 0.02. We also examine the parameterisation G 1 + z. Its short expansion age conflicts with the measured values of the expansion rate and the density in a flat Universe. In conclusion, supernova cosmology complements other experiments in limiting G variation. An important caveat is that it rests on the assumption that the same mass of 56Ni is burned to create the standard candle regardless of redshift. These two quantities, f and G, where f is the Chandrasekhar mass fraction burned, are degenerate. Constraining f variation alone requires more understanding of the SNIa mechanism.

Active flutter suppression using optical output feedback digital controllers

NASA Technical Reports Server (NTRS)

1982-01-01

A method for synthesizing digital active flutter suppression controllers using the concept of optimal output feedback is presented. A convergent algorithm is employed to determine constrained control law parameters that minimize an infinite time discrete quadratic performance index. Low order compensator dynamics are included in the control law and the compensator parameters are computed along with the output feedback gain as part of the optimization process. An input noise adjustment procedure is used to improve the stability margins of the digital active flutter controller. Sample rate variation, prefilter pole variation, control structure variation and gain scheduling are discussed. A digital control law which accommodates computation delay can stabilize the wing with reasonable rms performance and adequate stability margins.

Pristine extraterrestrial material with unprecedented nitrogen isotopic variation.

PubMed

Briani, Giacomo; Gounelle, Matthieu; Marrocchi, Yves; Mostefaoui, Smail; Leroux, Hugues; Quirico, Eric; Meibom, Anders

2009-06-30

Pristine meteoritic materials carry light element isotopic fractionations that constrain physiochemical conditions during solar system formation. Here we report the discovery of a unique xenolith in the metal-rich chondrite Isheyevo. Its fine-grained, highly pristine mineralogy has similarity with interplanetary dust particles (IDPs), but the volume of the xenolith is more than 30,000 times that of a typical IDP. Furthermore, an extreme continuum of N isotopic variation is present in this xenolith: from very light N isotopic composition (delta(15)N(AIR) = -310 +/- 20 per thousand), similar to that inferred for the solar nebula, to the heaviest ratios measured in any solar system material (delta(15)N(AIR) = 4,900 +/- 300 per thousand). At the same time, its hydrogen and carbon isotopic compositions exhibit very little variation. This object poses serious challenges for existing models for the origin of light element isotopic anomalies.

A sequential solution for anisotropic total variation image denoising with interval constraints

NASA Astrophysics Data System (ADS)

Xu, Jingyan; Noo, Frédéric

2017-09-01

We show that two problems involving the anisotropic total variation (TV) and interval constraints on the unknown variables admit, under some conditions, a simple sequential solution. Problem 1 is a constrained TV penalized image denoising problem; problem 2 is a constrained fused lasso signal approximator. The sequential solution entails finding first the solution to the unconstrained problem, and then applying a thresholding to satisfy the constraints. If the interval constraints are uniform, this sequential solution solves problem 1. If the interval constraints furthermore contain zero, the sequential solution solves problem 2. Here uniform interval constraints refer to all unknowns being constrained to the same interval. A typical example of application is image denoising in x-ray CT, where the image intensities are non-negative as they physically represent linear attenuation coefficient in the patient body. Our results are simple yet seem unknown; we establish them using the Karush-Kuhn-Tucker conditions for constrained convex optimization.

Determination of wave-function functionals: The constrained-search variational method

NASA Astrophysics Data System (ADS)

Pan, Xiao-Yin; Sahni, Viraht; Massa, Lou

2005-09-01

In a recent paper [Phys. Rev. Lett. 93, 130401 (2004)], we proposed the idea of expanding the space of variations in variational calculations of the energy by considering the approximate wave function ψ to be a functional of functions χ , ψ=ψ[χ] , rather than a function. A constrained search is first performed over all functions χ such that the wave-function functional ψ[χ] satisfies a physical constraint or leads to the known value of an observable. A rigorous upper bound to the energy is then obtained via the variational principle. In this paper we generalize the constrained-search variational method, applicable to both ground and excited states, to the determination of arbitrary Hermitian single-particle operators as applied to two-electron atomic and ionic systems. We construct analytical three-parameter ground-state functionals for the H- ion and the He atom through the constraint of normalization. We present the results for the total energy E , the expectations of the single-particle operators W=∑irin , n=-2,-1,1,2 , W=∑iδ(ri) , and W=∑iδ(ri-r) , the structure of the nonlocal Coulomb hole charge ρc(rr') , and the expectations of the two particle operators u2,u,1/u,1/u2 , where u=∣ri-rj∣ . The results for all the expectation values are remarkably accurate when compared with the 1078-parameter wave function of Pekeris, and other wave functions that are not functionals. We conclude by describing our current work on how the constrained-search variational method in conjunction with quantal density-functional theory is being applied to the many-electron case.

A brief survey of constrained mechanics and variational problems in terms of differential forms

NASA Technical Reports Server (NTRS)

Hermann, Robert

1994-01-01

There has been considerable interest recently in constrained mechanics and variational problems. This is in part due to applied interests (such as 'non-holonomic mechanics in robotics') and in other part due to the fact that several schools of 'pure' mathematics have found that this classical subject is of importance for what they are trying to do. I have made various attempts at developing these subjects since my Lincoln lab days of the late 1950's. In this Chapter, I will sketch a Unified point of view, using Cartan's approach with differential forms. This has the advantage from the C-O-R viewpoint being developed in this Volume that the extension from 'smooth' to 'generalized' data is very systematic and algebraic. (I will only deal with the 'smooth' point of view in this Chapter; I will develop the 'generalized function' material at a later point.) The material presented briefly here about Variational Calculus and Constrained Mechanics can be found in more detail in my books, 'Differential Geometry and the Calculus of Variations', 'Lie Algebras and Quantum Mechanics', and 'Geometry, Physics and Systems'.

Acoustic-radiation stress in solids. I - Theory

NASA Technical Reports Server (NTRS)

Cantrell, J. H., Jr.

1984-01-01

The general case of acoustic-radiation stress associated with quasi-compressional and quasi-shear waves propagating in infinite and semiinfinite lossless solids of arbitrary crystalline symmetry is studied. The Boussinesq radiation stress is defined and found to depend directly on an acoustic nonlinearity parameter which characterizes the radiation-induced static strain, a stress-generalized nonlinearity parameter which characterizes the stress nonlinearity, and the energy density of the propagating wave. Application of the Boltzmann-Ehrenfest principle of adiabatic invariance to a self-constrained system described by the nonlinear equations of motion allows the acoustic-radiation-induced static strain to be identified with a self-constrained variation in the time-averaged product of the internal energy density and displacement gradient. The time-averaged product is scaled by the acoustic nonlinearity parameter and represents the first-order nonlinearity in the virial theorem. Finally, the relationship between the Boussinesq and the Cauchy radiation stress is obtained in a closed three-dimensional form.

Missile Guidance Law Based on Robust Model Predictive Control Using Neural-Network Optimization.

PubMed

Li, Zhijun; Xia, Yuanqing; Su, Chun-Yi; Deng, Jun; Fu, Jun; He, Wei

2015-08-01

In this brief, the utilization of robust model-based predictive control is investigated for the problem of missile interception. Treating the target acceleration as a bounded disturbance, novel guidance law using model predictive control is developed by incorporating missile inside constraints. The combined model predictive approach could be transformed as a constrained quadratic programming (QP) problem, which may be solved using a linear variational inequality-based primal-dual neural network over a finite receding horizon. Online solutions to multiple parametric QP problems are used so that constrained optimal control decisions can be made in real time. Simulation studies are conducted to illustrate the effectiveness and performance of the proposed guidance control law for missile interception.

Postglacial migration supplements climate in determining plant species ranges in Europe

PubMed Central

Normand, Signe; Ricklefs, Robert E.; Skov, Flemming; Bladt, Jesper; Tackenberg, Oliver; Svenning, Jens-Christian

2011-01-01

The influence of dispersal limitation on species ranges remains controversial. Considering the dramatic impacts of the last glaciation in Europe, species might not have tracked climate changes through time and, as a consequence, their present-day ranges might be in disequilibrium with current climate. For 1016 European plant species, we assessed the relative importance of current climate and limited postglacial migration in determining species ranges using regression modelling and explanatory variables representing climate, and a novel species-specific hind-casting-based measure of accessibility to postglacial colonization. Climate was important for all species, while postglacial colonization also constrained the ranges of more than 50 per cent of the species. On average, climate explained five times more variation in species ranges than accessibility, but accessibility was the strongest determinant for one-sixth of the species. Accessibility was particularly important for species with limited long-distance dispersal ability, with southern glacial ranges, seed plants compared with ferns, and small-range species in southern Europe. In addition, accessibility explained one-third of the variation in species' disequilibrium with climate as measured by the realized/potential range size ratio computed with niche modelling. In conclusion, we show that although climate is the dominant broad-scale determinant of European plant species ranges, constrained dispersal plays an important supplementary role. PMID:21543356

Lidar Penetration Depth Observations for Constraining Cloud Longwave Feedbacks

NASA Astrophysics Data System (ADS)

Vaillant de Guelis, T.; Chepfer, H.; Noel, V.; Guzman, R.; Winker, D. M.; Kay, J. E.; Bonazzola, M.

2017-12-01

Satellite-borne active remote sensing Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations [CALIPSO; Winker et al., 2010] and CloudSat [Stephens et al., 2002] provide direct measurements of the cloud vertical distribution, with a very high vertical resolution. The penetration depth of the laser of the lidar Z_Opaque is directly linked to the LongWave (LW) Cloud Radiative Effect (CRE) at Top Of Atmosphere (TOA) [Vaillant de Guélis et al., in review]. In addition, this measurement is extremely stable in time making it an excellent observational candidate to verify and constrain the cloud LW feedback mechanism [Chepfer et al., 2014]. In this work, we present a method to decompose the variations of the LW CRE at TOA using cloud properties observed by lidar [GOCCP v3.0; Guzman et al., 2017]. We decompose these variations into contributions due to changes in five cloud properties: opaque cloud cover, opaque cloud altitude, thin cloud cover, thin cloud altitude, and thin cloud emissivity [Vaillant de Guélis et al., in review]. We apply this method, in the real world, to the CRE variations of CALIPSO 2008-2015 record, and, in climate model, to LMDZ6 and CESM simulations of the CRE variations of 2008-2015 period and of the CRE difference between a warm climate and the current climate. In climate model simulations, the same cloud properties as those observed by CALIOP are extracted from the CFMIP Observation Simulator Package (COSP) [Bodas-Salcedo et al., 2011] lidar simulator [Chepfer et al., 2008], which mimics the observations that would be performed by the lidar on board CALIPSO satellite. This method, when applied on multi-model simulations of current and future climate, could reveal the altitude of cloud opacity level observed by lidar as a strong constrain for cloud LW feedback, since the altitude feedback mechanism is physically explainable and the altitude of cloud opacity accurately observed by lidar.

Regularizing Unpredictable Variation: Evidence from a Natural Language Setting

ERIC Educational Resources Information Center

Hendricks, Alison Eisel; Miller, Karen; Jackson, Carrie N.

2018-01-01

While previous sociolinguistic research has demonstrated that children faithfully acquire probabilistic input constrained by sociolinguistic and linguistic factors (e.g., gender and socioeconomic status), research suggests children regularize inconsistent input-probabilistic input that is not sociolinguistically constrained (e.g., Hudson Kam &…

Joint seismic-infrasonic processing of recordings from a repeating source of atmospheric explosions.

PubMed

Gibbons, Steven J; Ringdal, Frode; Kvaerna, Tormod

2007-11-01

A database has been established of seismic and infrasonic recordings from more than 100 well-constrained surface explosions, conducted by the Finnish military to destroy old ammunition. The recorded seismic signals are essentially identical and indicate that the variation in source location and magnitude is negligible. In contrast, the infrasonic arrivals on both seismic and infrasound sensors exhibit significant variation both with regard to the number of detected phases, phase travel times, and phase amplitudes, which would be attributable to atmospheric factors. This data set provides an excellent database for studies in sound propagation, infrasound array detection, and direction estimation.

Tolerance to deer herbivory and resistance to insect herbivores in the common evening primrose (Oenothera biennis).

PubMed

Puentes, A; Johnson, M T J

2016-01-01

The evolution of plant defence in response to herbivory will depend on the fitness effects of damage, availability of genetic variation and potential ecological and genetic constraints on defence. Here, we examine the potential for evolution of tolerance to deer herbivory in Oenothera biennis while simultaneously considering resistance to natural insect herbivores. We examined (i) the effects of deer damage on fitness, (ii) the presence of genetic variation in tolerance and resistance, (iii) selection on tolerance, (iv) genetic correlations with resistance that could constrain evolution of tolerance and (v) plant traits that might predict defence. In a field experiment, we simulated deer damage occurring early and late in the season, recorded arthropod abundances, flowering phenology and measured growth rate and lifetime reproduction. Our study showed that deer herbivory has a negative effect on fitness, with effects being more pronounced for late-season damage. Selection acted to increase tolerance to deer damage, yet there was low and nonsignificant genetic variation in this trait. In contrast, there was substantial genetic variation in resistance to insect herbivores. Resistance was genetically uncorrelated with tolerance, whereas positive genetic correlations in resistance to insect herbivores suggest there exists diffuse selection on resistance traits. In addition, growth rate and flowering time did not predict variation in tolerance, but flowering phenology was genetically correlated with resistance. Our results suggest that deer damage has the potential to exert selection because browsing reduces plant fitness, but limited standing genetic variation in tolerance is expected to constrain adaptive evolution in O. biennis. © 2015 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2015 European Society For Evolutionary Biology.

Geometric constrained variational calculus. III: The second variation (Part II)

NASA Astrophysics Data System (ADS)

Massa, Enrico; Luria, Gianvittorio; Pagani, Enrico

2016-03-01

The problem of minimality for constrained variational calculus is analyzed within the class of piecewise differentiable extremaloids. A fully covariant representation of the second variation of the action functional based on a family of local gauge transformations of the original Lagrangian is proposed. The necessity of pursuing a local adaptation process, rather than the global one described in [1] is seen to depend on the value of certain scalar attributes of the extremaloid, here called the corners’ strengths. On this basis, both the necessary and the sufficient conditions for minimality are worked out. In the discussion, a crucial role is played by an analysis of the prolongability of the Jacobi fields across the corners. Eventually, in the appendix, an alternative approach to the concept of strength of a corner, more closely related to Pontryagin’s maximum principle, is presented.

Performance Analysis of BDS Medium-Long Baseline RTK Positioning Using an Empirical Troposphere Model.

PubMed

Shu, Bao; Liu, Hui; Xu, Longwei; Qian, Chuang; Gong, Xiaopeng; An, Xiangdong

2018-04-14

For GPS medium-long baseline real-time kinematic (RTK) positioning, the troposphere parameter is introduced along with coordinates, and the model is ill-conditioned due to its strong correlation with the height parameter. For BeiDou Navigation Satellite System (BDS), additional difficulties occur due to its special satellite constellation. In fact, relative zenith troposphere delay (RZTD) derived from high-precision empirical zenith troposphere models can be introduced. Thus, the model strength can be improved, which is also called the RZTD-constrained RTK model. In this contribution, we first analyze the factors affecting the precision of BDS medium-long baseline RTK; thereafter, 15 baselines ranging from 38 km to 167 km in different troposphere conditions are processed to assess the performance of RZTD-constrained RTK. Results show that the troposphere parameter is difficult to distinguish from the height component, even with long time filtering for BDS-only RTK. Due to the lack of variation in geometry for the BDS geostationary Earth orbit satellite, the long convergence time of ambiguity parameters may reduce the height precision of GPS/BDS-combined RTK in the initial period. When the RZTD-constrained model was used in BDS and GPS/BDS-combined situations compared with the traditional RTK, the standard deviation of the height component for the fixed solution was reduced by 52.4% and 34.0%, respectively.

Performance Analysis of BDS Medium-Long Baseline RTK Positioning Using an Empirical Troposphere Model

PubMed Central

Liu, Hui; Xu, Longwei; Qian, Chuang; Gong, Xiaopeng; An, Xiangdong

2018-01-01

For GPS medium-long baseline real-time kinematic (RTK) positioning, the troposphere parameter is introduced along with coordinates, and the model is ill-conditioned due to its strong correlation with the height parameter. For BeiDou Navigation Satellite System (BDS), additional difficulties occur due to its special satellite constellation. In fact, relative zenith troposphere delay (RZTD) derived from high-precision empirical zenith troposphere models can be introduced. Thus, the model strength can be improved, which is also called the RZTD-constrained RTK model. In this contribution, we first analyze the factors affecting the precision of BDS medium-long baseline RTK; thereafter, 15 baselines ranging from 38 km to 167 km in different troposphere conditions are processed to assess the performance of RZTD-constrained RTK. Results show that the troposphere parameter is difficult to distinguish from the height component, even with long time filtering for BDS-only RTK. Due to the lack of variation in geometry for the BDS geostationary Earth orbit satellite, the long convergence time of ambiguity parameters may reduce the height precision of GPS/BDS-combined RTK in the initial period. When the RZTD-constrained model was used in BDS and GPS/BDS-combined situations compared with the traditional RTK, the standard deviation of the height component for the fixed solution was reduced by 52.4% and 34.0%, respectively. PMID:29661999

He nui na ala e hiki aku ai: Factors Influencing Phonetic Variation in the Hawaiian Word "keia"

ERIC Educational Resources Information Center

Drager, Katie; Comstock, Bethany Kaleialohapau'ole Chun; Kneubuhl, Hina Puamohala

2017-01-01

Apart from a handful of studies (e.g., Kinney 1956), linguists know little about what variation exists in Hawaiian and what factors constrain the variation. In this paper, we present an analysis of phonetic variation in the word "keia," meaning "this," examining the social, linguistic, and probabilistic factors that constrain…

Constrained variational calculus for higher order classical field theories

NASA Astrophysics Data System (ADS)

Campos, Cédric M.; de León, Manuel; Martín de Diego, David

2010-11-01

We develop an intrinsic geometrical setting for higher order constrained field theories. As a main tool we use an appropriate generalization of the classical Skinner-Rusk formalism. Some examples of applications are studied, in particular to the geometrical description of optimal control theory for partial differential equations.

Buffer capacity, ecosystem feedbacks, and seawater chemistry under global change

NASA Astrophysics Data System (ADS)

Jury, C. P.; Thomas, F. I.; Atkinson, M. J.; Jokiel, P. L.; Onuma, M. A.; Kaku, N.; Toonen, R. J.

2013-12-01

Ocean acidification (OA) results in reduced seawater pH and aragonite saturation state (Ωarag), but also reduced seawater buffer capacity. As buffer capacity decreases, diel variation in seawater chemistry increases. However, a variety of ecosystem feedbacks can modulate changes in both average seawater chemistry and diel seawater chemistry variation. Here we model these effects for a coastal, reef flat ecosystem. We show that an increase in offshore pCO2 and temperature (to 900 μatm and +3°C) can increase diel pH variation by as much as a factor of 2.5 and can increase diel pCO2 variation by a factor of 4.6, depending on ecosystem feedbacks and seawater residence time. Importantly, these effects are different between day and night. With increasing seawater residence time and increasing feedback intensity, daytime seawater chemistry becomes more similar to present-day conditions while nighttime seawater chemistry becomes less similar to present-day conditions. Better constraining ecosystem feedbacks under global change will improve projections of coastal water chemistry, but this study shows the importance of considering changes in both average carbonate chemistry and diel chemistry variation for organisms and ecosystems. Further, we will discuss our recent work examining the effects of diel seawater chemistry variation on coral calcification rates.

The Holometer: An instrument to probe Planckian quantum geometry

DOE PAGES

Chou, Aaron; Glass, Henry; Gustafson, H. Richard; ...

2017-02-28

This paper describes the Fermilab Holometer, an instrument for measuring correlations of position variations over a four-dimensional volume of space-time. The apparatus consists of two co-located, but independent and isolated, 40 m power-recycled Michelson interferometers, whose outputs are cross-correlated to 25 MHz. The data are sensitive to correlations of differential position across the apparatus over a broad band of frequencies up to and exceeding the inverse light crossing time, 7.6 MHz. A noise model constrained by diagnostic and environmental data distinguishes among physical origins of measured correlations, and is used to verify shot-noise-limited performance. These features allow searches for exoticmore » quantum correlations that depart from classical trajectories at spacelike separations, with a strain noise power spectral density sensitivity smaller than the Planck time. As a result, the Holometer in current and future configurations is projected to provide precision tests of a wide class of models of quantum geometry at the Planck scale, beyond those already constrained by currently operating gravitational wave observatories.« less

An analysis of the transit times of TrES-1b

NASA Astrophysics Data System (ADS)

Steffen, Jason H.; Agol, Eric

2005-11-01

The presence of a second planet in a known, transiting-planet system will cause the time between transits to vary. These variations can be used to constrain the orbital elements and mass of the perturbing planet. We analyse the set of transit times of the TrES-1 system given in Charbonneau et al. We find no convincing evidence for a second planet in the TrES-1 system from those data. By further analysis, we constrain the mass that a perturbing planet could have as a function of the semi-major axis ratio of the two planets and the eccentricity of the perturbing planet. Near low-order, mean-motion resonances (within ~1 per cent fractional deviation), we find that a secondary planet must generally have a mass comparable to or less than the mass of the Earth - showing that these data are the first to have sensitivity to sub-Earth-mass planets. We compare the sensitivity of this technique to the mass of the perturbing planet with future, high-precision radial velocity measurements.

Simulation of X-ray absorption spectra with orthogonality constrained density functional theory.

PubMed

Derricotte, Wallace D; Evangelista, Francesco A

2015-06-14

Orthogonality constrained density functional theory (OCDFT) [F. A. Evangelista, P. Shushkov and J. C. Tully, J. Phys. Chem. A, 2013, 117, 7378] is a variational time-independent approach for the computation of electronic excited states. In this work we extend OCDFT to compute core-excited states and generalize the original formalism to determine multiple excited states. Benchmark computations on a set of 13 small molecules and 40 excited states show that unshifted OCDFT/B3LYP excitation energies have a mean absolute error of 1.0 eV. Contrary to time-dependent DFT, OCDFT excitation energies for first- and second-row elements are computed with near-uniform accuracy. OCDFT core excitation energies are insensitive to the choice of the functional and the amount of Hartree-Fock exchange. We show that OCDFT is a powerful tool for the assignment of X-ray absorption spectra of large molecules by simulating the gas-phase near-edge spectrum of adenine and thymine.

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

Diaz Cruz, J. Lorenzo

The standard Higgs mechanism employed in the Standard Model (SM) for electroweak symmetry breaking, relies on a homogenous Higgs vacuum expectation value (v.e.v.), i.e. a vacuum that does not depend on the position or the time coordinates. However, other non-homogeneous structures could also be considered, either at long or short distances. For instance, spatial variations of the Higgs v.e.v. on cosmological scales, would induce variations of the fundamental constants, and are severely constrained. Other possibilities, such as a discrete microscopic structure of the Higgs vacuum, or a confined Higgs mechanism associated with a strongly interacting Higgs sector, could be testedmore » and give some light on the electroweak-scale contributions to the cosmological constant.« less

Improving Grasp Skills Using Schema Structured Learning

NASA Technical Reports Server (NTRS)

Platt, Robert; Grupen, ROderic A.; Fagg, Andrew H.

2006-01-01

Abstract In the control-based approach to robotics, complex behavior is created by sequencing and combining control primitives. While it is desirable for the robot to autonomously learn the correct control sequence, searching through the large number of potential solutions can be time consuming. This paper constrains this search to variations of a generalized solution encoded in a framework known as an action schema. A new algorithm, SCHEMA STRUCTURED LEARNING, is proposed that repeatedly executes variations of the generalized solution in search of instantiations that satisfy action schema objectives. This approach is tested in a grasping task where Dexter, the UMass humanoid robot, learns which reaching and grasping controllers maximize the probability of grasp success.

Limits on the Time Evolution of Space Dimensions from Newton's Constant

NASA Astrophysics Data System (ADS)

Nasseri, Forough

Limits are imposed upon the possible rate of change of extra spatial dimensions in a decrumpling model Universe with time variable spatial dimensions (TVSD) by considering the time variation of (1+3)-dimensional Newton's constant. Previous studies on the time variation of (1+3)-dimensional Newton's constant in TVSD theory had not include the effects of the volume of the extra dimensions and the effects of the surface area of the unit sphere in D-space dimensions. Our main result is that the absolute value of the present rate of change of spatial dimensions to be less than about 10-14 yr-1. Our results would appear to provide a prima facie case for ruling the TVSD model out. We show that based on observational bounds on the present variation of Newton's constant, one would have to conclude that the spatial dimension of the Universe when the Universe was "at the Planck scale" to be less than or equal to 3.09. If the dimension of space when the Universe was "at the Planck scale" is constrained to be fractional and very close to 3, then the whole edifice of TVSD model loses credibility.

On the calculation of charge transfer transitions with standard density functionals using constrained variational density functional theory.

PubMed

Ziegler, Tom; Krykunov, Mykhaylo

2010-08-21

It is well known that time-dependent density functional theory (TD-DFT) based on standard gradient corrected functionals affords both a quantitative and qualitative incorrect picture of charge transfer transitions between two spatially separated regions. It is shown here that the well known failure can be traced back to the use of linear response theory. Further, it is demonstrated that the inclusion of higher order terms readily affords a qualitatively correct picture even for simple functionals based on the local density approximation. The inclusion of these terms is done within the framework of a newly developed variational approach to excitation energies called constrained variational density functional theory (CV-DFT). To second order [CV(2)-DFT] this theory is identical to adiabatic TD-DFT within the Tamm-Dancoff approximation. With inclusion of fourth order corrections [CV(4)-DFT] it affords a qualitative correct description of charge transfer transitions. It is finally demonstrated that the relaxation of the ground state Kohn-Sham orbitals to first order in response to the change in density on excitation together with CV(4)-DFT affords charge transfer excitations in good agreement with experiment. The new relaxed theory is termed R-CV(4)-DFT. The relaxed scheme represents an effective way in which to introduce double replacements into the description of single electron excitations, something that would otherwise require a frequency dependent kernel.

Transactions of the Eleventh Army Conference on Applied Mathematics and Computing

DTIC Science & Technology

1994-03-01

Hu, S., Laksmikantham, V. and Rama Mohana Rao , M., Nonlinear Variation of Parameters formula for Integro-differential equations of Volterra type, J. M...454-457, 1986 [31 N.K. Sinha and G.P. Rao , Identification of continuous-time systems, Boston: Kluwer Academic, 1991. [41 J. Ezzine and C.D. Johnson...Theory, pages 536-548, 1989. [15] J.A. O’Sullivan, P. Moulin, and D.L. Snyder. Cramer- rao bounds for constrained spectrum estimation with application

High regional climate sensitivity over continental China constrained by glacial-recent changes in temperature and the hydrological cycle.

PubMed

Eagle, Robert A; Risi, Camille; Mitchell, Jonathan L; Eiler, John M; Seibt, Ulrike; Neelin, J David; Li, Gaojun; Tripati, Aradhna K

2013-05-28

The East Asian monsoon is one of Earth's most significant climatic phenomena, and numerous paleoclimate archives have revealed that it exhibits variations on orbital and suborbital time scales. Quantitative constraints on the climate changes associated with these past variations are limited, yet are needed to constrain sensitivity of the region to changes in greenhouse gas levels. Here, we show central China is a region that experienced a much larger temperature change since the Last Glacial Maximum than typically simulated by climate models. We applied clumped isotope thermometry to carbonates from the central Chinese Loess Plateau to reconstruct temperature and water isotope shifts from the Last Glacial Maximum to present. We find a summertime temperature change of 6-7 °C that is reproduced by climate model simulations presented here. Proxy data reveal evidence for a shift to lighter isotopic composition of meteoric waters in glacial times, which is also captured by our model. Analysis of model outputs suggests that glacial cooling over continental China is significantly amplified by the influence of stationary waves, which, in turn, are enhanced by continental ice sheets. These results not only support high regional climate sensitivity in Central China but highlight the fundamental role of planetary-scale atmospheric dynamics in the sensitivity of regional climates to continental glaciation, changing greenhouse gas levels, and insolation.

Multiscale Dynamics of Aseismic Slip on Central San Andreas Fault

NASA Astrophysics Data System (ADS)

Khoshmanesh, M.; Shirzaei, M.

2018-03-01

Understanding the evolution of aseismic slip enables constraining the fault's seismic budget and provides insight into dynamics of creep. Inverting the time series of surface deformation measured along the Central San Andreas Fault obtained from interferometric synthetic aperture radar in combination with measurements of repeating earthquakes, we constrain the spatiotemporal distribution of creep during 1992-2010. We identify a new class of intermediate-term creep rate variations that evolve over decadal scale, releasing stress on the accelerating zone and loading adjacent decelerating patches. We further show that in short-term (<2 year period), creep avalanches, that is, isolated clusters of accelerated aseismic slip with velocities exceeding the long-term rate, govern the dynamics of creep. The statistical properties of these avalanches suggest existence of elevated pore pressure in the fault zone, consistent with laboratory experiments.

Historical precipitation predictably alters the shape and magnitude of microbial functional response to soil moisture.

PubMed

Averill, Colin; Waring, Bonnie G; Hawkes, Christine V

2016-05-01

Soil moisture constrains the activity of decomposer soil microorganisms, and in turn the rate at which soil carbon returns to the atmosphere. While increases in soil moisture are generally associated with increased microbial activity, historical climate may constrain current microbial responses to moisture. However, it is not known if variation in the shape and magnitude of microbial functional responses to soil moisture can be predicted from historical climate at regional scales. To address this problem, we measured soil enzyme activity at 12 sites across a broad climate gradient spanning 442-887 mm mean annual precipitation. Measurements were made eight times over 21 months to maximize sampling during different moisture conditions. We then fit saturating functions of enzyme activity to soil moisture and extracted half saturation and maximum activity parameter values from model fits. We found that 50% of the variation in maximum activity parameters across sites could be predicted by 30-year mean annual precipitation, an indicator of historical climate, and that the effect is independent of variation in temperature, soil texture, or soil carbon concentration. Based on this finding, we suggest that variation in the shape and magnitude of soil microbial response to soil moisture due to historical climate may be remarkably predictable at regional scales, and this approach may extend to other systems. If historical contingencies on microbial activities prove to be persistent in the face of environmental change, this approach also provides a framework for incorporating historical climate effects into biogeochemical models simulating future global change scenarios. © 2016 John Wiley & Sons Ltd.

What GRACE/GRACE-FO satellite gravity may tell about the atmosphere (and what not)

NASA Astrophysics Data System (ADS)

Eicker, Annette; Springer, Anne; Hense, Andreas; Panet, Isabelle; Kusche, Jürgen

2017-04-01

In this presentation we would like to discuss the present benefit and future potential of satellite gravity observations, as obtained from the satellite mission GRACE and its successor GRACE-Follow-On (GRACE-FO), for studying the atmospheric water cycle. In the first part of the presentation, we will show recent results of using GRACE to constrain atmospheric water budgets. GRACE-derived water storage changes (in combination with observed runoff) can be used to solve for the vertical water flux deficit of precipitation (P) minus evapotranspiration (E), which links the terrestrial and the atmospheric water balance equations. This relates gravity change to moisture flux divergence and water vapor change and thus provides, in principle, a link between GRACE/GRACE-FO and (area-averaged) GNSS integrated water vapor observations that may be exploited in the future. We will show that such an independent estimate of P minus E can be used to constrain land-atmosphere fluxes from monthly time scales to decadal trends and even provides meaningful flux information down to daily time steps. In the second part of the presentation, we would like to give an outlook towards the potential of using satellite gravity data directly for the estimation of atmospheric water mass changes. On the basis of ERA-Interim data, we provide a first assessment which suggests that an anticipated future double-pair gravity mission with enhanced temporal and spatial resolution would be sensitive to 'feeling' atmospheric water mass (water vapor) variations. However, whether these (faster) variations could be separated from dry air mass variations through modeling needs to be investigated. If possible, this would offer a completely new tool for validating atmospheric analyses and for improving engergy and mass budgets in models.

Cannibalism and activity rate in larval damselflies increase along a latitudinal gradient as a consequence of time constraints.

PubMed

Sniegula, Szymon; Golab, Maria J; Johansson, Frank

2017-07-14

Predation is ubiquitous in nature. One form of predation is cannibalism, which is affected by many factors such as size structure and resource density. However, cannibalism may also be influenced by abiotic factors such as seasonal time constraints. Since time constraints are greater at high latitudes, cannibalism could be stronger at such latitudes, but we know next to nothing about latitudinal variation in cannibalism. In this study, we examined cannibalism and activity in larvae of the damselfly Lestes sponsa along a latitudinal gradient across Europe. We did this by raising larvae from the egg stage at different temperatures and photoperiods corresponding to different latitudes. We found that the more seasonally time-constrained populations in northern latitudes and individuals subjected to greater seasonal time constraints exhibited a higher level of cannibalism. We also found that activity was higher at north latitude conditions, and thus correlated with cannibalism, suggesting that this behaviour mediates higher levels of cannibalism in time-constrained animals. Our results go counter to the classical latitude-predation pattern which predicts higher predation at lower latitudes, since we found that predation was stronger at higher latitudes. The differences in cannibalism might have implications for population dynamics along the latitudinal gradients, but further experiments are needed to explore this.

WE-FG-207B-05: Iterative Reconstruction Via Prior Image Constrained Total Generalized Variation for Spectral CT

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

Niu, S; Zhang, Y; Ma, J

Purpose: To investigate iterative reconstruction via prior image constrained total generalized variation (PICTGV) for spectral computed tomography (CT) using fewer projections while achieving greater image quality. Methods: The proposed PICTGV method is formulated as an optimization problem, which balances the data fidelity and prior image constrained total generalized variation of reconstructed images in one framework. The PICTGV method is based on structure correlations among images in the energy domain and high-quality images to guide the reconstruction of energy-specific images. In PICTGV method, the high-quality image is reconstructed from all detector-collected X-ray signals and is referred as the broad-spectrum image. Distinctmore » from the existing reconstruction methods applied on the images with first order derivative, the higher order derivative of the images is incorporated into the PICTGV method. An alternating optimization algorithm is used to minimize the PICTGV objective function. We evaluate the performance of PICTGV on noise and artifacts suppressing using phantom studies and compare the method with the conventional filtered back-projection method as well as TGV based method without prior image. Results: On the digital phantom, the proposed method outperforms the existing TGV method in terms of the noise reduction, artifacts suppression, and edge detail preservation. Compared to that obtained by the TGV based method without prior image, the relative root mean square error in the images reconstructed by the proposed method is reduced by over 20%. Conclusion: The authors propose an iterative reconstruction via prior image constrained total generalize variation for spectral CT. Also, we have developed an alternating optimization algorithm and numerically demonstrated the merits of our approach. Results show that the proposed PICTGV method outperforms the TGV method for spectral CT.« less

Constraining slip rates and spacings for active normal faults

NASA Astrophysics Data System (ADS)

Cowie, Patience A.; Roberts, Gerald P.

2001-12-01

Numerous observations of extensional provinces indicate that neighbouring faults commonly slip at different rates and, moreover, may be active over different time intervals. These published observations include variations in slip rate measured along-strike of a fault array or fault zone, as well as significant across-strike differences in the timing and rates of movement on faults that have a similar orientation with respect to the regional stress field. Here we review published examples from the western USA, the North Sea, and central Greece, and present new data from the Italian Apennines that support the idea that such variations are systematic and thus to some extent predictable. The basis for the prediction is that: (1) the way in which a fault grows is fundamentally controlled by the ratio of maximum displacement to length, and (2) the regional strain rate must remain approximately constant through time. We show how data on fault lengths and displacements can be used to model the observed patterns of long-term slip rate where measured values are sparse. Specifically, we estimate the magnitude of spatial variation in slip rate along-strike and relate it to the across-strike spacing between active faults.

Extremes in ecology: Avoiding the misleading effects of sampling variation in summary analyses

USGS Publications Warehouse

Link, W.A.; Sauer, J.R.

1996-01-01

Surveys such as the North American Breeding Bird Survey (BBS) produce large collections of parameter estimates. One's natural inclination when confronted with lists of parameter estimates is to look for the extreme values: in the BBS, these correspond to the species that appear to have the greatest changes in population size through time. Unfortunately, extreme estimates are liable to correspond to the most poorly estimated parameters. Consequently, the most extreme parameters may not match up with the most extreme parameter estimates. The ranking of parameter values on the basis of their estimates are a difficult statistical problem. We use data from the BBS and simulations to illustrate the potential misleading effects of sampling variation in rankings of parameters. We describe empirical Bayes and constrained empirical Bayes procedures which provide partial solutions to the problem of ranking in the presence of sampling variation.

Time-variable magma pressure at Kīlauea Volcano yields constraint on the volume and volatile content of shallow magma storage

NASA Astrophysics Data System (ADS)

Anderson, K. R.; Patrick, M. R.; Poland, M. P.; Miklius, A.

2015-12-01

Episodic depressurization-pressurization cycles of Kīlauea Volcano's shallow magma system cause variations in ground deformation, eruption rate, and surface height of the active summit lava lake. The mechanism responsible for these pressure-change cycles remains enigmatic, but associated monitoring signals often show a quasi-exponential temporal history that is consistent with a temporary reduction (or blockage) of supply to Kīlauea's shallow magma storage area. Regardless of their cause, the diverse signals produced by these deflation-inflation (DI) cycles offer an unrivaled opportunity to constrain properties of an active volcano's shallow magma reservoir and relation to its eruptive vents. We model transient behavior at Kīlauea Volcano using a simple mathematical model of an elastic reservoir that is coupled to magma flux through Kīlauea's East Rift Zone (ERZ) at a rate proportional to the difference in pressure between the summit reservoir and the ERZ eruptive vent (Newtonian flow). In this model, summit deflations and ERZ flux reductions are caused by a blockage in supply to the reservoir, while re-inflations occur as the system returns to a steady-state flux condition. The model naturally produces exponential variations in pressure and eruption rate which reasonably, albeit imperfectly, match observations during many of the transient events at Kīlauea. We constrain the model using a diverse range of observations including time-varying summit lava lake surface height and volume change, the temporal evolution of summit ground tilt, time-averaged eruption rate derived from TanDEM-X radar data, and height difference between the summit lava lake and the ERZ eruptive vent during brief eruptive pauses (Patrick et al., 2015). Formulating a Bayesian inverse and including independent prior constraint on magma density, host rock strength, and other properties of the system, we are able to place probabilistic constraints on the volume and volatile content of shallow magma storage, as well as properties of the ERZ conduit and influx of magma into Kīlauea's shallow magma reservoir. Reservoir influx parameters cannot in general be uniquely resolved, but reservoir volume and exsolved volatile content are well constrained; ERZ conduit radius may also be estimated given some simplifying assumptions.

High-Redshift SNe with Subaru and HST

NASA Astrophysics Data System (ADS)

Rubin, David; Suzuki, Nao; Regnault, Nicolas; Aldering, Gregory; Amanullah, Rahman; Antilogus, Pierre; Astier, Pierre; Barbary, Kyle; Betoule, Marc; Boone, Kyle Robert; Currie, Miles; Deustua, Susana; Doi, Mamoru; Fruchter, Andrew; Goobar, Ariel; Hayden, Brian; Hazenberg, Francois; Hook, Isobel; Huang, Xiaosheng; Jiang, Jian; Kato, Takahiro; Kim, Alex; Kowalski, Marek; Lidman, Chris; Linder, Eric; Maeda, Keiichi; Morokuma, Tomoki; Nordin, Jakob; Pain, Reynald; Perlmutter, Saul; Ruiz-Lapuente, Pilar; Sako, Masao; Myers Saunders, Clare; Spadafora, Anthony L.; Tanaka, Masaomi; Tominaga, Nozomu; Yasuda, Naoki; Yoshida, Naoki

2018-01-01

High-redshift type Ia supernovae are crucial for constraining any time variation in dark energy. Here, we present the first discoveries and light curves from the SUbaru Supernovae with Hubble Infrared (SUSHI) program, which combines high-redshift SN discoveries from the Subaru Strategic Program (SSP, as well as other Subaru time) with HST WFC3 IR followup. This program efficiently uses the wide field and high collecting area of Subaru Hyper Suprime-Cam for optical light curves, but still obtains a precision NIR color. We are on track to double the number of well-measured SNe Ia at z > 1.1, triggering on 23 SNe Ia in our first season.

The role of sexual selection and conflict in mediating among-population variation in mating strategies and sexually dimorphic traits in Sepsis punctum.

PubMed

Dmitriew, Caitlin; Blanckenhorn, Wolf U

2012-01-01

The black scavenger fly Sepsis punctum exhibits striking among-population variation in the direction and magnitude of sexual size dimorphism, modification to the male forelimb and pre-copulatory behaviour. In some populations, male-biased sexual size dimorphism is observed; in other, less dimorphic, populations males court prior to mating. Such variation in reproductive traits is of interest to evolutionary biologists because it has the potential to limit gene flow among populations, contributing to speciation. Here, we investigate whether large male body size and modified forefemur are associated with higher male mating success within populations, whether these traits are associated with higher mating success among populations, and if these traits carry viability costs that could constrain their response to sexual selection. Flies from five distinct populations were reared at high or low food, generating high and low quality males. The expression of body size, forelimb morphology and courtship rate were each greater at high food, but high food males experienced higher mating success or reduced latency to first copulation in only one of the populations. Among populations, overall mating success increased with the degree of male-bias in overall body size and forelimb modification, suggesting that these traits have evolved as a means of increasing male mating rate. The increased mating success observed in large-male populations raises the question of why variation in magnitude of dimorphism persists among populations. One reason may be that costs of producing a large size constrain the evolution of ever-larger males. We found no evidence that juvenile mortality under food stress was greater for large-male populations, but development time was considerably longer and may represent an important constraint in an ephemeral and competitive growth environment.

Determining the Ocean's Role on the Variable Gravity Field and Earth Rotation

NASA Technical Reports Server (NTRS)

Ponte, Rui M.; Frey, H. (Technical Monitor)

2000-01-01

A number of ocean models of different complexity have been used to study changes in the oceanic angular momentum (OAM) and mass fields and their relation to the variable Earth rotation and gravity field. Time scales examined range from seasonal to a few days. Results point to the importance of oceanic signals in driving polar motion, in particular the Chandler and annual wobbles. Results also show that oceanic signals have a measurable impact on length-of-day variations. Various circulation features and associated mass signals, including the North Pacific subtropical gyre, the equatorial currents, and the Antarctic Circumpolar Current play a significant role in oceanic angular momentum variability. The impact on OAM values of an optimization procedure that uses available data to constrain ocean model results was also tested for the first time. The optimization procedure yielded substantial changes, in OAM, related to adjustments in both motion and mass fields,as well as in the wind stress torques acting on the ocean. Constrained OAM values were found to yield noticeable improvements in the agreement with the observed Earth rotation parameters, particularly at the seasonal timescale.

APOSTLE: 11 TRANSIT OBSERVATIONS OF TrES-3b

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

Kundurthy, P.; Becker, A. C.; Agol, E.

2013-02-10

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 11 transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model, and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrES-3b are consistent with previous estimates to within the 2{sigma} confidence level. We improved the errors (by 10%-30%) on system parameters such as the orbital inclination (i {sub orb}), impact parameter (b), and stellar density ({rho}{sub *}) compared to previousmore » measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data shows no evidence for transit timing variations and our timing measurements are able to rule out super-Earth and gas giant companions in low-order mean motion resonance with TrES-3b.« less

Modeling slow-slip segmentation in Cascadia subduction zone constrained by tremor locations and gravity anomalies

NASA Astrophysics Data System (ADS)

Li, Duo; Liu, Yajing

2017-04-01

Along-strike segmentation of slow-slip events (SSEs) and nonvolcanic tremors in Cascadia may reflect heterogeneities of the subducting slab or overlying continental lithosphere. However, the nature behind this segmentation is not fully understood. We develop a 3-D model for episodic SSEs in northern and central Cascadia, incorporating both seismological and gravitational observations to constrain the heterogeneities in the megathrust fault properties. The 6 year automatically detected tremors are used to constrain the rate-state friction parameters. The effective normal stress at SSE depths is constrained by along-margin free-air and Bouguer gravity anomalies. The along-strike variation in the long-term plate convergence rate is also taken into consideration. Simulation results show five segments of ˜Mw6.0 SSEs spontaneously appear along the strike, correlated to the distribution of tremor epicenters. Modeled SSE recurrence intervals are equally comparable to GPS observations using both types of gravity anomaly constraints. However, the model constrained by free-air anomaly does a better job in reproducing the cumulative slip as well as more consistent surface displacements with GPS observations. The modeled along-strike segmentation represents the averaged slip release over many SSE cycles, rather than permanent barriers. Individual slow-slip events can still propagate across the boundaries, which may cause interactions between adjacent SSEs, as observed in time-dependent GPS inversions. In addition, the moment-duration scaling is sensitive to the selection of velocity criteria for determining when SSEs occur. Hence, the detection ability of the current GPS network should be considered in the interpretation of slow earthquake source parameter scaling relations.

Constraining the physical properties of compositionally distinctive surfaces on Mars from overlapping THEMIS observations

NASA Astrophysics Data System (ADS)

Ahern, A.; Rogers, D.

2017-12-01

Better constraints on the physical properties (e.g. grain size, rock abundance, cohesion, porosity and amount of induration) of Martian surface materials can lead to greater understanding of outcrop origin (e.g. via sedimentary, effusive volcanic, pyroclastic processes). Many outcrop surfaces on Mars likely contain near-surface (<3 cm) vertical heterogeneity in physical properties due to thin sediment cover, induration, and physical weathering, that can obscure measurement of the bulk thermal conductivity of the outcrop materials just below. Fortunately, vertical heterogeneity within near-surface materials can result in unique, and possibly predictable, diurnal and seasonal temperature patterns. The KRC thermal model has been utilized in a number of previous studies to predict thermal inertia of surface materials on Mars. Here we use KRC to model surface temperatures from overlapping Mars Odyssey THEMIS surface temperature observations that span multiple seasons and local times, in order to constrain both the nature of vertical heterogeneity and the underlying outcrop thermal inertia for various spectrally distinctive outcrops on Mars. We utilize spectral observations from TES and CRISM to constrain the particle size of the uppermost surface. For this presentation, we will focus specifically on chloride-bearing units in Terra Sirenum and Meridiani Planum, as well as mafic and feldspathic bedrock locations with distinct spectral properties, yet uncertain origins, in Noachis Terra and Nili Fossae. We find that many of these surfaces exhibit variations in apparent thermal inertia with season and local time that are consistent with low thermal inertia materials overlying higher thermal inertia substrates. Work is ongoing to compare surface temperature measurements with modeled two-layer scenarios in order to constrain the top layer thickness and bottom layer thermal inertia. The information will be used to better interpret the origins of these distinctive outcrops.

Coupling mammalian demography to climate through satellite time series of plant phenology

NASA Astrophysics Data System (ADS)

Stoner, D.; Sexton, J. O.; Nagol, J. R.; Ironside, K.; Choate, D.; Longshore, K.; Edwards, T., Jr.

2016-12-01

The seasonality of plant productivity governs the demography of primary and secondary consumers, and in arid ecosystems primary production is constrained by water availability. We relate the behavior, demography, and spatial distribution of large mammalian herbivores and their principal predator to remotely sensed indices of climate and vegetation across the western United States from 2000-2014. Terrain and plant community composition moderate the effects of climatological drought on primary productivity, resulting in spatial variation in ecosystem susceptibility to water stress. Herbivores track these patterns through habitat selection during key periods such as birthing and migration. Across a broad climatological gradient, timing of the start of growing season explains 75% of the variation in herbivore birth timing and 56% of the variation in neonatal survival rates. Initiation of autumn migration corresponds with the end of the growing season. Although indirectly coupled to primary production, carnivore home range size and population density are strongly correlated with plant productivity and growing-season length. Satellite measures of green reflectance during the peak of the growing season explain over 84% of the variation in carnivore home range size and 59% of the variation in density. Climate projections for the western United States predict warming temperatures and shifts in the timing and form of precipitation. Our analyses suggest that increased climatological variability will contribute to fluctuations in the composition and phenology of plant communities. These changes will propagate through consumer trophic levels, manifesting as increased home range area, shifts in the timing of migration, and greater volatility in large mammal populations. Combined with expansion and amplification of human land uses, these changes will likely have economic implications stemming from increased human-wildlife conflict and loss of ecosystem services.

An iterative algorithm for L1-TV constrained regularization in image restoration

NASA Astrophysics Data System (ADS)

Chen, K.; Loli Piccolomini, E.; Zama, F.

2015-11-01

We consider the problem of restoring blurred images affected by impulsive noise. The adopted method restores the images by solving a sequence of constrained minimization problems where the data fidelity function is the ℓ1 norm of the residual and the constraint, chosen as the image Total Variation, is automatically adapted to improve the quality of the restored images. Although this approach is general, we report here the case of vectorial images where the blurring model involves contributions from the different image channels (cross channel blur). A computationally convenient extension of the Total Variation function to vectorial images is used and the results reported show that this approach is efficient for recovering nearly optimal images.

Qatar-1: indications for possible transit timing variations

NASA Astrophysics Data System (ADS)

von Essen, C.; Schröter, S.; Agol, E.; Schmitt, J. H. M. M.

2013-07-01

Aims: Variations in the timing of transiting exoplanets provide a powerful tool for detecting additional planets in the system. Thus, the aim of this paper is to discuss the plausibility of transit timing variations (TTVs) on the Qatar-1 system by means of primary transit light curves analysis. Furthermore, we provide an interpretation of the timing variation. Methods: We observed Qatar-1 between March 2011 and October 2012 using the 1.2 m OLT telescope in Germany and the 0.6 m PTST telescope in Spain. We present 26 primary transits of the hot Jupiter Qatar-1b. In total, our light curves cover a baseline of 18 months. Results: We report on indications for possible long-term TTVs. Assuming that these TTVs are true, we present two different scenarios that could explain them. Our reported ~190 days TTV signal can be reproduced by either a weak perturber in resonance with Qatar-1b, or by a massive body in the brown dwarf regime. More observations and radial velocity monitoring are required to better constrain the perturber's characteristics. We also refine the ephemeris of Qatar-1b, which we find to be T0 = 2456157.42204 ± 0.0001 BJDTDB and P = 1.4200246 ± 0.0000007 days, and improve the system orbital parameters. Tables of the transit observations are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/555/A92

Relationship of D'' structure with the velocity variations near the inner-core boundary

NASA Astrophysics Data System (ADS)

Luo, Sheng-Nian; Ni, Sidao; Helmberger, Don

2002-06-01

Variations in regional differential times between PKiKP (i) and PKIKP (I) have been attributed to hemispheric P-velocity variations of about 1% in the upper 100 km of the inner core (referred to as HIC). The top of the inner core appears relatively fast beneath Asia where D'' is also fast. An alternative interpretation could be the lateral variation in P velocity at the lowermost outer core (HOC) producing the same differential times. To resolve this issue, we introduce the diffracted PKP phase near the B caustic (Bdiff) in the range of 139-145° epicenter distances, and the corresponding differential times between Bdiff and PKiKP and PKIKP as observed on broadband arrays. Due to the long-wavelength nature of Bdiff, we scaled the S-wave tomography model with k values (k ≡ dlnVs/dlnVp) to obtain large-scale P-wave velocity structure in the lower mantle as proposed by earlier studies. Waveform synthetics of Bdiff constructed with small k's predict complex waveforms not commonly observed, confirming the validity of large scaling factor k. With P-velocity in lower mantle constrained at large scale, the extra travel-time constraint imposed by Bdiff helps to resolve the HOC-HIC issue. Our preliminary results suggest k > 2 for the lowermost mantle and support HIC hypothesis. An important implication is that there appears to be a relationship of D'' velocity structures with the structures near the inner core boundary via core dynamics.

A MAP blind image deconvolution algorithm with bandwidth over-constrained

NASA Astrophysics Data System (ADS)

Ren, Zhilei; Liu, Jin; Liang, Yonghui; He, Yulong

2018-03-01

We demonstrate a maximum a posteriori (MAP) blind image deconvolution algorithm with bandwidth over-constrained and total variation (TV) regularization to recover a clear image from the AO corrected images. The point spread functions (PSFs) are estimated by bandwidth limited less than the cutoff frequency of the optical system. Our algorithm performs well in avoiding noise magnification. The performance is demonstrated on simulated data.

A Coordinated X-Ray and Optical Campaign of the Nearest Massive Eclipsing Binary, δ Orionis Aa. III. Analysis of Optical Photometric (MOST) and Spectroscopic (Ground-based) Variations

NASA Astrophysics Data System (ADS)

Pablo, Herbert; Richardson, Noel D.; Moffat, Anthony F. J.; Corcoran, Michael; Shenar, Tomer; Benvenuto, Omar; Fuller, Jim; Nazé, Yaël; Hoffman, Jennifer L.; Miroshnichenko, Anatoly; Maíz Apellániz, Jesús; Evans, Nancy; Eversberg, Thomas; Gayley, Ken; Gull, Ted; Hamaguchi, Kenji; Hamann, Wolf-Rainer; Henrichs, Huib; Hole, Tabetha; Ignace, Richard; Iping, Rosina; Lauer, Jennifer; Leutenegger, Maurice; Lomax, Jamie; Nichols, Joy; Oskinova, Lida; Owocki, Stan; Pollock, Andy; Russell, Christopher M. P.; Waldron, Wayne; Buil, Christian; Garrel, Thierry; Graham, Keith; Heathcote, Bernard; Lemoult, Thierry; Li, Dong; Mauclaire, Benjamin; Potter, Mike; Ribeiro, Jose; Matthews, Jaymie; Cameron, Chris; Guenther, David; Kuschnig, Rainer; Rowe, Jason; Rucinski, Slavek; Sasselov, Dimitar; Weiss, Werner

2015-08-01

We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system δ Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV P\\gt 400 years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations.

Constrained variation in Jastrow method at high density

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

Owen, J.C.; Bishop, R.F.; Irvine, J.M.

1976-11-01

A method is derived for constraining the correlation function in a Jastrow variational calculation which permits the truncation of the cluster expansion after two-body terms, and which permits exact minimization of the two-body cluster by functional variation. This method is compared with one previously proposed by Pandharipande and is found to be superior both theoretically and practically. The method is tested both on liquid /sup 3/He, by using the Lennard--Jones potential, and on the model system of neutrons treated as Boltzmann particles (''homework'' problem). Good agreement is found both with experiment and with other calculations involving the explicit evaluation ofmore » higher-order terms in the cluster expansion. The method is then applied to a more realistic model of a neutron gas up to a density of 4 neutrons per F/sup 3/, and is found to give ground-state energies considerably lower than those of Pandharipande. (AIP)« less

Dynamics of non-holonomic systems with stochastic transport

NASA Astrophysics Data System (ADS)

Holm, D. D.; Putkaradze, V.

2018-01-01

This paper formulates a variational approach for treating observational uncertainty and/or computational model errors as stochastic transport in dynamical systems governed by action principles under non-holonomic constraints. For this purpose, we derive, analyse and numerically study the example of an unbalanced spherical ball rolling under gravity along a stochastic path. Our approach uses the Hamilton-Pontryagin variational principle, constrained by a stochastic rolling condition, which we show is equivalent to the corresponding stochastic Lagrange-d'Alembert principle. In the example of the rolling ball, the stochasticity represents uncertainty in the observation and/or error in the computational simulation of the angular velocity of rolling. The influence of the stochasticity on the deterministically conserved quantities is investigated both analytically and numerically. Our approach applies to a wide variety of stochastic, non-holonomically constrained systems, because it preserves the mathematical properties inherited from the variational principle.

Pool desiccation and developmental thresholds in the common frog, Rana temporaria.

PubMed

Lind, Martin I; Persbo, Frida; Johansson, Frank

2008-05-07

The developmental threshold is the minimum size or condition that a developing organism must have reached in order for a life-history transition to occur. Although developmental thresholds have been observed for many organisms, inter-population variation among natural populations has not been examined. Since isolated populations can be subjected to strong divergent selection, population divergence in developmental thresholds can be predicted if environmental conditions favour fast or slow developmental time in different populations. Amphibian metamorphosis is a well-studied life-history transition, and using a common garden approach we compared the development time and the developmental threshold of metamorphosis in four island populations of the common frog Rana temporaria: two populations originating from islands with only temporary breeding pools and two from islands with permanent pools. As predicted, tadpoles from time-constrained temporary pools had a genetically shorter development time than those from permanent pools. Furthermore, the variation in development time among females from temporary pools was low, consistent with the action of selection on rapid development in this environment. However, there were no clear differences in the developmental thresholds between the populations, indicating that the main response to life in a temporary pool is to shorten the development time.

Using surface displacement derived from GRACE to constrain the water loading signal in cGPS measurements in the Amazon Basin

NASA Astrophysics Data System (ADS)

Jose, L.; Bennett, R. A.; Harig, C.

2017-12-01

Currently, cGPS data is well suited to track vertical changes in the Earth's surface. However, there are annual, semi-annual, and interannual signals within cGPS time series that are not well constrained. We hypothesize that these signals are primarily due to water loading. If this is the case, the conventional method of modeling cGPS data as an annual or semiannual sinusoid falls short, as such models cannot accurately capture all variations in surface displacement, especially those due to extreme hydrologic events. We believe that we can better correct the cGPS time series with another method we are developing wherein we use a time series of surface displacement derived from the GRACE geopotential field instead of a sinusoidal model to correct the data. Currently, our analysis is constrained to the Amazon Basin, where the signal due to water loading is large enough to appear in both the GRACE and cGPS measurements. The vertical signal from cGPS stations across the Amazon Basin show an apparent spatial correlation, which further supports our idea that these signals are due to a regional water loading signal. In our preliminary research, we used tsview for Matlab to find that the WRMS of the corrected cGPS time series can be reduced as much as 30% from the model corrected data to the GRACE corrected data. The Amazon, like many places around the world, has experienced extreme drought, in 2005, 2010, and recently in 2015. In addition to making the cGPS vertical signal more robust, the method we are developing has the potential to help us understand the effects of these weather events and track trends in water loading.

Le Chatelier's principle with multiple relaxation channels

NASA Astrophysics Data System (ADS)

Gilmore, R.; Levine, R. D.

1986-05-01

Le Chatelier's principle is discussed within the constrained variational approach to thermodynamics. The formulation is general enough to encompass systems not in thermal (or chemical) equilibrium. Particular attention is given to systems with multiple constraints which can be relaxed. The moderation of the initial perturbation increases as additional constraints are removed. This result is studied in particular when the (coupled) relaxation channels have widely different time scales. A series of inequalities is derived which describes the successive moderation as each successive relaxation channel opens up. These inequalities are interpreted within the metric-geometry representation of thermodynamics.

Fracture characterization from near-offset VSP inversion

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

Horne, S.; MacBeth, C.; Queen, J.

1997-01-01

A global optimization method incorporating a ray-tracing scheme is used to invert observations of shear-wave splitting from two near-offset VSPs recorded at the Conoco Borehole Test Facility, Kay County, Oklahoma. Inversion results suggest that the seismic anisotropy is due to a non-vertical fracture system. This interpretation is constrained by the VSP acquisition geometry for which two sources are employed along near diametrically opposite azimuths about the well heads. A correlation is noted between the time-delay variations between the fast and slow split shear waves and the sandstone formations.

Control of linear uncertain systems utilizing mismatched state observers

NASA Technical Reports Server (NTRS)

Goldstein, B.

1972-01-01

The control of linear continuous dynamical systems is investigated as a problem of limited state feedback control. The equations which describe the structure of an observer are developed constrained to time-invarient systems. The optimal control problem is formulated, accounting for the uncertainty in the design parameters. Expressions for bounds on closed loop stability are also developed. The results indicate that very little uncertainty may be tolerated before divergence occurs in the recursive computation algorithms, and the derived stability bound yields extremely conservative estimates of regions of allowable parameter variations.

The least-squares mixing models to generate fraction images derived from remote sensing multispectral data

NASA Technical Reports Server (NTRS)

Shimabukuro, Yosio Edemir; Smith, James A.

1991-01-01

Constrained-least-squares and weighted-least-squares mixing models for generating fraction images derived from remote sensing multispectral data are presented. An experiment considering three components within the pixels-eucalyptus, soil (understory), and shade-was performed. The generated fraction images for shade (shade image) derived from these two methods were compared by considering the performance and computer time. The derived shade images are related to the observed variation in forest structure, i.e., the fraction of inferred shade in the pixel is related to different eucalyptus ages.

Constraints on geomagnetic secular variation modeling from electromagnetism and fluid dynamics of the Earth's core

NASA Technical Reports Server (NTRS)

Benton, E. R.

1986-01-01

A spherical harmonic representation of the geomagnetic field and its secular variation for epoch 1980, designated GSFC(9/84), is derived and evaluated. At three epochs (1977.5, 1980.0, 1982.5) this model incorporates conservation of magnetic flux through five selected patches of area on the core/mantle boundary bounded by the zero contours of vertical magnetic field. These fifteen nonlinear constraints are included like data in an iterative least squares parameter estimation procedure that starts with the recently derived unconstrained field model GSFC (12/83). Convergence is approached within three iterations. The constrained model is evaluated by comparing its predictive capability outside the time span of its data, in terms of residuals at magnetic observatories, with that for the unconstrained model.

Investigating the physical properties of transiting hot Jupiters with the 1.5-m Kuiper Telescope

NASA Astrophysics Data System (ADS)

Turner, Jake D.; Leiter, Robin M.; Biddle, Lauren I.; Pearson, Kyle A.; Hardegree-Ullman, Kevin K.; Thompson, Robert M.; Teske, Johanna K.; Cates, Ian T.; Cook, Kendall L.; Berube, Michael P.; Nieberding, Megan N.; Jones, Christen K.; Raphael, Brandon; Wallace, Spencer; Watson, Zachary T.; Johnson, Robert E.

2017-12-01

We present new photometric data of 11 hot Jupiter transiting exoplanets (CoRoT-12b, HAT-P-5b, HAT-P-12b, HAT-P-33b, HAT-P-37b, WASP-2b, WASP-24b, WASP-60b, WASP-80b, WASP-103b and XO-3b) in order to update their planetary parameters and to constrain information about their atmospheres. These observations of CoRoT-12b, HAT-P-37b and WASP-60b are the first follow-up data since their discovery. Additionally, the first near-UV transits of WASP-80b and WASP-103b are presented. We compare the results of our analysis with previous work to search for transit timing variations (TTVs) and a wavelength dependence in the transit depth. TTVs may be evidence of a third body in the system, and variations in planetary radius with wavelength can help constrain the properties of the exoplanet's atmosphere. For WASP-103b and XO-3b, we find a possible variation in the transit depths which may be evidence of scattering in their atmospheres. The B-band transit depth of HAT-P-37b is found to be smaller than its near-IR transit depth and such a variation may indicate TiO/VO absorption. These variations are detected from 2-4.6σ, so follow-up observations are needed to confirm these results. Additionally, a flat spectrum across optical wavelengths is found for five of the planets (HAT-P-5b, HAT-P-12b, WASP-2b, WASP-24b and WASP-80b), suggestive that clouds may be present in their atmospheres. We calculate a refined orbital period and ephemeris for all the targets, which will help with future observations. No TTVs are seen in our analysis with the exception of WASP-80b and follow-up observations are needed to confirm this possible detection.

Towards high accuracy tests on the substellar IMF in young clusters. A survey in NGC 2024.

NASA Astrophysics Data System (ADS)

Da Rio, Nicola

2017-08-01

Measuring the Initial Mass Function in young clusters, and testing its universality, is a fundamental benchmark to constrain the physical processes and theoretical models of star formation. The shape and universality of the stellar IMF are well known. Our observational characterization of the substellar IMF, on the other hand, remains more uncertain, along with its possible environmental variations. Because of this, the physical processes that play a role in the formation of brown dwarfs are not fully constrained. In Cycle 22 we were awarded HST time to carry out the deepest spectro-photometric census of BDs in a young cluster: the Orion Nebula Cluster. Through deep WFC3/IR narrow band imaging, we are able to obtain Teff and A_V down to 15Mjup. Preliminary analysis limited to a portion of the total field of view allows us to classify several hundreds BDs, place them in the HRD and obtain, for an extinction limited sample, the complete and consistent IMF down to planetary masses. The substellar slope is consistent with the Galactic IMF but a rapid drop is found at the H-burning limit. We propose to carry out a nearly identical survey with HST in a younger, less massive nearby cluster: NGC2024 in the Flame Nebula. This will allow us to derive the complete census of the young population down to planetary masses, derive the IMF, enabling a consistent comparison with the results in the ONC. We will specifically look for statistically significant IMF variations with environmental properties (cluster mass, density) and investigate primordial mass segregation in the substellar regime. These results will significantly help to constrain the mechanisms involved in BD formation.

Heritability, covariation and natural selection on 24 traits of common evening primrose (Oenothera biennis) from a field experiment.

PubMed

Johnson, M T J; Agrawal, A A; Maron, J L; Salminen, J-P

2009-06-01

This study explored genetic variation and co-variation in multiple functional plant traits. Our goal was to characterize selection, heritabilities and genetic correlations among different types of traits to gain insight into the evolutionary ecology of plant populations and their interactions with insect herbivores. In a field experiment, we detected significant heritable variation for each of 24 traits of Oenothera biennis and extensive genetic covariance among traits. Traits with diverse functions formed several distinct groups that exhibited positive genetic covariation with each other. Genetic variation in life-history traits and secondary chemistry together explained a large proportion of variation in herbivory (r(2) = 0.73). At the same time, selection acted on lifetime biomass, life-history traits and two secondary compounds of O. biennis, explaining over 95% of the variation in relative fitness among genotypes. The combination of genetic covariances and directional selection acting on multiple traits suggests that adaptive evolution of particular traits is constrained, and that correlated evolution of groups of traits will occur, which is expected to drive the evolution of increased herbivore susceptibility. As a whole, our study indicates that an examination of genetic variation and covariation among many different types of traits can provide greater insight into the evolutionary ecology of plant populations and plant-herbivore interactions.

A Variational Assimilation Method for Satellite and Conventional Data: Model 2 (version 1)

NASA Technical Reports Server (NTRS)

Achtemeier, Gary L.

1991-01-01

The Model II variational data assimilation model is the second of the four variational models designed to blend diverse meteorological data into a dynamically constrained data set. Model II differs from Model I in that it includes the thermodynamic equation as the fifth dynamical constraint. Thus, Model II includes all five of the primative equations that govern atmospheric flow for a dry atmosphere.

Geodynamic contributions to global climatic change

NASA Technical Reports Server (NTRS)

Bills, Bruce G.

1992-01-01

Orbital and rotational variations perturb the latitudinal and seasonal pattern of incident solar radiation, producing major climatic change on time scales of 10(exp 4)-10(exp 6) years. The orbital variations are oblivious to internal structure and processes, but the rotational variations are not. A program of investigation whose objective would be to explore and quantify three aspects of orbital, rotational, and climatic interactions is described. An important premise of this investigation is the synergism between geodynamics and paleoclimate. Better geophysical models of precessional dynamics are needed in order to accurately reconstruct the radiative input to climate models. Some of the paleoclimate proxy records contain information relevant to solid Earth processes, on time scales which are difficult to constrain otherwise. Specific mechanisms which will be addressed include: (1) climatic consequences of deglacial polar motion; and (2) precessional and climatic consequences of glacially induced perturbations in the gravitational oblateness and partial decoupling of the mantle and core. The approach entails constructing theoretical models of the rotational, deformational, radiative, and climatic response of the Earth to known orbital perturbations, and comparing these with extensive records of paleoclimate proxy data. Several of the mechanisms of interest may participate in previously unrecognized feed-back loops in the climate dynamics system. A new algorithm for estimating climatically diagnostic locations and seasons from the paleoclimate time series is proposed.

Precision Timing of PSR J0437-4715: An Accurate Pulsar Distance, a High Pulsar Mass, and a Limit on the Variation of Newton's Gravitational Constant

NASA Astrophysics Data System (ADS)

Verbiest, J. P. W.; Bailes, M.; van Straten, W.; Hobbs, G. B.; Edwards, R. T.; Manchester, R. N.; Bhat, N. D. R.; Sarkissian, J. M.; Jacoby, B. A.; Kulkarni, S. R.

2008-05-01

Analysis of 10 years of high-precision timing data on the millisecond pulsar PSR J0437-4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, dot Pb , determined at the 1.5% level of precision (Dk = 157.0 +/- 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 solar system ephemerides. The precise measurement of dot Pb allows a limit on the variation of Newton's gravitational constant, |Ġ/G| <= 23 × 10-12 yr-1. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a⊙/c| <= 1.5 × 10-18 s-1 at 95% confidence, and derive a pulsar mass, mpsr = 1.76 +/- 0.20 M⊙, one of the highest estimates so far obtained.

Transit timing variations for planets co-orbiting in the horseshoe regime

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

Vokrouhlický, David; Nesvorný, David, E-mail: vokrouhl@cesnet.cz, E-mail: davidn@boulder.swri.edu

2014-08-10

Although not yet detected, pairs of exoplanets in 1:1 mean motion resonance probably exist. Low eccentricity, near-planar orbits, which in the comoving frame follow horseshoe trajectories, are one of the possible stable configurations. Here we study transit timing variations (TTVs) produced by mutual gravitational interaction of planets in this orbital architecture, with the goal to develop methods that can be used to recognize this case in observational data. In particular, we use a semi-analytic model to derive parametric constraints that should facilitate data analysis. We show that characteristic traits of the TTVs can directly constrain the (1) ratio of planetarymore » masses and (2) their total mass (divided by that of the central star) as a function of the minimum angular separation as seen from the star. In an ideal case, when transits of both planets are observed and well characterized, the minimum angular separation can also be inferred from the data. As a result, parameters derived from the observed transit timing series alone can directly provide both planetary masses scaled to the central star mass.« less

Perspective: tipping the scales: search for drifting constants from molecular spectra.

PubMed

Jansen, Paul; Bethlem, Hendrick L; Ubachs, Wim

2014-01-07

Transitions in atoms and molecules provide an ideal test ground for constraining or detecting a possible variation of the fundamental constants of nature. In this perspective, we review molecular species that are of specific interest in the search for a drifting proton-to-electron mass ratio μ. In particular, we outline the procedures that are used to calculate the sensitivity coefficients for transitions in these molecules and discuss current searches. These methods have led to a rate of change in μ bounded to 6 × 10(-14)/yr from a laboratory experiment performed in the present epoch. On a cosmological time scale, the variation is limited to ∣Δμ∕μ∣ < 10(-5) for look-back times of 10-12× 10(9) years and to ∣Δμ∕μ∣ < 10(-7) for look-back times of 7× 10(9) years. The last result, obtained from high-redshift observation of methanol, translates into μ̇/μ=(1.4±1.4)×10(-17)/yr if a linear rate of change is assumed.

Variational method based on Retinex with double-norm hybrid constraints for uneven illumination correction

NASA Astrophysics Data System (ADS)

Li, Shuo; Wang, Hui; Wang, Liyong; Yu, Xiangzhou; Yang, Le

2018-01-01

The uneven illumination phenomenon reduces the quality of remote sensing image and causes interference in the subsequent processing and applications. A variational method based on Retinex with double-norm hybrid constraints for uneven illumination correction is proposed. The L1 norm and the L2 norm are adopted to constrain the textures and details of reflectance image and the smoothness of the illumination image, respectively. The problem of separating the illumination image from the reflectance image is transformed into the optimal solution of the variational model. In order to accelerate the solution, the split Bregman method is used to decompose the variational model into three subproblems, which are calculated by alternate iteration. Two groups of experiments are implemented on two synthetic images and three real remote sensing images. Compared with the variational Retinex method with single-norm constraint and the Mask method, the proposed method performs better in both visual evaluation and quantitative measurements. The proposed method can effectively eliminate the uneven illumination while maintaining the textures and details of the remote sensing image. Moreover, the proposed method using split Bregman method is more than 10 times faster than the method with the steepest descent method.

Heat resistance throughout ontogeny: body size constrains thermal tolerance.

PubMed

Klockmann, Michael; Günter, Franziska; Fischer, Klaus

2017-02-01

Heat tolerance is a trait of paramount ecological importance and may determine a species' ability to cope with ongoing climate change. Although critical thermal limits have consequently received substantial attention in recent years, their potential variation throughout ontogeny remained largely neglected. We investigate whether such neglect may bias conclusions regarding a species' sensitivity to climate change. Using a tropical butterfly, we found that developmental stages clearly differed in heat tolerance. It was highest in pupae followed by larvae, adults and finally eggs and hatchlings. Strikingly, most of the variation found in thermal tolerance was explained by differences in body mass, which may thus impose a severe constraint on adaptive variation in stress tolerance. Furthermore, temperature acclimation was beneficial by increasing heat knock-down time and therefore immediate survival under heat stress, but it affected reproduction negatively. Extreme temperatures strongly reduced survival and subsequent reproductive success even in our highly plastic model organism, exemplifying the potentially dramatic impact of extreme weather events on biodiversity. We argue that predictions regarding a species' fate under changing environmental conditions should consider variation in thermal tolerance throughout ontogeny, variation in body mass and acclimation responses as important predictors of stress tolerance. © 2016 John Wiley & Sons Ltd.

Reconstructing spatial and temporal patterns of paleoglaciation across Central Asia

NASA Astrophysics Data System (ADS)

Stroeven, Arjen P.

2014-05-01

Understanding the behaviour of mountain glaciers and ice caps, the evolution of mountain landscapes, and testing global climate models all require well-constrained information on past spatial and temporal patterns of glacier change. Particularly important are transitional regions that have high spatial and temporal variation in glacier activity and that can provide a sensitive record of past climate change. Central Asia is an extreme continental location with glaciers that have responded sensitively to variations in major regional climate systems. As an international team, we are reconstructing glacial histories of several areas of the Tibetan Plateau as well as along the Tian Shan, Altai and Kunlun Mountains. Building on previous work, we are using remote sensing-based geomorphological mapping augmented with field observations to map out glacial landforms and the maximum distributions of erratics. We then use cosmogenic nuclide Be-10 and Al-26, optically stimulated luminescence, and electron spin resonance dating of moraines and other landforms to compare dating techniques and to constrain the ages of defined extents of paleo-glaciers and ice caps. Comparing consistently dated glacial histories across central Asia provides an opportunity to examine shifts in the dominance patterns of climate systems over time in the region. Results to date show significant variations in the timing and extent of glaciation, including areas in the southeast Tibetan Plateau and Tian Shan with extensive valley and small polythermal ice cap glaciation during the global last glacial maximum in contrast to areas in central and northeast Tibetan Plateau that had very limited valley glacier expansion then. Initial numerical modelling attempting to simulate mapped and dated paleoglacial extents indicates that relatively limited cooling is sufficient to produce observed past expansions of glaciers across the Tibetan Plateau, and predicts complex basal thermal regimes in some locations that match patterns of past glacial erosion inferred from landform patterns and ages. Future modelling will examine glacier behaviour along major mountain ranges across central Asia.

An Efficient Augmented Lagrangian Method with Applications to Total Variation Minimization

DTIC Science & Technology

2012-08-17

the classic augmented Lagrangian multiplier method, we propose, analyze and test an algorithm for solving a class of equality-constrained non-smooth...method, we propose, analyze and test an algorithm for solving a class of equality-constrained non-smooth optimization problems (chie y but not...significantly outperforming several state-of-the-art solvers on most tested problems. The resulting MATLAB solver, called TVAL3, has been posted online [23]. 2

Implications of ecological energetics and biophysical and developmental constraints for life history variation in dinosaurs

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

Dunham, A.E.; Overall, K.L.; Forster, C.A.

1988-01-01

There has been much recent speculation concerning the nature of life history variation in dinosaurs (Case, 1978; Bakker, 1986; Horner, 1982, 1984a). The purpose of this paper is to review the data on dinosaur life histories and to examine the nature and magnitude of the demographic and physiological factors that must have constrained life history variation in this group. 145 refs., 8 figs., 3 tabs.

Constrained and Unconstrained Variational Finite Element Formulation of Solutions to a Stress Wave Problem - a Numerical Comparison.

DTIC Science & Technology

1982-10-01

Element Unconstrained Variational Formulations," Innovativ’e Numerical Analysis For the Applied Engineering Science, R. P. Shaw, et at, Fitor...Initial Boundary Value of Gun Dynamics Solved by Finite Element Unconstrained Variational Formulations," Innovative Numerical Analysis For the Applied ... Engineering Science, R. P. Shaw, et al, Editors, University Press of Virginia, Charlottesville, pp. 733-741, 1980. 2 J. J. Wu, "Solutions to Initial

Mantle viscosity structure constrained by joint inversions of seismic velocities and density

NASA Astrophysics Data System (ADS)

Rudolph, M. L.; Moulik, P.; Lekic, V.

2017-12-01

The viscosity structure of Earth's deep mantle affects the thermal evolution of Earth, the ascent of mantle upwellings, sinking of subducted oceanic lithosphere, and the mixing of compositional heterogeneities in the mantle. Modeling the long-wavelength dynamic geoid allows us to constrain the radial viscosity profile of the mantle. Typically, in inversions for the mantle viscosity structure, wavespeed variations are mapped into density variations using a constant- or depth-dependent scaling factor. Here, we use a newly developed joint model of anisotropic Vs, Vp, density and transition zone topographies to generate a suite of solutions for the mantle viscosity structure directly from the seismologically constrained density structure. The density structure used to drive our forward models includes contributions from both thermal and compositional variations, including important contributions from compositionally dense material in the Large Low Velocity Provinces at the base of the mantle. These compositional variations have been neglected in the forward models used in most previous inversions and have the potential to significantly affect large-scale flow and thus the inferred viscosity structure. We use a transdimensional, hierarchical, Bayesian approach to solve the inverse problem, and our solutions for viscosity structure include an increase in viscosity below the base of the transition zone, in the shallow lower mantle. Using geoid dynamic response functions and an analysis of the correlation between the observed geoid and mantle structure, we demonstrate the underlying reason for this inference. Finally, we present a new family of solutions in which the data uncertainty is accounted for using covariance matrices associated with the mantle structure models.

Intelligent Sampling of Hazardous Particle Populations in Resource-Constrained Environments

NASA Astrophysics Data System (ADS)

McCollough, J. P.; Quinn, J. M.; Starks, M. J.; Johnston, W. R.

2017-10-01

Sampling of anomaly-causing space environment drivers is necessary for both real-time operations and satellite design efforts, and optimizing measurement sampling helps minimize resource demands. Relating these measurements to spacecraft anomalies requires the ability to resolve spatial and temporal variability in the energetic charged particle hazard of interest. Here we describe a method for sampling particle fluxes informed by magnetospheric phenomenology so that, along a given trajectory, the variations from both temporal dynamics and spatial structure are adequately captured while minimizing oversampling. We describe the coordinates, sampling method, and specific regions and parameters employed. We compare resulting sampling cadences with data from spacecraft spanning the regions of interest during a geomagnetically active period, showing that the algorithm retains the gross features necessary to characterize environmental impacts on space systems in diverse orbital regimes while greatly reducing the amount of sampling required. This enables sufficient environmental specification within a resource-constrained context, such as limited telemetry bandwidth, processing requirements, and timeliness.

A Model for Engaging Students in a Research Experience Involving Variational Techniques, Mathematica, and Descent Methods.

ERIC Educational Resources Information Center

Mahavier, W. Ted

2002-01-01

Describes a two-semester numerical methods course that serves as a research experience for undergraduate students without requiring external funding or the modification of current curriculum. Uses an engineering problem to introduce students to constrained optimization via a variation of the traditional isoperimetric problem of finding the curve…

Spatial scale of deformation constrained by combinations of InSAR and GPS observations in Southern California

NASA Astrophysics Data System (ADS)

Lohman, R. B.; Scott, C. P.

2014-12-01

Efforts to understand the buildup and release of strain within the Earth's crust often rely on well-characterized observations of ground deformation, over time scales that include interseismic periods, earthquakes, and transient deformation episodes. Constraints on current rates of surface deformation in 1-, 2- or 3-dimensions can be obtained by examining sets of GPS and Interferometric Synthetic Aperture Radar (InSAR) observations, both alone and in combination. Contributions to the observed signal often include motion along faults, seasonal cycles of subsidence and recharge associated with aquifers, anthropogenic extraction of hydrocarbons, and variations in atmospheric water vapor and ionospheric properties. Here we examine methods for extracting time-varying ground deformation signals from combinations of InSAR and GPS data, real and synthetic, applied to Southern California. We show that two methods for combining the data through removal of a GPS-constrained function (a plane, and filtering) from the InSAR result in a clear tradeoff between the contribution from the two datatypes at diffferent spatial scales. We also show that the contribution to the secular rates at GPS sites from seasonal signals is large enough to be a significant error in this estimation process, and should be accounted for.

The effects of sample size on population genomic analyses--implications for the tests of neutrality.

PubMed

Subramanian, Sankar

2016-02-20

One of the fundamental measures of molecular genetic variation is the Watterson's estimator (θ), which is based on the number of segregating sites. The estimation of θ is unbiased only under neutrality and constant population growth. It is well known that the estimation of θ is biased when these assumptions are violated. However, the effects of sample size in modulating the bias was not well appreciated. We examined this issue in detail based on large-scale exome data and robust simulations. Our investigation revealed that sample size appreciably influences θ estimation and this effect was much higher for constrained genomic regions than that of neutral regions. For instance, θ estimated for synonymous sites using 512 human exomes was 1.9 times higher than that obtained using 16 exomes. However, this difference was 2.5 times for the nonsynonymous sites of the same data. We observed a positive correlation between the rate of increase in θ estimates (with respect to the sample size) and the magnitude of selection pressure. For example, θ estimated for the nonsynonymous sites of highly constrained genes (dN/dS < 0.1) using 512 exomes was 3.6 times higher than that estimated using 16 exomes. In contrast this difference was only 2 times for the less constrained genes (dN/dS > 0.9). The results of this study reveal the extent of underestimation owing to small sample sizes and thus emphasize the importance of sample size in estimating a number of population genomic parameters. Our results have serious implications for neutrality tests such as Tajima D, Fu-Li D and those based on the McDonald and Kreitman test: Neutrality Index and the fraction of adaptive substitutions. For instance, use of 16 exomes produced 2.4 times higher proportion of adaptive substitutions compared to that obtained using 512 exomes (24% vs 10 %).

The exit-time problem for a Markov jump process

NASA Astrophysics Data System (ADS)

Burch, N.; D'Elia, M.; Lehoucq, R. B.

2014-12-01

The purpose of this paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developed nonlocal vector calculus. This calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.

Constraining DALECv2 using multiple data streams and ecological constraints: analysis and application

DOE PAGES

Delahaies, Sylvain; Roulstone, Ian; Nichols, Nancy

2017-07-10

We use a variational method to assimilate multiple data streams into the terrestrial ecosystem carbon cycle model DALECv2 (Data Assimilation Linked Ecosystem Carbon). Ecological and dynamical constraints have recently been introduced to constrain unresolved components of this otherwise ill-posed problem. We recast these constraints as a multivariate Gaussian distribution to incorporate them into the variational framework and we demonstrate their advantage through a linear analysis. By using an adjoint method we study a linear approximation of the inverse problem: firstly we perform a sensitivity analysis of the different outputs under consideration, and secondly we use the concept of resolution matricesmore » to diagnose the nature of the ill-posedness and evaluate regularisation strategies. We then study the non-linear problem with an application to real data. Finally, we propose a modification to the model: introducing a spin-up period provides us with a built-in formulation of some ecological constraints which facilitates the variational approach.« less

Constraining DALECv2 using multiple data streams and ecological constraints: analysis and application

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

Delahaies, Sylvain; Roulstone, Ian; Nichols, Nancy

We use a variational method to assimilate multiple data streams into the terrestrial ecosystem carbon cycle model DALECv2 (Data Assimilation Linked Ecosystem Carbon). Ecological and dynamical constraints have recently been introduced to constrain unresolved components of this otherwise ill-posed problem. We recast these constraints as a multivariate Gaussian distribution to incorporate them into the variational framework and we demonstrate their advantage through a linear analysis. By using an adjoint method we study a linear approximation of the inverse problem: firstly we perform a sensitivity analysis of the different outputs under consideration, and secondly we use the concept of resolution matricesmore » to diagnose the nature of the ill-posedness and evaluate regularisation strategies. We then study the non-linear problem with an application to real data. Finally, we propose a modification to the model: introducing a spin-up period provides us with a built-in formulation of some ecological constraints which facilitates the variational approach.« less

Geometric constrained variational calculus I: Piecewise smooth extremals

NASA Astrophysics Data System (ADS)

Massa, Enrico; Bruno, Danilo; Luria, Gianvittorio; Pagani, Enrico

2015-05-01

A geometric setup for constrained variational calculus is presented. The analysis deals with the study of the extremals of an action functional defined on piecewise differentiable curves, subject to differentiable, non-holonomic constraints. Special attention is paid to the tensorial aspects of the theory. As far as the kinematical foundations are concerned, a fully covariant scheme is developed through the introduction of the concept of infinitesimal control. The standard classification of the extremals into normal and abnormal ones is discussed, pointing out the existence of an algebraic algorithm assigning to each admissible curve a corresponding abnormality index, related to the co-rank of a suitable linear map. Attention is then shifted to the study of the first variation of the action functional. The analysis includes a revisitation of Pontryagin's equations and of the Lagrange multipliers method, as well as a reformulation of Pontryagin's algorithm in Hamiltonian terms. The analysis is completed by a general result, concerning the existence of finite deformations with fixed endpoints.

Thermodynamics constrains allometric scaling of optimal development time in insects.

PubMed

Dillon, Michael E; Frazier, Melanie R

2013-01-01

Development time is a critical life-history trait that has profound effects on organism fitness and on population growth rates. For ectotherms, development time is strongly influenced by temperature and is predicted to scale with body mass to the quarter power based on 1) the ontogenetic growth model of the metabolic theory of ecology which describes a bioenergetic balance between tissue maintenance and growth given the scaling relationship between metabolism and body size, and 2) numerous studies, primarily of vertebrate endotherms, that largely support this prediction. However, few studies have investigated the allometry of development time among invertebrates, including insects. Abundant data on development of diverse insects provides an ideal opportunity to better understand the scaling of development time in this ecologically and economically important group. Insects develop more quickly at warmer temperatures until reaching a minimum development time at some optimal temperature, after which development slows. We evaluated the allometry of insect development time by compiling estimates of minimum development time and optimal developmental temperature for 361 insect species from 16 orders with body mass varying over nearly 6 orders of magnitude. Allometric scaling exponents varied with the statistical approach: standardized major axis regression supported the predicted quarter-power scaling relationship, but ordinary and phylogenetic generalized least squares did not. Regardless of the statistical approach, body size alone explained less than 28% of the variation in development time. Models that also included optimal temperature explained over 50% of the variation in development time. Warm-adapted insects developed more quickly, regardless of body size, supporting the "hotter is better" hypothesis that posits that ectotherms have a limited ability to evolutionarily compensate for the depressing effects of low temperatures on rates of biological processes. The remaining unexplained variation in development time likely reflects additional ecological and evolutionary differences among insect species.

CONSTRAINTS ON SPATIAL VARIATIONS IN THE FINE-STRUCTURE CONSTANT FROM PLANCK

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

O'Bryan, Jon; Smidt, Joseph; De Bernardis, Francesco

2015-01-01

We use the cosmic microwave background (CMB) anisotropy data from Planck to constrain the spatial fluctuations of the fine-structure constant α at a redshift of 1100. We use a quadratic estimator to measure the four-point correlation function of the CMB temperature anisotropies and extract the angular power spectrum fine-structure constant spatial variations projected along the line of sight at the last scattering surface. At tens of degree angular scales and above, we constrain the fractional rms fluctuations of the fine-structure constant to be (δα/α){sub rms} < 3.4 × 10{sup –3} at the 68% confidence level. We find no evidence formore » a spatially varying α at a redshift of 10{sup 3}.« less

A survey of volcano deformation in the central Andes using InSAR: Evidence for deep, slow inflation

NASA Astrophysics Data System (ADS)

Pritchard, M. E.; Simons, M.

2001-12-01

We use interferometric synthetic aperture radar (InSAR) to survey about 50 volcanos of the central Andes (15-27o S) for deformation during the 1992-2000 time interval. Because of the remote location of these volcanos, the activity of most are poorly constrained. Using the ERS-1/2 C-band radars (5.6 cm), we observe good interferometric correlation south of about 21o S, but poor correlation north of that latitude, especially in southern Peru. This variation is presumably related to regional climate variations. Our survey reveals broad (10's of km), roughly axisymmetric deformation at 2 volcanic centers with no previously documented deformation. At Uturuncu volcano, in southwestern Bolivia, the deformation rate can be constrained with radar data from several satellite tracks and is about 1 cm/year between 1992 and 2000. We find a second source of volcanic deformation located between Lastarria and Cordon del Azufre volcanos near the Chile/Argentina border. There is less radar data to constrain the deformation in this area, but the rate is also about 1 cm/yr between 1996 and 2000. While the spatial character of the deformation field appears to be affected by atmosphere at both locations, we do not think that the entire signal is atmospheric, because the signal is observed in several interferograms and nearby edifices do not show similar patterns. The deformation signal appears to be time-variable, although it is difficult to determine whether this is due to real variations in the deformation source or atmospheric effects. We model the deformation with both a uniform point-source source of inflation, and a tri-axial point-source ellipsoid, and compare both elastic half-space and layered-space models. We also explore the effects of local topography upon the deformation field using the method of Williams and Wadge (1998). We invert for source parameters using the global search Neighborhood Algorithm of Sambridge (1998). Preliminary results indicate that the sources at both Uturuncu and Lastarria/Cordon del Azufre volcanos are model-dependent, but are generally greater than 10 km deep. This depth suggests a potential relationship between the deformation source at Uturuncu and the large Altiplano-Puna Magmatic Complex that has been imaged seismically (e.g. Chmielowski et al., 1999), although the deformation at Lastarria/Cordon del Azufre lies outside the region of lowest seismic velocities (Yuan et al., 2000).

Regional variation in acute care length of stay after orthopaedic surgery total joint replacement surgery and hip fracture surgery.

PubMed

Fitzgerald, John D; Weng, Haoling H; Soohoo, Nelson F; Ettner, Susan L

2013-01-01

To examine change in regional variations in acute care length of stay (LOS) after orthopedic surgery following expiration of the New York (NY) State exemption to the Prospective Payment System and implementation of the Medicare Short Stay Transfer Policy. Time series analyses were conducted to evaluate change in LOS across regions after policy implementations. Small area analyses were conducted to examine residual variation in LOS. The dataset included A 100% sample of fee-for-service Medicare patients undergoing surgical repair for hip fracture or elective joint replacement surgery between 1996 and 2001. Data files from Centers for Medicare and Medicaid Services 1996-2001 Medicare Provider Analysis and Review file, 1999 Provider of Service file, and data from the 2000 United States Census were used for analysis. In 1996, LOS in NY after orthopedic procedures was much longer than the remainder of the country. After policy changes, LOS fell. However, significant residual variation in LOS persisted. This residual variation was likely partly explained by differences variation in regional managed care market penetration, patient management practices and unmeasured characteristics associated with the hospital location. NY hospitals responded to changes in reimbursement policy, reducing variation in LOS. However, even after 5 years of financial pressure to constrain costs, other factors still have a strong impact on delivery of patient care.

Biogeography of time partitioning in mammals.

PubMed

Bennie, Jonathan J; Duffy, James P; Inger, Richard; Gaston, Kevin J

2014-09-23

Many animals regulate their activity over a 24-h sleep-wake cycle, concentrating their peak periods of activity to coincide with the hours of daylight, darkness, or twilight, or using different periods of light and darkness in more complex ways. These behavioral differences, which are in themselves functional traits, are associated with suites of physiological and morphological adaptations with implications for the ecological roles of species. The biogeography of diel time partitioning is, however, poorly understood. Here, we document basic biogeographic patterns of time partitioning by mammals and ecologically relevant large-scale patterns of natural variation in "illuminated activity time" constrained by temperature, and we determine how well the first of these are predicted by the second. Although the majority of mammals are nocturnal, the distributions of diurnal and crepuscular species richness are strongly associated with the availability of biologically useful daylight and twilight, respectively. Cathemerality is associated with relatively long hours of daylight and twilight in the northern Holarctic region, whereas the proportion of nocturnal species is highest in arid regions and lowest at extreme high altitudes. Although thermal constraints on activity have been identified as key to the distributions of organisms, constraints due to functional adaptation to the light environment are less well studied. Global patterns in diversity are constrained by the availability of the temporal niche; disruption of these constraints by the spread of artificial lighting and anthropogenic climate change, and the potential effects on time partitioning, are likely to be critical influences on species' future distributions.

A Test of Carbon and Oxygen Stable Isotope Ratio Process Models in Tree Rings.

NASA Astrophysics Data System (ADS)

Roden, J. S.; Farquhar, G. D.

2008-12-01

Stable isotopes ratios of carbon and oxygen in tree ring cellulose have been used to infer environmental change. Process-based models have been developed to clarify the potential of historic tree ring records for meaningful paleoclimatic reconstructions. However, isotopic variation can be influenced by multiple environmental factors making simplistic interpretations problematic. Recently, the dual isotope approach, where the variation in one stable isotope ratio (e.g. oxygen) is used to constrain the interpretation of variation in another (e.g. carbon), has been shown to have the potential to de-convolute isotopic analysis. However, this approach requires further testing to determine its applicability for paleo-reconstructions using tree-ring time series. We present a study where the information needed to parameterize mechanistic models for both carbon and oxygen stable isotope ratios were collected in controlled environment chambers for two species (Pinus radiata and Eucalyptus globulus). The seedlings were exposed to treatments designed to modify leaf temperature, transpiration rates, stomatal conductance and photosynthetic capacity. Both species were grown for over 100 days under two humidity regimes that differed by 20%. Stomatal conductance was significantly different between species and for seedlings under drought conditions but not between other treatments or humidity regimes. The treatments produced large differences in transpiration rate and photosynthesis. Treatments that effected photosynthetic rates but not stomatal conductance influenced carbon isotope discrimination more than those that influenced primarily conductance. The various treatments produced a range in oxygen isotope ratios of 7 ‰. Process models predicted greater oxygen isotope enrichment in tree ring cellulose than observed. The oxygen isotope ratios of bulk leaf water were reasonably well predicted by current steady-state models. However, the fractional difference between models that predict bulk leaf water versus the site of evaporation did not increase with transpiration rates. In conclusion, although the dual isotope approach may better constrain interpretation of isotopic variation, more work is required before its predictive power can be applied to tree-ring archives.

Radio to gamma-ray variability study of blazar S5 0716+714

DOE PAGES

Rani, B.; Krichbaum, T. P.; Fuhrmann, L.; ...

2013-03-13

In this paper, we present the results of a series of radio, optical, X-ray, and γ-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multifrequency observations were obtained using several ground- and space-based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend on a time scale of ~350 days. Episodes of fast variability recur on time scales of ~60-70 days. The intense and simultaneous activity at optical and γ-ray frequencies favors the synchrotron self-Compton mechanism for the production of the high-energy emission. Twomore » major low-peaking radio flares were observed during this high optical/γ-ray activity period. The radio flares are characterized by a rising and a decaying stage and agrees with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield robust and self-consistent lower limits of δ ≥ 20 and equipartition magnetic field B eq ≥ 0.36 G. Causality arguments constrain the size of emission region θ ≤ 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and γ-rays. Theoptical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Finally, here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.« less

Body Size Adaptations to Altitudinal Climatic Variation in Neotropical Grasshoppers of the Genus Sphenarium (Orthoptera: Pyrgomorphidae)

PubMed Central

2015-01-01

Altitudinal clines in body size can result from the effects of natural and sexual selection on growth rates and developing times in seasonal environments. Short growing and reproductive seasons constrain the body size that adults can attain and their reproductive success. Little is known about the effects of altitudinal climatic variation on the diversification of Neotropical insects. In central Mexico, in addition to altitude, highly heterogeneous topography generates diverse climates that can occur even at the same latitude. Altitudinal variation and heterogeneous topography open an opportunity to test the relative impact of climatic variation on body size adaptations. In this study, we investigated the relationship between altitudinal climatic variation and body size, and the divergence rates of sexual size dimorphism (SSD) in Neotropical grasshoppers of the genus Sphenarium using a phylogenetic comparative approach. In order to distinguish the relative impact of natural and sexual selection on the diversification of the group, we also tracked the altitudinal distribution of the species and trends of both body size and SSD on the phylogeny of Sphenarium. The correlative evidence suggests no relationship between altitude and body size. However, larger species were associated with places having a warmer winter season in which the temporal window for development and reproduction can be longer. Nonetheless, the largest species were also associated with highly seasonal environments. Moreover, large body size and high levels of SSD have evolved independently several times throughout the history of the group and male body size has experienced a greater evolutionary divergence than females. These lines of evidence suggest that natural selection, associated with seasonality and sexual selection, on maturation time and body size could have enhanced the diversification of this insect group. PMID:26684616

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

Saide, Pablo E.; Peterson, David A.; de Silva, Arlindo

We couple airborne, ground-based, and satellite observations; conduct regional simulations; and develop and apply an inversion technique to constrain hourly smoke emissions from the Rim Fire, the third largest observed in California, USA. Emissions constrained with multiplatform data show notable nocturnal enhancements (sometimes over a factor of 20), correlate better with daily burned area data, and are a factor of 2–4 higher than a priori estimates, highlighting the need for improved characterization of diurnal profiles and day-to-day variability when modeling extreme fires. Constraining only with satellite data results in smaller enhancements mainly due to missing retrievals near the emissions source,more » suggesting that top-down emission estimates for these events could be underestimated and a multiplatform approach is required to resolve them. Predictions driven by emissions constrained with multiplatform data present significant variations in downwind air quality and in aerosol feedback on meteorology, emphasizing the need for improved emissions estimates during exceptional events.« less

Degree-constrained multicast routing for multimedia communications

NASA Astrophysics Data System (ADS)

Wang, Yanlin; Sun, Yugeng; Li, Guidan

2005-02-01

Multicast services have been increasingly used by many multimedia applications. As one of the key techniques to support multimedia applications, the rational and effective multicast routing algorithms are very important to networks performance. When switch nodes in networks have different multicast capability, multicast routing problem is modeled as the degree-constrained Steiner problem. We presented two heuristic algorithms, named BMSTA and BSPTA, for the degree-constrained case in multimedia communications. Both algorithms are used to generate degree-constrained multicast trees with bandwidth and end to end delay bound. Simulations over random networks were carried out to compare the performance of the two proposed algorithms. Experimental results show that the proposed algorithms have advantages in traffic load balancing, which can avoid link blocking and enhance networks performance efficiently. BMSTA has better ability in finding unsaturated links and (or) unsaturated nodes to generate multicast trees than BSPTA. The performance of BMSTA is affected by the variation of degree constraints.

What Do District Health Managers in Ghana Use Their Working Time for? A Case Study of Three Districts.

PubMed

Bonenberger, Marc; Aikins, Moses; Akweongo, Patricia; Bosch-Capblanch, Xavier; Wyss, Kaspar

2015-01-01

Ineffective district health management potentially impacts on health system performance and service delivery. However, little is known about district health managing practices and time allocation in resource-constrained health systems. Therefore, a time use study was conducted in order to understand current time use practices of district health managers in Ghana. All 21 district health managers working in three districts of the Eastern Region were included in the study and followed for a period of three months. Daily retrospective interviews about their time use were conducted, covering 1182 person-days of observation. Total time use of the sample population was assessed as well as time use stratified by managerial position. Differences of time use over time were also evaluated. District health managers used most of their working time for data management (16.6%), attending workshops (12.3%), financial management (8.7%), training of staff (7.1%), drug and supply management (5.0%), and travelling (9.6%). The study found significant variations of time use across the managerial cadres as well as high weekly variations of time use impulsed mainly by a national vertical program. District health managers in Ghana use substantial amounts of their working time in only few activities and vertical programs greatly influence their time use. Our findings suggest that efficiency gains are possible for district health managers. However, these are unlikely to be achieved without improvements within the general health system, as inefficiencies seem to be largely caused by external factors.

Modeling Wide-Angle Seismic Data from the Hi-CLIMB Experiment in Tibet

NASA Astrophysics Data System (ADS)

Nowack, R. L.; Griffin, J. D.; Tseng, T.; Chen, W.

2009-12-01

Using data from local and regional events recorded by the Hi-CLIMB array in Tibet, we utilize seismic attributes, including arrival times, Hilbert amplitudes and pulse frequencies, to constrain structures of seismic wave speed and attenuation in the crust and the upper mantle in western China. We construct more than 30 high-quality, regional seismic profiles, and select 14 of these, which show excellent crustal and Pn arrivals, for further analysis. Travel-times from events at regional distances constrain large-scale velocity structures, and four close-in events provide further details on crustal structure. We use the 3-D ray tracer, CRT, to model the travel-times. Initial results indicate that the Moho beneath the Lhasa terrane of southern Tibet is over 73 km deep with a high Pn speed of about 8.2 km/s. In contrast, the Qiangtang terrane farther north shows a thinner crust, by up to 10 km, and a low Pn speed of 7.8-7.9 km/s. Preliminary estimates of upper mantle velocity gradients are between .003 and .004 km/s per km, consistent with previous results by Phillips et al. (2007). We also use P to SV conversions from teleseismic earthquakes to independently constrain variations in speeds of Pn and depths of the Moho. For instance, amplitudes of the SsPmP phase, when its last reflection off the Moho is near-critical, are particularly sensitive to the contrast in seismic wave speeds across the crust-mantle interface; and results from these additional data are consistent with those from modeling of travel-times. Additional seismic attributes, extracted from wave-trains containing Pn and major crustal phases, are being compared with results of numerical modeling based on the spectral element method and asymptotic calculations in laterally varying media, where both lateral and vertical gradients in seismic wave speeds can strongly affect Pn amplitudes and pulse frequencies.

Hyperspherical lowest-order constrained-variational approximation to resonant Bose-Einstein condensates

NASA Astrophysics Data System (ADS)

Sze, M. W. C.; Sykes, A. G.; Blume, D.; Bohn, J. L.

2018-03-01

We study the ground-state properties of a system of N harmonically trapped bosons of mass m interacting with two-body contact interactions, from small to large scattering lengths. This is accomplished in a hyperspherical coordinate system that is flexible enough to describe both the overall scale of the gas and two-body correlations. By adapting the lowest-order constrained-variational method, we are able to semiquantitatively attain Bose-Einstein condensate ground-state energies even for gases with infinite scattering length. In the large-particle-number limit, our method provides analytical estimates for the energy per particle E0/N ≈2.5 N1 /3ℏ ω and two-body contact C2/N ≈16 N1 /6√{m ω /ℏ } for a Bose gas on resonance, where ω is the trap frequency.

The exit-time problem for a Markov jump process

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

Burch, N.; D'Elia, Marta; Lehoucq, Richard B.

2014-12-15

The purpose of our paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developedmore » nonlocal vector calculus. Furthermore, this calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.« less

A Long-term study of three rotating radio transients

NASA Astrophysics Data System (ADS)

Bhattacharyya, B.; Lyne, A. G.; Stappers, B. W.; Weltevrede, P.; Keane, E. F.; McLaughlin, M. A.; Kramer, M.; Jordan, C.; Bassa, C.

2018-04-01

We present the longest-term timing study so far of three Rotating Radio Transients (RRATs) - J1819-1458, J1840-1419 and J1913+1330 - performed using the Lovell, Parkes and Green Bank telescopes over the past decade. We study long-term and short-term variations of the pulse emission rate from these RRATs and report a marginal indication of a long-term increase in pulse detection rate over time for PSR J1819-1458 and J1913+1330. For PSR J1913+1330, we also observe a two orders of magnitude variation in the observed pulse detection rates across individual epochs, which may constrain the models explaining the origin of RRAT pulses. PSR J1913+1330 is also observed to exhibit a weak persistent emission mode. We investigate the post-glitch timing properties of J1819-1458 (the only RRAT for which glitches are observed) and discuss the implications for possible glitch models. Its post-glitch over-recovery of the frequency derivative is magnetar-like and similar behaviour is only observed for two other pulsars, both of which have relatively high magnetic field strengths. Following the over-recovery we also observe that some fraction of the pre-glitch frequency derivative is gradually recovered.

A multi-centennial time series of well-constrained ΔR values for the Irish Sea derived using absolutely-dated shell samples from the mollusc Arctica islandica

NASA Astrophysics Data System (ADS)

Butler, P. G.; Scourse, J. D.; Richardson, C. A.; Wanamaker, A. D., Jr.

2009-04-01

Determinations of the local correction (ΔR) to the globally averaged marine radiocarbon reservoir age are often isolated in space and time, derived from heterogeneous sources and constrained by significant uncertainties. Although time series of ΔR at single sites can be obtained from sediment cores, these are subject to multiple uncertainties related to sedimentation rates, bioturbation and interspecific variations in the source of radiocarbon in the analysed samples. Coral records provide better resolution, but these are available only for tropical locations. It is shown here that it is possible to use the shell of the long-lived bivalve mollusc Arctica islandica as a source of high resolution time series of absolutely-dated marine radiocarbon determinations for the shelf seas surrounding the North Atlantic ocean. Annual growth increments in the shell can be crossdated and chronologies can be constructed in a precise analogue with the use of tree-rings. Because the calendar dates of the samples are known, ΔR can be determined with high precision and accuracy and because all the samples are from the same species, the time series of ΔR values possesses a high degree of internal consistency. Presented here is a multi-centennial (AD 1593 - AD 1933) time series of 31 ΔR values for a site in the Irish Sea close to the Isle of Man. The mean value of ΔR (-62 14C yrs) does not change significantly during this period but increased variability is apparent before AD 1750.

Digging into the corona: A modeling framework trained with Sun-grazing comet observations

NASA Astrophysics Data System (ADS)

Jia, Y. D.; Pesnell, W. D.; Bryans, P.; Downs, C.; Liu, W.; Schwartz, S. J.

2017-12-01

Images of comets diving into the low corona have been captured a few times in the past decade. Structures visible at various wavelengths during these encounters indicate a strong variation of the ambient conditions of the corona. We combine three numerical models: a global coronal model, a particle transportation model, and a cometary plasma interaction model into one framework to model the interaction of such Sun-grazing comets with plasma in the low corona. In our framework, cometary vapors are ionized via multiple channels and then captured by the coronal magnetic field. In seconds, these ions are further ionized into their highest charge state, which is revealed by certain coronal emission lines. Constrained by observations, we apply our framework to trace back to the local conditions of the ambient corona, and their spatial/time variation over a broad range of scales. Once trained by multiple stages of the comet's journey in the low corona, we illustrate how this framework can leverage these unique observations to probe the structure of the solar corona and solar wind.

Variations in Temperature at the Base of the Lithosphere Beneath the Archean Superior Province, Canada

NASA Astrophysics Data System (ADS)

Mareschal, J.; Jaupart, C. P.

2013-12-01

Most of the variations in surface heat flux in stable continents are caused by variations in crustal heat production, with an almost uniform heat flux at the base of the crust ( 15+/-3 mW/m2). Such relatively small differences in Moho heat flux cannot be resolved by heat flow data alone, but they lead to important lateral variations in lithospheric temperatures and thicknesses. In order to better constrain temperatures in the lower lithosphere, we have combined surface heat flow and heat production data from the southern Superior Province in Canada with vertical shear wave velocity profiles obtained from surface wave inversion. We use the Monte-Carlo method to generate lithospheric temperature profiles from which shear wave velocity can be calculated for a given mantle composition. We eliminate thermal models which yield lithospheric and sub-lithospheric velocities that do not fit the shear wave velocity profile. Surface heat flux being constrained, the free parameters of the thermal model are: the mantle heat flux, the mantle heat production, the crustal differentiation index (ratio of surface to bulk crustal heat production) and the temperature of the mantle isentrope. Two conclusions emerge from this study. One is that, for some profiles, the vertical variations in shear wave velocities cannot be accounted for by temperature alone but also require compositional changes within the lithosphere. The second is that there are long wavelength horizontal variations in mantle temperatures (~80-100K) at the base of the lithosphere and in the mantle below

Relationships between contourite deposition, climate and slope instability: new insights from the Demerara Plateau (French Guyana)

NASA Astrophysics Data System (ADS)

Tallobre, C.; Bassetti, M. A.; Loncke, L.; Giresse, P.; Bayon, G.; Buscail, R.

2015-12-01

A Contourite Depositional System (CDS) has been described at the Demerara Plateau (DP) based on seismic investigations, but little is known about the mechanisms of associated sediment deposition and its interaction with past deep ocean circulation patterns (e.g. bottom current velocity) and bottom morphology related to ancient event of slope instability. The new seismic and bathymetric data recently acquired allow describing in details the CDS on the DP. Erosional and syn-sedimentary features on the seafloor (comet tail, « longitudinal waves », contourite drifts and moats) have been observed, helping to constrain the sedimentary processes at the origin of the CDS construction. Also, the recovery and multi-proxy analysis of sediment cores allows the characterization of sedimentary environments and possible relation with climate forcing. These sediment cores are characterized by the presence of several beds rich in glauconite grains. Glauconite can form at the sediment/water interface by winnowing effect that prevent sediment deposition and increase the residence time at the seafloor. Under strong winnowing conditions, glauconite grains can develop at several stages of maturity. We observed that the residence time and hence the maturity of glauconite is reflected by the color changes (light to dark green), the presence of crack on grains, the formation of (secondary) glauconite lamellae and decrease of grain porosity. A chronological framework (based on radiocarbon dates and δ18O variations) of contourite sequences at the studied location indicates correlation with grain-size parameters (sortable silt) and allows one to further constrain their dynamics through time. The combination of these proxies allows us to estimate and understand the evolution and the impact of the bottom current on sedimentation on the DP during the last 80 ky. These results show the potentiality of the glauconite study to estimate the relative variation of bottom current velocity at margins.

Io's Heat Flow: A Model Including "Warm" Polar Regions

NASA Astrophysics Data System (ADS)

Veeder, G. J.; Matson, D. L.; Johnson, T. V.; Davies, A. G.; Blaney, D. L.

2002-12-01

Some 90 percent of Io's surface is thermally "passive" material. It is separate from the sites of active volcanic eruptions. Though "passive", its thermal behavior continues to be a challenge for modelers. The usual approach is to take albedo, average daytime temperature, temperature as a function of time of day, etc., and attempt to match these constraints with a uniform surface with a single value of thermal inertia. Io is a case where even globally averaged observations are inconsistent with a single-thermal-inertia model approach. The Veeder et al. (1994) model for "passive" thermal emission addressed seven constraints derived from a decade of ground-based, global observations - average albedo plus infrared fluxes at three separate wavelengths (4.8, 8.7, and 20 microns) for both daytime and eclipsed conditions. This model has only two components - a unit of infinite thermal inertia and a unit of zero thermal inertia. The free parameters are the areal coverage ratio of the two units and their relative albedos (constrained to match the known average albedo). This two-parameter model agreed with the global radiometric data and also predicted significantly higher non-volcanic nighttime temperatures than traditional ("lunar-like") single thermal inertia models. Recent observations from the Galileo infrared radiometer show relatively uniform minimum-night-time temperatures. In particular, they show little variation with either latitude or time of night (Spencer et al., 2000; Rathbun et al., 2002). Additionally, detailed analyses of Io's scattering properties and reflectance variations have led to the interesting conclusion that Io's albedo at regional scales varies little with latitude (Simonelli, et al., 2001). This effectively adds four new observational constraints - lack of albedo variation with latitude, average minimum nighttime temperature and lack of variation of temperature with either latitude or longitude. We have made the fewest modifications necessary for the Veeder et al. model to match these new constrains - we added two model parameters to characterize the volcanically heated high-latitude units. These are the latitude above which the unit exists and its nighttime temperature. The resulting four-parameter model is the first that encompasses all of the available observations of Io's thermal emission and that quantitatively satisfies all eleven observational constraints. While no model is unique, this model is significant because it is the first to accommodate widespread polar regions that are relatively "warm". This work was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under contract to NASA.

Seismic velocity and crustal thickness inversions: Moon and Mars

NASA Astrophysics Data System (ADS)

Drilleau, Melanie; Blanchette-Guertin, Jean-François; Kawamura, Taichi; Lognonné, Philippe; Wieczorek, Mark

2017-04-01

We present results from new inversions of seismic data arrival times acquired by the Apollo active and passive experiments. Markov chain Monte Carlo inversions are used to constrain (i) 1-D lunar crustal and upper mantle velocity models and (ii) 3-D lateral crustal thickness models under the Apollo stations and the artificial and natural impact sites. A full 3-D model of the lunar crustal thickness is then obtained using the GRAIL gravimetric data, anchored by the crustal thicknesses under each Apollo station and impact site. To avoid the use of any seismic reference model, a Bayesian inversion technique is implemented. The advantage of such an approach is to obtain robust probability density functions of interior structure parameters governed by uncertainties on the seismic data arrival times. 1-D seismic velocities are parameterized using C1-Bézier curves, which allow the exploration of both smoothly varying models and first-order discontinuities. The parameters of the inversion include the seismic velocities of P and S waves as a function of depth, the thickness of the crust under each Apollo station and impact epicentre. The forward problem consists in a ray tracing method enabling both the relocation of the natural impact epicenters, and the computation of time corrections associated to the surface topography and the crustal thickness variations under the stations and impact sites. The results show geology-related differences between the different sites, which are due to contrasts in megaregolith thickness and to shallow subsurface composition and structure. Some of the finer structural elements might be difficult to constrain and might fall within the uncertainties of the dataset. However, we use the more precise LROC-located epicentral locations for the lunar modules and Saturn-IV upper stage artificial impacts, reducing some of the uncertainties observed in past studies. In the framework of the NASA InSight/SEIS mission to Mars, the method developed in this study will be used to constrain the Martian crustal thickness as soon as the first data will be available (late 2018). For Insight, impacts will be located by MRO data differential analysis, which provide a known location enabling the direct inversion of all differential travel times with respect to P arrival time. We have performed resolution tests to investigate to what extend impact events might help us to constrain the Martian crustal thickness. Due to the high flexibility of the Bayesian algorithm, the interior model will be refined each time a new event will be detected.

Wide variation in hospital and physician payment rates evidence of provider market power.

PubMed

Ginsburg, Paul B

2010-11-01

Wide variation in private insurer payment rates to hospitals and physicians across and within local markets suggests that some providers, particularly hospitals, have significant market power to negotiate higher-than-competitive prices, according to a new study by the Center for Studying Health System Change (HSC). Looking across eight health care markets--Cleveland; Indianapolis; Los Angeles; Miami; Milwaukee; Richmond, Va.; San Francisco; and rural Wisconsin--average inpatient hospital payment rates of four large national insurers ranged from 147 percent of Medicare in Miami to 210 percent in San Francisco. In extreme cases, some hospitals command almost five times what Medicare pays for inpatient services and more than seven times what Medicare pays for outpatient care. Variation within markets was just as dramatic. For example, the hospital with prices at the 25th percentile of Los Angeles hospitals received 84 percent of Medicare rates for inpatient care, while the hospital with prices at the 75th percentile received 184 percent of Medicare rates. The highest-priced Los Angeles hospital with substantial inpatient claims volume received 418 percent of Medicare. While not as pronounced, significant variation in physician payment rates also exists across and within markets and by specialty. Few would characterize the variation in hospital and physician payment rates found in this study to be consistent with a highly competitive market. Purchasers and public policy makers can address provider market power, or the ability to negotiate higher-than-competitive prices, through two distinct approaches. One is to pursue market approaches to strengthen competitive forces, while the other is to constrain payment rates through regulation.

Auto-biometric for M-mode echocardiography

NASA Astrophysics Data System (ADS)

Zhang, Wei; Park, Jinhyong; Zhou, S. Kevin

2010-03-01

In this paper we present a system for fast and accurate detection of anatomical structures (calipers) in M-mode images. The task is challenging because of dramatic variations in their appearances. We propose to solve the problem in a progressive manner, which ensures both robustness and efficiency. It first obtains rough caliper localization using the intensity profile image. Then run a constrained search for accurate caliper positions. Markov Random Field (MRF) and warping image detectors are used for jointly considering appearance information and the geometric relationship between calipers. Extensive experiments show that our system achieves more accurate results and uses less time in comparison with previously reported work.

On the Interaction Between Gravity Waves and Atmospheric Thermal Tides

NASA Astrophysics Data System (ADS)

Agner, Ryan Matthew

Gravity waves and thermal tides are two of the most important dynamical features of the atmosphere. They are both generated in the lower atmosphere and propagate upward transporting energy and momentum to the upper atmosphere. This dissertation focuses on the interaction of these waves in the Mesosphere and Lower Thermosphere (MLT) region of the atmosphere using both observational data and Global Circulation Model (GCMs). The first part of this work focuses on observations of gravity wave interactions with the tides using both LIDAR data at the Star Fire Optical Range (SOR, 35?N, 106.5?W) and a meteor radar data at the Andes LIDAR Observatory (ALO, 30.3?S, 70.7?W). At SOR, the gravity waves are shown to enhance or damp the amplitude of the diurnal variations dependent on altitude while the phase is always delayed. The results compare well with previous mechanistic model results and with the Japanese Atmospheric General circulation model for Upper Atmosphere Research (JAGUAR) high resolution global circulation model. The meteor radar observed the GWs to almost always enhance the tidal amplitudes and either delay or advance the phase depending on the altitude. When compared to previous radar results from the same meteor radar when it was located in Maui, Hawaii, the Chile results are very similar while the LIDAR results show significant differences. This is because of several instrument biases when calculating GW momentum fluxes that is not significant when determining the winds. The radar needs to perform large amounts of all-sky averaging across many weeks, while the LIDAR directly detects waves in a small section of sky. The second part of this work focuses on gravity wave parameterization scheme effects on the tides in GCMs. The Specified Dynamics Whole Atmosphere Community Climate Model (SD-WACCM) and the extended Canadian Middle Atmosphere Model (eCMAM) are used for this analysis. The gravity wave parameterization schemes in the eCMAM (Hines scheme) have been shown to enhance the tidal amplitudes compared to observations while the parameterization scheme in SD-WACCM (Lindzen scheme) overdamps the tides. It is shown here that the Hines scheme assumption that only small scale gravity waves force the atmosphere do not create enough drag to properly constrain the tidal amplitudes. The Lindzen scheme produces too much drag because all wave scales are assumed to be saturated thus continuing to provide forcing on the atmosphere above the breaking altitude. The final part of this work investigates GWs, tides and their interactions on a local time scale instead of a global scale in the two GCMs. The local time GWs in eCMAM are found to have a strong seasonal dependence, with the majority of the forcings at the winter pole at latitudes where the diurnal variations are weak limiting their interactions. In SD-WACCM, the largest local GW forcings are located at mid latitudes near where the diurnal variations peak causing them to dampen the diurnal amplitudes. On a local time level the diurnal variations may be a summation of many tidal modes. The analysis reveals that in eCMAM the DW1 tidal mode is by far the dominant mode accounting for the local time variations. The high amount of modulation of GWs by the DW1 tidal winds does not allow it to be properly constrained, causing it to dominate the local time diurnal variations. Similarly, the DW1 projection of GW forcing is dominant over all other other modes and contributes the most to the local time diurnal GW variations. The local time wind variations in SD-WACCM are in uenced by several tidal modes because the DW1 tide is of compatible amplitudes to other modes. This is because of the increased damping on the tide by the GWs. It is also found that the local GW diurnal variations have significant contributions from all tidal modes due to the time and location of the forcing being dependent only on the tropospheric source regions and not the at altitude tidal winds.

Ground-state densities from the Rayleigh-Ritz variation principle and from density-functional theory.

PubMed

Kvaal, Simen; Helgaker, Trygve

2015-11-14

The relationship between the densities of ground-state wave functions (i.e., the minimizers of the Rayleigh-Ritz variation principle) and the ground-state densities in density-functional theory (i.e., the minimizers of the Hohenberg-Kohn variation principle) is studied within the framework of convex conjugation, in a generic setting covering molecular systems, solid-state systems, and more. Having introduced admissible density functionals as functionals that produce the exact ground-state energy for a given external potential by minimizing over densities in the Hohenberg-Kohn variation principle, necessary and sufficient conditions on such functionals are established to ensure that the Rayleigh-Ritz ground-state densities and the Hohenberg-Kohn ground-state densities are identical. We apply the results to molecular systems in the Born-Oppenheimer approximation. For any given potential v ∈ L(3/2)(ℝ(3)) + L(∞)(ℝ(3)), we establish a one-to-one correspondence between the mixed ground-state densities of the Rayleigh-Ritz variation principle and the mixed ground-state densities of the Hohenberg-Kohn variation principle when the Lieb density-matrix constrained-search universal density functional is taken as the admissible functional. A similar one-to-one correspondence is established between the pure ground-state densities of the Rayleigh-Ritz variation principle and the pure ground-state densities obtained using the Hohenberg-Kohn variation principle with the Levy-Lieb pure-state constrained-search functional. In other words, all physical ground-state densities (pure or mixed) are recovered with these functionals and no false densities (i.e., minimizing densities that are not physical) exist. The importance of topology (i.e., choice of Banach space of densities and potentials) is emphasized and illustrated. The relevance of these results for current-density-functional theory is examined.

Monitoring temporal opacity fluctuations of large structures with muon radiography: a calibration experiment using a water tower

PubMed Central

Jourde, Kevin; Gibert, Dominique; Marteau, Jacques; de Bremond d’Ars, Jean; Gardien, Serge; Girerd, Claude; Ianigro, Jean-Christophe

2016-01-01

Usage of secondary cosmic muons to image the geological structures density distribution significantly developed during the past ten years. Recent applications demonstrate the method interest to monitor magma ascent and volcanic gas movements inside volcanoes. Muon radiography could be used to monitor density variations in aquifers and the critical zone in the near surface. However, the time resolution achievable by muon radiography monitoring remains poorly studied. It is biased by fluctuation sources exterior to the target, and statistically affected by the limited number of particles detected during the experiment. The present study documents these two issues within a simple and well constrained experimental context: a water tower. We use the data to discuss the influence of atmospheric variability that perturbs the signal, and propose correction formulas to extract the muon flux variations related to the water level changes. Statistical developments establish the feasibility domain of muon radiography monitoring as a function of target thickness (i.e. opacity). Objects with a thickness comprised between ≈50 ± 30 m water equivalent correspond to the best time resolution. Thinner objects have a degraded time resolution that strongly depends on the zenith angle, whereas thicker objects (like volcanoes) time resolution does not. PMID:26971718

A SEMI-LAGRANGIAN TWO-LEVEL PRECONDITIONED NEWTON-KRYLOV SOLVER FOR CONSTRAINED DIFFEOMORPHIC IMAGE REGISTRATION.

PubMed

Mang, Andreas; Biros, George

2017-01-01

We propose an efficient numerical algorithm for the solution of diffeomorphic image registration problems. We use a variational formulation constrained by a partial differential equation (PDE), where the constraints are a scalar transport equation. We use a pseudospectral discretization in space and second-order accurate semi-Lagrangian time stepping scheme for the transport equations. We solve for a stationary velocity field using a preconditioned, globalized, matrix-free Newton-Krylov scheme. We propose and test a two-level Hessian preconditioner. We consider two strategies for inverting the preconditioner on the coarse grid: a nested preconditioned conjugate gradient method (exact solve) and a nested Chebyshev iterative method (inexact solve) with a fixed number of iterations. We test the performance of our solver in different synthetic and real-world two-dimensional application scenarios. We study grid convergence and computational efficiency of our new scheme. We compare the performance of our solver against our initial implementation that uses the same spatial discretization but a standard, explicit, second-order Runge-Kutta scheme for the numerical time integration of the transport equations and a single-level preconditioner. Our improved scheme delivers significant speedups over our original implementation. As a highlight, we observe a 20 × speedup for a two dimensional, real world multi-subject medical image registration problem.

Study of Plasma Waves Observed onboard Rosetta in the 67P/ChuryumovGerasimenko Comet Environment Using High Time Resolution Density Data Inferred from RPC-MIP and RPC-LAP Cross-calibration

NASA Astrophysics Data System (ADS)

Breuillard, H.; Henri, P.; Vallières, X.; Eriksson, A. I.; Odelstad, E.; Johansson, F. L.; Richter, I.; Goetz, C.; Wattieaux, G.; Tsurutani, B.; Hajra, R.; Le Contel, O.

2017-12-01

During two years, the groundbreaking ESA/Rosetta mission was able to escort comet 67P where previous cometary missions were only limited to flybys. This enabled for the first time to make in-situ measurements of the evolution of a comet's plasma environment. The density and temperature measured by Rosetta are derived from RPC-Mutual Impedance Probe (MIP) and RPC-Langmuir Probe (LAP). On one hand, low time resolution electron density are calculated using the plasma frequency extracted from the MIP mutual impedance spectra. On the other hand, high time resolution density fluctuations are estimated from the spacecraft potential measured by LAP. In this study, using a simple spacecraft charging model, we perform a cross-calibration of MIP plasma density and LAP spacecraft potential variations to obtain high time resolution measurements of the electron density. These results are also used to constrain the electron temperature. Then we make use of these new dataset, together with RPC-MAG magnetic field measurements, to investigate for the first time the compressibility and the correlations between plasma and magnetic field variations, for both singing comet waves and steepened waves observed, respectively during low and high cometary outgassing activity, in the plasma environment of comet 67P.

VLBI observations of the 0957 + 561 gravitational lens system

NASA Technical Reports Server (NTRS)

Gorenstein, M. V.; Falco, E. E.; Shapiro, I. I.; Bartel, N.; Bonometti, R. J.; Cohen, N. L.; Rogers, A. E. E.; Marcaide, J. M.; Clark, T. A.

1988-01-01

A series of VLBI observations of the gravitational lens system 0957 + 561 at a wavelength of 13 cm has yielded the positions of the A and B images, the relative magnification of their largest discernible radio structures, and the time variability of their smallest discernible radio structures. These observations have also allowed upper limits to be placed on the flux density of an expected third image. The positions and relative magnification of the A and B images provide new information with which to constrain models of the lens that forms the images. The detection of variations in the flux densities of the cores of A and B suggests that observations at shorter wavelengths may reveal superluminal motion, which may in turn provide a means to measure the relative time delay.

New Late Neolithic (c. 7000-5000 BC) archeointensity data from Syria. Reconstructing 9000 years of archeomagnetic field intensity variations in the Middle East

NASA Astrophysics Data System (ADS)

Gallet, Yves; Molist Montaña, Miquel; Genevey, Agnès; Clop García, Xavier; Thébault, Erwan; Gómez Bach, Anna; Le Goff, Maxime; Robert, Béatrice; Nachasova, Inga

2015-01-01

We present new archeomagnetic intensity data from two Late Neolithic archeological sites (Tell Halula and Tell Masaïkh) in Syria. These data, from 24 groups of potsherds encompassing 15 different time levels, are obtained using the Triaxe experimental protocol, which takes into account both the thermoremanent magnetization anisotropy and cooling rate effects on intensity determinations. They allow us to recover the geomagnetic intensity variations in the Middle East, between ∼7000 BC and ∼5000 BC, i.e. during the so-called pre-Halaf, proto-Halaf, Halaf and Halaf-Ubaid Transitional cultural phases. The data are compared with previous archeointensity results of similar ages from Northern Iraq (Yarim Tepe II and Tell Sotto) and Bulgaria. We find that previous dating of the Iraqi material was in error. When corrected, all northern Mesopotamian data show a relatively good consistency and also reasonably match with the Bulgarian archeointensity dataset. Using a compilation of available data, we construct a geomagnetic field intensity variation curve for the Middle East encompassing the past 9000 years, which makes it presently the longest known regional archeomagnetic intensity record. We further use this compilation to constrain variations in dipole field moment over most of the Holocene. In particular, we discuss the possibility that a significant dipole moment maximum occurred during the third millennium BC, which cannot easily be identified in available time-varying global geomagnetic field reconstructions.

Constraining response output on conjunctive fixed-ratio 1 fixed-time reinforcement schedules: Effects on the postreinforcement pause.

PubMed

Lopez, F; Pereira, C

1985-03-01

Two experiments used response-restriction procedures in order to test the independence of the factors determining response rate and the factors determining the size of the postreinforcement pause on interval schedules. Responding was restricted by response-produced blackout or by retracting the lever. In Experiment 1 with a Conjunctive FR 1 FT schedule, the blackout procedure reduced the postreinforcement pause more than the lever-retraction procedure did, and both procedures produced shorter pauses than did the schedule without response restriction. In Experiment 2 the interreinforcement interval was also manipulated, and the size of the pause was an increasing function of the interreinforcement interval, but the rate of increase was lower than that produced by fixed interval schedules of comparable interval durations. The assumption of functional independence of the postreinforcement pause and terminal rate in fixed interval schedules is questioned since data suggest that pause reductions resulted from constraining variation in response number compared to equivalent periodic schedules in which response number was allowed to vary. Copyright © 1985. Published by Elsevier B.V.

Optimal Control of Evolution Mixed Variational Inclusions

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

Alduncin, Gonzalo, E-mail: alduncin@geofisica.unam.mx

2013-12-15

Optimal control problems of primal and dual evolution mixed variational inclusions, in reflexive Banach spaces, are studied. The solvability analysis of the mixed state systems is established via duality principles. The optimality analysis is performed in terms of perturbation conjugate duality methods, and proximation penalty-duality algorithms to mixed optimality conditions are further presented. Applications to nonlinear diffusion constrained problems as well as quasistatic elastoviscoplastic bilateral contact problems exemplify the theory.

A model for seasonal changes in GPS positions and seismic wave speeds due to thermoelastic and hydrologic variations

USGS Publications Warehouse

Tsai, V.C.

2011-01-01

It is known that GPS time series contain a seasonal variation that is not due to tectonic motions, and it has recently been shown that crustal seismic velocities may also vary seasonally. In order to explain these changes, a number of hypotheses have been given, among which thermoelastic and hydrology-induced stresses and strains are leading candidates. Unfortunately, though, since a general framework does not exist for understanding such seasonal variations, it is currently not possible to quickly evaluate the plausibility of these hypotheses. To fill this gap in the literature, I generalize a two-dimensional thermoelastic strain model to provide an analytic solution for the displacements and wave speed changes due to either thermoelastic stresses or hydrologic loading, which consists of poroelastic stresses and purely elastic stresses. The thermoelastic model assumes a periodic surface temperature, and the hydrologic models similarly assume a periodic near-surface water load. Since all three models are two-dimensional and periodic, they are expected to only approximate any realistic scenario; but the models nonetheless provide a quantitative framework for estimating the effects of thermoelastic and hydrologic variations. Quantitative comparison between the models and observations is further complicated by the large uncertainty in some of the relevant parameters. Despite this uncertainty, though, I find that maximum realistic thermoelastic effects are unlikely to explain a large fraction of the observed annual variation in a typical GPS displacement time series or of the observed annual variations in seismic wave speeds in southern California. Hydrologic loading, on the other hand, may be able to explain a larger fraction of both the annual variations in displacements and seismic wave speeds. Neither model is likely to explain all of the seismic wave speed variations inferred from observations. However, more definitive conclusions cannot be made until the model parameters are better constrained. Copyright ?? 2011 by the American Geophysical Union.

Variations of System IV period of the sulfur ions in the Io torus for the volcanic event in 2015 observed by the HISAKI/EXCEED

NASA Astrophysics Data System (ADS)

Ryo, A.; Misawa, H.; Tsuchiya, F.; Kagitani, M.; Yoshioka, K.; Yoshikawa, I.; Kimura, T.; Murakami, G.; Yamazaki, A.

2017-12-01

Previous ground-based and probe observations of Io plasma torus (IPT) in various wavelengths have detected a periodic variation which is longer than System III period (9.925 h). It has been called System IV period. Although various ideas to explain the origin of System IV period are proposed , little progress has been made. A previous study reported that the System IV period was 10.07 h during the Cassini's Jupiter flyby in 2000. This period was shorter than the typical period of 10.21h (Brown 1995). The Io's volcanic event just before Cassini's flyby might shorten the System IV period (Steffl et al., 2006). However, the certainty of this inference is not clear from the Cassini observation. Here, we analyzed time variations in intensities of EUV emissions from IPT obtained by the HISAKI satellite to understand the mechanism responsible for the System IV period and the influence of Io's volcanic event on IPT. The data used in this study was obtained from Dec. 2014 to the middle of May 2015. During this period, the Io's volcanic event from Jan. to Mar. 2015 was reported from the observation of Iogenic sodium emission (Yoneda et al., 2015). To find variations of the System IV period, the temporal variation of System III longitude at peak EUV intensity was derived by fitting light curves of three ion species (S II 76.5 nm + 126 nm, S III 68 nm, and S IV 65.7 nm + 140.5 nm). The System IV periods of S II before and after the Io's volcanic event were 10.16±0.008 h and 10.03±0.006 h, respectively. On the other hands, the System IV period of S II was 9.95±0.003 h during the Io's volcanic event. This is the first observational evidence which shows that the System IV period has shortened during the Io's volcanic event, which is consistent with the result from Cassini observation. Origin of the System IV period has been discussed with sub-corotation of plasma in IPT. Since the sub-corotation occurs due to mass loading of newly picked-up ions into IPT, it is expected that the System IV period becomes long during Io's volcanic event. However, the result shows the opposite feature and will give important information to constrain the origin of the System IV period. To verify above consideration, we will derive time variation in plasma parameters in IPT from the HISAKI data with plasma diagnosis analysis. This analysis will enable us to constrain the origin of the System IV period.

Systematic variations of argon diffusion in feldspars and implications for thermochronometry

DOE PAGES

Cassata, William S.; Renne, Paul R.

2013-03-07

Coupled information about the time-dependent production and temperature-dependent diffusion of radiogenic argon in feldspars can be used to constrain the thermal evolution attending a host of Earth and planetary processes. To better assess the accuracy of thermal models, an understanding of the mechanisms and pathways by which argon diffuses in feldspars is desirable. Here we present step-heating Ar diffusion experiments conducted on feldspars with diverse compositions, structural states, and microstructural characteristics. The experiments reveal systematic variations in diffusive behavior that appear closely related to these variables, with apparent closure temperatures for 0.1–1 mm grains of ~200–400 °C (assuming a 10more » °C/Ma cooling rate). Given such variability, there is no broadly applicable set of diffusion parameters that can be utilized in feldspar thermal modeling; sample-specific data are required. Diffusion experiments conducted on oriented cleavage flakes do not reveal directionally-dependent diffusive anisotropy to within the resolution limits of our approach (approximately a factor of 2). Additional experiments aimed at constraining the physical significance of the diffusion domain are presented and indicate that unaltered feldspar crystals with or without coherent exsolution lamellae diffuse at the grain-scale, whereas feldspars containing hydrothermal alteration and/or incoherent sub-grain intergrowths do not. Arrhenius plots for argon diffusion in plagioclase and alkali feldspars appear to reflect a confluence of intrinsic diffusion kinetics and structural transitions that occur during incremental heating experiments. These structural transitions, along with sub-grain domain size variations, cause deviations from linearity (i.e., upward and downward curvature) on Arrhenius plots. An atomistic model for Arrhenius behavior is proposed that incorporates the variable lattice deformations of different feldspars in response to heating and compression. Furthermore, the resulting implications for accurately extrapolating laboratory-derived diffusion parameters to natural settings and over geologic time are discussed. We find that considerable inaccuracies may exist in published thermal histories obtained using multiple diffusion domain (MDD) models fit to Arrhenius plots for exsolved alkali feldspar, where the inferred Ar partial retention zones may be spuriously hot.« less

Soil thermal dynamics, snow cover, and frozen depth under five temperature treatments in an ombrotrophic bog: Constrained forecast with data assimilation: Forecast With Data Assimilation

DOE PAGES

Huang, Yuanyuan; Jiang, Jiang; Ma, Shuang; ...

2017-08-18

We report that accurate simulation of soil thermal dynamics is essential for realistic prediction of soil biogeochemical responses to climate change. To facilitate ecological forecasting at the Spruce and Peatland Responses Under Climatic and Environmental change site, we incorporated a soil temperature module into a Terrestrial ECOsystem (TECO) model by accounting for surface energy budget, snow dynamics, and heat transfer among soil layers and during freeze-thaw events. We conditioned TECO with detailed soil temperature and snow depth observations through data assimilation before the model was used for forecasting. The constrained model reproduced variations in observed temperature from different soil layers,more » the magnitude of snow depth, the timing of snowfall and snowmelt, and the range of frozen depth. The conditioned TECO forecasted probabilistic distributions of soil temperature dynamics in six soil layers, snow, and frozen depths under temperature treatments of +0.0, +2.25, +4.5, +6.75, and +9.0°C. Air warming caused stronger elevation in soil temperature during summer than winter due to winter snow and ice. And soil temperature increased more in shallow soil layers in summer in response to air warming. Whole ecosystem warming (peat + air warmings) generally reduced snow and frozen depths. The accuracy of forecasted snow and frozen depths relied on the precision of weather forcing. Uncertainty is smaller for forecasting soil temperature but large for snow and frozen depths. Lastly, timely and effective soil thermal forecast, constrained through data assimilation that combines process-based understanding and detailed observations, provides boundary conditions for better predictions of future biogeochemical cycles.« less

Soil thermal dynamics, snow cover, and frozen depth under five temperature treatments in an ombrotrophic bog: Constrained forecast with data assimilation: Forecast With Data Assimilation

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

Huang, Yuanyuan; Jiang, Jiang; Ma, Shuang

We report that accurate simulation of soil thermal dynamics is essential for realistic prediction of soil biogeochemical responses to climate change. To facilitate ecological forecasting at the Spruce and Peatland Responses Under Climatic and Environmental change site, we incorporated a soil temperature module into a Terrestrial ECOsystem (TECO) model by accounting for surface energy budget, snow dynamics, and heat transfer among soil layers and during freeze-thaw events. We conditioned TECO with detailed soil temperature and snow depth observations through data assimilation before the model was used for forecasting. The constrained model reproduced variations in observed temperature from different soil layers,more » the magnitude of snow depth, the timing of snowfall and snowmelt, and the range of frozen depth. The conditioned TECO forecasted probabilistic distributions of soil temperature dynamics in six soil layers, snow, and frozen depths under temperature treatments of +0.0, +2.25, +4.5, +6.75, and +9.0°C. Air warming caused stronger elevation in soil temperature during summer than winter due to winter snow and ice. And soil temperature increased more in shallow soil layers in summer in response to air warming. Whole ecosystem warming (peat + air warmings) generally reduced snow and frozen depths. The accuracy of forecasted snow and frozen depths relied on the precision of weather forcing. Uncertainty is smaller for forecasting soil temperature but large for snow and frozen depths. Lastly, timely and effective soil thermal forecast, constrained through data assimilation that combines process-based understanding and detailed observations, provides boundary conditions for better predictions of future biogeochemical cycles.« less

Orbital time scale records of Asian eolian dust from the Sea of Japan since the early Pliocene

NASA Astrophysics Data System (ADS)

Zhang, Wenfang; De Vleeschouwer, David; Shen, Ji; Zhang, Zeke; Zeng, Lin

2018-05-01

A high-resolution potassium content record of sediments from Integrated Ocean Drilling Program (IODP) Site U1422 drilled in the northern Sea of Japan was employed to yield an astronomical timescale for this ∼205 m thick sedimentary archive. The K content was quantified using the natural gamma radiation (NGR) data routinely measured on DV JOIDES Resolution. The U1422 K (wt.%) series shows an increasing trend with time, which parallels the growth of North Hemisphere ice sheets since ∼4 Ma, as revealed by the global benthic foraminifer's oxygen isotope stack (LR04). We propose that K content variations reflect changes in the relative contributions of Asian eolian dust and volcanic weathering products, in response to changes in global ice volume. Using the shipboard age model, constrained by palaeomagnetism and nanofossils datums, we tuned the U1422 K content variations to the LR04 stack. Our tuned age model for site U1422 goes back to 3.9 Ma. The studied record indicates that Asian dust generation is enhanced during periods of global cooling. The synchronous variations between the U1422 K record and the LR04 stack also suggest that global cooling played a dominant role in promoting the central Asian aridification since ∼3.9 Ma, while the uplift of Tibetan Plateau plays a secondary role, considering many disagreements still exists between the timing and amplitude of the Qinghai-Tibetan uplift and the evolution of central Asian aridity since the early Pliocene.

What Do District Health Managers in Ghana Use Their Working Time for? A Case Study of Three Districts

PubMed Central

Bonenberger, Marc; Aikins, Moses; Akweongo, Patricia; Bosch-Capblanch, Xavier; Wyss, Kaspar

2015-01-01

Background Ineffective district health management potentially impacts on health system performance and service delivery. However, little is known about district health managing practices and time allocation in resource-constrained health systems. Therefore, a time use study was conducted in order to understand current time use practices of district health managers in Ghana. Methods All 21 district health managers working in three districts of the Eastern Region were included in the study and followed for a period of three months. Daily retrospective interviews about their time use were conducted, covering 1182 person-days of observation. Total time use of the sample population was assessed as well as time use stratified by managerial position. Differences of time use over time were also evaluated. Results District health managers used most of their working time for data management (16.6%), attending workshops (12.3%), financial management (8.7%), training of staff (7.1%), drug and supply management (5.0%), and travelling (9.6%). The study found significant variations of time use across the managerial cadres as well as high weekly variations of time use impulsed mainly by a national vertical program. Conclusions District health managers in Ghana use substantial amounts of their working time in only few activities and vertical programs greatly influence their time use. Our findings suggest that efficiency gains are possible for district health managers. However, these are unlikely to be achieved without improvements within the general health system, as inefficiencies seem to be largely caused by external factors. PMID:26068907

A model of the spatial and temporal variation of the Uranus thermal structure

NASA Technical Reports Server (NTRS)

Bezard, B.; Gautier, D.

1986-01-01

Seasonal variability of the temperature structure of Uranus is modeled for all latitudes in the .0004 to 2 bar pressure range in anticipation of the Voyager encounter in January 1986. Atmospheric heating in the model results on the one hand from an internal heat source and, on the other hand, from absorption of solar energy by methane and by non-conservative aerosols located between the 0.5 and 2 bar levels. Various cases for the behavior of the internal heat flux are investigated, such as constant with latitude or constrained to yield a time-averaged thermal emission independent of latitude. Meridional transport of heat in the stably stratified atmosphere is not taken into account. The results indicate that the Voyager encounter time, very small north-south temperature asymmetry should be expected. Moreover, the northern hemisphere, although not illuminated, should emit as much energy (within one percent) as the southern hemisphere at this date. At a given latitude, extreme temperatures are reached at the equinoxes. At the poles, seasonal amplitudes of about 10 K in the upper stratosphere and 6 K at the 0.6 bar level are predicted, and the variation with time of the emission to space is found to be at most 20 percent. The atmosphere of Uranus appears to be characterized by very long radiative response times (mainly due to its cold temperature) which inhibit the large seasonal variations that one could otherwise expect in view of the high obliquity of the planet and its long orbital period.

THE SPECTRAL-TIMING PROPERTIES OF UPPER AND LOWER kHz QPOs

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

Peille, Philippe; Barret, Didier; Uttley, Phil, E-mail: philippe.peille@irap.omp.eu

2015-10-01

Soft lags from the emission of the lower kilohertz quasi-periodic oscillations (kHz QPOs) of neutron star low-mass X-ray binaries have been reported from 4U1608-522 and 4U1636-536. Those lags hold prospects for constraining the origin of the QPO emission. In this paper, we investigate the spectral-timing properties of both the lower and upper kHz QPOs from the neutron star binary 4U1728-34, using the entire Rossi X-Ray Timing Explorer archive on this source. We show that the lag-energy spectra of the two QPOs are systematically different: while the lower kHz QPO shows soft lags, the upper kHz QPO shows either a flatmore » lag-energy spectrum or hard variations lagging softer variations. This suggests two different QPO-generation mechanisms. We also performed the first spectral deconvolution of the covariance spectra of both kHz QPOs. The QPO spectra are consistent with Comptonized blackbody emission, similar to the one found in the time-averaged spectrum, but with a higher seed-photon temperature, suggesting that a more compact inner region of the Comptonization layer (boundary/spreading layer, corona) is responsible for the QPO emission. Considering our results together with other recent findings, this leads us to the hypothesis that the lower kHz QPO signal is generated by coherent oscillations of the compact boundary layer region itself. The upper kHz QPO signal may then be linked to less-coherent accretion-rate variations produced in the inner accretion disk, and is then detected when they reach the boundary layer.« less

Timing Measurements and Their Implications for Four Binary Millisecond Pulsars

NASA Astrophysics Data System (ADS)

Bell, J. F.; Bailes, M.; Manchester, R. N.; Lyne, A. G.; Camilo, F.; Sandhu, J. S.

1997-04-01

We present timing observations of four millisecond pulsars, using data obtained over three years at the Australia Telescope National Facility (ATNF) Parkes and Nuffield Radio Astronomy Laboratory (NRAL) Jodrell Bank radio telescopes. Astrometric, spin and binary parameters are updated, and substantially improved for three pulsars, PSRs J0613-0200, J1045-4509 and J1643-1224. We have measured the time variation of the projected semimajor axis of the PSR J0437-4715 orbit due to its proper motion, and use it to constrain the inclination of the orbit and the mass of the companion. Some evidence is found for changes in the dispersion measures of PSRs J1045-4509 and J1643-1224. Limits are placed on the existence of planetary mass companions, ruling out companions with masses and orbits similar to the terrestrial planets of the Solar system for eight pulsars.

The rotation of Titan and Ganymede

NASA Astrophysics Data System (ADS)

Van Hoolst, Tim; Coyette, Alexis; Baland, Rose-Marie; Trinh, Antony

2016-10-01

The rotation rates of Titan and Ganymede, the largest satellites of Saturn and Jupiter, are on average equal to their orbital mean motion. Here we discuss small deviations from the average rotation for both satellites and evaluate the polar motion of Titan induced by its surface fluid layers. We examine different causes at various time scales and assess possible consequences and the potential of using librations and polar motion as probes of the interior structure of the satellites.The rotation rate of Titan and Ganymede cannot be constant on the orbital time scale as a result of the gravitational torque of the central planet acting on the satellites. Titan is moreover expected to show significant polar motion and additional variations in the rotation rate due to angular momentum exchange with the atmosphere, mainly at seasonal periods. Observational evidence for deviations from the synchronous state has been reported several times for Titan but is unfortunately inconclusive. The measurements of the rotation variations are based on determinations of the shift in position of Cassini radar images taken during different flybys. The ESA JUICE (JUpiter ICy moons Explorer) mission will measure the rotation variations of Ganymede during its orbital phase around the satellite starting in 2032.We report on different theoretical aspects of the librations and polar motion. We consider the influence of the rheology of the ice shell and take into account Cassini measurements of the external gravitational field and of the topography of Titan and similar Galileo data about Ganymede. We also evaluate the librations and polar motion induced by Titan's hydrocarbon seas and use the most recent results of Titan's atmosphere dynamics. We finally evaluate the potential of rotation variations to constrain the satellite's interior structure, in particular its ice shell and ocean.

Variation in Seed Germination of 134 Common Species on the Eastern Tibetan Plateau: Phylogenetic, Life History and Environmental Correlates

PubMed Central

Xu, Jing; Li, Wenlong; Zhang, Chunhui; Liu, Wei; Du, Guozhen

2014-01-01

Seed germination is a crucial stage in the life history of a species because it represents the pathway from adult to offspring, and it can affect the distribution and abundance of species in communities. In this study, we examined the effects of phylogenetic, life history and environmental factors on seed germination of 134 common species from an alpine/subalpine meadow on the eastern Tibetan Plateau. In one-way ANOVAs, phylogenetic groups (at or above order) explained 13.0% and 25.9% of the variance in germination percentage and mean germination time, respectively; life history attributes, such as seed size, dispersal mode, explained 3.7%, 2.1% of the variance in germination percentage and 6.3%, 8.7% of the variance in mean germination time, respectively; the environmental factors temperature and habitat explained 4.7%, 1.0% of the variance in germination percentage and 13.5%, 1.7% of the variance in mean germination time, respectively. Our results demonstrated that elevated temperature would lead to a significant increase in germination percentage and an accelerated germination. Multi-factorial ANOVAs showed that the three major factors contributing to differences in germination percentage and mean germination time in this alpine/subalpine meadow were phylogenetic attributes, temperature and seed size (explained 10.5%, 4.7% and 1.4% of the variance in germination percentage independently, respectively; and explained 14.9%, 13.5% and 2.7% of the variance in mean germination time independently, respectively). In addition, there were strong associations between phylogenetic group and life history attributes, and between life history attributes and environmental factors. Therefore, germination variation are constrained mainly by phylogenetic inertia in a community, and seed germination variation correlated with phylogeny is also associated with life history attributes, suggesting a role of niche adaptation in the conservation of germination variation within lineages. Meanwhile, selection can maintain the association between germination behavior and the environmental conditions within a lineage. PMID:24893308

Ecological and sampling constraints on defining landscape fire severity

USGS Publications Warehouse

Key, C.H.

2006-01-01

Ecological definition and detection of fire severity are influenced by factors of spatial resolution and timing. Resolution determines the aggregation of effects within a sampling unit or pixel (alpha variation), hence limiting the discernible ecological responses, and controlling the spatial patchiness of responses distributed throughout a burn (beta variation). As resolution decreases, alpha variation increases, extracting beta variation and complexity from the spatial model of the whole burn. Seasonal timing impacts the quality of radiometric data in terms of transmittance, sun angle, and potential contrast between responses within burns. Detection sensitivity candegrade toward the end of many fire seasons when low sun angles, vegetation senescence, incomplete burning, hazy conditions, or snow are common. Thus, a need exists to supersede many rapid response applications when remote sensing conditions improve. Lag timing, or timesince fire, notably shapes the ecological character of severity through first-order effects that only emerge with time after fire, including delayed survivorship and mortality. Survivorship diminishes the detected magnitude of severity, as burned vegetation remains viable and resprouts, though at first it may appear completely charred or consumed above ground. Conversely, delayed mortality increases the severity estimate when apparently healthy vegetation is in fact damaged by heat to the extent that it dies over time. Both responses dependon fire behavior and various species-specific adaptations to fire that are unique to the pre-firecomposition of each burned area. Both responses can lead initially to either over- or underestimating severity. Based on such implications, three sampling intervals for short-term burn severity are identified; rapid, initial, and extended assessment, sampled within about two weeks, two months, and depending on the ecotype, from three months to one year after fire, respectively. Spatial and temporal conditions of sampling strategies constrain data quality and ecological information obtained about fire severity. Though commonly overlooked, such considerations determine the objectives and hypotheses that are appropriate for each application, and are especially important when building comparative studies or long-term reference databases on fire severity.

Risk of herbivore attack and heritability of ontogenetic trajectories in plant defense.

PubMed

Ochoa-López, Sofía; Rebollo, Roberto; Barton, Kasey E; Fornoni, Juan; Boege, Karina

2018-06-01

Ontogeny has been identified as a main source of variation in the expression of plant phenotypes. However, there is limited information on the mechanisms behind the evolution of ontogenetic trajectories in plant defense. We explored if risk of attack, herbivore damage, heritability, and phenotypic plasticity can promote or constrain the evolutionary potential of ontogenetic trajectories in three defensive traits. We exposed 20 genotypes of Turnera velutina to contrasting environments (shadehouse and field plots), and measured the cyanogenic potential, trichome density, and sugar content in extrafloral nectar in seedlings, juveniles and reproductive plants. We also assessed risk of attack through oviposition preferences, and quantified herbivore damage in the field. We estimated genetic variance, broad sense heritability, and evolvability of the defensive traits at each ontogenetic stage, and of the ontogenetic trajectories themselves. For plants growing in the shadehouse, we found genetic variation and broad sense heritability for cyanogenic potential in seedlings, and for trichome density at all ontogenetic stages. Genetic variation and heritability of ontogenetic trajectories was detected for trichome density only. These genetic pre-requisites for evolution, however, were not detected in the field, suggesting that environmental variation and phenotypic plastic responses mask any heritable variation. Finally, ontogenetic trajectories were found to be plastic, differing between shadehouse and field conditions for the same genetic families. Overall, we provide support for the idea that changes in herbivore pressure can be a mechanism behind the evolution of ontogenetic trajectories. This evolutionary potential, however, can be constrained by phenotypic plasticity expressed in heterogeneous environments.

Stable isotope series from elephant ivory reveal lifetime histories of a true dietary generalist

PubMed Central

Codron, Jacqueline; Codron, Daryl; Sponheimer, Matt; Kirkman, Kevin; Duffy, Kevin J.; Raubenheimer, Erich J.; Mélice, Jean-Luc; Grant, Rina; Clauss, Marcus; Lee-Thorp, Julia A.

2012-01-01

Longitudinal studies have revealed how variation in resource use within consumer populations can impact their dynamics and functional significance in communities. Here, we investigate multi-decadal diet variations within individuals of a keystone megaherbivore species, the African elephant (Loxodonta africana), using serial stable isotope analysis of tusks from the Kruger National Park, South Africa. These records, representing the longest continuous diet histories documented for any extant species, reveal extensive seasonal and annual variations in isotopic—and hence dietary—niches of individuals, but little variation between them. Lack of niche distinction across individuals contrasts several recent studies, which found relatively high levels of individual niche specialization in various taxa. Our result is consistent with theory that individual mammal herbivores are nutritionally constrained to maintain broad diet niches. Individual diet specialization would also be a costly strategy for large-bodied taxa foraging over wide areas in spatio-temporally heterogeneous environments. High levels of within-individual diet variability occurred within and across seasons, and persisted despite an overall increase in inferred C4 grass consumption through the twentieth century. We suggest that switching between C3 browsing and C4 grazing over extended time scales facilitates elephant survival through environmental change, and could even allow recovery of overused resources. PMID:22337695

Stable isotope series from elephant ivory reveal lifetime histories of a true dietary generalist.

PubMed

Codron, Jacqueline; Codron, Daryl; Sponheimer, Matt; Kirkman, Kevin; Duffy, Kevin J; Raubenheimer, Erich J; Mélice, Jean-Luc; Grant, Rina; Clauss, Marcus; Lee-Thorp, Julia A

2012-06-22

Longitudinal studies have revealed how variation in resource use within consumer populations can impact their dynamics and functional significance in communities. Here, we investigate multi-decadal diet variations within individuals of a keystone megaherbivore species, the African elephant (Loxodonta africana), using serial stable isotope analysis of tusks from the Kruger National Park, South Africa. These records, representing the longest continuous diet histories documented for any extant species, reveal extensive seasonal and annual variations in isotopic--and hence dietary--niches of individuals, but little variation between them. Lack of niche distinction across individuals contrasts several recent studies, which found relatively high levels of individual niche specialization in various taxa. Our result is consistent with theory that individual mammal herbivores are nutritionally constrained to maintain broad diet niches. Individual diet specialization would also be a costly strategy for large-bodied taxa foraging over wide areas in spatio-temporally heterogeneous environments. High levels of within-individual diet variability occurred within and across seasons, and persisted despite an overall increase in inferred C(4) grass consumption through the twentieth century. We suggest that switching between C(3) browsing and C(4) grazing over extended time scales facilitates elephant survival through environmental change, and could even allow recovery of overused resources.

Minimal complexity control law synthesis

NASA Technical Reports Server (NTRS)

Bernstein, Dennis S.; Haddad, Wassim M.; Nett, Carl N.

1989-01-01

A paradigm for control law design for modern engineering systems is proposed: Minimize control law complexity subject to the achievement of a specified accuracy in the face of a specified level of uncertainty. Correspondingly, the overall goal is to make progress towards the development of a control law design methodology which supports this paradigm. Researchers achieve this goal by developing a general theory of optimal constrained-structure dynamic output feedback compensation, where here constrained-structure means that the dynamic-structure (e.g., dynamic order, pole locations, zero locations, etc.) of the output feedback compensation is constrained in some way. By applying this theory in an innovative fashion, where here the indicated iteration occurs over the choice of the compensator dynamic-structure, the paradigm stated above can, in principle, be realized. The optimal constrained-structure dynamic output feedback problem is formulated in general terms. An elegant method for reducing optimal constrained-structure dynamic output feedback problems to optimal static output feedback problems is then developed. This reduction procedure makes use of star products, linear fractional transformations, and linear fractional decompositions, and yields as a byproduct a complete characterization of the class of optimal constrained-structure dynamic output feedback problems which can be reduced to optimal static output feedback problems. Issues such as operational/physical constraints, operating-point variations, and processor throughput/memory limitations are considered, and it is shown how anti-windup/bumpless transfer, gain-scheduling, and digital processor implementation can be facilitated by constraining the controller dynamic-structure in an appropriate fashion.

An additional pulsating mode (7.35 mHz) and pulsations timing variations of PG 1613+426

NASA Astrophysics Data System (ADS)

Otani, Tomomi; Oswalt, Terry D.; Majewski, Patrice; Jordan, Riley; Amaral, Marc; Moss, Adam

2017-12-01

We present the detection of an additional pulsation mode (7.35 mHz) of a subdwarf B star, PG 1613+426, and periodic Observed minus Calculated (O-C) variations for two existing pulsations. PG 1613+426 is near the hot end of the sdB instability strip. One pulsation mode (6.94 mHz) was detected so far by Bonanno et al. (2002) and another pulsation mode candidate (7.05 mHz) was proposed with a confidence level above 90% by Kuassivi and Ferlet (2005). To constrain sdB star evolutional scenarios, this star was monitored in 2010, 2011, 2015, and 2017 as a part of a project for finding companions to sdB stars using the pulsation timing method. The photometric analysis of those data shows an additional 7.35 mHz pulsation mode as well as the previously detected 6.93 mHz mode. However the 7.05 mHz mode was not detected. Nightly amplitude changes of 7.35 mHz mode were observed in the 2011 data, however the 2017 data did not show nightly amplitude shifts. O-C variations were detected in both 6.93 mHz and 7.35 mHz pulsations, indicating that PG 1613+426 may have a low mass companion star. However, more observations are needed to confirm it.

An Alcohol Test for Drifting Constants

NASA Astrophysics Data System (ADS)

Jansen, P.; Bagdonaite, J.; Ubachs, W.; Bethlem, H. L.; Kleiner, I.; Xu, L.-H.

2013-06-01

The Standard Model of physics is built on the fundamental constants of nature, however without providing an explanation for their values, nor requiring their constancy over space and time. Molecular spectroscopy can address this issue. Recently, we found that microwave transitions in methanol are extremely sensitive to a variation of the proton-to-electron mass ratio μ, due to a fortuitous interplay between classically forbidden internal rotation and rotation of the molecule as a whole. In this talk, we will explain the origin of this effect and how the sensitivity coefficients in methanol are calculated. In addition, we set a limit on a possible cosmological variation of μ by comparing transitions in methanol observed in the early Universe with those measured in the laboratory. Based on radio-astronomical observations of PKS1830-211, we deduce a constraint of Δμ/μ=(0.0± 1.0)× 10^{-7} at redshift z = 0.89, corresponding to a look-back time of 7 billion years. While this limit is more constraining and systematically more robust than previous ones, the methanol method opens a new search territory for probing μ-variation on cosmological timescales. P. Jansen, L.-H. Xu, I. Kleiner, W. Ubachs, and H.L. Bethlem Phys. Rev. Lett. {106}(100801) 2011. J. Bagdonaite, P. Jansen, C. Henkel, H.L. Bethlem, K.M. Menten, and W. Ubachs Science {339}(46) 2013.

Multi-Compartment T2 Relaxometry Using a Spatially Constrained Multi-Gaussian Model

PubMed Central

Raj, Ashish; Pandya, Sneha; Shen, Xiaobo; LoCastro, Eve; Nguyen, Thanh D.; Gauthier, Susan A.

2014-01-01

The brain’s myelin content can be mapped by T2-relaxometry, which resolves multiple differentially relaxing T2 pools from multi-echo MRI. Unfortunately, the conventional fitting procedure is a hard and numerically ill-posed problem. Consequently, the T2 distributions and myelin maps become very sensitive to noise and are frequently difficult to interpret diagnostically. Although regularization can improve stability, it is generally not adequate, particularly at relatively low signal to noise ratio (SNR) of around 100–200. The purpose of this study was to obtain a fitting algorithm which is able to overcome these difficulties and generate usable myelin maps from noisy acquisitions in a realistic scan time. To this end, we restrict the T2 distribution to only 3 distinct resolvable tissue compartments, modeled as Gaussians: myelin water, intra/extra-cellular water and a slow relaxing cerebrospinal fluid compartment. We also impose spatial smoothness expectation that volume fractions and T2 relaxation times of tissue compartments change smoothly within coherent brain regions. The method greatly improves robustness to noise, reduces spatial variations, improves definition of white matter fibers, and enhances detection of demyelinating lesions. Due to efficient design, the additional spatial aspect does not cause an increase in processing time. The proposed method was applied to fast spiral acquisitions on which conventional fitting gives uninterpretable results. While these fast acquisitions suffer from noise and inhomogeneity artifacts, our preliminary results indicate the potential of spatially constrained 3-pool T2 relaxometry. PMID:24896833

Measurements of Kepler Planet Masses and Eccentricities from Transit Timing Variations: Analytic and N-body Results

NASA Astrophysics Data System (ADS)

Hadden, Sam; Lithwick, Yoram

2015-12-01

Several Kepler planets reside in multi-planet systems where gravitational interactions result in transit timing variations (TTVs) that provide exquisitely sensitive probes of their masses of and orbits. Measuring these planets' masses and orbits constrains their bulk compositions and can provide clues about their formation. However, inverting TTV measurements in order to infer planet properties can be challenging: it involves fitting a nonlinear model with a large number of parameters to noisy data, often with significant degeneracies between parameters. I present results from two complementary approaches to TTV inversion: Markov chain Monte Carlo simulations that use N-body integrations to compute transit times and a simplified analytic model for computing the TTVs of planets near mean motion resonances. The analytic model allows for straightforward interpretations of N-body results and provides an independent estimate of parameter uncertainties that can be compared to MCMC results which may be sensitive to factors such as priors. We have conducted extensive MCMC simulations along with analytic fits to model the TTVs of dozens of Kepler multi-planet systems. We find that the bulk of these sub-Jovian planets have low densities that necessitate significant gaseous envelopes. We also find that the planets' eccentricities are generally small but often definitively non-zero.

Discovery of pulsed OH maser emission stimulated by a pulsar.

PubMed

Weisberg, Joel M; Johnston, Simon; Koribalski, Bärbel; Stanimirovic, Snezana

2005-07-01

Stimulated emission of radiation has not been directly observed in astrophysical situations up to this time. Here we demonstrate that photons from pulsar B1641-45 stimulate pulses of excess 1720-megahertz line emission in an interstellar hydroxyl (OH) cloud. As this stimulated emission is driven by the pulsar, it varies on a few-millisecond time scale, which is orders of magnitude shorter than the quickest OH maser variations previously detected. Our 1612-megahertz spectra are inverted copies of the 1720-megahertz spectra. This "conjugate line" phenomenon enables us to constrain the properties of the interstellar OH line-producing gas. We also show that pulsar signals undergo significantly deeper OH absorption than do other background sources, which confirms earlier tentative findings that OH clouds are clumpier on small scales than are neutral hydrogen clouds.

Extreme-scale motions in turbulent plane Couette flows

NASA Astrophysics Data System (ADS)

Lee, Myoungkyu; Moser, Robert D.

2018-05-01

We study the size of large-scale motions in turbulent plane Couette flows at moderate Reynolds number up to $Re_\\tau$ = 500. Direct numerical simulation domains were as large as $100\\pi\\delta\\times2\\delta\\times5\\pi\\delta$, where $\\delta$ is half the distance between the walls. The results indicate that there are structures with streamwise extent, as measured by the wavelength, as long as 78$\\delta$ and at least 310$\\delta$ at $Re_\\tau$ = 220 and 500, respectively. The presence of these very long structures is apparent in the spectra of all three velocity components and the Reynolds stress. In DNS using a smaller domain, the large structures are constrained, eliminating the streamwise variations present in the larger domain. Effects of a smaller domain are also present in the mean velocity and the streamwise velocity variance in the outer flow.

Semismooth Newton method for gradient constrained minimization problem

NASA Astrophysics Data System (ADS)

Anyyeva, Serbiniyaz; Kunisch, Karl

2012-08-01

In this paper we treat a gradient constrained minimization problem, particular case of which is the elasto-plastic torsion problem. In order to get the numerical approximation to the solution we have developed an algorithm in an infinite dimensional space framework using the concept of the generalized (Newton) differentiation. Regularization was done in order to approximate the problem with the unconstrained minimization problem and to make the pointwise maximum function Newton differentiable. Using semismooth Newton method, continuation method was developed in function space. For the numerical implementation the variational equations at Newton steps are discretized using finite elements method.

Use of constrained optimization in the conceptual design of a medium-range subsonic transport

NASA Technical Reports Server (NTRS)

Sliwa, S. M.

1980-01-01

Constrained parameter optimization was used to perform the optimal conceptual design of a medium range transport configuration. The impact of choosing a given performance index was studied, and the required income for a 15 percent return on investment was proposed as a figure of merit. A number of design constants and constraint functions were systematically varied to document the sensitivities of the optimal design to a variety of economic and technological assumptions. A comparison was made for each of the parameter variations between the baseline configuration and the optimally redesigned configuration.

Variational Trajectory Optimization Tool Set: Technical description and user's manual

NASA Technical Reports Server (NTRS)

Bless, Robert R.; Queen, Eric M.; Cavanaugh, Michael D.; Wetzel, Todd A.; Moerder, Daniel D.

1993-01-01

The algorithms that comprise the Variational Trajectory Optimization Tool Set (VTOTS) package are briefly described. The VTOTS is a software package for solving nonlinear constrained optimal control problems from a wide range of engineering and scientific disciplines. The VTOTS package was specifically designed to minimize the amount of user programming; in fact, for problems that may be expressed in terms of analytical functions, the user needs only to define the problem in terms of symbolic variables. This version of the VTOTS does not support tabular data; thus, problems must be expressed in terms of analytical functions. The VTOTS package consists of two methods for solving nonlinear optimal control problems: a time-domain finite-element algorithm and a multiple shooting algorithm. These two algorithms, under the VTOTS package, may be run independently or jointly. The finite-element algorithm generates approximate solutions, whereas the shooting algorithm provides a more accurate solution to the optimization problem. A user's manual, some examples with results, and a brief description of the individual subroutines are included.

Global solar wind variations over the last four centuries.

PubMed

Owens, M J; Lockwood, M; Riley, P

2017-01-31

The most recent "grand minimum" of solar activity, the Maunder minimum (MM, 1650-1710), is of great interest both for understanding the solar dynamo and providing insight into possible future heliospheric conditions. Here, we use nearly 30 years of output from a data-constrained magnetohydrodynamic model of the solar corona to calibrate heliospheric reconstructions based solely on sunspot observations. Using these empirical relations, we produce the first quantitative estimate of global solar wind variations over the last 400 years. Relative to the modern era, the MM shows a factor 2 reduction in near-Earth heliospheric magnetic field strength and solar wind speed, and up to a factor 4 increase in solar wind Mach number. Thus solar wind energy input into the Earth's magnetosphere was reduced, resulting in a more Jupiter-like system, in agreement with the dearth of auroral reports from the time. The global heliosphere was both smaller and more symmetric under MM conditions, which has implications for the interpretation of cosmogenic radionuclide data and resulting total solar irradiance estimates during grand minima.

Interactive effects of genotype and food quality on consumer growth rate and elemental content.

PubMed

Prater, Clay; Wagner, Nicole D; Frost, Paul C

2017-05-01

Consumer body stoichiometry is a key trait that links organismal physiology to population and ecosystem-level dynamics. However, as elemental composition has traditionally been considered to be constrained within a species, the ecological and evolutionary factors shaping consumer elemental composition have not been clearly resolved. To this end, we examined the causes and extent of variation in the body phosphorus (P) content and the expression of P-linked traits, mass specific growth rate (MSGR), and P use efficiency (PUE) of the keystone aquatic consumer Daphnia using lake surveys and common garden experiments. While daphnid body %P was relatively constrained in field assemblages sampled across an environmental P gradient, unique genotypes isolated from these lakes showed highly variable phenotypic responses when raised across dietary P gradients in the laboratory. Specifically, we observed substantial inter- and intra-specific variation and differences in daphnid responses within and among our study lakes. While variation in Daphnia body %P was mostly due to plastic phenotypic changes, we documented considerable genetic differences in daphnid MSGR and PUE, and relationships between MSGR and body P content were highly variable among genotypes. Overall, our study found that consumer responses to food quality may differ considerably among genotypes and that relationships between organismal life-history traits and body stoichiometry may be strongly influenced by genetic and environmental variation in natural assemblages. © 2017 by the Ecological Society of America.

80-Myr history of buoyancy and volcanic fluxes along the trails of the Walvis and St. Helena hotspots (South Atlantic)

NASA Astrophysics Data System (ADS)

Vidal, V.; Adam, C.; Escartin, J.

2007-12-01

Walvis and St.~Helena are the only long-lived hotspot chains in the South Atlantic. Therefore, their characterization is important to constrain the processes associated with mantle plume formation, their temporal evolution, and the interaction with plate and mantle dynamics in the region. We study the temporal evolution of plume buoyancy and magma production rate along both hotspot chains, which are constrained from the swell and volume of volcanic materials emplaced along the chain. The regional depth anomaly is calculated by correcting the 2' bathymetry grid of Smith & Sandwell (1997) for thermal subsidence and sediment loading. We separate the topography associated with volcanism and the swell surrounding the hotspot chains using the MiFil filtering method (Adam et al., 2005). We then estimate the temporal variations associated with both parameters by computing volumes along the hotspot tracks. Neither Walvis nor St.~Helena show a 'classical' hotspot behavior. We find that two plumes are at the origin of the St.~Helena chain. This study also shows a swell associated with the Circe seamount, supporting the existence of a hotspot NW of the St.~Helena trail. The variation in swell and volcanic fluxes suggests temporal variability in the plume behavior at time scales of 10-20~m.y. and 5~m.y., which may be related to oscillations and instabilities of the plume conduit, respectively. Cumulative fluxes in the area are largest for Walvis and weakest for Circe, and all are significantly lower than that reported for the Hawai'i hotspot.

Numeric stratigraphic modeling: Testing sequence Numeric stratigraphic modeling: Testing sequence stratigraphic concepts using high resolution geologic examples

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

Armentrout, J.M.; Smith-Rouch, L.S.; Bowman, S.A.

1996-08-01

Numeric simulations based on integrated data sets enhance our understanding of depositional geometry and facilitate quantification of depositional processes. Numeric values tested against well-constrained geologic data sets can then be used in iterations testing each variable, and in predicting lithofacies distributions under various depositional scenarios using the principles of sequence stratigraphic analysis. The stratigraphic modeling software provides a broad spectrum of techniques for modeling and testing elements of the petroleum system. Using well-constrained geologic examples, variations in depositional geometry and lithofacies distributions between different tectonic settings (passive vs. active margin) and climate regimes (hothouse vs. icehouse) can provide insight tomore » potential source rock and reservoir rock distribution, maturation timing, migration pathways, and trap formation. Two data sets are used to illustrate such variations: both include a seismic reflection profile calibrated by multiple wells. The first is a Pennsylvanian mixed carbonate-siliciclastic system in the Paradox basin, and the second a Pliocene-Pleistocene siliciclastic system in the Gulf of Mexico. Numeric simulations result in geometry and facies distributions consistent with those interpreted using the integrated stratigraphic analysis of the calibrated seismic profiles. An exception occurs in the Gulf of Mexico study where the simulated sediment thickness from 3.8 to 1.6 Ma within an upper slope minibasin was less than that mapped using a regional seismic grid. Regional depositional patterns demonstrate that this extra thickness was probably sourced from out of the plane of the modeled transect, illustrating the necessity for three-dimensional constraints on two-dimensional modeling.« less

A variational approach to parameter estimation in ordinary differential equations.

PubMed

Kaschek, Daniel; Timmer, Jens

2012-08-14

Ordinary differential equations are widely-used in the field of systems biology and chemical engineering to model chemical reaction networks. Numerous techniques have been developed to estimate parameters like rate constants, initial conditions or steady state concentrations from time-resolved data. In contrast to this countable set of parameters, the estimation of entire courses of network components corresponds to an innumerable set of parameters. The approach presented in this work is able to deal with course estimation for extrinsic system inputs or intrinsic reactants, both not being constrained by the reaction network itself. Our method is based on variational calculus which is carried out analytically to derive an augmented system of differential equations including the unconstrained components as ordinary state variables. Finally, conventional parameter estimation is applied to the augmented system resulting in a combined estimation of courses and parameters. The combined estimation approach takes the uncertainty in input courses correctly into account. This leads to precise parameter estimates and correct confidence intervals. In particular this implies that small motifs of large reaction networks can be analysed independently of the rest. By the use of variational methods, elements from control theory and statistics are combined allowing for future transfer of methods between the two fields.

Analyzing endocrine system conservation and evolution.

PubMed

Bonett, Ronald M

2016-08-01

Analyzing variation in rates of evolution can provide important insights into the factors that constrain trait evolution, as well as those that promote diversification. Metazoan endocrine systems exhibit apparent variation in evolutionary rates of their constituent components at multiple levels, yet relatively few studies have quantified these patterns and analyzed them in a phylogenetic context. This may be in part due to historical and current data limitations for many endocrine components and taxonomic groups. However, recent technological advancements such as high-throughput sequencing provide the opportunity to collect large-scale comparative data sets for even non-model species. Such ventures will produce a fertile data landscape for evolutionary analyses of nucleic acid and amino acid based endocrine components. Here I summarize evolutionary rate analyses that can be applied to categorical and continuous endocrine traits, and also those for nucleic acid and protein-based components. I emphasize analyses that could be used to test whether other variables (e.g., ecology, ontogenetic timing of expression, etc.) are related to patterns of rate variation and endocrine component diversification. The application of phylogenetic-based rate analyses to comparative endocrine data will greatly enhance our understanding of the factors that have shaped endocrine system evolution. Copyright © 2016 Elsevier Inc. All rights reserved.

Life cycle stage and water depth affect flooding-induced adventitious root formation in the terrestrial species Solanum dulcamara

PubMed Central

Zhang, Qian; Visser, Eric J. W.; de Kroon, Hans; Huber, Heidrun

2015-01-01

Background and Aims Flooding can occur at any stage of the life cycle of a plant, but often adaptive responses of plants are only studied at a single developmental stage. It may be anticipated that juvenile plants may respond differently from mature plants, as the amount of stored resources may differ and morphological changes can be constrained. Moreover, different water depths may require different strategies to cope with the flooding stress, the expression of which may also depend on developmental stage. This study investigated whether flooding-induced adventitious root formation and plant growth were affected by flooding depth in Solanum dulcamara plants at different developmental stages. Methods Juvenile plants without pre-formed adventitious root primordia and mature plants with primordia were subjected to shallow flooding or deep flooding for 5 weeks. Plant growth and the timing of adventitious root formation were monitored during the flooding treatments. Key Results Adventitious root formation in response to shallow flooding was significantly constrained in juvenile S. dulcamara plants compared with mature plants, and was delayed by deep flooding compared with shallow flooding. Complete submergence suppressed adventitious root formation until up to 2 weeks after shoots restored contact with the atmosphere. Independent of developmental stage, a strong positive correlation was found between adventitious root formation and total biomass accumulation during shallow flooding. Conclusions The potential to deploy an escape strategy (i.e. adventitious root formation) may change throughout a plant’s life cycle, and is largely dependent on flooding depth. Adaptive responses at a given stage of the life cycle thus do not necessarily predict how the plant responds to flooding in another growth stage. As variation in adventitious root formation also correlates with finally attained biomass, this variation may form the basis for variation in resistance to shallow flooding among plants. PMID:26105188

Life cycle stage and water depth affect flooding-induced adventitious root formation in the terrestrial species Solanum dulcamara.

PubMed

Zhang, Qian; Visser, Eric J W; de Kroon, Hans; Huber, Heidrun

2015-08-01

Flooding can occur at any stage of the life cycle of a plant, but often adaptive responses of plants are only studied at a single developmental stage. It may be anticipated that juvenile plants may respond differently from mature plants, as the amount of stored resources may differ and morphological changes can be constrained. Moreover, different water depths may require different strategies to cope with the flooding stress, the expression of which may also depend on developmental stage. This study investigated whether flooding-induced adventitious root formation and plant growth were affected by flooding depth in Solanum dulcamara plants at different developmental stages. Juvenile plants without pre-formed adventitious root primordia and mature plants with primordia were subjected to shallow flooding or deep flooding for 5 weeks. Plant growth and the timing of adventitious root formation were monitored during the flooding treatments. Adventitious root formation in response to shallow flooding was significantly constrained in juvenile S. dulcamara plants compared with mature plants, and was delayed by deep flooding compared with shallow flooding. Complete submergence suppressed adventitious root formation until up to 2 weeks after shoots restored contact with the atmosphere. Independent of developmental stage, a strong positive correlation was found between adventitious root formation and total biomass accumulation during shallow flooding. The potential to deploy an escape strategy (i.e. adventitious root formation) may change throughout a plant's life cycle, and is largely dependent on flooding depth. Adaptive responses at a given stage of the life cycle thus do not necessarily predict how the plant responds to flooding in another growth stage. As variation in adventitious root formation also correlates with finally attained biomass, this variation may form the basis for variation in resistance to shallow flooding among plants. © The Author 2015. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For Permissions, please email: journals.permissions@oup.com.

The Evolution of the Seismic-Aseismic Transition During the Earthquake Cycle: Constraints from the Time-Dependent Depth Distribution of Aftershocks

NASA Astrophysics Data System (ADS)

Rolandone, F.; Bürgmann, R.; Nadeau, R.; Freed, A.

2003-12-01

We have demonstrated that in the aftermath of large earthquakes, the depth extent of aftershocks shows an immediate deepening from pre-earthquake levels, followed by a time-dependent postseismic shallowing. We use these seismic data to constrain the variation of the depth of the seismic-aseismic transition with time throughout the earthquake cycle. Most studies of the seismic-aseismic transition have focussed on the effect of temperature and/or lithology on the transition either from brittle faulting to viscous flow or from unstable to stable sliding. They have shown that the maximum depth of seismic activity is well correlated with the spatial variations of these two parameters. However, little has been done to examine how the maximum depth of seismogenic faulting varies locally, at the scale of a fault segment, during the course of the earthquake cycle. Geologic and laboratory observations indicate that the depth of the seismic-aseismic transition should vary with strain rate and thus change with time throughout the earthquake cycle. We quantify the time-dependent variations in the depth of seismicity on various strike-slip faults in California before and after large earthquakes. We specifically investigate (1) the deepening of the aftershocks relative to the background seismicity, (2) the time constant of the postseismic shallowing of the deepest earthquakes, and (3) the correlation of the time-dependent pattern with the coseismic slip distribution and the expected stress increase. Together with geodetic measurements, these seismological observations form the basis for developing more sophisticated models for the mechanical evolution of strike-slip shear zones during the earthquake cycle. We develop non-linear viscoelastic models, for which the brittle-ductile transition is not fixed, but varies with assumed temperature and calculated stress gradients. We use them to place constraints on strain rate at depth, on time-dependent rheology, and on the partitioning of deformation between brittle faulting and distributed viscous flow associated with the earthquake cycle.

Time functions of deep earthquakes from broadband and short-period stacks

USGS Publications Warehouse

Houston, H.; Benz, H.M.; Vidale, J.E.

1998-01-01

To constrain dynamic source properties of deep earthquakes, we have systematically constructed broadband time functions of deep earthquakes by stacking and scaling teleseismic P waves from U.S. National Seismic Network, TERRAscope, and Berkeley Digital Seismic Network broadband stations. We examined 42 earthquakes with depths from 100 to 660 km that occurred between July 1, 1992 and July 31, 1995. To directly compare time functions, or to group them by size, depth, or region, it is essential to scale them to remove the effect of moment, which varies by more than 3 orders of magnitude for these events. For each event we also computed short-period stacks of P waves recorded by west coast regional arrays. The comparison of broadband with short-period stacks yields a considerable advantage, enabling more reliable measurement of event duration. A more accurate estimate of the duration better constrains the scaling procedure to remove the effect of moment, producing scaled time functions with both correct timing and amplitude. We find only subtle differences in the broadband time-function shape with moment, indicating successful scaling and minimal effects of attenuation at the periods considered here. The average shape of the envelopes of the short-period stacks is very similar to the average broadband time function. The main variations seen with depth are (1) a mild decrease in duration with increasing depth, (2) greater asymmetry in the time functions of intermediate events compared to deep ones, and (3) unexpected complexity and late moment release for events between 350 and 550 km, with seven of the eight events in that depth interval displaying markedly more complicated time functions with more moment release late in the rupture than most events above or below. The first two results are broadly consistent with our previous studies, while the third is reported here for the first time. The greater complexity between 350 and 550 km suggests greater heterogeneity in the failure process in that depth range. Copyright 1998 by the American Geophysical Union.

Seismic velocity anisotropy and heterogeneity beneath the Mantle Electromagnetic and Tomography Experiment (MELT) region of the East Pacific Rise from analysis of P and S body waves

USGS Publications Warehouse

Hammond, W.C.; Toomey, D.R.

2003-01-01

We use teleseismic P and S delay times and shear wave splitting measurements to constrain isotropic and anisotropic heterogeneity in the mantle beneath the southern East Pacific Rise (SEPR). The data comprise 462 P and S delay times and 18 shear wave splitting observations recorded during the Mantle Electromagnetic and Tomography (MELT) Experiment. We estimate the mantle melt content (F) and temperature (T) variation from the isotropic velocity variation. Our results indicate that the maximum variation in F beneath our array is between zero and ???1.2%, and maximum variation in T is between zero and ???100 K. We favor an explanation having partial contributions from both T and F. We approximate the seismic anisotropy of the upper mantle with hexagonal symmetry, consistent with the assumption of two dimensionality of mantle flow. Our new tomographic technique uses a nonlinear inversion of P and slow S polarization delay times to simultaneously solve for coupled VP and VS heterogeneity throughout the model and for the magnitude of anisotropy within discrete domains. The domain dimensions and the dip of the anisotropy are fixed for each inversion but are varied in a grid search, obtaining the misfit of the models to the body wave delay data and to split times of vertically propagating S waves. The data misfit and the isotropic heterogeneity are sensitive to domain dimensions and dip of anisotropy. In a region centered beneath the SEPR the best average dip of the hexagonal symmetry axis is horizontal or dipping shallowly (<30??) west. Given the resolution of our data, a subaxial region characterized by vertically aligned symmetry axes may exist but is limited to be <80 km deep. We infer that the mantle flow beneath the SEPR is consistent with shallow asthenospheric return flow from the direction of the South Pacific superswell.

Constraining Earth System Models in the Tropics with Multiple Satellite Observations

NASA Astrophysics Data System (ADS)

Shi, M.; Liu, J.; Saatchi, S. S.; Chan, S.; Yu, Y.; Zhao, M.

2016-12-01

Because of the impacts of cloud and atmospheric aerosol on spectral observations and the saturation of spectral observations over dense forests, the current spectral observations (e.g., Moderate Resolution Imaging Spectroradiometer) have large uncertainties in the tropics. Nevertheless, the backscatter observations from the SeaWinds Scatterometer onboard QuikSCAT (QSCAT) are sensitive to the variations of canopy water content and structure of forest canopy, and are not affected by clouds and atmospheric aerosols. In addition, the lack of sensitivity of the Soil Moisture Active Passive (SMAP) Level 1C brightness temperature (TB) to soil moisture under dense forest canopies (e.g., forests in tropics) makes the SMAP TB data a direct indicator of canopy properties. In this study, we use a variety of new satellite observations, including the QSCAT backscatter observations, the Gravity Recovery and Climate Experiment (GRACE) satellite's observed temporal gravity field variations, and the SMAP Level 1C TB, to constrain the carbon (C) cycle simulated by the Community Land Model version 4.5 BGC (CLM4.5) for the 2005 Amazonia drought and 2015 El Nino. Our results show that the leaf C pool size simulated by CLM4.5 decreases dramatically in southwest Amazonia in the 2005 drought, and recovers slowly afterward (after about 3 years). This result is consistent with the long-term C-recovery after the 2005 Amazonia drought observed by QSCAT. The slow C pool recovery is associated with large fire disturbance and the slow water storage recovery simulated by CLM4.5 and observed by GRACE. We will also discuss the impact of the 2015 El Nino on the tropical C dynamics constrained by SMAP Level 1C data. This study represents an innovative way of using satellite microwave observations to constrain C cycle in an Earth system model.

Stoichiometry of a semi-aquatic plethodontid salamander: Intraspecific variation due to location, size and diet

EPA Science Inventory

Ecological stoichiometry provides a framework to investigate an organism's relationship to nutrient cycles. An organism's stoichiometry is thought to constrain its contribution to nutrient cycles (recycling or storage), and to limit its growth and reproduction. Factors that influ...

Bulk diffusion in a kinetically constrained lattice gas

NASA Astrophysics Data System (ADS)

Arita, Chikashi; Krapivsky, P. L.; Mallick, Kirone

2018-03-01

In the hydrodynamic regime, the evolution of a stochastic lattice gas with symmetric hopping rules is described by a diffusion equation with density-dependent diffusion coefficient encapsulating all microscopic details of the dynamics. This diffusion coefficient is, in principle, determined by a Green-Kubo formula. In practice, even when the equilibrium properties of a lattice gas are analytically known, the diffusion coefficient cannot be computed except when a lattice gas additionally satisfies the gradient condition. We develop a procedure to systematically obtain analytical approximations for the diffusion coefficient for non-gradient lattice gases with known equilibrium. The method relies on a variational formula found by Varadhan and Spohn which is a version of the Green-Kubo formula particularly suitable for diffusive lattice gases. Restricting the variational formula to finite-dimensional sub-spaces allows one to perform the minimization and gives upper bounds for the diffusion coefficient. We apply this approach to a kinetically constrained non-gradient lattice gas in two dimensions, viz. to the Kob-Andersen model on the square lattice.

Inferring lateral density variations in Great Geneva Basin, western Switzerland from wells and gravity data.

NASA Astrophysics Data System (ADS)

Carrier, Aurore; Lupi, Matteo; Clerc, Nicolas; Rusillon, Elme; Do Couto, Damien

2017-04-01

In the framework of sustainable energy development Switzerland supports the growth of renewable energies. SIG (Services Industriels de Genève) and the Canton of Geneva intend to develop the use of hydrothermal energy in western Switzerland. As a Mesozoïc-formed sedimentary basin, the Great Geneva Basin (GGB) shares geological and petrophysical similarities with the Munich area (Baviera, Germany) and Paris Basin (France). The latter already provide significant amounts of geothermal energy for district heating. The prospection phase has been launched in 2014 by SIG and aims at identifying relevant geological units and defining their geometries. Lower Cretaceous and Tertiary geological units have first been targeted as potential layers. At the depth we find these units (and according to the normal geothermal gradient), low enthalpy geothermal resources are rather expected. In this framework, our study aims at constraining and refining lateral and vertical heterogeneities of Quaternary to Cretaceous sedimentary layers in GGB. Linear velocity law is inverted at wells and then interpolated to the whole basin for each geological layer. Using time pickings from available data and Quaternary information from previous studies time to depth conversion is performed. Thickness map of every geological unit is then produced. Tertiary thickness ranges from 0 m at the NW border of the GGB at the foothill of the Jura Mountains to 3000 m in the SE of the GGB at the border with the French Alps. These observations are consistent with field and well observations. The produced thickness map will be used as a geometry support for gravity data inversion and then density lateral variations estimation. Unconstrained, and a priori constrained inversion has been performed in GGB using Gauss-Newton algorithms. Velocity versus density relationships will then enable to refine velocity law interpolation. Our procedure allowed us to reduce the uncertainty of key target formation and represents an important step towards the development of geothermal energy in the Great Geneva Basin.

Analysis of Gaspra lightcurves using Galileo shape and photometric models

NASA Technical Reports Server (NTRS)

Simonelli, Damon P.; Veverka, J.; Thomas, P. C.; Helfenstein, P.; Belton, M. J. S.

1995-01-01

Galileo-based models for the shape of 951 Gaspra and the global-average photometric behavior of its surface have been used to model a representative subset of the asteroid's telescopic lightcurves. Fitting the synthetic lightcurves to the observed timing of lightcurve extrema, and knowing the orientation of Gaspra's axes at the time of the Galileo flyby, leads to a sidereal rotation period for the asteroid of 7.042024 +/- 0.000020 hr, a slight change from the period reported by Magnusson et al. (1992). Initially, the shapes, amplitudes, and absolute photometry of the synthetic and observed lightcurves agree with each other to within 0.05-0.1 mag. Small modifications to the Gaspra shape model on sides of the asteroid poorly imaged by Galileo (changes of 700 m or less in the southern hemisphere at longitudes 90 deg-270 deg W) reduce the typical discrepancies to approximately 0.05 mag in lightcurve shape and less than 0.03 mag in absolute photometry. The result demonstrates that Earth-based lightcurves can be used to refine the shape of a spacecraft-imaged irregular object in areas that are poorly constrained by the spacecraft observations. The consistency and phase-angle dependence of the Galileo-based model for Gaspra photometry, supports the accuracy of the absolute calibration of the Galileo SSI camera, and confirms the Earth-based determination of the V-filter geometric albedo of the asteroid (0.22 +/- 0.03; Tholen et al., submitted for publication). Remaining discrepancies between the synthetic and observed lightcurves show no indication of systematic latitudinal variations in albedo and also cannot be explained entirely by isolated albedo spots. These discrepancies are most likely caused by (1) small, remaining, hard-to-constrain errors in the Gaspra shape model and/or (2) moderate variations in macroscopic roughness across the asteroid's surface, in particular making longitudes 130 deg to 300 deg W moderately rougher than the opposite hemisphere.

Quantum dynamics by the constrained adiabatic trajectory method

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

Leclerc, A.; Jolicard, G.; Guerin, S.

2011-03-15

We develop the constrained adiabatic trajectory method (CATM), which allows one to solve the time-dependent Schroedinger equation constraining the dynamics to a single Floquet eigenstate, as if it were adiabatic. This constrained Floquet state (CFS) is determined from the Hamiltonian modified by an artificial time-dependent absorbing potential whose forms are derived according to the initial conditions. The main advantage of this technique for practical implementation is that the CFS is easy to determine even for large systems since its corresponding eigenvalue is well isolated from the others through its imaginary part. The properties and limitations of the CATM are exploredmore » through simple examples.« less

Coupled Long-Term Simulation of Reach-Scale Water and Heat Fluxes Across the River-Groundwater Interface for Retrieving Hyporheic Residence Times and Temperature Dynamics

NASA Astrophysics Data System (ADS)

Munz, Matthias; Oswald, Sascha E.; Schmidt, Christian

2017-11-01

Flow patterns in conjunction with seasonal and diurnal temperature variations control ecological and biogeochemical conditions in hyporheic sediments. In particular, hyporheic temperatures have a great impact on many temperature-sensitive microbial processes. In this study, we used 3-D coupled water flow and heat transport simulations applying the HydroGeoSphere code in combination with high-resolution observations of hydraulic heads and temperatures to quantify reach-scale water and heat flux across the river-groundwater interface and hyporheic temperature dynamics of a lowland gravel bed river. The model was calibrated in order to constrain estimates of the most sensitive model parameters. The magnitude and variations of the simulated temperatures matched the observed ones, with an average mean absolute error of 0.7°C and an average Nash Sutcliffe efficiency of 0.87. Our results indicate that nonsubmerged streambed structures such as gravel bars cause substantial thermal heterogeneity within the saturated sediment at the reach scale. Individual hyporheic flow path temperatures strongly depend on the flow path residence time, flow path depth, river, and groundwater temperature. Variations in individual hyporheic flow path temperatures were up to 7.9°C, significantly higher than the daily average (2.8°C), but still lower than the average seasonal hyporheic temperature difference (19.2°C). The distribution between flow path temperatures and residence times follows a power law relationship with exponent of about 0.37. Based on this empirical relation, we further estimated the influence of hyporheic flow path residence time and temperature on oxygen consumption which was found to partly increase by up to 29% in simulations.

What is the correct cost functional for variational data assimilation?

NASA Astrophysics Data System (ADS)

Bröcker, Jochen

2018-03-01

Variational approaches to data assimilation, and weakly constrained four dimensional variation (WC-4DVar) in particular, are important in the geosciences but also in other communities (often under different names). The cost functions and the resulting optimal trajectories may have a probabilistic interpretation, for instance by linking data assimilation with maximum aposteriori (MAP) estimation. This is possible in particular if the unknown trajectory is modelled as the solution of a stochastic differential equation (SDE), as is increasingly the case in weather forecasting and climate modelling. In this situation, the MAP estimator (or "most probable path" of the SDE) is obtained by minimising the Onsager-Machlup functional. Although this fact is well known, there seems to be some confusion in the literature, with the energy (or "least squares") functional sometimes been claimed to yield the most probable path. The first aim of this paper is to address this confusion and show that the energy functional does not, in general, provide the most probable path. The second aim is to discuss the implications in practice. Although the mentioned results pertain to stochastic models in continuous time, they do have consequences in practice where SDE's are approximated by discrete time schemes. It turns out that using an approximation to the SDE and calculating its most probable path does not necessarily yield a good approximation to the most probable path of the SDE proper. This suggest that even in discrete time, a version of the Onsager-Machlup functional should be used, rather than the energy functional, at least if the solution is to be interpreted as a MAP estimator.

TIME-RESOLVED SPECTROSCOPY OF THE POLAR EU CANCRI IN THE OPEN CLUSTER MESSIER 67

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

Williams, Kurtis A.; Howell, Steve B.; Liebert, James

2013-05-15

We present time-resolved spectroscopic and polarimetric observations of the AM Her system EU Cnc. EU Cnc is located near the core of the old open cluster Messier 67; new proper motion measurements indicate that EU Cnc is indeed a member of the star cluster, and this system therefore is useful to constrain the formation and evolution of magnetic cataclysmic variables. The spectra exhibit two-component emission features with independent radial velocity variations as well as time-variable cyclotron emission indicating a magnetic field strength of 41 MG. The period of the radial velocity and cyclotron hump variations are consistent with the previouslymore » known photometric period, and the spectroscopic flux variations are consistent in amplitude with previous photometric amplitude measurements. The secondary star is also detected in the spectrum. We also present polarimetric imaging measurements of EU Cnc that show a clear detection of polarization, and the degree of polarization drops below our detection threshold at phases when the cyclotron emission features are fading or not evident. The combined data are all consistent with the interpretation that EU Cnc is a low-state polar in the cluster Messier 67. The mass function of the system gives an estimate of the accretor mass of M{sub WD} {>=} 0.68 M{sub Sun} with M{sub WD} Almost-Equal-To 0.83 M{sub Sun} for an average inclination. We are thus able to place a lower limit on the progenitor mass of the accreting white dwarf of {>=}1.43 M{sub Sun }.« less

The Mid-Latitude Positive Bay and the MPB Index of Substorm Activity

NASA Astrophysics Data System (ADS)

McPherron, Robert L.; Chu, Xiangning

2017-03-01

Substorms are a major source of magnetic activity. At substorm expansion phase onset a westward current flows through the expanding aurora. This current is the ionospheric closure of the substorm current wedge produced by diversion of tail current along magnetic field lines. At low latitudes the field-aligned currents create a systematic pattern in the north (X) and east (Y) components of the surface magnetic field. The rise and decay in X is called a midlatitude positive bay whose start is a proxy for expansion onset. In this paper we describe a new index called the midlatitude positive bay index (MPB) which monitors the power in the substorm perturbations of X and Y. The index is obtained by removing the main field, storm time variations, and the solar quiet (Sq) variation from the measured field. These are estimated with spline fits and principal component analysis. The residuals of X and Y are high pass filtered to eliminate variations with period longer than 3 hours. The sum of squares of the X and Y power is determined at each of 35 midlatitude stations. The average power in night time stations is the MPB index. The index series is standardized and intervals above a fixed threshold are taken as possible bay signatures. Post processing constrains these to have reasonable values of rise time, strength, and duration. Minima in the index before and after the peak are taken as the start and end of the bay. The MPB and AL indices can be used to identify quiet intervals in the magnetic field.

Transit Timing Variation analysis with Kepler light curves of KOI 227 and Kepler 93b

NASA Astrophysics Data System (ADS)

Dulz, Shannon; Reed, Mike

2017-01-01

By searching for transit signals in approximately 150,000 stars, NASA’s Kepler Space telescope found thousands of exoplanets over its primary mission from 2009 to 2013 (Tenenbaum et al. 2014, ApJS, 211, 6). Yet, a detailed follow-up examination of Kepler light curves may contribute more evidence on system dynamics and planetary atmospheres of these objects. Kepler’s continuous observing of these systems over the mission duration produced light curves of sufficient duration to allow for the search for transit timing variations. Transit timing variations over the course of many orbits may indicate a precessing orbit or the existence of a non-transiting third body such as another exoplanet. Flux contributions of the planet just prior to secondary eclipse may provide a measurement of bond albedo from the day-side of the transiting planet. Any asymmetries of the transit shape may indicate thermal asymmetries which can measure upper atmosphere motion of the planet. These two factors can constrain atmospheric models of close orbiting exoplanets. We first establish our procedure with the well-documented TTV system, KOI 227 (Nesvorny et al. 2014, ApJ, 790, 31). Using the test case of KOI 227, we analyze Kepler-93b for TTVs and day-side flux contributions. Kepler-93b is likely a rocky planet with R = 1.50 ± 0.03 Earth Radii and M = 2.59 ± 2.0 Earth Masses (Marcy et al. 2014, ApJS, 210, 20). This research is funded by a NASA EPSCoR grant.

Finite Time Control Design for Bilateral Teleoperation System With Position Synchronization Error Constrained.

PubMed

Yang, Yana; Hua, Changchun; Guan, Xinping

2016-03-01

Due to the cognitive limitations of the human operator and lack of complete information about the remote environment, the work performance of such teleoperation systems cannot be guaranteed in most cases. However, some practical tasks conducted by the teleoperation system require high performances, such as tele-surgery needs satisfactory high speed and more precision control results to guarantee patient' health status. To obtain some satisfactory performances, the error constrained control is employed by applying the barrier Lyapunov function (BLF). With the constrained synchronization errors, some high performances, such as, high convergence speed, small overshoot, and an arbitrarily predefined small residual constrained synchronization error can be achieved simultaneously. Nevertheless, like many classical control schemes only the asymptotic/exponential convergence, i.e., the synchronization errors converge to zero as time goes infinity can be achieved with the error constrained control. It is clear that finite time convergence is more desirable. To obtain a finite-time synchronization performance, the terminal sliding mode (TSM)-based finite time control method is developed for teleoperation system with position error constrained in this paper. First, a new nonsingular fast terminal sliding mode (NFTSM) surface with new transformed synchronization errors is proposed. Second, adaptive neural network system is applied for dealing with the system uncertainties and the external disturbances. Third, the BLF is applied to prove the stability and the nonviolation of the synchronization errors constraints. Finally, some comparisons are conducted in simulation and experiment results are also presented to show the effectiveness of the proposed method.

Automated design of minimum drag light aircraft fuselages and nacelles

NASA Technical Reports Server (NTRS)

Smetana, F. O.; Fox, S. R.; Karlin, B. E.

1982-01-01

The constrained minimization algorithm of Vanderplaats is applied to the problem of designing minimum drag faired bodies such as fuselages and nacelles. Body drag is computed by a variation of the Hess-Smith code. This variation includes a boundary layer computation. The encased payload provides arbitrary geometric constraints, specified a priori by the designer, below which the fairing cannot shrink. The optimization may include engine cooling air flows entering and exhausting through specific port locations on the body.

Recent Changes in Global Photosynthesis and Terrestrial Ecosystem Respiration Constrained From Multiple Observations

NASA Astrophysics Data System (ADS)

Li, Wei; Ciais, Philippe; Wang, Yilong; Yin, Yi; Peng, Shushi; Zhu, Zaichun; Bastos, Ana; Yue, Chao; Ballantyne, Ashley P.; Broquet, Grégoire; Canadell, Josep G.; Cescatti, Alessandro; Chen, Chi; Cooper, Leila; Friedlingstein, Pierre; Le Quéré, Corinne; Myneni, Ranga B.; Piao, Shilong

2018-01-01

To assess global carbon cycle variability, we decompose the net land carbon sink into the sum of gross primary productivity (GPP), terrestrial ecosystem respiration (TER), and fire emissions and apply a Bayesian framework to constrain these fluxes between 1980 and 2014. The constrained GPP and TER fluxes show an increasing trend of only half of the prior trend simulated by models. From the optimization, we infer that TER increased in parallel with GPP from 1980 to 1990, but then stalled during the cooler periods, in 1990-1994 coincident with the Pinatubo eruption, and during the recent warming hiatus period. After each of these TER stalling periods, TER is found to increase faster than GPP, explaining a relative reduction of the net land sink. These results shed light on decadal variations of GPP and TER and suggest that they exhibit different responses to temperature anomalies over the last 35 years.

Multi-Sensor Constrained Time Varying Emissions Estimation of Black Carbon: Attributing Urban and Fire Sources Globally

NASA Astrophysics Data System (ADS)

Cohen, J. B.

2015-12-01

The short lifetime and heterogeneous distribution of Black Carbon (BC) in the atmosphere leads to complex impacts on radiative forcing, climate, and health, and complicates analysis of its atmospheric processing and emissions. Two recent papers have estimated the global and regional emissions of BC using advanced statistical and computational methods. One used a Kalman Filter, including data from AERONET, NOAA, and other ground-based sources, to estimate global emissions of 17.8+/-5.6 Tg BC/year (with the increase attributable to East Asia, South Asia, Southeast Asia, and Eastern Europe - all regions which have had rapid urban, industrial, and economic expansion). The second additionally used remotely sensed measurements from MISR and a variance maximizing technique, uniquely quantifying fire and urban sources in Southeast Asia, as well as their large year-to-year variability over the past 12 years, leading to increases from 10% to 150%. These new emissions products, when run through our state-of-the art modelling system of chemistry, physics, transport, removal, radiation, and climate, match 140 ground stations and satellites better in both an absolute and a temporal sense. New work now further includes trace species measurements from OMI, which are used with the variance maximizing technique to constrain the types of emissions sources. Furthermore, land-use change and fire estimation products from MODIS are also included, which provide other constraints on the temporal and spatial nature of the variations of intermittent sources like fires or new permanent sources like expanded urbanization. This talk will introduce a new, top-down constrained, weekly varying BC emissions dataset, show that it produces a better fit with observations, and draw conclusions about the sources and impacts from urbanization one hand, and fires on another hand. Results specific to the Southeast and East Asia will demonstrate inter- and intra-annual variations, such as the function of the wet and dry seasons. Further, the impacts of missing data due to cloud coverage and of long-range transport from highly polluted areas to relatively clean downwind areas will be demonstrated. More general results will also be discussed in relation to the global anthropogenic aerosol distribution.

Constraining the models' response of tropical low clouds to SST forcings using CALIPSO observations

NASA Astrophysics Data System (ADS)

Cesana, G.; Del Genio, A. D.; Ackerman, A. S.; Brient, F.; Fridlind, A. M.; Kelley, M.; Elsaesser, G.

2017-12-01

Low-cloud response to a warmer climate is still pointed out as being the largest source of uncertainty in the last generation of climate models. To date there is no consensus among the models on whether the tropical low cloudiness would increase or decrease in a warmer climate. In addition, it has been shown that - depending on their climate sensitivity - the models either predict deeper or shallower low clouds. Recently, several relationships between inter-model characteristics of the present-day climate and future climate changes have been highlighted. These so-called emergent constraints aim to target relevant model improvements and to constrain models' projections based on current climate observations. Here we propose to use - for the first time - 10 years of CALIPSO cloud statistics to assess the ability of the models to represent the vertical structure of tropical low clouds for abnormally warm SST. We use a simulator approach to compare observations and simulations and focus on the low-layered clouds (i.e. z < 3.2km) as well the more detailed level perspective of clouds (40 levels from 0 to 19km). Results show that in most models an increase of the SST leads to a decrease of the low-layer cloud fraction. Vertically, the clouds deepen namely by decreasing the cloud fraction in the lowest levels and increasing it around the top of the boundary-layer. This feature is coincident with an increase of the high-level cloud fraction (z > 6.5km). Although the models' spread is large, the multi-model mean captures the observed variations but with a smaller amplitude. We then employ the GISS model to investigate how changes in cloud parameterizations affect the response of low clouds to warmer SSTs on the one hand; and how they affect the variations of the model's cloud profiles with respect to environmental parameters on the other hand. Finally, we use CALIPSO observations to constrain the model by determining i) what set of parameters allows reproducing the observed relationships and ii) what are the consequences on the cloud feedbacks. These results point toward process-oriented constraints of low-cloud responses to surface warming and environmental parameters.

Direct and indirect selection on flowering time, water-use efficiency (WUE, δ 13C), and WUE plasticity to drought in Arabidopsis thaliana

PubMed Central

Kenney, Amanda M; McKay, John K; Richards, James H; Juenger, Thomas E

2014-01-01

Flowering time and water-use efficiency (WUE) are two ecological traits that are important for plant drought response. To understand the evolutionary significance of natural genetic variation in flowering time, WUE, and WUE plasticity to drought in Arabidopsis thaliana, we addressed the following questions: (1) How are ecophysiological traits genetically correlated within and between different soil moisture environments? (2) Does terminal drought select for early flowering and drought escape? (3) Is WUE plasticity to drought adaptive and/or costly? We measured a suite of ecophysiological and reproductive traits on 234 spring flowering accessions of A. thaliana grown in well-watered and season-ending soil drying treatments, and quantified patterns of genetic variation, correlation, and selection within each treatment. WUE and flowering time were consistently positively genetically correlated. WUE was correlated with WUE plasticity, but the direction changed between treatments. Selection generally favored early flowering and low WUE, with drought favoring earlier flowering significantly more than well-watered conditions. Selection for lower WUE was marginally stronger under drought. There were no net fitness costs of WUE plasticity. WUE plasticity (per se) was globally neutral, but locally favored under drought. Strong genetic correlation between WUE and flowering time may facilitate the evolution of drought escape, or constrain independent evolution of these traits. Terminal drought favored drought escape in these spring flowering accessions of A. thaliana. WUE plasticity may be favored over completely fixed development in environments with periodic drought. PMID:25512847

Secular Orbit and Spin Variations of Asteroid (16) Psyche

NASA Astrophysics Data System (ADS)

Bills, B. G.; Park, R. S.; Scott, B.

2016-12-01

The obliquity, or angular separation between spin and orbit poles, of asteroid (16) Psyche is currently 95 degrees. We are interested in knowing how much that angular separation varies, on time scales of 104 to 106 years. To answer that question, we have done several related analyses. On short time scales, the orbital element variations of Psyche are dominated by perturbations from Jupiter. Jupiter's dominance has two basic causes: first is the large mass and relatively close position of Jupiter, and second is a 19:8 mean motion resonance. Jupiter completes 8 orbits in 94.9009 years, while Psyche takes 94.9107 years to complete 19 orbits. As a result of this, all of the orbital elements of Psyche exhibit significant periodic variations, with a 94.9 year period dominating. There are also significant variations at the synodic period, which is 8.628 years, or 1/11 of the resonant period. Over a 1000 year time span, centered on the present, the eccentricity varies from 0.133 to 0.140, and the inclination varies from 2.961 to 3.229 degrees. On longer time scales, the orbital elements of Psyche vary considerably more than that, due to secular perturbations from the planets. The secular variations are modeled as the response of interacting mass rings, rather than point masses. Again, Jupiter is the main perturbing influence on Psyche. The eccentricity and inclination both oscillate, with dominant periods of 18.667 kyr. The range of values seen over a million year time span, is 0.057 to 0.147 for eccentricity, and 0.384 to 4.777 degrees for inclination. Using a recent shape model, and assumption of uniform density, to constrain relevant moments of inertia, we estimate the spin pole precession rate parameter to be 8.53 arcsec/year. The current spin pole is at ecliptic {lon, lat} = { 32, -7} deg, whereas the orbit pole is at {lon, lat} = {60.47, 86.91} deg. The current obliquity is thus 94.3 degree. Using nominal values of the input parameters, the recovered spin pole trajectory is such that, over a million year time span, centered on the present, the minimum and maximum values of obliquity are 92.36 and 98.56 deg. The obliquity oscillates with dominant periods of 18.45 and 48.40 kyr.

Variation in handshape and orientation in British Sign Language: The case of the ‘1’ hand configuration

PubMed Central

Fenlon, Jordan; Schembri, Adam; Rentelis, Ramas; Cormier, Kearsy

2013-01-01

This paper investigates phonological variation in British Sign Language (BSL) signs produced with a ‘1’ hand configuration in citation form. Multivariate analyses of 2084 tokens reveals that handshape variation in these signs is constrained by linguistic factors (e.g., the preceding and following phonological environment, grammatical category, indexicality, lexical frequency). The only significant social factor was region. For the subset of signs where orientation was also investigated, only grammatical function was important (the surrounding phonological environment and social factors were not significant). The implications for an understanding of pointing signs in signed languages are discussed. PMID:23805018

Working Memory in Children: A Time-Constrained Functioning Similar to Adults

ERIC Educational Resources Information Center

Portrat, Sophie; Camos, Valerie; Barrouillet, Pierre

2009-01-01

Within the time-based resource-sharing (TBRS) model, we tested a new conception of the relationships between processing and storage in which the core mechanisms of working memory (WM) are time constrained. However, our previous studies were restricted to adults. The current study aimed at demonstrating that these mechanisms are present and…

Numerical methods for the inverse problem of density functional theory

DOE PAGES

Jensen, Daniel S.; Wasserman, Adam

2017-07-17

Here, the inverse problem of Kohn–Sham density functional theory (DFT) is often solved in an effort to benchmark and design approximate exchange-correlation potentials. The forward and inverse problems of DFT rely on the same equations but the numerical methods for solving each problem are substantially different. We examine both problems in this tutorial with a special emphasis on the algorithms and error analysis needed for solving the inverse problem. Two inversion methods based on partial differential equation constrained optimization and constrained variational ideas are introduced. We compare and contrast several different inversion methods applied to one-dimensional finite and periodic modelmore » systems.« less

Numerical methods for the inverse problem of density functional theory

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

Jensen, Daniel S.; Wasserman, Adam

Here, the inverse problem of Kohn–Sham density functional theory (DFT) is often solved in an effort to benchmark and design approximate exchange-correlation potentials. The forward and inverse problems of DFT rely on the same equations but the numerical methods for solving each problem are substantially different. We examine both problems in this tutorial with a special emphasis on the algorithms and error analysis needed for solving the inverse problem. Two inversion methods based on partial differential equation constrained optimization and constrained variational ideas are introduced. We compare and contrast several different inversion methods applied to one-dimensional finite and periodic modelmore » systems.« less

On the linear relation between the mean and the standard deviation of a response time distribution.

PubMed

Wagenmakers, Eric-Jan; Brown, Scott

2007-07-01

Although it is generally accepted that the spread of a response time (RT) distribution increases with the mean, the precise nature of this relation remains relatively unexplored. The authors show that in several descriptive RT distributions, the standard deviation increases linearly with the mean. Results from a wide range of tasks from different experimental paradigms support a linear relation between RT mean and RT standard deviation. Both R. Ratcliff's (1978) diffusion model and G. D. Logan's (1988) instance theory of automatization provide explanations for this linear relation. The authors identify and discuss 3 specific boundary conditions for the linear law to hold. The law constrains RT models and supports the use of the coefficient of variation to (a) compare variability while controlling for differences in baseline speed of processing and (b) assess whether changes in performance with practice are due to quantitative speedup or qualitative reorganization. Copyright 2007 APA.

Comment on “Symplectic integration of magnetic systems”: A proof that the Boris algorithm is not variational

DOE PAGES

Ellison, C. L.; Burby, J. W.; Qin, H.

2015-11-01

One popular technique for the numerical time advance of charged particles interacting with electric and magnetic fields according to the Lorentz force law [1], [2], [3] and [4] is the Boris algorithm. Its popularity stems from simple implementation, rapid iteration, and excellent long-term numerical fidelity [1] and [5]. Excellent long-term behavior strongly suggests the numerical dynamics exhibit conservation laws analogous to those governing the continuous Lorentz force system [6]. Moreover, without conserved quantities to constrain the numerical dynamics, algorithms typically dissipate or accumulate important observables such as energy and momentum over long periods of simulated time [6]. Identification of themore » conservative properties of an algorithm is important for establishing rigorous expectations on the long-term behavior; energy-preserving, symplectic, and volume-preserving methods each have particular implications for the qualitative numerical behavior [6], [7], [8], [9], [10] and [11].« less

On optimal infinite impulse response edge detection filters

NASA Technical Reports Server (NTRS)

Sarkar, Sudeep; Boyer, Kim L.

1991-01-01

The authors outline the design of an optimal, computationally efficient, infinite impulse response edge detection filter. The optimal filter is computed based on Canny's high signal to noise ratio, good localization criteria, and a criterion on the spurious response of the filter to noise. An expression for the width of the filter, which is appropriate for infinite-length filters, is incorporated directly in the expression for spurious responses. The three criteria are maximized using the variational method and nonlinear constrained optimization. The optimal filter parameters are tabulated for various values of the filter performance criteria. A complete methodology for implementing the optimal filter using approximating recursive digital filtering is presented. The approximating recursive digital filter is separable into two linear filters operating in two orthogonal directions. The implementation is very simple and computationally efficient, has a constant time of execution for different sizes of the operator, and is readily amenable to real-time hardware implementation.

A chance-constrained stochastic approach to intermodal container routing problems.

PubMed

Zhao, Yi; Liu, Ronghui; Zhang, Xi; Whiteing, Anthony

2018-01-01

We consider a container routing problem with stochastic time variables in a sea-rail intermodal transportation system. The problem is formulated as a binary integer chance-constrained programming model including stochastic travel times and stochastic transfer time, with the objective of minimising the expected total cost. Two chance constraints are proposed to ensure that the container service satisfies ship fulfilment and cargo on-time delivery with pre-specified probabilities. A hybrid heuristic algorithm is employed to solve the binary integer chance-constrained programming model. Two case studies are conducted to demonstrate the feasibility of the proposed model and to analyse the impact of stochastic variables and chance-constraints on the optimal solution and total cost.

A chance-constrained stochastic approach to intermodal container routing problems

PubMed Central

Zhao, Yi; Zhang, Xi; Whiteing, Anthony

2018-01-01

We consider a container routing problem with stochastic time variables in a sea-rail intermodal transportation system. The problem is formulated as a binary integer chance-constrained programming model including stochastic travel times and stochastic transfer time, with the objective of minimising the expected total cost. Two chance constraints are proposed to ensure that the container service satisfies ship fulfilment and cargo on-time delivery with pre-specified probabilities. A hybrid heuristic algorithm is employed to solve the binary integer chance-constrained programming model. Two case studies are conducted to demonstrate the feasibility of the proposed model and to analyse the impact of stochastic variables and chance-constraints on the optimal solution and total cost. PMID:29438389

Dark energy equation of state parameter and its evolution at low redshift

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

Tripathi, Ashutosh; Sangwan, Archana; Jassal, H.K., E-mail: ashutosh_tripathi@fudan.edu.cn, E-mail: archanakumari@iisermohali.ac.in, E-mail: hkjassal@iisermohali.ac.in

In this paper, we constrain dark energy models using a compendium of observations at low redshifts. We consider the dark energy as a barotropic fluid, with the equation of state a constant as well the case where dark energy equation of state is a function of time. The observations considered here are Supernova Type Ia data, Baryon Acoustic Oscillation data and Hubble parameter measurements. We compare constraints obtained from these data and also do a combined analysis. The combined observational constraints put strong limits on variation of dark energy density with redshift. For varying dark energy models, the range ofmore » parameters preferred by the supernova type Ia data is in tension with the other low redshift distance measurements.« less

An ensemble constrained variational analysis of atmospheric forcing data and its application to evaluate clouds in CAM5: Ensemble 3DCVA and Its Application

DOE PAGES

Tang, Shuaiqi; Zhang, Minghua; Xie, Shaocheng

2016-01-05

Large-scale atmospheric forcing data can greatly impact the simulations of atmospheric process models including Large Eddy Simulations (LES), Cloud Resolving Models (CRMs) and Single-Column Models (SCMs), and impact the development of physical parameterizations in global climate models. This study describes the development of an ensemble variationally constrained objective analysis of atmospheric large-scale forcing data and its application to evaluate the cloud biases in the Community Atmospheric Model (CAM5). Sensitivities of the variational objective analysis to background data, error covariance matrix and constraint variables are described and used to quantify the uncertainties in the large-scale forcing data. Application of the ensemblemore » forcing in the CAM5 SCM during March 2000 intensive operational period (IOP) at the Southern Great Plains (SGP) of the Atmospheric Radiation Measurement (ARM) program shows systematic biases in the model simulations that cannot be explained by the uncertainty of large-scale forcing data, which points to the deficiencies of physical parameterizations. The SCM is shown to overestimate high clouds and underestimate low clouds. These biases are found to also exist in the global simulation of CAM5 when it is compared with satellite data.« less

An ensemble constrained variational analysis of atmospheric forcing data and its application to evaluate clouds in CAM5: Ensemble 3DCVA and Its Application

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

Tang, Shuaiqi; Zhang, Minghua; Xie, Shaocheng

Large-scale atmospheric forcing data can greatly impact the simulations of atmospheric process models including Large Eddy Simulations (LES), Cloud Resolving Models (CRMs) and Single-Column Models (SCMs), and impact the development of physical parameterizations in global climate models. This study describes the development of an ensemble variationally constrained objective analysis of atmospheric large-scale forcing data and its application to evaluate the cloud biases in the Community Atmospheric Model (CAM5). Sensitivities of the variational objective analysis to background data, error covariance matrix and constraint variables are described and used to quantify the uncertainties in the large-scale forcing data. Application of the ensemblemore » forcing in the CAM5 SCM during March 2000 intensive operational period (IOP) at the Southern Great Plains (SGP) of the Atmospheric Radiation Measurement (ARM) program shows systematic biases in the model simulations that cannot be explained by the uncertainty of large-scale forcing data, which points to the deficiencies of physical parameterizations. The SCM is shown to overestimate high clouds and underestimate low clouds. These biases are found to also exist in the global simulation of CAM5 when it is compared with satellite data.« less

Temporal Variations in 234U/238U Activity Ratios in the Lower Mississippi River due to Changes in Source Tributary Discharges

NASA Astrophysics Data System (ADS)

Grzymko, T. J.; Marcantonio, F.; McKee, B. A.; Stewart, C. M.

2004-12-01

The world's 25 largest river systems contribute nearly 50% of all freshwater to the global ocean and carry large quantities of dissolved trace metals annually. Trace metal concentrations in these systems show large variances on seasonal time scales. In order to constrain the causes of these variations, consistent sampling on sub-seasonal time intervals is essential. Here, we focus on the Mississippi River, the seventh largest river in the world in terms of freshwater discharge and the third largest in terms of drainage basin area. Biweekly sampling of the lower Mississippi River at New Orleans was performed from January 2003 to August 2004. Uranium concentrations and 234U/238U activity ratios were measured for the dissolved component (<0.2 μ m-fraction) of river water. Over the course of this study, dissolved U activity ratios spanned a range of about 25%, from 1.23 to 1.60. Dissolved U concentrations ranged from 0.28 to 1.06 ppb. The relationship between concentrations, activity ratios, and lower river discharge is complicated, and no clear pattern is observed on both biweekly and seasonal timescales. However, there does seem to be a relationship between the larger seasonal trends in the lower Mississippi River and variations in the discharge of its upstream tributaries. To constrain this relationship, we have sampled water from the Missouri River, the upper Mississippi River above the confluence with the Missouri, the Ohio River, and the Arkansas River in February, April, and August of 2004. For the upstream samples measured thus far, the highest dissolved uranium concentrations are observed for the Missouri River at 2.02 ppb, while the lowest are found in the Ohio River at 0.38 ppb. Dissolved 234U/238U activity ratios are as unique for each tributary and vary from 1.36 in the Ohio River to 1.51 in the Missouri River. A preliminary mass balance analysis reveals that the lower river uranium activity ratios are controlled simply by the quantity and isotope signature of the waters discharged from the upstream tributaries. A discussion of the implications of this work for global ocean budgets of uranium will be presented.

Evaluating the Utility of Adjoint-based Inverse Modeling with Aircraft and Surface Measurements during ARCTAS-CARB to Constrain Wildfire Emissions of Black Carbon

NASA Astrophysics Data System (ADS)

Henze, D. K.; Guerrette, J.; Bousserez, N.

2016-12-01

Wildfires contribute significantly to regional haze events globally, and they are potentially becoming more commonplace with increasing droughts due to climate change. Aerosol emissions from wildfires are highly uncertain, with global annual totals varying by a factor of 2 to 3 and regional rates varying by up to a factor of 10. At the high resolution required to predict PM2.5 exposure events, this variance is attributable to differences in methodology, differing land cover datasets, spatial variation in fire locations, and limited understanding of fast transient fire behavior. Here we apply an adjoint-based online chemical inverse modeling tool, WRFDA-Chem, to constrain black carbon aerosol (BC) emissions from fires during the 2008 ARCTAS-CARB field campaign. We identify several weaknesses in the prior diurnal distribution of emissions, including a missing early morning emission peak associated with local, persistent, large-scale forest fires. On 22 June, 2008, aircraft observations are able to reduce the spread between FINNv1.0 and QFEDv2.4r8 from ×3.5 to ×2.1. On 23 and 24 June, the spread is reduced from ×3.4 to ×1.4. Using posterior error estimates, we found that emission variance improvements are limited to a small footprint surrounding the measurements. Relative BB emission variances are reduced by up to 35% near aircraft flight paths and up to 60% near IMPROVE surface sites. Due to the spatial variation of observations on multiple days, and the heterogeneous biomass burning errors on daily scales, cross-validation was not successful. Future high-resolution measurements need to be carefully planned to characterize biomass burning emission errors and control for day-to-day variation. In general, the 4D-Var inversion framework would benefit from reduced wall-time. For the problem presented, incremental 4D-Var requires 20 hours on 96 cores to reach practical optimization convergence and generate the posterior covariance matrix for a 24-hour assimilation window. We will present initial computational comparisons with a recently developed method to parallelize those calculations, which will reduce wall-time by a factor of 5 or more for all WRFDA 4D-Var applications.

Spatiotemporal variations of radar glacier zones in the Karakoram Mountains

NASA Astrophysics Data System (ADS)

Lund, Jewell

2017-04-01

Glaciers of the Karakoram Mountains are important climate indicators for densely populated South Central Asia. Glacial meltwater is a significant source of runoff in the Indus River Basin, upon which 60 million people rely for food security, economy and hydropower in Pakistan and India. Contrary to worldwide and also Himalayan trends, Karakoram glaciers have recently been verified in near balance, with some glaciers even gaining mass or surging. This 'Karakoram anomaly' is of interest, and many hypotheses exist for its causes. Complex climatology, coupled with the challenges of field study in this region, illicit notable uncertainties both in observation and prediction of glacial status. Constraining spatiotemporal variations in glacial mass balance will elucidate the extent and possible longevity of this anomaly, and its implications for water resources, as climate continues to change. Depending on snowpack conditions during image acquisition, different snow and ice zones on a glacier are identifiable in synthetic aperture radar (SAR) images. The identification and monitoring of radar glacier zones over time can provide indicators of relative glacial mass balance to compliment field studies in a region with sparse field measurement. We will present spatiotemporal evolution of basic radar glacier zones such as wet snow, bare ice, percolation, and firn for glaciers feeding into the Upper Indus Basin. We will incorporate both ascending and descending passes of Sentinel-1 series C -band sensors, and possibly ALOS-2 PALSAR-2 L-band images. We may also explore the impacts of these variations on timing and intensity of runoff.

Phonologically Driven Variability: The Case of Determiners

ERIC Educational Resources Information Center

Bürki, Audrey; Laganaro, Marina; Alario, F.-Xavier

2014-01-01

Speakers usually produce words in connected speech. In such contexts, the form in which many words are uttered is influenced by the phonological properties of neighboring words. The current article examines the representations and processes underlying the production of phonologically constrained word form variations. For this purpose, we consider…

Watershed Land Use and Seasonal Variation Constrain the Influence of Riparian Canopy Cover on Stream Ecosystem Metabolism

EPA Science Inventory

While watershed and local scale controls on stream metabolism have been independently investigated, little is known about how controls exerted at these different scales interact to determine stream metabolic rates, or how these interactions vary across seasons. To address this ...

Comparison of ocean mass content change from direct and inversion based approaches

NASA Astrophysics Data System (ADS)

Uebbing, Bernd; Kusche, Jürgen; Rietbroek, Roelof

2017-04-01

The GRACE satellite mission provides an indispensable tool for measuring oceanic mass variations. Such time series are essential to separate global mean sea level rise in thermosteric and mass driven contributions, and thus to constrain ocean heat content and (deep) ocean warming when viewed together with altimetry and Argo data. However, published estimates over the GRACE era differ, not only depending on the time window considered. Here, we will look into sources of such differences with direct and inverse approaches. Deriving ocean mass time series requires several processing steps; choosing a GRACE (and altimetry and Argo) product, data coverage, masks and filters to be applied in either spatial or spectral domain, corrections related to spatial leakage, GIA and geocenter motion need to be accounted for. In this study, we quantify the effects of individual processing choices and assumptions of the direct and inversion based approaches to derive ocean mass content change. Furthermore, we compile the different estimates from existing literature and sources, to highlight the differences.

Simultaneous travel time tomography for updating both velocity and reflector geometry in triangular/tetrahedral cell model

NASA Astrophysics Data System (ADS)

Bai, Chao-ying; He, Lei-yu; Li, Xing-wang; Sun, Jia-yu

2018-05-01

To conduct forward and simultaneous inversion in a complex geological model, including an irregular topography (or irregular reflector or velocity anomaly), we in this paper combined our previous multiphase arrival tracking method (referred as triangular shortest-path method, TSPM) in triangular (2D) or tetrahedral (3D) cell model and a linearized inversion solver (referred to as damped minimum norms and constrained least squares problem solved using the conjugate gradient method, DMNCLS-CG) to formulate a simultaneous travel time inversion method for updating both velocity and reflector geometry by using multiphase arrival times. In the triangular/tetrahedral cells, we deduced the partial derivative of velocity variation with respective to the depth change of reflector. The numerical simulation results show that the computational accuracy can be tuned to a high precision in forward modeling and the irregular velocity anomaly and reflector geometry can be accurately captured in the simultaneous inversion, because the triangular/tetrahedral cell can be easily used to stitch the irregular topography or subsurface interface.

Simultaneous travel time tomography for updating both velocity and reflector geometry in triangular/tetrahedral cell model

NASA Astrophysics Data System (ADS)

Bai, Chao-ying; He, Lei-yu; Li, Xing-wang; Sun, Jia-yu

2017-12-01

To conduct forward and simultaneous inversion in a complex geological model, including an irregular topography (or irregular reflector or velocity anomaly), we in this paper combined our previous multiphase arrival tracking method (referred as triangular shortest-path method, TSPM) in triangular (2D) or tetrahedral (3D) cell model and a linearized inversion solver (referred to as damped minimum norms and constrained least squares problem solved using the conjugate gradient method, DMNCLS-CG) to formulate a simultaneous travel time inversion method for updating both velocity and reflector geometry by using multiphase arrival times. In the triangular/tetrahedral cells, we deduced the partial derivative of velocity variation with respective to the depth change of reflector. The numerical simulation results show that the computational accuracy can be tuned to a high precision in forward modeling and the irregular velocity anomaly and reflector geometry can be accurately captured in the simultaneous inversion, because the triangular/tetrahedral cell can be easily used to stitch the irregular topography or subsurface interface.

Distributed Soil Moisture Estimation in a Mountainous Semiarid Basin: Constraining Soil Parameter Uncertainty through Field Studies

NASA Astrophysics Data System (ADS)

Yatheendradas, S.; Vivoni, E.

2007-12-01

A common practice in distributed hydrological modeling is to assign soil hydraulic properties based on coarse textural datasets. For semiarid regions with poor soil information, the performance of a model can be severely constrained due to the high model sensitivity to near-surface soil characteristics. Neglecting the uncertainty in soil hydraulic properties, their spatial variation and their naturally-occurring horizonation can potentially affect the modeled hydrological response. In this study, we investigate such effects using the TIN-based Real-time Integrated Basin Simulator (tRIBS) applied to the mid-sized (100 km2) Sierra Los Locos watershed in northern Sonora, Mexico. The Sierra Los Locos basin is characterized by complex mountainous terrain leading to topographic organization of soil characteristics and ecosystem distributions. We focus on simulations during the 2004 North American Monsoon Experiment (NAME) when intensive soil moisture measurements and aircraft- based soil moisture retrievals are available in the basin. Our experiments focus on soil moisture comparisons at the point, topographic transect and basin scales using a range of different soil characterizations. We compare the distributed soil moisture estimates obtained using (1) a deterministic simulation based on soil texture from coarse soil maps, (2) a set of ensemble simulations that capture soil parameter uncertainty and their spatial distribution, and (3) a set of simulations that conditions the ensemble on recent soil profile measurements. Uncertainties considered in near-surface soil characterization provide insights into their influence on the modeled uncertainty, into the value of soil profile observations, and into effective use of on-going field observations for constraining the soil moisture response uncertainty.

Combined constraints on the structure and physical properties of the East Antarctic lithosphere from geology and geophysics.

NASA Astrophysics Data System (ADS)

Reading, A. M.; Staal, T.; Halpin, J.; Whittaker, J. M.; Morse, P. E.

2017-12-01

The lithosphere of East Antarctica is one of the least explored regions of the planet, yet it is gaining in importance in global scientific research. Continental heat flux density and 3D glacial isostatic adjustment studies, for example, rely on a good knowledge of the deep structure in constraining model inputs.In this contribution, we use a multidisciplinary approach to constrain lithospheric domains. To seismic tomography models, we add constraints from magnetic studies and also new geological constraints. Geological knowledge exists around the periphery of East Antarctica and is reinforced in the knowledge of plate tectonic reconstructions. The subglacial geology of the Antarctic hinterland is largely unknown but the plate reconstructions allow the well-posed extrapolation of major terranes into the interior of the continent, guided by the seismic tomography and magnetic images. We find that the northern boundary of the lithospheric domain centred on the Gamburtsev Subglacial Mountains has a possible trend that runs south of the Lambert Glacier region, turning coastward through Wilkes Land. Other periphery-to-interior connections are less well constrained and the possibility of lithospheric domains that are entirely sub-glacial is high. We develop this framework to include a probabilistic method of handling alternate models and quantifiable uncertainties. We also show first results in using a Bayesian approach to predicting lithospheric boundaries from multivariate data.Within the newly constrained domains, we constrain heat flux (density) as the sum of basal heat flux and upper crustal heat flux. The basal heat flux is constrained by geophysical methods while the upper crustal heat flux is constrained by geology or predicted geology. In addition to heat flux constraints, we also consider the variations in friction experienced by moving ice sheets due to varying geology.

Historical trace metal accumulation in the sediments of an urbanized region of the Lake Champlain watershed, Burlington, Vermont

USGS Publications Warehouse

Mecray, E.L.; King, J.W.; Appleby, P.G.; Hunt, A.S.

2001-01-01

This study documents the history of pollution inputs in the Burlington region of Lake Champlain, Vermont using measurements of anthropogenic metals (Cu, Zn, Cr, Pb, Cd, and Ag) in four age-dated sediment cores. Sediments record a history of contamination in a region and can be used to assess the changing threat to biota over time and to evaluate the effectiveness of discharge regulations on anthropogenic inputs. Grain size, magnetic susceptibility, radiometric dating and pollen stratigraphy were combined with trace metal data to provide an assessment of the history of contamination over the last 350 yr in the Burlington region of Lake Champlain. Magnetic susceptibility was initially used to identify land-use history for each site because it is a proxy indicator of soil erosion. Historical trends in metal inputs in the Burlington region from the seventeenth through the twentieth centuries are reflected in downcore variations in metal concentrations and accumulation rates. Metal concentrations increase above background values in the early to mid nineteenth century. The metal input rate to the sediments increases around 1920 and maximum concentrations and accumulation rates are observed in the late 1960s. Decreases in concentration and accumulation rate between 1970 and the present are observed, for most metals. The observed trends are primarily a function of variations in anthropogenic inputs and not variations in sediment grain size. Grain size data were used to remove texture variations from the metal profiles and results show trends in the anthropogenic metal signals remain. Radiometric dating and pollen stratigraphy provide well-constrained dates for the sediments thereby allowing the metal profiles to be interpreted in terms of land-use history.This study documents the history of pollution inputs in the Burlington region of Lake Champlain, Vermont using measurements of anthropogenic metals (Cu, Zn, Cr, Pb, Cd, and Ag) in four age-dated sediment cores. Sediments record a history of contamination in a region and can be used to assess the changing threat to biota over time and to evaluate the effectiveness of discharge regulations on anthropogenic inputs. Grain size, magnetic susceptibility, radiometric dating and pollen stratigraphy were combined with trace metal data to provide an assessment of the history of contamination over the last 350 yr in the Burlington region of Lake Champlain. Magnetic susceptibility was initially used to identify land-use history for each site because it is a proxy indicator of soil erosion. Historical trends in metal inputs in the Burlington region from the seventeenth through the twentieth centuries are reflected in downcore variations in metal concentrations and accumulation rates. Metal concentrations increase above background values in the early to mid nineteenth century. The metal input rate to the sediments increases around 1920 and maximum concentrations and accumulation rates are observed in the late 1960s. Decreases in concentration and accumulation rate between 1970 and the present are observed for most metals. The observed trends are primarily a function of variations in anthropogenic inputs and not variations in sediment grain size. Grain size data were used to remove texture variations from the metal profiles and results show trends in the anthropogenic metal signals remain. Radiometric dating and pollen stratigraphy provide well-constrained dates for the sediments thereby allowing the metal profiles to be interpreted in terms of land-use history.

The role of residence time in diagnostic models of global carbon storage capacity: model decomposition based on a traceable scheme.

PubMed

Yizhao, Chen; Jianyang, Xia; Zhengguo, Sun; Jianlong, Li; Yiqi, Luo; Chengcheng, Gang; Zhaoqi, Wang

2015-11-06

As a key factor that determines carbon storage capacity, residence time (τE) is not well constrained in terrestrial biosphere models. This factor is recognized as an important source of model uncertainty. In this study, to understand how τE influences terrestrial carbon storage prediction in diagnostic models, we introduced a model decomposition scheme in the Boreal Ecosystem Productivity Simulator (BEPS) and then compared it with a prognostic model. The result showed that τE ranged from 32.7 to 158.2 years. The baseline residence time (τ'E) was stable for each biome, ranging from 12 to 53.7 years for forest biomes and 4.2 to 5.3 years for non-forest biomes. The spatiotemporal variations in τE were mainly determined by the environmental scalar (ξ). By comparing models, we found that the BEPS uses a more detailed pool construction but rougher parameterization for carbon allocation and decomposition. With respect to ξ comparison, the global difference in the temperature scalar (ξt) averaged 0.045, whereas the moisture scalar (ξw) had a much larger variation, with an average of 0.312. We propose that further evaluations and improvements in τ'E and ξw predictions are essential to reduce the uncertainties in predicting carbon storage by the BEPS and similar diagnostic models.

The role of residence time in diagnostic models of global carbon storage capacity: model decomposition based on a traceable scheme

PubMed Central

Yizhao, Chen; Jianyang, Xia; Zhengguo, Sun; Jianlong, Li; Yiqi, Luo; Chengcheng, Gang; Zhaoqi, Wang

2015-01-01

As a key factor that determines carbon storage capacity, residence time (τE) is not well constrained in terrestrial biosphere models. This factor is recognized as an important source of model uncertainty. In this study, to understand how τE influences terrestrial carbon storage prediction in diagnostic models, we introduced a model decomposition scheme in the Boreal Ecosystem Productivity Simulator (BEPS) and then compared it with a prognostic model. The result showed that τE ranged from 32.7 to 158.2 years. The baseline residence time (τ′E) was stable for each biome, ranging from 12 to 53.7 years for forest biomes and 4.2 to 5.3 years for non-forest biomes. The spatiotemporal variations in τE were mainly determined by the environmental scalar (ξ). By comparing models, we found that the BEPS uses a more detailed pool construction but rougher parameterization for carbon allocation and decomposition. With respect to ξ comparison, the global difference in the temperature scalar (ξt) averaged 0.045, whereas the moisture scalar (ξw) had a much larger variation, with an average of 0.312. We propose that further evaluations and improvements in τ′E and ξw predictions are essential to reduce the uncertainties in predicting carbon storage by the BEPS and similar diagnostic models. PMID:26541245

The Evolution of Quasar C IV and Si IV Broad Absorption Lines over Multi-year Timescales

NASA Astrophysics Data System (ADS)

Gibson, Robert R.; Brandt, W. N.; Gallagher, S. C.; Hewett, Paul C.; Schneider, Donald P.

2010-04-01

We investigate the variability of C IV λ1549 broad absorption line (BAL) troughs over rest-frame timescales of up to ≈7 yr in 14 quasars at redshifts z >~ 2.1. For nine sources at sufficiently high redshift, we also compare the C IV and Si IV λ1400 absorption variation. We compare shorter and longer term variability using spectra from up to four different epochs per source and find complex patterns of variation in the sample overall. The scatter in the change of absorption equivalent width (EW), ΔEW, increases with the time between observations. BALs do not, in general, strengthen or weaken monotonically, and variation observed over shorter (lsimmonths) timescales is not predictive of multi-year variation. We find no evidence for asymmetry in the distribution of ΔEW that would indicate that BALs form and decay on different timescales, and we constrain the typical BAL lifetime to be gsim30 yr. The BAL absorption for one source, LBQS 0022+0150, has weakened and may now be classified as a mini-BAL. Another source, 1235+1453, shows evidence of variable, blue continuum emission that is relatively unabsorbed by the BAL outflow. C IV and Si IV BAL shape changes are related in at least some sources. Given their high velocities, BAL outflows apparently traverse large spatial regions and may interact with parsec-scale structures such as an obscuring torus. Assuming BAL outflows are launched from a rotating accretion disk, notable azimuthal symmetry is required in the outflow to explain the relatively small changes observed in velocity structure over times up to 7 yr.

Orbital variability of the PSR J2051-0827 binary system

NASA Astrophysics Data System (ADS)

Doroshenko, O.; Löhmer, O.; Kramer, M.; Jessner, A.; Wielebinski, R.; Lyne, A. G.; Lange, Ch.

2001-11-01

We have carried out high-precision timing measurements of the binary millisecond pulsar PSR J2051-0827 with the Effelsberg 100-m radio telescope of the Max-Planck-Institut für Radioastronomie and with the Lovell 76-m radio telescope at Jodrell Bank. The 6.5-yrs radio timing measurements have revealed a significant secular variation of the projected semi-major axis of the pulsar at a rate of dot xequiv d(a1 sin i)/dt = (-0.23+/- 0.03)x 10-12, which is probably caused by the Newtonian spin-orbit coupling in this binary system leading to a precession of the orbital plane. The required misalignment of the spin and orbital angular momenta of the companion are evidence for an asymmetric supernova explosion. We have also confirmed that the orbital period is currently decreasing at a rate of dot Pb=(-15.5 +/- 0.8)x 10-12 s s-1 and have measured second and third orbital period derivatives d2Pb/dt2=(+2.1 +/- 0.3)x 10-20: s-1 and d3Pb/dt3 =(3.6 +/- 0.6)x 10-28: s-2, which indicate a quasi-cyclic orbital period variation similar to those found in another eclipsing pulsar system, PSR B1957+20. The observed variation of the orbital parameters constrains the maximal value of the companion radius to Rc: max ~ 0.06: Rsun and implies that the companion is underfilling its Roche lobe by 50%. The derived variation in the quadrupole moment of the companion is probably caused by tidal dissipation similar to the mechanism proposed for PSR B1957+20. We conclude that the companion is at least partially non-degenerate, convective and magnetically active.

Regulation of water flux through tropical forest canopy trees: do universal rules apply?

PubMed

Meinzer, F C; Goldstein, G; Andrade, J L

2001-01-01

Tropical moist forests are notable for their richness in tree species. The presence of such a diverse tree flora presents potential problems for scaling up estimates of water use from individual trees to entire stands and for drawing generalizations about physiological regulation of water use in tropical trees. We measured sapwood area or sap flow, or both, in 27 co-occurring canopy species in a Panamanian forest to determine the extent to which relationships between tree size, sapwood area and sap flow were species-specific, or whether they were constrained by universal functional relationships between tree size, conducting xylem area, and water use. For the 24 species in which active xylem area was estimated over a range of size classes, diameter at breast height (DBH) accounted for 98% of the variation in sapwood area and 67% of the variation in sapwood depth when data for all species were combined. The DBH alone also accounted for > or = 90% of the variation in both maximum and total daily sap flux density in the outermost 2 cm of sapwood for all species taken together. Maximum sap flux density measured near the base of the tree occurred at about 1,400 h in the largest trees and 1,130 h in the smallest trees studied, and DBH accounted for 93% of the variation in the time of day at which maximum sap flow occurred. The shared relationship between tree size and time of maximum sap flow at the base of the tree suggests that a common relationship between diurnal stem water storage capacity and tree size existed. These results are consistent with a recent hypothesis that allometric scaling of plant vascular systems, and therefore water use, is universal.

NoMelt Experiment: High-resolution constraints on Pacific upper mantle fabric inferred from radial and azimuthal anisotropy

NASA Astrophysics Data System (ADS)

Russell, J. B.; Gaherty, J. B.; Lin, P. P.; Lizarralde, D.; Collins, J. A.; Hirth, G.; Evans, R. L.

2017-12-01

Observations of seismic anisotropy in the ocean basins are important for constraining deformation and melting processes in the upper mantle. The NoMelt OBS array was deployed on relatively pristine, 70 Ma seafloor in the central Pacific with the aim of constraining upper mantle circulation and the evolution of the lithosphere-asthenosphere system. Surface-waves traversing the array provide a unique opportunity to estimate a comprehensive set of anisotropic parameters. Azimuthal variations in Rayleigh-wave velocity over a period band of 15-180 s suggest strong anisotropic fabric both in the lithosphere and deep in the asthenosphere. High-frequency ambient noise (4-10 s) provides constraints on average VSV and VSH as well as azimuthal variations in both VS and VP in the upper ˜10 km of the mantle. Our best fitting models require radial anisotropy in the uppermost mantle with VSH > VSV by 3 - 7% and as much as 2% radial anisotropy in the crust. Additionally, we find a strong azimuthal dependence for Rayleigh- and Love-wave velocities, with Rayleigh 2θ fast direction parallel to the fossil spreading direction (FSD) and Love 2θ and 4θ fast directions shifted 90º and 45º from the FSD, respectively. These are some of the first direct observations of the Love 2θ and 4θ azimuthal signal, which allows us to directly invert for anisotropic terms G, B, and E in the uppermost Pacific lithosphere, for the first time. Together, these observations of radial and azimuthal anisotropy provide a comprehensive picture of oceanic mantle fabric and are consistent with horizontal alignment of olivine with the a-axis parallel to fossil spreading and having an orthorhombic or hexagonal symmetry.

Frequency Drives Lexical Access in Reading but not in Speaking: The Frequency-Lag Hypothesis

PubMed Central

Gollan, Tamar H.; Slattery, Timothy J.; Goldenberg, Diane; van Assche, Eva; Duyck, Wouter; Rayner, Keith

2010-01-01

To contrast mechanisms of lexical access in production versus comprehension we compared the effects of word-frequency (high, low), context (none, low-constraining, high-constraining), and level of English proficiency (monolinguals, Spanish-English bilinguals, Dutch-English bilinguals), on picture naming, lexical decision, and eye fixation times. Semantic constraint effects were larger in production than in reading. Frequency effects were larger in production than in reading without constraining context, but larger in reading than in production with constraining context. Bilingual disadvantages were modulated by frequency in production but not in eye fixation times, were not smaller in low-constraining context, and were reduced by high-constraining context only in production and only at the lowest level of English proficiency. These results challenge existing accounts of bilingual disadvantages, and reveal fundamentally different processes during lexical access across modalities, entailing a primarily semantically driven search in production, but a frequency driven search in comprehension. The apparently more interactive process in production than comprehension could simply reflect a greater number of frequency-sensitive processing stages in production. PMID:21219080

Using GIA observables to constrain the thermal contribution to lateral variations in mantle viscosity

NASA Astrophysics Data System (ADS)

Wu, P.; Wang, H.; van der Wal, W.; Shum, C.; Lee, H.; Braun, A.

2009-05-01

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, observations of the glacial isostatic adjustment (GIA) process are used to constrain the thermal contribution to lateral variations in upper and lower mantle viscosities. The Coupled Laplace-Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic earth with self-gravitating oceans, induced by either the ICE-5G or ICE-4G deglaciation models. GIA observations include global historic relative sea level data, GPS uplift rates in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide. The lateral viscosity perturbations are inferred from the high resolution seismic tomography model of Grand (2002) by using a conversion relation that takes into account both anelastic and anharmonic effects (Karato 2008). To determine the contribution of thermal effects in the upper and lower mantle, the scaling factor b is also inserted into the conversion relation: For b = 1, lateral velocity variations are caused by thermal effects alone; while b < 1 indicates a decreasing contribution of thermal effects; eventually when b = 0, there is no lateral viscosity variations exist and the Earth is laterally homogeneous. The value of b in the upper mantle is b1 while that in the lower mantle is b2. The lateral viscosity variations computed this way are superposed on a reference model that is able to give a reasonably good fit to the GIA observations. The parameter space for (b1, b2) is then searched to find the combination that yields the best improvement in fitting the GIA data in Laurentide, Fennoscandia or globally.

A 3000-year record of ground-rupturing earthquakes along the central North Anatolian fault near Lake Ladik, Turkey

USGS Publications Warehouse

Fraser, J.; Pigati, J.S.; Hubert-Ferrari, A.; Vanneste, K.; Avsar, U.; Altinok, S.

2009-01-01

The North Anatolian fault (NAF) is a ???1500 km long, arcuate, dextral strike-slip fault zone in northern Turkey that extends from the Karliova triple junction to the Aegean Sea. East of Bolu, the fault zone exhibits evidence of a sequence of large (Mw >7) earthquakes that occurred during the twentieth century that displayed a migrating earthquake sequence from east to west. Prolonged human occupation in this region provides an extensive, but not exhaustive, historical record of large earthquakes prior to the twentieth century that covers much of the last 2000 yr. In this study, we extend our knowledge of rupture events in the region by evaluating the stratigraphy and chronology of sediments exposed in a paleoseismic trench across a splay of the NAF at Destek, ???6:5 km east of Lake Ladik (40.868?? N, 36.121?? E). The trenched fault strand forms an uphill-facing scarp and associated sediment trap below a small catchment area. The trench exposed a narrow fault zone that has juxtaposed a sequence of weakly defined paleosols interbedded with colluvium against highly fractured bedrock. We mapped magnetic susceptibility variations on the trench walls and found evidence for multiple visually unrecognized colluvial wedges. This technique was also used to constrain a predominantly dip-slip style of displacement on this fault splay. Sediments exposed in the trench were dated using both charcoal and terrestrial gastropod shells to constrain the timing of the earthquake events. While the gastropod shells consistently yielded 14 C ages that were too old (by ???900 yr), we obtained highly reliable 14 C ages from the charcoal by dating multiple components of the sample material. Our radiocarbon chronology constrains the timing of seven large earthquakes over the past 3000 yr prior to the 1943 Tosya earthquake, including event ages of (2?? error): A.D. 1437-1788, A.D. 1034-1321, A.D. 549-719, A.D. 17-585 (1-3 events), 35 B.C.-A.D. 28, 700-392 B.C., 912-596 B.C. Our results indicate an average interevent time of 385 166?? yr (1??).

A variational approach to dynamics of flexible multibody systems

NASA Technical Reports Server (NTRS)

Wu, Shih-Chin; Haug, Edward J.; Kim, Sung-Soo

1989-01-01

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body references frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed.

Phenomenological constraints on accretion of non-annihilating dark matter on the PSR B1257+12 pulsar from orbital dynamics of its planets

NASA Astrophysics Data System (ADS)

Iorio, Lorenzo

2010-11-01

We analytically compute the effects that a pulsar's mass variation, whatever its physical origin may be, has on the standard Keplerian changes ΔτKep in the times of arrival of its pulses due to potential test particle companions, and on their orbital dynamics over long time scales. We apply our results to the planetary system of the PSR B1257+12 pulsar, located in the Galaxy at ~ 600 pc from us, to phenomenologically constrain a putative accretion of non-annihilating dark matter on the hosting neutron star. By comparing our prediction for Δτdot M/M to the root-mean-square accuracy of the timing residuals δ(Δτ) = 3.0μs we find for the mass variation rate dot M/M <= 1.3 × 10-6 yr-1. Actually, considerations related to the pulsar's lifetime, of the order of Δt ~ 0.8 Gyr, and to the currently accepted picture of the formation of its planets point toward a tighter constrain on the mass accretion rate, i.e. dot M/M <= 10-9 yr-1. Otherwise, the planets would have formed at about 300-700 au from PSR B1257+12, i.e. too far with respect to the expected extension of 1-2 au of the part of the protoplanetary disk containing the solid constituents from which they likely originated. In fact, an even smaller upper limit, dot M/M <= 10-11 yr-1, would likely be more realistic to avoid certain technical inconsistencies with the quality of the fit of the timing data, performed by keeping the standard value M = 1.4Modot fixed for the neutron star's mass. Anyway, the entire pulsar data set should be re-processed by explicitly modeling the mass variation rate and solving for it. Model-dependent theoretical predictions for the pulsar's mass accretion, in the framework of the mirror matter scenario, yield a mass increment rate of about 10-16 yr-1 for a value of the density of mirror matter ρdm as large as 10-17 g cm-3 = 5.6 × 106 GeV cm-3. Such a rate corresponds to a fractional mass variation of ΔM/M ~ 10-7 over the pulsar's lifetime. It would imply a formation of a black hole from the accreted dark matter inner core for values of the dark matter particle's mass mdm larger than 3 × 103 Gev, which are, thus, excluded since PSR B1257+12 is actually not such a kind of compact object. Instead, by assuming ρdm ~ 10-24 g cm-3 = 0.56 GeV cm-3, the mass accretion rate would be dot M/M ~ 10-23 yr-1, with a fractional mass variation of the order of ΔM/M ~ 10-14. It rules out mdm >= 8 × 106 Gev. Extreme values ρdm = 1.8 × 10-13 g cm-3 = 1011 GeV cm-3 for non-annihilating dark matter in central spike may yield the constraint dot M/M <= 10-11 yr-1 over Δt = 0.8 Gyr, it rules out mdm >= 12 Gev.

Cosmology of a Friedmann-Lamaître-Robertson-Walker 3-brane, late-time cosmic acceleration, and the cosmic coincidence.

PubMed

Doolin, Ciaran; Neupane, Ishwaree P

2013-04-05

A late epoch cosmic acceleration may be naturally entangled with cosmic coincidence--the observation that at the onset of acceleration the vacuum energy density fraction nearly coincides with the matter density fraction. In this Letter we show that this is indeed the case with the cosmology of a Friedmann-Lamaître-Robertson-Walker (FLRW) 3-brane in a five-dimensional anti-de Sitter spacetime. We derive the four-dimensional effective action on a FLRW 3-brane, from which we obtain a mass-reduction formula, namely, M(P)(2) = ρ(b)/|Λ(5)|, where M(P) is the effective (normalized) Planck mass, Λ(5) is the five-dimensional cosmological constant, and ρ(b) is the sum of the 3-brane tension V and the matter density ρ. Although the range of variation in ρ(b) is strongly constrained, the big bang nucleosynthesis bound on the time variation of the effective Newton constant G(N) = (8πM(P)(2))(-1) is satisfied when the ratio V/ρ ≳ O(10(2)) on cosmological scales. The same bound leads to an effective equation of state close to -1 at late epochs in accordance with astrophysical and cosmological observations.

Adaptive form-finding method for form-fixed spatial network structures

NASA Astrophysics Data System (ADS)

Lan, Cheng; Tu, Xi; Xue, Junqing; Briseghella, Bruno; Zordan, Tobia

2018-02-01

An effective form-finding method for form-fixed spatial network structures is presented in this paper. The adaptive form-finding method is introduced along with the example of designing an ellipsoidal network dome with bar length variations being as small as possible. A typical spherical geodesic network is selected as an initial state, having bar lengths in a limit group number. Next, this network is transformed into the ellipsoidal shape as desired by applying compressions on bars according to the bar length variations caused by transformation. Afterwards, the dynamic relaxation method is employed to explicitly integrate the node positions by applying residual forces. During the form-finding process, the boundary condition of constraining nodes on the ellipsoid surface is innovatively considered as reactions on the normal direction of the surface at node positions, which are balanced with the components of the nodal forces in a reverse direction induced by compressions on bars. The node positions are also corrected according to the fixed-form condition in each explicit iteration step. In the serial results of time history, the optimal solution is found from a time history of states by properly choosing convergence criteria, and the presented form-finding procedure is proved to be applicable for form-fixed problems.

Neuroendocrine mechanisms underlying behavioral stability: implications for the evolutionary origin of personality.

PubMed

Duckworth, Renée A

2015-12-01

Personality traits are behaviors that show limited flexibility over time and across contexts, and thus understanding their origin requires an understanding of what limits behavioral flexibility. Here, I suggest that insight into the evolutionary origin of personality traits requires determining the relative importance of selection and constraint in producing limits to behavioral flexibility. Natural selection as the primary cause of limits to behavioral flexibility assumes that the default state of behavior is one of high flexibility and predicts that personality variation arises through evolution of buffering mechanisms to stabilize behavioral expression, whereas the constraint hypothesis assumes that the default state is one of limited flexibility and predicts that the neuroendocrine components that underlie personality variation are those most constrained in flexibility. Using recent work on the neurobiology of sensitive periods and maternal programming of offspring behavior, I show that some of the most stable aspects of the neuroendocrine system are structural components and maternally induced epigenetic effects. Evidence of numerous constraints to changes in structural features of the neuroendocrine system and far fewer constraints to flexibility of epigenetic systems suggests that structural constraints play a primary role in the origin of behavioral stability and that epigenetic programming may be more important in generating adaptive variation among individuals. © 2015 New York Academy of Sciences.

High-precision photometry of WASP-12 b transits

NASA Astrophysics Data System (ADS)

Maciejewski, G.; Errmann, R.; Raetz, St.; Seeliger, M.; Spaleniak, I.; Neuhäuser, R.

2011-04-01

Aims: The transiting extrasolar planet WASP-12 b was found to be one of the most intensely irradiated exoplanets. It is unexpectedly bloated and is losing mass that may accrete into the host star. Our aim was to refine the parameters of this intriguing system and search for signs of transit timing variations. Methods: We gathered high-precision light curves for two transits of WASP-12 b. Assuming various limb-darkening laws, we generated best-fitting models and redetermined the parameters of the system. Error estimates were derived by the prayer-bead method and Monte Carlo simulations. Results: System parameters obtained by us are found to agree with previous studies within one sigma. Use of the non-linear limb-darkening laws results in the best-fitting models. With two new mid-transit times, the ephemeris was refined to BJDTDB = (2 454 508.97682 ± 0.00020) + (1.09142245 ± 0.00000033)E. Interestingly, indications of transit timing variation are detected at the level of 3.4 sigma. This signal can be induced by an additional planet in the system. Simplified numerical simulations show that a perturber could be a terrestrial-type planet if both planets are in a low-order orbital resonance. However, we emphasise that further observations are needed to confirm variation and to constrain properties of the perturber. Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA), operated jointly by the Max-Planck Institut für Astronomie and the Instituto de Astrofisica de Andalucia (CSIC).Photometric data are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/528/A65

Evaluation of the evolving stress field of the Yellowstone volcanic plateau, 1988 to 2010, from earthquake first-motion inversions

NASA Astrophysics Data System (ADS)

Russo, E.; Waite, G. P.; Tibaldi, A.

2017-03-01

Although the last rhyolite eruption occurred around 70 ka ago, the silicic Yellowstone volcanic field is still considered active due to high hydrothermal and seismic activity and possible recent magma intrusions. Geodetic measurements document complex deformation patterns in crustal strain and seismic activity likewise reveal spatial and temporal variations in the stress field. We use earthquake data recorded between 1988 and 2010 to investigate these variations and their possible causes in more detail. Earthquake relocations and a set of 369 well-constrained, double-couple, focal mechanism solutions were computed. Events were grouped according to location and time to investigate trends in faulting. The majority of the events have normal-faulting solutions, subordinate strike-slip kinematics, and very rarely, reverse motions. The dominant direction of extension throughout the 0.64 Ma Yellowstone caldera is nearly ENE, consistent with the perpendicular direction of alignments of volcanic vents within the caldera, but our study also reveals spatial and temporal variations. Stress-field solutions for different areas and time periods were calculated from earthquake focal mechanism inversion. A well-resolved rotation of σ3 was found, from NNE-SSW near the Hebgen Lake fault zone, to ENE-WSW near Norris Junction. In particular, the σ3 direction changed throughout the years around Norris Geyser Basin, from being ENE-WSW, as calculated in the study by Waite and Smith (2004), to NNE-SSW, while the other σ3 directions are mostly unchanged over time. The presence of ;chocolate tablet; structures, with two sets of nearly perpendicular normal faults, was identified in many stages of the deformation history both in the Norris Geyser Basin area and inside the caldera.

Coupling Poisson rectangular pulse and multiplicative microcanonical random cascade models to generate sub-daily precipitation timeseries

NASA Astrophysics Data System (ADS)

Pohle, Ina; Niebisch, Michael; Müller, Hannes; Schümberg, Sabine; Zha, Tingting; Maurer, Thomas; Hinz, Christoph

2018-07-01

To simulate the impacts of within-storm rainfall variabilities on fast hydrological processes, long precipitation time series with high temporal resolution are required. Due to limited availability of observed data such time series are typically obtained from stochastic models. However, most existing rainfall models are limited in their ability to conserve rainfall event statistics which are relevant for hydrological processes. Poisson rectangular pulse models are widely applied to generate long time series of alternating precipitation events durations and mean intensities as well as interstorm period durations. Multiplicative microcanonical random cascade (MRC) models are used to disaggregate precipitation time series from coarse to fine temporal resolution. To overcome the inconsistencies between the temporal structure of the Poisson rectangular pulse model and the MRC model, we developed a new coupling approach by introducing two modifications to the MRC model. These modifications comprise (a) a modified cascade model ("constrained cascade") which preserves the event durations generated by the Poisson rectangular model by constraining the first and last interval of a precipitation event to contain precipitation and (b) continuous sigmoid functions of the multiplicative weights to consider the scale-dependency in the disaggregation of precipitation events of different durations. The constrained cascade model was evaluated in its ability to disaggregate observed precipitation events in comparison to existing MRC models. For that, we used a 20-year record of hourly precipitation at six stations across Germany. The constrained cascade model showed a pronounced better agreement with the observed data in terms of both the temporal pattern of the precipitation time series (e.g. the dry and wet spell durations and autocorrelations) and event characteristics (e.g. intra-event intermittency and intensity fluctuation within events). The constrained cascade model also slightly outperformed the other MRC models with respect to the intensity-frequency relationship. To assess the performance of the coupled Poisson rectangular pulse and constrained cascade model, precipitation events were stochastically generated by the Poisson rectangular pulse model and then disaggregated by the constrained cascade model. We found that the coupled model performs satisfactorily in terms of the temporal pattern of the precipitation time series, event characteristics and the intensity-frequency relationship.

Rupture process of a multiple main shock sequence: analysis of teleseismic, local and field observations of the Tennant Creek, Australia, earthquakes of January 22, 1988

USGS Publications Warehouse

Choy, G.L.; Bowman, J.R.

1990-01-01

On January 22, 1988, three large intraplate earthquakes (with MS 6.3, 6.4 and 6.7) occurred within a 12-hour period near Tennant Creek, Australia. Broadband displacement and velocity records of body waves from teleseismically recorded data are analyzed to determine source mechanisms, depths, and complexity of rupture of each of the three main shocks. Hypocenters of an additional 150 foreshocks and aftershocks constrained by local arrival time data and field observations of surface rupture are used to complement the source characteristics of the main shocks. The interpretation of the combined data sets suggests that the overall rupture process involved unusually complicated stress release. Rupture characteristics suggest that substantial slow slip occurred on each of the three fault interfaces that was not accompanied by major energy release. Variation of focal depth and the strong increase of moment and radiated energy with each main shock imply that lateral variations of strength were more important than vertical gradients of shear stress in controlling the progression of rupture. -from Authors

Archaeomagnetic results from southern Italy and their bearing on geomagnetic secular variation

NASA Astrophysics Data System (ADS)

Evans, M. E.; Hoye, G. S.

2005-07-01

Archaeodirectional results from kilns and other baked structures in southern Italy are presented. They are generally compatible with the much larger data sets from France and Bulgaria. In particular, a summary of all the results associated with the well-known eruption of Vesuvius that destroyed Pompeii ( n = 9, D = 355°, I = 58°, α95 = 1.5°) provides a reliable archaeomagnetic anchor point supporting the French and Bulgarian master curves. It is extremely well-constrained in time and it comprises independent studies carried out in four different countries. Furthermore, it is derived from a diverse set of features agreement amongst which argues strongly against significant perturbations due to magnetic refraction, structural disturbance, or depositional shallowing. In terms of geomagnetic secular variation, we interpret the western European archaeomagnetic data summarized here in terms of an open loop caused by westward drift, followed by an inclination low spanning the first few centuries CE representing the signal of a static flux pulse that reaches a maximum magnetic moment of a few percent of the earth's main central dipole.

Global solar wind variations over the last four centuries

PubMed Central

Owens, M. J.; Lockwood, M.; Riley, P.

2017-01-01

The most recent “grand minimum” of solar activity, the Maunder minimum (MM, 1650–1710), is of great interest both for understanding the solar dynamo and providing insight into possible future heliospheric conditions. Here, we use nearly 30 years of output from a data-constrained magnetohydrodynamic model of the solar corona to calibrate heliospheric reconstructions based solely on sunspot observations. Using these empirical relations, we produce the first quantitative estimate of global solar wind variations over the last 400 years. Relative to the modern era, the MM shows a factor 2 reduction in near-Earth heliospheric magnetic field strength and solar wind speed, and up to a factor 4 increase in solar wind Mach number. Thus solar wind energy input into the Earth’s magnetosphere was reduced, resulting in a more Jupiter-like system, in agreement with the dearth of auroral reports from the time. The global heliosphere was both smaller and more symmetric under MM conditions, which has implications for the interpretation of cosmogenic radionuclide data and resulting total solar irradiance estimates during grand minima. PMID:28139769

SCM-Forcing Data

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

Xie, Shaocheng; Tang, Shuaiqi; Zhang, Yunyan

2016-07-01

Single-Column Model (SCM) Forcing Data are derived from the ARM facility observational data using the constrained variational analysis approach (Zhang and Lin 1997 and Zhang et al., 2001). The resulting products include both the large-scale forcing terms and the evaluation fields, which can be used for driving the SCMs and Cloud Resolving Models (CRMs) and validating model simulations.

Effect of water management variation on As and Cd accumulation or rice grain

USDA-ARS?s Scientific Manuscript database

Because of the current interest in As accumulation in rice, we examined rice grain As and Cd levels from a field test evaluating the effect of irrigation water management. The original study was conducted to test water saving production methods because limitation on water supply is constraining prod...

Seed availability constrains plant species sorting along a soil fertility gradient

Treesearch

Bryan L. Foster; Erin J. Questad; Cathy D. Collins; Cheryl A. Murphy; Timothy L. Dickson; Val H. Smith

2011-01-01

1. Spatial variation in species composition within and among communities may be caused by deterministic, niche-based species sorting in response to underlying environmental heterogeneity as well as by stochastic factors such as dispersal limitation and variable species pools. An important goal in ecology is to reconcile deterministic and stochastic perspectives of...

A chance constraint estimation approach to optimizing resource management under uncertainty

Treesearch

Michael Bevers

2007-01-01

Chance-constrained optimization is an important method for managing risk arising from random variations in natural resource systems, but the probabilistic formulations often pose mathematical programming problems that cannot be solved with exact methods. A heuristic estimation method for these problems is presented that combines a formulation for order statistic...

Research Review: A Neuroscience Framework for Pediatric Anxiety Disorders

ERIC Educational Resources Information Center

Pine, Daniel S.

2007-01-01

Across a range of mammalian species, early developmental variations in fear-related behaviors constrain patterns of anxious behavior throughout life. Individual differences in anxiety among rodents and non-human primates have been shown to reflect early-life influences of genes and the environment on brain circuitry. However, in humans, the manner…

Dew-worms in white nights: High latitude light constrains earthworm (Lumbricus terrestris) behaviour at the soil surface

USDA-ARS?s Scientific Manuscript database

Soil is an effective barrier to light penetration that limits the direct influence of light on belowground organisms. Variation in aboveground light conditions, however, is important to soil-dwelling animals that are periodically active on the soil surface. A prime example is the earthworm Lumbricus...

Genome-Wide Associations for Multiple Pest Resistances in a Northwestern United States Elite Spring Wheat Panel

USDA-ARS?s Scientific Manuscript database

Northern areas of the western United States are one of the most productive wheat growing regions in the United States. Increasing productivity through breeding is hindered by several biotic stresses which slow and constrain targeted yield improvement. In order to understand genetic variation for str...

Audible Conversational Affordances and Constraints of Verbalizing Professional Noticing during Prospective Teacher Lesson Study

ERIC Educational Resources Information Center

Amador, Julie M.; Carter, Ingrid S.

2018-01-01

Lesson Study, which assumes multiple variations, is a collaborative method that provides an opportunity for teachers to notice professionally. We analyzed how audible conversational components afforded and constrained the verbalization of professional noticing of a lesson study team. The team, comprised of six prospective teachers, a classroom…

Detection of a Sharp Structural Boundary in Lowermost Mantle Beneath Alaska by Core Phase PKPbc-df Differential Travel Times - Observation from the Anomalous South Sandwich Islands to Alaska Path

NASA Astrophysics Data System (ADS)

Xin, L.; Kawakatsu, H.; Takeuchi, N.

2017-12-01

Differential travel time residuals of PKPbc and PKPdf for the path from South Sandwich Islands (SSI) to Alaska are usually used to constrain anisotropy of the western hemisphere of the Earth's inner-core. For this polar path, it has been found that PKPbc-df differential residuals are generally anomalously larger than data that sample other regions, and also show strong lateral variation. Due to sparse distribution of seismic stations in Alaska in early times, previous researches have been unable to propose a good model to explain this particular data set. Using data recorded by the current dense stations in Alaska for SSI earthquakes, we reexamine the anomalous behavior of core phase PKPbc-df differential travel times and try to explain the origin. The data sample the inner-core for the polar paths, as well as the lowermost mantle beneath Alaska. Our major observations are: (1) fractional travel time residuals of PKPbc-df increase rapidly within 2° (up to 1%). (2) A clear shift of the residual pattern could be seen for earthquakes with different locations. (3) The residual shows systematic lateral variation: at northern part, no steep increase of residual can be seen. A sharp lateral structural boundary with a P-wave velocity contrast of about 3% at lowermost mantle beneath East Alaska is invoked to explain the steep increase of the observed residuals. By combining the effects of a uniformly anisotropic inner-core and the heterogeneity, the observed residual patterns could be well reproduced. This high velocity anomaly might be related with an ancient subducted slab. Lateral variation of the PKPbc-df residuals suggests that the heterogeneity layer is not laterally continuous and may terminate beneath Northeastern Alaska. We also conclude that core phases may be strongly affected by heterogeneities at lowermost mantle, and should be carefully treated if they are used to infer the inner-core structure.

Satellite Gravity Drilling the Earth

NASA Technical Reports Server (NTRS)

vonFrese, R. R. B.; Potts, L. V.; Leftwich, T. E.; Kim, H. R.; Han, S.-H.; Taylor, P. T.; Ashgharzadeh, M. F.

2005-01-01

Analysis of satellite-measured gravity and topography can provide crust-to-core mass variation models for new insi@t on the geologic evolution of the Earth. The internal structure of the Earth is mostly constrained by seismic observations and geochemical considerations. We suggest that these constraints may be augmented by gravity drilling that interprets satellite altitude free-air gravity observations for boundary undulations of the internal density layers related to mass flow. The approach involves separating the free-air anomalies into terrain-correlated and -decorrelated components based on the correlation spectrum between the anomalies and the gravity effects of the terrain. The terrain-decorrelated gravity anomalies are largely devoid of the long wavelength interfering effects of the terrain gravity and thus provide enhanced constraints for modeling mass variations of the mantle and core. For the Earth, subcrustal interpretations of the terrain-decorrelated anomalies are constrained by radially stratified densities inferred from seismic observations. These anomalies, with frequencies that clearly decrease as the density contrasts deepen, facilitate mapping mass flow patterns related to the thermodynamic state and evolution of the Earth's interior.

The Pendulum: From Constrained Fall to the Concept of Potential

ERIC Educational Resources Information Center

Bevilacqua, Fabio; Falomo, Lidia; Fregonese, Lucio; Giannetto, Enrico; Giudice, Franco; Mascheretti, Paolo

2006-01-01

Kuhn underlined the relevance of Galileo's gestalt switch in the interpretation of a swinging body from constrained fall to time metre. But the new interpretation did not eliminate the older one. The constrained fall, both in the motion of pendulums and along inclined planes, led Galileo to the law of free fall. Experimenting with physical…

Patterns of Diel Variation in Nitrate Concentrations in the Potomac River

NASA Astrophysics Data System (ADS)

Burns, D. A.; Miller, M. P.; Pellerin, B. A.; Capel, P. D.

2015-12-01

The Potomac River is the second largest source of nitrogen to Chesapeake Bay, where reducing nutrient loads has been a focus of efforts to improve estuarine trophic status. Two years of high frequency sensor measurements of nitrate (NO3-) concentrations in the Upper Potomac River at the Little Falls gage were analyzed to quantify seasonal variation in the magnitude and timing of the apparent loss of NO3- from the water column that results from diel-driven processes. In addition to broad seasonal and flow-driven variation in NO3- concentrations, clear diel patterns were evident in the river, especially during low flow conditions that follow stormflow by several days. Diel variation was about 0.01 mg N/L in winter and 0.02 to 0.03 mg N/L in summer with intermediate values during spring and fall. This variation was equivalent to <1% of the mean daily NO3- concentration in winter and about 4% in summer; however, variation >10% occurred during some summer days. Maximum diel concentrations occurred during mid- to late-morning in most seasons, with the most repeatable patterns in summer and wider variation in timing during fall and winter. Diel NO3- loss diminished loads by about 0.6% in winter and 1.3% in summer, and diel-driven processes were minor compared to estimates of total in-stream NO3- loss that averaged about one-third of the inferred groundwater NO3- contribution to the river network. The magnitude of diel NO3- variation was more strongly related to metrics based on water temperature and discharge than to metrics based on photosynthetically active radiation. Despite the fairly low diminishment of NO3- loads attributable to diel variation, estimates of diel NO3- uptake were fairly high compared to published values from smaller streams and rivers. The diel NO3- patterns observed in the Potomac River are consistent with photosynthesis of periphyton as a principal driver which may be linked to denitrification through the release of labile carbon. The extent to which these diel patterns are related to measures of aquatic metabolism are unknown as is the role of dispersion in obscuring diel patterns. Improvements to these diel estimates will require additional measures such as dissolved oxygen and ammonium, and the use of a second upstream measurement station to better constrain NO3- uptake values.

Rapid variations in fluid chemistry constrain hydrothermal phase separation at the Main Endeavour Field

NASA Astrophysics Data System (ADS)

Love, Brooke; Lilley, Marvin; Butterfield, David; Olson, Eric; Larson, Benjamin

2017-02-01

Previous work at the Main Endeavour Field (MEF) has shown that chloride concentration in high-temperature vent fluids has not exceeded 510 mmol/kg (94% of seawater), which is consistent with brine condensation and loss at depth, followed by upward flow of a vapor phase toward the seafloor. Magmatic and seismic events have been shown to affect fluid temperature and composition and these effects help narrow the possibilities for sub-surface processes. However, chloride-temperature data alone are insufficient to determine details of phase separation in the upflow zone. Here we use variation in chloride and gas content in a set of fluid samples collected over several days from one sulfide chimney structure in the MEF to constrain processes of mixing and phase separation. The combination of gas (primarily magmatic CO2 and seawater-derived Ar) and chloride data, indicate that neither variation in the amount of brine lost, nor mixing of the vapor phase produced at depth with variable quantities of (i) brine or (ii) altered gas rich seawater that has not undergone phase separation, can explain the co-variation of gas and chloride content. The gas-chloride data require additional phase separation of the ascending vapor-like fluid. Mixing and gas partitioning calculations show that near-critical temperature and pressure conditions can produce the fluid compositions observed at Sully vent as a vapor-liquid conjugate pair or as vapor-liquid pair with some remixing, and that the gas partition coefficients implied agree with theoretically predicted values.