Non-constant discounting in finite horizon: The free terminal time case

This paper derives the HJB (Hamilton-Jacobi-Bellman) equation for sophisticated agents in a finite horizon dynamic optimization problem with non-constant discounting in a continuous setting, by using a dynamic programming approach. A simple example is used in order to illustrate the applicability of this HJB equation, by suggesting a method for constructing the subgame perfect equilibrium solution to the problem.Conditions for the observational equivalence with an associated problem with constantdiscounting are analyzed. Special attention is paid to the case of free terminal time. Strotz¿s model (an eating cake problem of a nonrenewable resource with non-constant discounting) is revisited.

Time consistent Pareto solutions in common access resource gameswith asymmetric players

In the analysis of equilibrium policies in a di erential game, if agents have different time preference rates, the cooperative (Pareto optimum) solution obtained by applying the Pontryagin's Maximum Principle becomes time inconsistent. In this work we derive a set of dynamic programming equations (in discrete and continuous time) whose solutions are time consistent equilibrium rules for N-player cooperative di erential games in which agents di er in their instantaneous utility functions and also in their discount rates of time preference. The results are applied to the study of a cake-eating problem describing the management of a common property exhaustible natural resource. The extension of the results to a simple common property renewable natural resource model in in nite horizon is also discussed.

A territorial model of commuting in Catalonia, 1986-1996

Commuting consists in the fact that an important fraction of workers in developed countries do not reside close to their workplaces but at long distances from them, so they have to travel to their jobs and then back home daily. Although most workers hold a job in the same municipality where they live or in a neighbouring one, an important fraction of workers face long daily trips to get to their workplace and then back home.Even if we divide Catalonia (Spain) in small aggregations of municipalities, trying to make them as close to local labour markets as possible, we will find out that some of them have a positive commuting balance, attracting many workers from other areas and providing local jobs for almost all their resident workers. On the other side, other zones seem to be mostly residential, so an important fraction of their resident workers hold jobs in different local labour markets. Which variables influence an area¿s role as an attraction pole or a residential zone? In previous papers (Artís et al, 1998a, 2000; Romaní, 1999) we have brought out the main individual variables that influence commuting by analysing a sample of Catalan workers and their commuting decisions. In this paper we perform an analysis of the territorial variables that influence commuting, using data for aggregate commuting flows in Catalonia from the 1991 and 1996 Spanish Population Censuses.These variables influence commuting in two different ways: a zone with a dense, welldeveloped economical structure will have a high density of jobs. Work demand cannot be fulfilled with resident workers, so it spills over local boundaries. On the other side, this economical activity has a series of side-effects like pollution, congestion or high land prices which make these areas less desirable to live in. Workers who can afford it may prefer to live in less populated, less congested zones, where they can find cheaper land, larger homes and a better quality of life. The penalty of this decision is an increased commuting time. Our aim in this paper is to highlight the influence of local economical structure and amenities endowment in the workplace-residence location decision. A place-to-place logit commuting models is estimated for 1991 and 1996 in order to find the economical and amenities variables with higher influence in commuting decisions. From these models, we can outline a first approximation to the evolution of these variables in the 1986-1996 period. Data have been obtained from aggregate flow travel-matrix from the 1986, 1991 and 1996 Spanish Population Censuses

Political economy of municipal water service privatization in Spain: a duration model analysis

This empirical work applies a duration model to the study of factors determining privatization of local water services. I assess how factors determining privatization decision evolve as time goes by. A sample of 133 Spanish municipalities during the six terms of office taken place during the 1980-2002 period is analyzed. A dynamic neighboring effect is hypothesized and successfully tested. In a first stage, private water supply firms may try to expand to regions where there is no service privatized, in order to spread over this region after having being installed thanks to its scale advantages. Other factors influencing privatization decision evolve during the two decades under study, from the priority to fix old infrastructures to the concern about service efficiency. Some complementary results regarding political and budgetary factors are also obtained

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

A joint back calculation model for the imputation of the date of HIV infection in a prevalent cohort.

In studies of the natural history of HIV-1 infection, the time scale of primary interest is the time since infection. Unfortunately, this time is very often unknown for HIV infection and using the follow-up time instead of the time since infection is likely to provide biased results because of onset confounding. Laboratory markers such as the CD4 T-cell count carry important information concerning disease progression and can be used to predict the unknown date of infection. Previous work on this topic has made use of only one CD4 measurement or based the imputation on incident patients only. However, because of considerable intrinsic variability in CD4 levels and because incident cases are different from prevalent cases, back calculation based on only one CD4 determination per person or on characteristics of the incident sub-cohort may provide unreliable results. Therefore, we propose a methodology based on the repeated individual CD4 T-cells marker measurements that use both incident and prevalent cases to impute the unknown date of infection. Our approach uses joint modelling of the time since infection, the CD4 time path and the drop-out process. This methodology has been applied to estimate the CD4 slope and impute the unknown date of infection in HIV patients from the Swiss HIV Cohort Study. A procedure based on the comparison of different slope estimates is proposed to assess the goodness of fit of the imputation. Results of simulation studies indicated that the imputation procedure worked well, despite the intrinsic high volatility of the CD4 marker.

