Using spatially distributed parameters and multi-response objective functions to solve parameterization of complex applications of semi-distributed hydrological models

Application of semi-distributed hydrological models to large, heterogeneous watersheds deals with several problems. On one hand, the spatial and temporal variability in catchment features should be adequately represented in the model parameterization, while maintaining the model complexity in an acceptable level to take advantage of state-of-the-art calibration techniques. On the other hand, model complexity enhances uncertainty in adjusted model parameter values, therefore increasing uncertainty in the water routing across the watershed. This is critical for water quality applications, where not only streamflow, but also a reliable estimation of the surface versus subsurface contributions to the runoff is needed. In this study, we show how a regularized inversion procedure combined with a multiobjective function calibration strategy successfully solves the parameterization of a complex application of a water quality-oriented hydrological model. The final value of several optimized parameters showed significant and consistentdifferences across geological and landscape features. Although the number of optimized parameters was significantly increased by the spatial and temporal discretization of adjustable parameters, the uncertainty in water routing results remained at reasonable values. In addition, a stepwise numerical analysis showed that the effects on calibration performance due to inclusion of different data types in the objective function could be inextricably linked. Thus caution should be taken when adding or removing data from an aggregated objective function.

The Role of model uncertainty and learning in the U.S. postwar policy response to oil prices

This paper studies optimal monetary policy in a framework that explicitly accounts for policymakers' uncertainty about the channels of transmission of oil prices into the economy. More specfically, I examine the robust response to the real price of oil that US monetary authorities would have been recommended to implement in the period 1970 2009; had they used the approach proposed by Cogley and Sargent (2005b) to incorporate model uncertainty and learning into policy decisions. In this context, I investigate the extent to which regulator' changing beliefs over different models of the economy play a role in the policy selection process. The main conclusion of this work is that, in the specific environment under analysis, one of the underlying models dominates the optimal interest rate response to oil prices. This result persists even when alternative assumptions on the model's priors change the pattern of the relative posterior probabilities, and can thus be attributed to the presence of model uncertainty itself.

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

From SNPs to pathways: integration of functional effect of sequence variations on models of cell signalling pathways

Background: Single nucleotide polymorphisms (SNPs) are the most frequent type of sequence variation between individuals, and represent a promising tool for finding genetic determinants of complex diseases and understanding the differences in drug response. In this regard, it is of particular interest to study the effect of non-synonymous SNPs in the context of biological networks such as cell signalling pathways. UniProt provides curated information about the functional and phenotypic effects of sequence variation, including SNPs, as well as on mutations of protein sequences. However, no strategy has been developed to integrate this information with biological networks, with the ultimate goal of studying the impact of the functional effect of SNPs in the structure and dynamics of biological networks. Results: First, we identified the different challenges posed by the integration of the phenotypic effect of sequence variants and mutations with biological networks. Second, we developed a strategy for the combination of data extracted from public resources, such as UniProt, NCBI dbSNP, Reactome and BioModels. We generated attribute files containing phenotypic and genotypic annotations to the nodes of biological networks, which can be imported into network visualization tools such as Cytoscape. These resources allow the mapping and visualization of mutations and natural variations of human proteins and their phenotypic effect on biological networks (e.g. signalling pathways, protein-protein interaction networks, dynamic models). Finally, an example on the use of the sequence variation data in the dynamics of a network model is presented. Conclusion: In this paper we present a general strategy for the integration of pathway and sequence variation data for visualization, analysis and modelling purposes, including the study of the functional impact of protein sequence variations on the dynamics of signalling pathways. This is of particular interest when the SNP or mutation is known to be associated to disease. We expect that this approach will help in the study of the functional impact of disease-associated SNPs on the behaviour of cell signalling pathways, which ultimately will lead to a better understanding of the mechanisms underlying complex diseases.

Business Cycles, Unemployment Insurance, and the Calibration of Matching Models

This paper points out an empirical puzzle that arises when an RBC economy with a job matching function is used to model unemployment. The standard model can generate sufficiently large cyclical fluctuations in unemployment, or a sufficiently small response of unemployment to labor market policies, but it cannot do both. Variable search and separation, finite UI benefit duration, efficiency wages, and capital all fail to resolve this puzzle. However, both sticky wages and match-specific productivity shocks help the model reproduce the stylized facts: both make the firm's flow of surplus more procyclical, thus making hiring more procyclical too.

Business cycles, unemployment insurance and the calibration of matching models

