Age and geochemistry of mantle peridotites and diorite dykes from the Baldissero body: Insights into the Paleozoic-Mesozoic evolution of the Southern Alps

Trace element and isotopic data obtained for mantle spinel Iherzolites and diorite dykes from the Baldissero massif (Ivrea-Verbano Zone, Western Italy) provide new, valuable constraints on the petrologic and geodynamic evolution of the Southern Alps in Paleozoic to Mesozoic times. Whole rock and mineral chemistry indicates that Baldissero Iherzolites can be regarded as refractory mantle residues following limited melt extraction. In particular, the Light Rare Earth Elements (LREE)-depleted and fractionated compositions of whole rock and clinopyroxene closely match modelling results for refractory residues after low degrees (similar to 4-5%) of near-fractional melting of depleted mantle, possibly under garnet-facies conditions. Following this, the peridotite sequence experienced subsolidus re-equilibration at lithospheric spinel-facies conditions and intrusion of several generations of dykes. However, Iherzolites far from dykes show very modest metasomatic changes, as evidenced by the crystallisation of accessory titanian pargasite and the occurrence of very slight enrichments in highly incompatible trace elements (e.g. Nb). The Re-Os data for Iherzolites far from the dykes yield a 376 Ma (Upper Devonian) model age that is considered to record a partial melting event related to the Variscan orogenic cycle s.l. Dioritic dykes cutting the mantle sequence have whole rock, clinopyroxene and plagioclase characterised by high radiogenic Nd and low radiogenic Sr, which point to a depleted to slightly enriched mantle source. Whole rock and mafic phases of diorites have high Mg# values that positively correlate with the incompatible trace element concentrations. The peridotite at the dyke contact is enriched in orthopyroxene, iron and incompatible trace elements with respect to the Iherzolites far from dykes. Numerical simulations indicate that the geochemical characteristics of the diorites can be explained by flow of a hydrous, silica-saturated melt accompanied by reaction with the ambient peridotite and fractional crystallisation. The composition of the more primitive melts calculated in equilibrium with the diorite minerals show tholeiitic to transitional affinity. Internal Sm-Nd, three-point isochrons obtained for two dykes suggest an Upper Triassic-Lower Jurassic emplacement age (from 204 31 to 198 29 Ma). Mesozoic igneous events are unknown in the southern Ivrea-Verbano Zone (IVZ), but the intrusion of hydrous melts, mostly silica-saturated, have been well documented in the Finero region, i.e. the northernmost part of IVZ and Triassic magmatism with calc-alkaline to shoshonitic affinity is abundant throughout the Central-Eastern Alps. The geochemical and chronological features of the Baldissero diorites shed new light on the geodynamic evolution of the Southern Alps before the opening of the Jurassic Tethys. (C) 2010 Elsevier B.V. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Some restriction tests in a new class of regression models for proportions

The main purpose of this work is to study the behaviour of Skovgaard`s [Skovgaard, I.M., 2001. Likelihood asymptotics. Scandinavian journal of Statistics 28, 3-32] adjusted likelihood ratio statistic in testing simple hypothesis in a new class of regression models proposed here. The proposed class of regression models considers Dirichlet distributed observations, and the parameters that index the Dirichlet distributions are related to covariates and unknown regression coefficients. This class is useful for modelling data consisting of multivariate positive observations summing to one and generalizes the beta regression model described in Vasconcellos and Cribari-Neto [Vasconcellos, K.L.P., Cribari-Neto, F., 2005. Improved maximum likelihood estimation in a new class of beta regression models. Brazilian journal of Probability and Statistics 19,13-31]. We show that, for our model, Skovgaard`s adjusted likelihood ratio statistics have a simple compact form that can be easily implemented in standard statistical software. The adjusted statistic is approximately chi-squared distributed with a high degree of accuracy. Some numerical simulations show that the modified test is more reliable in finite samples than the usual likelihood ratio procedure. An empirical application is also presented and discussed. (C) 2009 Elsevier B.V. All rights reserved.