Intensity correlation function of dye lasers: Short-time behavior

We propose an equation to calculate the intensity correlation function of a dye-laser model with a pump parameter subject to finite-bandwidth fluctuations. The equation is valid, in the weak-noise limit, for all times. It incorporates novel non-Markovian features. Results are given for the short-time behavior of the correlation function. It exhibits a characteristic initial plateau. Our findings are supported by a numerical simulation of the model.

Modeling changes in organic carbon stocks for distinct soils in southeastern brazil after four eucalyptus rotations using the century model

Soil organic matter (SOM) plays an important role in carbon (C) cycle and soil quality. Considering the complexity of factors that control SOM cycling and the long time it usually takes to observe changes in SOM stocks, modeling constitutes a very important tool to understand SOM cycling in forest soils. The following hypotheses were tested: (i) soil organic carbon (SOC) stocks would be higher after several rotations of eucalyptus than in low-productivity pastures; (ii) SOC values simulated by the Century model would describe the data better than the mean of observations. So, the aims of the current study were: (i) to evaluate the SOM dynamics using the Century model to simulate the changes of C stocks for two eucalyptus chronosequences in the Rio Doce Valley, Minas Gerais State, Brazil; and (ii) to compare the C stocks simulated by Century with the C stocks measured in soils of different Orders and regions of the Rio Doce Valley growing eucalyptus. In Belo Oriente (BO), short-rotation eucalyptus plantations had been cultivated for 4.0; 13.0, 22.0, 32.0 and 34.0 years, at a lower elevation and in a warmer climate, while in Virginópolis (VG), these time periods were 8.0, 19.0 and 33.0 years, at a higher elevation and in a milder climate. Soil samples were collected from the 0-20 cm layer to estimate C stocks. Results indicate that the C stocks simulated by the Century model decreased after 37 years of poorly managed pastures in areas previously covered by native forest in the regions of BO and VG. The substitution of poorly managed pastures by eucalyptus in the early 1970´s led to an average increase of C of 0.28 and 0.42 t ha-1 year-1 in BO and VG, respectively. The measured C stocks under eucalyptus in distinct soil Orders and independent regions with variable edapho-climate conditions were not far from the values estimated by the Century model (root mean square error - RMSE = 20.9; model efficiency - EF = 0.29) despite the opposite result obtained with the statistical procedure to test the identity of analytical methods. Only for lower soil C stocks, the model over-estimated the C stock in the 0-20 cm layer. Thus, the Century model is highly promising to detect changes in C stocks in distinct soil orders under eucalyptus, as well as to indicate the impact of harvest residue management on SOM in future rotations.

Relaxation time of processes driven by multiplicative noise

We consider systems described by nonlinear stochastic differential equations with multiplicative noise. We study the relaxation time of the steady-state correlation function as a function of noise parameters. We consider the white- and nonwhite-noise case for a prototype model for which numerical data are available. We discuss the validity of analytical approximation schemes. For the white-noise case we discuss the results of a projector-operator technique. This discussion allows us to give a generalization of the method to the non-white-noise case. Within this generalization, we account for the growth of the relaxation time as a function of the correlation time of the noise. This behavior is traced back to the existence of a non-Markovian term in the equation for the correlation function.

Electrical conductivity of soil irrigated with swine wastewater estimated by time-domain reflectometry¹

Wastewater application to soil is an alternative for fertilization and water reuse. However, particular care must be taken with this practice, since successive wastewater applications can cause soil salinization. Time-domain reflectometry (TDR) allows for the simultaneous and continuous monitoring of both soil water content and apparent electrical conductivity and thus for the indirect measurement of the electrical conductivity of the soil solution. This study aimed to evaluate the suitability of TDR for the indirect determination of the electrical conductivity (ECse) of the saturated soil extract by using an empirical equation for the apparatus TDR Trase 6050X1. Disturbed soil samples saturated with swine wastewater were used, at soil proportions of 0, 0.45, 0.90, 1.80, 2.70, and 3.60 m³ m-3. The probes were equipped with three handmade 0.20 cm long rods. The fit of the empirical model that associated the TDR measured values of electrical conductivity (EC TDR) to ECse was excellent, indicating this approach as suitable for the determination of electrical conductivity of the soil solution.

Follicular lymphoma at relapse after rituximab containing regimens: comparison of time to event intervals prior to and after (90) Y-ibritumomab-tiuxetan.