This paper theoretically and empirically documents a puzzle that arises when an RBC economy with a job matching function is used to model unemployment. The standard model can generate sufficiently large cyclical fluctuations in unemployment, or a sufficiently small response of unemployment to labor market policies, but it cannot do both. Variable search and separation, finite UI benefit duration, efficiency wages, and capital all fail to resolve this puzzle. However, either sticky wages or match-specific productivity shocks can improve the model's performance by making the firm's flow of surplus more procyclical, which makes hiring more procyclical too.

Periodic Solutions in Low-Dimensional Climatic Models

Classic climatic models use constitutive laws without any response time. A more realistic approach to the natural processes governing climate dynamics must introduce response time for heat and radiation fluxes. Extended irreversible thermodynamics (EIT) is a good thermodynamical framework for introducing nonclassical constitutive laws. In the present study EIT has been used to analyze a Budyko–Sellers one-dimensional energybalance model developed by G. R. North. The results present self-sustained periodic oscillations when the response time is greater than a critical value. The high-frequency (few kiloyears) damped and nondamped oscillations obtained can be related to abrupt climatic changes without any variation in the external forcing of the system

Second differential of the entropy as a criterion for the stability in low-dimensional climate models

The second differential of the entropy is used for analysing the stability of a thermodynamic climatic model. A delay time for the heat flux is introduced whereby it becomes an independent variable. Two different expressions for the second differential of the entropy are used: one follows classical irreversible thermodynamics theory; the second is related to the introduction of response time and is due to the extended irreversible thermodynamics theory. the second differential of the classical entropy leads to unstable solutions for high values of delay times. the extended expression always implies stable states for an ice-free earth. When the ice-albedo feedback is included, a discontinuous distribution of stable states is found for high response times. Following the thermodynamic analysis of the model, the maximum rates of entropy production at the steady state are obtained. A latitudinally isothermal earth produces the extremum in global entropy production. the material contribution to entropy production (by which we mean the production of entropy by material transport of heat) is a maximum when the latitudinal distribution of temperatures becomes less homogeneous than present values

Les set vides del paper: formats de presentació i models de distribució en línia dels treballs publicats en revistes científiques digitals

Study of the publication models and the means of accessing scientific literature in the current environment of digital communication and the web. The text introduces the concept of journal article as a well-defined and stable unit within the publishing world, and as a nucleus on which professional and scholarly communication has been based since its beginnings in the 17th century. The transformation of scientific communication that the digital world has enabled is analysed. Descriptions are provided of some of the practices undertaken by authors, research organisations, publishers and library-related institutions as a response to the new possibilities being unveiled for articles, both as products as well as for their creation and distribution processes. These transformations affect the very nature of articles as a minimal unit -both unique and stable- of scientific communication. The article concludes by noting that under varying documentary forms of publisher aggregation and bibliographic control -sometimes simultaneously and, even, apparently contradictory- there flourishes a more pluralistic type of scientific communication. This pluralism offers: more possibilities for communication among authors; fewer levels of intermediaries such as agents that intervene and provide added value to the products; greater availability for users both economically speaking and from the point of view of access; and greater interaction and wealth of contents, thanks to the new hypertext and multimedia possibilities.

Nonlinear inverse dynamic models of gas sensing systems based on chemical sensor arrays for quantitative measurements

Gas sensing systems based on low-cost chemical sensor arrays are gaining interest for the analysis of multicomponent gas mixtures. These sensors show different problems, e.g., nonlinearities and slow time-response, which can be partially solved by digital signal processing. Our approach is based on building a nonlinear inverse dynamic system. Results for different identification techniques, including artificial neural networks and Wiener series, are compared in terms of measurement accuracy.

A time-domain method for the analysis of thermal impedance response preserving the convolution form

The study of the thermal behavior of complex packages as multichip modules (MCM¿s) is usually carried out by measuring the so-called thermal impedance response, that is: the transient temperature after a power step. From the analysis of this signal, the thermal frequency response can be estimated, and consequently, compact thermal models may be extracted. We present a method to obtain an estimate of the time constant distribution underlying the observed transient. The method is based on an iterative deconvolution that produces an approximation to the time constant spectrum while preserving a convenient convolution form. This method is applied to the obtained thermal response of a microstructure as analyzed by finite element method as well as to the measured thermal response of a transistor array integrated circuit (IC) in a SMD package.

Impulse response recovery of linear systems through weighted cumulant slice