Localization in splitting of matter waves

In this paper we present an analysis of how matter waves, guided as propagating modes in potential structures, are split under adiabatic conditions. The description is formulated in terms of localized states obtained through a unitary transformation acting on the mode functions. The mathematical framework results in coupled propagation equations that are decoupled in the asymptotic regions as well before as after the split. The resulting states have the advantage of describing propagation in situations, for instance matter-wave interferometers, where local perturbations make the transverse modes of the guiding potential unsuitable as a basis. The different regimes of validity of adiabatic propagation schemes based on localized versus delocalized basis states are also outlined. Nontrivial dynamics for superposition states propagating through split potential structures is investigated through numerical simulations. For superposition states the influence of longitudinal wave-packet extension on the localization is investigated and shown to be accurately described in quantitative terms using the adiabatic formulations presented here.

Coherence dynamics of two-mode condensates in asymmetric potentials

Detection of weak forces with an accuracy beyond the standard quantum limit holds promise both for fundamental research and for technological applications. Schemes involving ultracold atoms for such measurements are now considered to be prime candidates for increased sensitivity. In this paper we use a combination of analytical and numerical techniques to investigate the possible subshot-noise estimation of applied force fields through detection of coherence dynamics of Bose-condensed atoms in asymmetric double-well traps. Following a semiclassical description of the system dynamics and fringe visibility, we present numerical simulations of the full quantum dynamics that demonstrate the dynamical production of phase squeezing beyond the standard quantum limit. Nonlinear interactions are found to limit the achievable amount to a finite value determined by the external weak force.

Formulação LTSn para problemas de transporte sem simetria azimutal e problemas dependentes do tempo

Neste trabalho, estendemos, de forma analítica, a formulação LTSN à problemas de transporte unidimensionais sem simetria azimutal. Para este problema, também apresentamos a solução com dependência contínua na variável angular, a partir da qual é estabelecido um método iterativo de solução da equação de transporte unidimensional. Também discutimos como a formulação LTSN é aplicada na resolução de problemas de transporte unidimensionais dependentes do tempo, tanto de forma aproximada pela inversão numérica do fluxo transformado na variável tempo, bem como analiticamente, pela aplicação do método LTSNnas equações nodais. Simulações numéricas e comparações com resultados disponíveis na literatura são apresentadas.

Estudo do comportamento de um solo residual melhorado atraves de tecnicas mecanicas e fisico-quimicas e sua aplicacao a analise de fundacoes superficiais

O presente trabalho tem como objetivo estudar o comportamento de um solo residual melhorado através do uso de técnicas mecânicas de compactação e da adição de cimento. Complementarmente foram realizadas análises numéricas destes materiais tratados quando utilizados como base de fundações superficiais. O programa experimental incluiu a retirada de amostras intactas e de material amolgado para a execução de ensaios triaxiais saturados drenados com medida interna de deformações, a fim de estudar o comportamento do solo natural e do solo tratado, quer por compactação, quer por adição de cimento e compactação. Além disto, tais ensaios são determinantes na obtenção de parâmetros constitutivos para a realização de simulações numéricas. O Método dos Elementos Finitos foi utilizado para simular o comportamento carga versus recalque de placas assentes sobre o solo natural e sobre camadas de solo melhorado. O modelo Hiperbólico foi empregado na análise numérica para modelar o comportamento tensãodeformação dos materiais. Os resultados das simulações dos ensaios de placa sobre camadas de solo melhorado demonstraram que houve um aumento significativo da capacidade de suporte, além de uma redução considerável dos recalques, quando comparados ao comportamento do solo natural.

Two additions to Lucas's "inflation and welfare"

