Data Assimilation and Inverse Methods in Terms of a Probabilistic Formulation

The weak-constraint inverse for nonlinear dynamical models is discussed and derived in terms of a probabilistic formulation. The well-known result that for Gaussian error statistics the minimum of the weak-constraint inverse is equal to the maximum-likelihood estimate is rederived. Then several methods based on ensemble statistics that can be used to find the smoother (as opposed to the filter) solution are introduced and compared to traditional methods. A strong point of the new methods is that they avoid the integration of adjoint equations, which is a complex task for real oceanographic or atmospheric applications. they also avoid iterative searches in a Hilbert space, and error estimates can be obtained without much additional computational effort. the feasibility of the new methods is illustrated in a two-layer quasigeostrophic model.

Rossby wave propagation on potential vorticity fronts with finite width

The horizontal gradient of potential vorticity (PV) across the tropopause typically declines with lead time in global numerical weather forecasts and tends towards a steady value dependent on model resolution. This paper examines how spreading the tropopause PV contrast over a broader frontal zone affects the propagation of Rossby waves. The approach taken is to analyse Rossby waves on a PV front of finite width in a simple single-layer model. The dispersion relation for linear Rossby waves on a PV front of infinitesimal width is well known; here an approximate correction is derived for the case of a finite width front, valid in the limit that the front is narrow compared to the zonal wavelength. Broadening the front causes a decrease in both the jet speed and the ability of waves to propagate upstream. The contribution of these changes to Rossby wave phase speeds cancel at leading order. At second order the decrease in jet speed dominates, meaning phase speeds are slower on broader PV fronts. This asymptotic phase speed result is shown to hold for a wide class of single-layer dynamics with a varying range of PV inversion operators. The phase speed dependence on frontal width is verified by numerical simulations and also shown to be robust at finite wave amplitude, and estimates are made for the error in Rossby wave propagation speeds due to the PV gradient error present in numerical weather forecast models.

Resonant Wave Interactions in the Equatorial Waveguide

Weakly nonlinear interactions among equatorial waves have been explored in this paper using the adiabatic version of the equatorial beta-plane primitive equations in isobaric coordinates. Assuming rigid lid vertical boundary conditions, the conditions imposed at the surface and at the top of the troposphere were expanded in a Taylor series around two isobaric surfaces in an approach similar to that used in the theory of surface-gravity waves in deep water and capillary-gravity waves. By adopting the asymptotic method of multiple time scales, the equatorial Rossby, mixed Rossby-gravity, inertio-gravity, and Kelvin waves, as well as their vertical structures, were obtained as leading-order solutions. These waves were shown to interact resonantly in a triad configuration at the O(epsilon) approximation. The resonant triads whose wave components satisfy a resonance condition for their vertical structures were found to have the most significant interactions, although this condition is not excluding, unlike the resonant conditions for the zonal wavenumbers and meridional modes. Thus, the analysis has focused on such resonant triads. In general, it was found that for these resonant triads satisfying the resonance condition in the vertical direction, the wave with the highest absolute frequency always acts as an energy source (or sink) for the remaining triad components, as usually occurs in several other physical problems in fluid dynamics. In addition, the zonally symmetric geostrophic modes act as catalyst modes for the energy exchanges between two dispersive waves in a resonant triad. The integration of the reduced asymptotic equations for a single resonant triad shows that, for the initial mode amplitudes characterizing realistic magnitudes of atmospheric flow perturbations, the modes in general exchange energy on low-frequency (intraseasonal and/or even longer) time scales, with the interaction period being dependent upon the initial mode amplitudes. Potential future applications of the present theory to the real atmosphere with the inclusion of diabatic forcing, dissipation, and a more realistic background state are also discussed.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Two essays on education and governament policies

Whether human capital increases or decreases wage uncertainty is an open ques- tion from an empirical standpoint. Yet, most policy prescriptions regarding human capital formation are based on models that impose riskiness on this type of invest- ment. We slightly deviate from the rest of the literature by allowing for non-linear income taxes in a two period model. This enables us to derive prescriptions that are robust to the risk characteristics of human capital: savings should be discouraged, human capital investments encouraged and both types of investment driven to an e¢ cient level from an aggregate perspective. These prescriptions are also robust to what choices are observed, even though the policy instruments used to implement them are not.

A transformação do modelo de gestão no futebol II