Identifiability of the so-called ω-slice algorithm is proven for ARMA linear systems. Although proofs were developed in the past for the simpler cases of MA and AR models, they were not extendible to general exponential linear systems. The results presented in this paper demonstrate a unique feature of the ω-slice method, which is unbiasedness and consistency when order is overdetermined, regardless of the IIR or FIR nature of the underlying system, and numerical robustness.

Quantitative developmental response to the length of exposure to long photoperiod in wheat and barley

The present study was designed to analyse the effect of the length of exposure to a long photoperiod imposed c. 3 weeks after sowing in spring wheat (cv. UQ189) and barley (cv. Arapiles) to (i) establish whether the response to the number of cycles of exposure is quantitative or qualitative, (ii) determine the existence of a commitment to particular stages well before the stage has been observable, and (iii) study the interrelationships between the effects on final leaf number and phyllochron when the stimulus is provided several days after seedling emergence. Both wheat and barley seemed to respond quantitatively to the number of long-day cycles they were exposed to. However, wheat showed a requirement of approximately 4 long-day cycles to be able to produce a significant response in time to heading. The barley cultivar used in the study was responsive to the minimum length of exposure. The response to extended photoperiod cycles during the stem elongation phase was due to the ‘ memory’ photoperiod effects being related, in the case of wheat, to the fact that the pre-terminal spikelet appearance phase saturated its photoperiod response well before that stage was reached. Therefore, the commitment to the terminal spikelet appearance in wheat may be reached well before this stage could be recognized. As the response in duration to heading exceeded that of the final leaf number, and the stem elongation phase responded to memory effects of photoperiod, the phyllochron of both cereals was responsive to the treatments accelerating the average phyllochron when exposed to longer periods of long days. The response in average phyllochron was due to a switch from bi-linear to linear models of leaf number v. time when the conditions were increasingly inductive, with the phyllochron of the initial (6–8) leaves being similar for all treatments (within each species), and from then on increased.

Disentangling the formation of contrasting tree line physiognomies combining model selection and Bayesian parameterization for simulation models.

Alpine tree-line ecotones are characterized by marked changes at small spatial scales that may result in a variety of physiognomies. A set of alternative individual-based models was tested with data from four contrasting Pinus uncinata ecotones in the central Spanish Pyrenees to reveal the minimal subset of processes required for tree-line formation. A Bayesian approach combined with Markov chain Monte Carlo methods was employed to obtain the posterior distribution of model parameters, allowing the use of model selection procedures. The main features of real tree lines emerged only in models considering nonlinear responses in individual rates of growth or mortality with respect to the altitudinal gradient. Variation in tree-line physiognomy reflected mainly changes in the relative importance of these nonlinear responses, while other processes, such as dispersal limitation and facilitation, played a secondary role. Different nonlinear responses also determined the presence or absence of krummholz, in agreement with recent findings highlighting a different response of diffuse and abrupt or krummholz tree lines to climate change. The method presented here can be widely applied in individual-based simulation models and will turn model selection and evaluation in this type of models into a more transparent, effective, and efficient exercise.

Ecological and physiological variance in T-cell mediated immune response in Cory's shearwaters

T-cell mediated immune response (CMI) hasbeen widely studied in relation to individual andfitness components in birds. However, few studieshave simultaneously examined individual and socialfactors and habitat-mediated variance in theimmunity of chicks and adults from the samepopulation and in the same breeding season. Weinvestigated ecological and physiological variancein CMI of male and female nestlings and adults in abreeding population of Cory's Shearwaters(Calonectrisdiomedea) in theMediterranean Sea. Explanatory variables includedindividual traits (body condition, carbon andnitrogen stable isotope ratios, plasma totalproteins, triglycerides, uric acid, osmolarity,β-hydroxy-butyrate, erythrocyte meancorpuscular diameter, hematocrit, andhemoglobin) and burrow traits(temperature, isolation, and physicalstructure). During incubation, immune responseof adult males was significantly greater than thatof females. Nestlings exhibited a lower immuneresponse than adults. Ecological and physiologicalfactors affecting immune response differed betweenadults and nestlings. General linear models showedthat immune response in adult males was positivelyassociated with burrow isolation, suggesting thatmales breeding at higher densities suffer immunesystem suppression. In contrast, immune response inchicks was positively associated with bodycondition and plasma triglyceride levels.Therefore, adult immune response appears to beassociated with social stress, whereas a trade-offbetween immune function and fasting capability mayexist for nestlings. Our results, and those fromprevious studies, provide support for anasymmetrical influence of ecological andphysiological factors on the health of differentage and sex groups within a population, and for theimportance of simultaneously considering individualand population characteristics in intraspecificstudies of immune response.