This work adds to Lucas (2000) by providing analytical solutions to two problems that are solved only numerically by the author. The first part uses a theorem in control theory (Arrow' s sufficiency theorem) to provide sufficiency conditions to characterize the optimum in a shopping-time problem where the value function need not be concave. In the original paper the optimality of the first-order condition is characterized only by means of a numerical analysis. The second part of the paper provides a closed-form solution to the general-equilibrium expression of the welfare costs of inflation when the money demand is double logarithmic. This closed-form solution allows for the precise calculation of the difference between the general-equilibrium and Bailey's partial-equilibrium estimates of the welfare losses due to inflation. Again, in Lucas's original paper, the solution to the general-equilibrium-case underlying nonlinear differential equation is done only numerically, and the posterior assertion that the general-equilibrium welfare figures cannot be distinguished from those derived using Bailey's formula rely only on numerical simulations as well.

The private memory of aggregate shocks

Em economias caracterizadas por choques agregados e privados, mostramos que a alocação ótima restrita pode depender de forma não-trivial dos choques agregados. Usando versões dos modelos de Atkeson e Lucas (1992) e Mirrlees (1971) de dois períodos, é mostrado que a alocação ótima apresenta memória com relação aos choques agregados mesmo eles sendo i.i.d. e independentes dos choques individuais, quando esses últimos choques não são totalmente persistentes. O fato de os choques terem efeitos persistentes na alocação mesmo sendo informação pública, foi primeiramente apresentado em Phelan (1994). Nossas simulações numéricas indicam que esse não é um resultado pontual: existe uma relação contínua entre persistência de tipos privados e memória do choque agregado.

Simple calvo economies with heterogeneous pricing

A motivação para este trabalho vem dos principais resultados de Carvalho e Schwartzman (2008), onde a heterogeneidade surge a partir de diferentes regras de ajuste de preço entre os setores. Os momentos setoriais da duração da rigidez nominal são su cientes para explicar certos efeitos monetários. Uma vez que concordamos que a heterogeneidade é relevante para o estudo da rigidez de preços, como poderíamos escrever um modelo com o menor número possível de setores, embora com um mínimo de heterogeneidade su ciente para produzir qualquer impacto monetário desejado, ou ainda, qualquer três momentos da duração? Para responder a esta questão, este artigo se restringe a estudar modelos com hazard-constante e considera que o efeito acumulado e a dinâmica de curto-prazo da política monetária são boas formas de se resumir grandes economias heterogêneas. Mostramos que dois setores são su cientes para resumir os efeitos acumulados de choques monetários, e economias com 3 setores são boas aproximações para a dinâmica destes efeitos. Exercícios numéricos para a dinâmica de curto prazo de uma economia com rigidez de informação mostram que aproximar 500 setores usando apenas 3 produz erros inferiores a 3%. Ou seja, se um choque monetário reduz o produto em 5%, a economia aproximada produzirá um impacto entre 4,85% e 5,15%. O mesmo vale para a dinâmica produzida por choques de nível de moeda em uma economia com rigidez de preços. Para choques na taxa de crescimento da moeda, a erro máximo por conta da aproximação é de 2,4%.

A note on learning chaotic sunspot equilibrium

In this paper we prove convergence to chaotic sunspot equilibrium through two learning rules used in the bounded rationality literature. The rst one shows the convergence of the actual dynamics generated by simple adaptive learning rules to a probability distribution that is close to the stationary measure of the sunspot equilibrium; since this stationary measure is absolutely continuous it results in a robust convergence to the stochastic equilibrium. The second one is based on the E-stability criterion for testing stability of rational expectations equilibrium, we show that the conditional probability distribution de ned by the sunspot equilibrium is expectational stable under a reasonable updating rule of this parameter. We also report some numerical simulations of the processes proposed.

A roughly smooth optimal consumption path: smoothing the rough annuity puzzle

A inconsistência entre a teoria e o comportamento empírico dos agentes no que tange ao mercado privado de pensões tem se mostrado um dos mais resistentes puzzles presentes na literatura econômica. Em modelos de otimização intertemporal de consumo e poupança sob incerteza em relação ao tempo de vida dos agentes, anuidades são ativos dominantes, anulando ou restringindo fortemente a demanda por ativos cujos retornos não estão relacionados à probabilidade de sobrevivência. Na prática, entretanto, consumidores são extremamente céticos em relação às anuidades. Em oposição ao seguro contra longevidade oferecido pelas anuidades, direitos sobre esses ativos - essencialmente ilíquidos - cessam no caso de morte do titular. Nesse sentido, choques não seguráveis de liquidez e a presença de bequest motives foram consideravelmente explorados como possíveis determinantes da baixa demanda verificada. Apesar dos esforços, o puzzle persiste. Este trabalho amplia a dominância teórica das anuidades sobre ativos não contingentes em mercados incompletos; total na ausência de bequest motives, e parcial, quando os agentes se preocupam com possíveis herdeiros. Em linha com a literatura, simulações numéricas atestam que uma parcela considerável do portfolio ótimo dos agentes seria constituída de anuidades mesmo diante de choques de liquidez, bequest motives, e preços não atuarialmente justos. Em relação a um aspecto relativamente negligenciado pela academia, mostramos que o tempo ótimo de conversão de poupança em anuidades está diretamente relacionado à curva salarial dos agentes. Finalmente, indicamos que, caso as preferências dos agentes sejam tais que o nível de consumo ótimo decaia com a idade, a demanda por anuidades torna-se bastante sensível ao sobrepreço (em relação àquele atuarialmente justo) praticado pela indústria, chegando a níveis bem mais compatíveis com a realidade empírica.

Essays on economic regulation

It is known that stock prices of public listed regulated companies react to price revisions by the regulator, and that the information conveyed by this price reaction might be used by the regulator on the contract design. This paper builds on Laffont and Tirole's (1986) regulation model with observable costs to better understand the effects the inclusion of the stock market can have on the regulator-regulated firm relationship. Our numerical simulations show that the inclusion of the market induce more powerful incentive schemes, with higher cost-reducing efforts, smaller informational rent by the firms and higher overall social welfare. In particular, we find that when the regulator is committed, the presence of short-term investors can make the first-best contract feasible, and that in the non-commitment case the market affects the firm's strategy by making it reveal more information about its cost than it normally would.

Mid-auction information acquisition

This paper studies a model of a sequential auction where bidders are allowed to acquire further information about their valuations of the object in the middle of the auction. It is shown that, in any equilibrium where the distribution of the final price is atornless, a bidder's best response has a simple characterization. In particular, the optimal information acquisition point is the same, regardless of the other bidders' actions. This makes it natural to focus on symmetric, undominated equilibria, as in the Vickrey auction. An existence theorem for such a class of equilibria is presented. The paper also presents some results and numerical simulations that compare this sequential auction with the one-shot auction. 8equential auctions typically yield more expected revenue for the seller than their one-shot counterparts. 80 the possibility of mid-auction information acquisition can provide an explanation for why sequential procedures are more often adopted.

Multistep stochastic mirror descent for risk-averse convex stochastic programs based on extended polyhedral risk measures

We consider risk-averse convex stochastic programs expressed in terms of extended polyhedral risk measures. We derive computable con dence intervals on the optimal value of such stochastic programs using the Robust Stochastic Approximation and the Stochastic Mirror Descent (SMD) algorithms. When the objective functions are uniformly convex, we also propose a multistep extension of the Stochastic Mirror Descent algorithm and obtain con dence intervals on both the optimal values and optimal solutions. Numerical simulations show that our con dence intervals are much less conservative and are quicker to compute than previously obtained con dence intervals for SMD and that the multistep Stochastic Mirror Descent algorithm can obtain a good approximate solution much quicker than its nonmultistep counterpart. Our con dence intervals are also more reliable than asymptotic con dence intervals when the sample size is not much larger than the problem size.