Este relatório corresponde à primeira fase do nosso segundo projeto de pesquisa intitulado "A transformação do Modelo de Gestão no Futebol" e encaminhado ao Núcleo de Pesquisas e Publicações em Junho de 2000. É a fase de pesquisa de campo em que objetivo principal é a condução de uma pesquisa exploratória sobre o processo de transformação do modelo de gestão estratégica nos clubes de futebol que se transformaram ou estão se transformando "organizacional e/ou institucionalmente" em empresas. Tendo em vista esse intuito exploratório desta pesquisa de campo, dois estudos de caso foram conduzidos - C. R. Flamengo e A. D. São Caetano - e nos permitiram entender melhor: (1) as razões para a parceria entre clube e investidor, (2) o modelo de parceria adotado e principalmente (3) as principais questões que devem ser consideradas para a concepção de um modelo de gestão profissional em um clube de futebol.

A model of capital accumulation and rent-seeking

A simple model incorporating rent-seeking into the standard neoclassical model of capital accumulation is presented. It embodies the idea that the performance of an economy depends on the efficiency of its institutions. It is shown that welfare is positively affected by the institutional efficiency, although output is not necessarily so. It is also shown that an economy with a monopolistic rent-seeker performs better than one with a competitive rent-seeking industry.

Model selection, estimation and forecasting in VAR models with short-run and long-run restrictions

We study the joint determination of the lag length, the dimension of the cointegrating space and the rank of the matrix of short-run parameters of a vector autoregressive (VAR) model using model selection criteria. We consider model selection criteria which have data-dependent penalties as well as the traditional ones. We suggest a new two-step model selection procedure which is a hybrid of traditional criteria and criteria with data-dependant penalties and we prove its consistency. Our Monte Carlo simulations measure the improvements in forecasting accuracy that can arise from the joint determination of lag-length and rank using our proposed procedure, relative to an unrestricted VAR or a cointegrated VAR estimated by the commonly used procedure of selecting the lag-length only and then testing for cointegration. Two empirical applications forecasting Brazilian inflation and U.S. macroeconomic aggregates growth rates respectively show the usefulness of the model-selection strategy proposed here. The gains in different measures of forecasting accuracy are substantial, especially for short horizons.

Model selection, estimation and forecasting in VAR models with short-run and long-run restrictions

We study the joint determination of the lag length, the dimension of the cointegrating space and the rank of the matrix of short-run parameters of a vector autoregressive (VAR) model using model selection criteria. We consider model selection criteria which have data-dependent penalties as well as the traditional ones. We suggest a new two-step model selection procedure which is a hybrid of traditional criteria and criteria with data-dependant penalties and we prove its consistency. Our Monte Carlo simulations measure the improvements in forecasting accuracy that can arise from the joint determination of lag-length and rank using our proposed procedure, relative to an unrestricted VAR or a cointegrated VAR estimated by the commonly used procedure of selecting the lag-length only and then testing for cointegration. Two empirical applications forecasting Brazilian in ation and U.S. macroeconomic aggregates growth rates respectively show the usefulness of the model-selection strategy proposed here. The gains in di¤erent measures of forecasting accuracy are substantial, especially for short horizons.

Model selection, estimation and forecasting in VAR models with short-run and long-run restrictions

We study the joint determination of the lag length, the dimension of the cointegrating space and the rank of the matrix of short-run parameters of a vector autoregressive (VAR) model using model selection criteria. We suggest a new two-step model selection procedure which is a hybrid of traditional criteria and criteria with data-dependant penalties and we prove its consistency. A Monte Carlo study explores the finite sample performance of this procedure and evaluates the forecasting accuracy of models selected by this procedure. Two empirical applications confirm the usefulness of the model selection procedure proposed here for forecasting.

A model to predict R134a refrigerant leakage through the radial clearance of rolling piston compressors

In this paper a non-isothermal two-phase model for oil-R134a refrigerant mixture flow is presented to predict the R134a leakage through the radial clearance of rolling piston compressors. The flow is divided in a liquid single-phase region and in a two-phase region, in which the homogeneous model is used to simulate the flow. The refrigerant leakage is determined using the mixture mass flow rate and the refrigerant mass fraction variation along the flow. The results are obtained for inlet pressures varying from 200 to 700 kPa, inlet temperatures ranging from 40 to 60 degrees C, and minimal clearances between 10 and 60 mu m. The results are firstly compared to existing isothermal model data, showing that there is a significant difference between the leakage flow rates predicted by isothermal and non-isothermal models. Finally, a useful general equation for compressor designers is proposed to calculate the refrigerant leakage for a large range of operation conditions. (C) 2012 Elsevier Ltd and IIR. All rights reserved.