Radioimmunotherapies with Zevalin® (RIT-Z) showed encouraging results in patients with relapsed/refractory follicular lymphoma (FL), leading frequently to failure-free intervals longer than those achieved by the last previous therapy. We compared time-to-event variables obtained before and after RIT-Z in patients with relapsed FL, previously exposed to rituximab. All patients with relapsed non-transformed, non-refractory, non-rituximab-naïve FL who have been treated with RIT-Z in two different centres in Europe were included. Staging and response were assessed by contrast-enhanced CT in all patients; PET/CT was performed according to local availability. Event-free survival (EFS) and time to next treatment (TTNT) following the last previous therapy and after RIT-Z were compared. Pre-therapy characteristics were tested in univariate analyses for prediction of outcomes. A description of the patterns of relapse was also provided. Among 70 patients treated, only 16 fulfilled the inclusion criteria. They were treated with a median of 3 prior lines of chemo-immunotherapies, including a median of 2 rituximab-containing regimens; 6 patients had undergone myeloablative chemotherapy with autologous stem cell rescue (ASCT). Overall response rates were 10 (62%) CR/CRu, 3 (19%) PR and 3 (19%) PD; response rates were similar in patients with prior ASCT. After RIT-Z only few patients obtained EFS and TTNT longer than after the last previous therapy. All four patients receiving rituximab maintenance were without progression 12 months after RIT-Z. Relapses occurred in both previously and newly involved sites; a significant association was found between the number of pathologic sites involved prior to RIT-Z and subsequent TTNT. Despite the excellent response rate, the duration of response was shorter than the previous one confirming the known trend of relapses to occur earlier after subsequent treatments. Rituximab maintenance after RIT-Z showed encouraging results in terms of prolonging EFS, warranting further studies. Copyright © 2010 John Wiley & Sons, Ltd.

Competition between beta-subunits of cardiac L-type calcium channels

Cardiac L-type Ca (CaV1.2) channels are composed of a pore forming CaV1.2-α1 subunit and auxiliary β- and α2δ-subunits. β-subunits are important not only for surface expression of the channel pore but also for modulation of channel gating properties. Different β-subunits differentially modulate channel activity (Hullin et al., PLOSone, 2007) and thus L-type Ca2+ channel gating is altered when β-subunit expression pattern is changed. In human heart failure increased activity of single ventricular L-type Ca2+-channels is associated with an increased expression of β2-subunits. Interestingly, induction of β2-subunit over-expression in hearts of transgenic mice resembled this heart failure phenotype of hyperactive single L-type Ca2+-channel channels (Beetz et al., Cardiovasc Res. 2009). We hypothesised that competition of less stimulating β-subunits (e.g. β1) with β-subunits causing strong channel stimulation (e.g. β2) might be a means to treat dysfunctional L-type Ca2+-channel activity. To test this hypothesis, we performed whole-cell and single-channel measurements employing recombinant CaV1.2 channels expressed in HEK293 cells together with both β- and β1a2b-subunits. Whole-cell analysis revealed no differences of maximum L-type Ca2+-current densities [pA/pF] with coexpression of either β1a-subunits (-52±3.8), β2b-subunits (-61.5±6.6) or the mixtures of β- and β1a2b-subunits with the plasmid transfection ratio of 2:1 (-60.2±1.6) and 1:1 (-56.7±2.6) respectively. However, steady state inactivation kinetics differed between particular β-subunit and the relative amount of β-subunit presence in the mixtures (β1a1a-subunit (-41.1±1.0), β2b-subunits (-35.1±1.1), mixture 2:1 (-40.3±1.5), and mixture 1:1 (-38.4±2.0); [mV]; p<0.05, students t-test). Using a novel single-channel analysis, switching of gating between β1-like and β2-like modes was monitored on a minute time-scale when both β-subunits were co-expressed in the same cells, but the larger amount of β1a-subunits is required for the effective switching of gating. Our results indicate a model of mutually exclusive binding and effective competition between several β-subunits suggesting that hyperactive channel gating mediated e.g. by β2-subunits can be normalized by β1-subunits. Therefore, competitive replacement between different L-type Ca2+-channel β-subunits might serve as a novel therapeutic strategy for e.g. heart failure.

1/c expansion of a separable model of direct-interaction type

The properties of a proposed model of N point particles in direct interaction are considered in the limit of small velocities. It is shown that, in this limit, time correlations cancel out and that Newtonian dynamics is recovered for the system in a natural way.

Issue of time in generally covariant theories and the Komar-Bergmann approach to observables in general relativity

Diffeomorphism-induced symmetry transformations and time evolution are distinct operations in generally covariant theories formulated in phase space. Time is not frozen. Diffeomorphism invariants are consequently not necessarily constants of the motion. Time-dependent invariants arise through the choice of an intrinsic time, or equivalently through the imposition of time-dependent gauge fixation conditions. One example of such a time-dependent gauge fixing is the Komar-Bergmann use of Weyl curvature scalars in general relativity. An analogous gauge fixing is also imposed for the relativistic free particle and the resulting complete set time-dependent invariants for this exactly solvable model are displayed. In contrast with the free particle case, we show that gauge invariants that are simultaneously constants of motion cannot exist in general relativity. They vary with intrinsic time.

The nonlinear relaxation time and quasideterministic-theory approaches to characterize the decay of unestable states

We present the relationship between nonlinear-relaxation-time (NLRT) and quasideterministic approaches to characterize the decay of an unstable state. The universal character of the NLRT is established. The theoretical results are applied to study the dynamical relaxation of the Landau model in one and n variables and also a laser model.