Artificial neural networks: A novel approach to analysing the nutritional ecology of a blowfly species, Chrysomya megacephala

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Pharmacological evidence for beta(2)-adrenoceptor in right atria from stressed female rats

The purpose of the present study was to demonstrate a physiological response to TA2005, a potent beta(2)-adrenoceptor (beta(2)-AR) selective agonist, in right atria isolated from stressed female rats under the influence of the estrous cycle. We obtained concentration-response curves to the agonist in the presence and in the absence of selective antagonists in right atria isolated from female rats submitted to three daily foot-shock sessions (30 min duration, 120 pulses of 1.0 mA, 1.0 s, applied at random intervals of 5-25 s) and sacrificed at estrus or diestrus. Our results showed that the pD(2), values of TA2005 were not influenced by estrous cycle phase or foot-shock stress. However, in right atria from stressed rats sacrificed during diestrus, the concentration-response curve to TA2005 was biphasic, with a response being obtained at concentrations of 0.1 nM, whereas during estrus no response was observed at doses lower than 3 nM. ICI118,551, a beta(2)-AR antagonist, abolished the response to nanomolar concentrations of TA2005 in right atria from stressed rats at diestrus, with no changes in a,agonist pot values in right atria from control rats (7.47 +/- 0.09, p > 0.05) but a 3-fold decrease in pD(2), values of TA2005 in right atria from foot shock stressed rats (7.90 +/- 0.07, p less than or equal to 0.05). Concentration-response curves to TA2005 in the presence of ICI118,551 were best fitted by a one-site model equation. The beta(1)-AR antagonist, CGP20712A, shifted to the right only the second part of the concentration-response curves to the agonist, unmasking the putative Pz-AR-mediated response to the agonist in tissues isolated from stressed rats at diestrus. Under this condition, concentration-response curves to the agonist were best fitted by a two-site model equation. pot and maximum response of TA2005 interaction with beta(1)- and putative Pz-adrenoceptor components were calculated. Schild analyses gave a pK(B) value for CGP20712A that was typical for the interaction with beta(1)-AR in each experimental group. pK(B), values for ICI118,551 could not be obtained in stressed rats sacrificed at diestrus since Schild plot slopes were lower than 1.0. In right atria from control rats, ICI118,551 pK(B), values were similar to reported values for the interaction of the antagonist with beta(1)-AR. These results confirm that a heterogenous PAR population mediating the chronotropic response to catecholamines can be demonstrated in right atria from foot shock stressed female rats sacrificed at diestrus. The stress-induced response seems to be mediated by the beta(2)-AR subtype. Right atria from rats sacrificed during estrus are protected against stress-induced alterations on the homogeneity of PAR population.

Heat transfer in fixed bed: a model non-linearity approach

Studies have been carried out on the heat transfer in a packed bed of glass beads percolated by air at moderate flow rates. Rigorous statistic analysis of the experimental data was carried out and the traditional two parameter model was used to represent them. The parameters estimated were the effective radial thermal conductivity, k, and the wall coefficient, h, through the least squares method. The results were evaluated as to the boundary bed inlet temperature, T-o, number of terms of the solution series and number of experimental points used in the estimate. Results indicated that a small difference in T-o was sufficient to promote great modifications in the estimated parameters and in the statistical properties of the model. The use of replicas at points of high parametric information of the model improved the results, although analysis of the residuals has resulted in the rejection of this alternative. In order to evaluate cion-linearity of the model, Bates and Watts (1988) curvature measurements and the Box (1971) biases of the coefficients were calculated. The intrinsic curvatures of the model (IN) tend to be concentrated at low bed heights and those due to parameter effects (PE) are spread all over the bed. The Box biases indicated both parameters as responsible for the curvatures PE, h being somewhat more problematic. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

CHIRAL SCHWINGER MODEL WITH A PODOLSKY TERM AT FINITE-TEMPERATURE

We study an exactly solvable two-dimensional model which mimics the basic features of the standard model. This model combines chiral coupling with an infrared behavior which resembles low energy QCD. This is done by adding a Podolsky higher-order derivative term in the gauge field to the Lagrangian of the usual chiral Schwinger model. We adopt a finite temperature regularization procedure in order to calculate the non-trivial fermionic Jacobian and obtain the photon and fermion propagators, first at zero temperature and then at finite temperature in the imaginary and real time formalisms. Both singular and non-singular cases, corresponding to the choice of the regularization parameter, are treated. In the nonsingular case there is a tachyonic mode as usual in a higher order derivative theory, however in the singular case there is no tachyonic excitation in the spectrum.