PICES-GLOBEC International Program on Climate Change and Carrying Capacity: Report of the 1999 MONITOR and REX Workshops, and the 2000 MODEL Workshop on Lower Trophic Level Modelling

Table of Contents [pdf, 0.11 Mb] Executive Summary [pdf, 0.07 Mb] MODEL Task Team Workshop Report Final Report of the International Workshop to Develop a Prototype Lower Trophic Level Ecosystem Model for Comparison of Different Marine Ecosystems in the North Pacific [pdf, 11.64 Mb] Report of the 1999 MONITOR Task Team Workshop [pdf, 0.32 Mb] Report of the 1999 REX Task Team Workshop Herring and Euphausiid population dynamics Douglas E. Hay and Bruce McCarter Spatial, temporal and life-stage variation in herring diets in British Columbia [pdf, 0.10 Mb] Augustus J. Paul and J. M. Paul Over winter changes in herring from Prince William Sound, Alaska [pdf, 0.08 Mb] N. G. Chupisheva Qualitative texture characteristic of herring (Clupea pallasi pallasi) pre-larvae developed from the natural and artificial spawning-grounds in Severnaya Bay (Peter the Great Bay) [pdf, 0.07 Mb] Gordon A. McFarlane, Richard J. Beamish and Jake SchweigertPacific herring: Common factors have opposite impacts in adjacent ecosystems [pdf, 0.15 Mb] Tokimasa Kobayashi, Keizou Yabuki, Masayoshi Sasaki and Jun-Ichi Kodama Long-term fluctuation of the catch of Pacific herring in Northern Japan [pdf, 0.39 Mb] Jacqueline M. O’Connell Holocene fish remains from Saanich Inlet, British Columbia, Canada [pdf, 0.40 Mb] Elsa R. Ivshina and Irina Y. Bragina On relationship between crustacean zooplankton (Euphausiidae and Copepods) and Sakhalin-Hokkaido herring (Tatar Strait, Sea of Japan) [pdf, 0.14 Mb] Stein Kaartvbeedt Fish predation on krill and krill antipredator behaviour [pdf, 0.08 Mb] Nikolai I. Naumenko Euphausiids and western Bering Sea herring feeding [pdf, 0.07 Mb] David M. Checkley, Jr. Interactions Between Fish and Euphausiids and Potential Relations to Climate and Recruitment [pdf, 0.08 Mb] Vladimir I. Radchenko and Elena P. Dulepova Shall we expect the Korf-Karaginsky herring migrations into the offshore western Bering Sea? [pdf, 0.75 Mb] Young Shil Kang Euphausiids in the Korean waters and its relationship with major fish resources [pdf, 0.29 Mb] William T. Peterson, Leah Feinberg and Julie Keister Ecological Zonation of euphausiids off central Oregon [pdf, 0.11 Mb] Scott M. Rumsey Environmentally forced variability in larval development and stage-structure: Implications for the recruitment of Euphausia pacifica (Hansen) in the Southern California Bight [pdf, 3.26 Mb] Scott M. Rumsey Inverse modelling of developmental parameters in Euphausia pacifica: The relative importance of spawning history and environmental forcing to larval stage-frequency distributions [pdf, 98.79 Mb] Michio J. Kishi, Hitoshi Motono & Kohji Asahi An ecosystem model with zooplankton vertical migration focused on Oyashio region [pdf, 33.32 Mb] PICES-GLOBEC Implementation Panel on Climate Change and Carrying Capacity Program Executive Committee and Task Team List [pdf, 0.05 Mb] (Document pdf contains 142 pages)

A rational fraction polynomials model to study vertical dynamic wheel-rail interaction

This paper presents a model designed to study vertical interactions between wheel and rail when the wheel moves over a rail welding. The model focuses on the spatial domain, and is drawn up in a simple fashion from track receptances. The paper obtains the receptances from a full track model in the frequency domain already developed by the authors, which includes deformation of the rail section and propagation of bending, elongation and torsional waves along an infinite track. Transformation between domains was secured by applying a modified rational fraction polynomials method. This obtains a track model with very few degrees of freedom, and thus with minimum time consumption for integration, with a good match to the original model over a sufficiently broad range of frequencies. Wheel-rail interaction is modelled on a non-linear Hertzian spring, and consideration is given to parametric excitation caused by the wheel moving over a sleeper, since this is a moving wheel model and not a moving irregularity model. The model is used to study the dynamic loads and displacements emerging at the wheel-rail contact passing over a welding defect at different speeds.

Experimental Studies on Model Scramjet in a Mach 6 High Enthalpy Free-Jet Wind Tunnel

A side-wall compression scramjet model with different combustor geometries has been tested in a propulsion tunnel that typically provides the testing flow with Mach number of 5.8, total temperature of 1800K, total pressure of 4.5MPa and mass flow rate of 4kg/s. This kerosene-fueled scramjet model consists of a side-wall compression inlet, a combustor and a thrust nozzle. A strut was used to increase the contraction ratio and to inject fuels, as well as a mixing enhancement device. Several wall cavities were also employed for flame-holding. In order to shorten the ignition delay time of the kerosene fuel, a little amount of hydrogen was used as a pilot flame. The pressure along the combustor has an evident raise after ignition occurred. Consequently thrust was observed during the fuel-on period. However, the thrust was still less than the drag of the scramjet model. For this reason, the drag variation produced by different strut and cavities was tested. Typical results showed that the cavities do not influence the drag so much, but the length of the strut does.

Discontinuous growth mechanisms of mode II crack under high-speed impact conditions

A recoverable plate impact testing technology has been used for studying the growth mechanisms of mode II crack. The results show that interactions of microcracks ahead of a crack tip cause the crack growth unsteadily. Failure mode transitions of materials were observed. Based on the observations, a discontinuous crack growth model was established. Analysis shows that the shear crack grows unsteady as the growth speed is between the Rayleigh wave speed c(R) and the shear wave speed c(s); however, when the growth speed approaches root 2c(s), the crack grows steadily. The transient microcrack growth makes the main crack speed to jump from subsonic to intersonic and the steady growth of all the sub-cracks leads the main crack to grow stably at an intersonic speed.

Determination Of Interfacial Properties Between Metal Film And Ceramic Substrate With An Adhesive Layer

Peel test measurements and inverse analysis to determine the interfacial mechanical parameters for the metal film/ceramic system are performed, considering that there exist an epoxy interface layer between film and ceramic. In the present investigation, Al films with a series of thicknesses between 20 and 250 mu m and three peel angles of 90, 135 and 180 degrees are considered. A finite element model with the cohesive zone elements is used to simulate the peel test process. The finite element results are taken as the training data of a neural network in the inverse analysis. The interfacial cohesive energy and the separation strength can be determined based on the inverse analysis and peel experimental result. (C) 2008 Elsevier Ltd. All rights reserved.

Efficient learning of decomposable models with a bounded clique size

The learning of probability distributions from data is a ubiquitous problem in the fields of Statistics and Artificial Intelligence. During the last decades several learning algorithms have been proposed to learn probability distributions based on decomposable models due to their advantageous theoretical properties. Some of these algorithms can be used to search for a maximum likelihood decomposable model with a given maximum clique size, k, which controls the complexity of the model. Unfortunately, the problem of learning a maximum likelihood decomposable model given a maximum clique size is NP-hard for k > 2. In this work, we propose a family of algorithms which approximates this problem with a computational complexity of O(k · n^2 log n) in the worst case, where n is the number of implied random variables. The structures of the decomposable models that solve the maximum likelihood problem are called maximal k-order decomposable graphs. Our proposals, called fractal trees, construct a sequence of maximal i-order decomposable graphs, for i = 2, ..., k, in k − 1 steps. At each step, the algorithms follow a divide-and-conquer strategy based on the particular features of this type of structures. Additionally, we propose a prune-and-graft procedure which transforms a maximal k-order decomposable graph into another one, increasing its likelihood. We have implemented two particular fractal tree algorithms called parallel fractal tree and sequential fractal tree. These algorithms can be considered a natural extension of Chow and Liu’s algorithm, from k = 2 to arbitrary values of k. Both algorithms have been compared against other efficient approaches in artificial and real domains, and they have shown a competitive behavior to deal with the maximum likelihood problem. Due to their low computational complexity they are especially recommended to deal with high dimensional domains.

Asymptotic scaling in the two-dimensional O(3) Nonlinear sigma model: a Monte Carlo study on parallel computers

We investigate the 2d O(3) model with the standard action by Monte Carlo simulation at couplings β up to 2.05. We measure the energy density, mass gap and susceptibility of the model, and gather high statistics on lattices of size L ≤ 1024 using the Floating Point Systems T-series vector hypercube and the Thinking Machines Corp.'s Connection Machine 2. Asymptotic scaling does not appear to set in for this action, even at β = 2.10, where the correlation length is 420. We observe a 20% difference between our estimate m/Λ^─_(Ms) = 3.52(6) at this β and the recent exact analytical result . We use the overrelaxation algorithm interleaved with Metropolis updates and show that decorrelation time scales with the correlation length and the number of overrelaxation steps per sweep. We determine its effective dynamical critical exponent to be z' = 1.079(10); thus critical slowing down is reduced significantly for this local algorithm that is vectorizable and parallelizable.

We also use the cluster Monte Carlo algorithms, which are non-local Monte Carlo update schemes which can greatly increase the efficiency of computer simulations of spin models. The major computational task in these algorithms is connected component labeling, to identify clusters of connected sites on a lattice. We have devised some new SIMD component labeling algorithms, and implemented them on the Connection Machine. We investigate their performance when applied to the cluster update of the two dimensional Ising spin model.

Finally we use a Monte Carlo Renormalization Group method to directly measure the couplings of block Hamiltonians at different blocking levels. For the usual averaging block transformation we confirm the renormalized trajectory (RT) observed by Okawa. For another improved probabilistic block transformation we find the RT, showing that it is much closer to the Standard Action. We then use this block transformation to obtain the discrete β-function of the model which we compare to the perturbative result. We do not see convergence, except when using a rescaled coupling β_E to effectively resum the series. For the latter case we see agreement for m/ Λ^─_(Ms) at , β = 2.14, 2.26, 2.38 and 2.50. To three loops m/Λ^─_(Ms) = 3.047(35) at β = 2.50, which is very close to the exact value m/ Λ^─_(Ms) = 2.943. Our last point at β = 2.62 disagrees with this estimate however.

Studies of molecular bonding, interactions and decomposition reactions on the (001) surface of ruthenium

The interactions of N₂, formic acid and acetone on the Ru(001) surface are studied using thermal desorption mass spectrometry (TDMS), electron energy loss spectroscopy (EELS), and computer modeling.

Low energy electron diffraction (LEED), EELS and TDMS were used to study chemisorption of N₂ on Ru(001). Adsorption at 75 K produces two desorption states. Adsorption at 95 K fills only the higher energy desorption state and produces a (√3 x √3)R30° LEED pattern. EEL spectra indicate both desorption states are populated by N₂ molecules bonded "on-top" of Ru atoms.

Monte Carlo simulation results are presented on Ru(001) using a kinetic lattice gas model with precursor mediated adsorption, desorption and migration. The model gives good agreement with experimental data. The island growth rate was computed using the same model and is well fit by R(t)^m - R(t₀)^m = At, with m approximately 8. The island size was determined from the width of the superlattice diffraction feature.

The techniques, algorithms and computer programs used for simulations are documented. Coordinate schemes for indexing sites on a 2-D hexagonal lattice, programs for simulation of adsorption and desorption, techniques for analysis of ordering, and computer graphics routines are discussed.

The adsorption of formic acid on Ru(001) has been studied by EELS and TDMS. Large exposures produce a molecular multilayer species. A monodentate formate, bidentate formate, and a hydroxyl species are stable intermediates in formic acid decomposition. The monodentate formate species is converted to the bidentate species by heating. Formic acid decomposition products are CO₂, CO, H₂, H₂O and oxygen adatoms. The ratio of desorbed CO with respect to CO₂ increases both with slower heating rates and with lower coverages.

The existence of two different forms of adsorbed acetone, side-on, bonded through the oxygen and acyl carbon, and end-on, bonded through the oxygen, have been verified by EELS. On Pt(111), only the end-on species is observed. On dean Ru(001) and p(2 x 2)O precovered Ru(001), both forms coexist. The side-on species is dominant on clean Ru(001), while O stabilizes the end-on form. The end-on form desorbs molecularly. Bonding geometry stability is explained by surface Lewis acidity and by comparison to organometallic coordination complexes.

I. Application of multi-Regge theory to production processes. II. High energy model for proton-proton scattering

This dissertation consists of two parts. The first part presents an explicit procedure for applying multi-Regge theory to production processes. As an illustrative example, the case of three body final states is developed in detail, both with respect to kinematics and multi-Regge dynamics. Next, the experimental consistency of the multi-Regge hypothesis is tested in a specific high energy reaction; the hypothesis is shown to provide a good qualitative fit to the data. In addition, the results demonstrate a severe suppression of double Pomeranchon exchange, and show the coupling of two "Reggeons" to an external particle to be strongly damped as the particle's mass increases. Finally, with the use of two body Regge parameters, order of magnitude estimates of the multi-Regge cross section for various reactions are given.

The second part presents a diffraction model for high energy proton-proton scattering. This model developed by Chou and Yang assumes high energy elastic scattering results from absorption of the incident wave into the many available inelastic channels, with the absorption proportional to the amount of interpenetrating hadronic matter. The assumption that the hadronic matter distribution is proportional to the charge distribution relates the scattering amplitude for pp scattering to the proton form factor. The Chou-Yang model with the empirical proton form factor as input is then applied to calculate a high energy, fixed momentum transfer limit for the scattering cross section, This limiting cross section exhibits the same "dip" or "break" structure indicated in present experiments, but falls significantly below them in magnitude. Finally, possible spin dependence is introduced through a weak spin-orbit type term which gives rather good agreement with pp polarization data.

Local and Global Externalities, Environmental Policies and Growth

33 p.

Baryon number violation beyond the standard model

This thesis describes simple extensions of the standard model with new sources of baryon number violation but no proton decay. The motivation for constructing such theories comes from the shortcomings of the standard model to explain the generation of baryon asymmetry in the universe, and from the absence of experimental evidence for proton decay. However, lack of any direct evidence for baryon number violation in general puts strong bounds on the naturalness of some of those models and favors theories with suppressed baryon number violation below the TeV scale. The initial part of the thesis concentrates on investigating models containing new scalars responsible for baryon number breaking. A model with new color sextet scalars is analyzed in more detail. Apart from generating cosmological baryon number, it gives nontrivial predictions for the neutron-antineutron oscillations, the electric dipole moment of the neutron, and neutral meson mixing. The second model discussed in the thesis contains a new scalar leptoquark. Although this model predicts mainly lepton flavor violation and a nonzero electric dipole moment of the electron, it includes, in its original form, baryon number violating nonrenormalizable dimension-five operators triggering proton decay. Imposing an appropriate discrete symmetry forbids such operators. Finally, a supersymmetric model with gauged baryon and lepton numbers is proposed. It provides a natural explanation for proton stability and predicts lepton number violating processes below the supersymmetry breaking scale, which can be tested at the Large Hadron Collider. The dark matter candidate in this model carries baryon number and can be searched for in direct detection experiments as well. The thesis is completed by constructing and briefly discussing a minimal extension of the standard model with gauged baryon, lepton, and flavor symmetries.

High pressure states in condensed matter: I. High pressure behavior of the iron-sulphur system with applications to the earth's core. II. Empirical equation of state for organic compounds at high pressures

Part I:

The earth's core is generally accepted to be composed primarily of iron, with an admixture of other elements. Because the outer core is observed not to transmit shear waves at seismic frequencies, it is known to be liquid or primarily liquid. A new equation of state is presented for liquid iron, in the form of parameters for the 4th order Birch-Murnaghan and Mie-Grüneisen equations of state. The parameters were constrained by a set of values for numerous properties compiled from the literature. A detailed theoretical model is used to constrain the P-T behavior of the heat capacity, based on recent advances in the understanding of the interatomic potentials for transition metals. At the reference pressure of 10⁵ Pa and temperature of 1811 K (the normal melting point of Fe), the parameters are: ρ = 7037 kg/m³, K_S0 = 110 GPa, K_S' = 4.53, K_S" = -.0337 GPa-1, and γ = 2.8, with γ α ρ^-1.17. Comparison of the properties predicted by this model with the earth model PREM indicates that the outer core is 8 to 10 % less dense than pure liquid Fe at the same conditions. The inner core is also found to be 3 to 5% less dense than pure liquid Fe, supporting the idea of a partially molten inner core. The density deficit of the outer core implies that the elements dissolved in the liquid Fe are predominantly of lower atomic weight than Fe. Of the candidate light elements favored by researchers, only sulfur readily dissolves into Fe at low pressure, which means that this element was almost certainly concentrated in the core at early times. New melting data are presented for FeS and FeS₂ which indicate that the FeS₂ is the S-hearing liquidus solid phase at inner core pressures. Consideration of the requirement that the inner core boundary be observable by seismological means and the freezing behavior of solutions leads to the possibility that the outer core may contain a significant fraction of solid material. It is found that convection in the outer core is not hindered if the solid particles are entrained in the fluid flow. This model for a core of Fe and S admits temperatures in the range 3450K to 4200K at the top of the core. An all liquid Fe-S outer core would require a temperature of about 4900 K at the top of the core.

Part II.

The abundance of uses for organic compounds in the modern world results in many applications in which these materials are subjected to high pressures. This leads to the desire to be able to describe the behavior of these materials under such conditions. Unfortunately, the number of compounds is much greater than the number of experimental data available for many of the important properties. In the past, one approach that has worked well is the calculation of appropriate properties by summing the contributions from the organic functional groups making up molecules of the compounds in question. A new set of group contributions for the molar volume, volume thermal expansivity, heat capacity, and the Rao function is presented for functional groups containing C, H, and O. This set is, in most cases, limited in application to low molecular liquids. A new technique for the calculation of the pressure derivative of the bulk modulus is also presented. Comparison with data indicates that the presented technique works very well for most low molecular hydrocarbon liquids and somewhat less well for oxygen-bearing compounds. A similar comparison of previous results for polymers indicates that the existing tabulations of group contributions for this class of materials is in need of revision. There is also evidence that the Rao function contributions for polymers and low molecular compounds are somewhat different.

Programming chemical kinetics: engineering dynamic reaction networks with DNA strand displacement

Over the last century, the silicon revolution has enabled us to build faster, smaller and more sophisticated computers. Today, these computers control phones, cars, satellites, assembly lines, and other electromechanical devices. Just as electrical wiring controls electromechanical devices, living organisms employ "chemical wiring" to make decisions about their environment and control physical processes. Currently, the big difference between these two substrates is that while we have the abstractions, design principles, verification and fabrication techniques in place for programming with silicon, we have no comparable understanding or expertise for programming chemistry.

In this thesis we take a small step towards the goal of learning how to systematically engineer prescribed non-equilibrium dynamical behaviors in chemical systems. We use the formalism of chemical reaction networks (CRNs), combined with mass-action kinetics, as our programming language for specifying dynamical behaviors. Leveraging the tools of nucleic acid nanotechnology (introduced in Chapter 1), we employ synthetic DNA molecules as our molecular architecture and toehold-mediated DNA strand displacement as our reaction primitive.

Abstraction, modular design and systematic fabrication can work only with well-understood and quantitatively characterized tools. Therefore, we embark on a detailed study of the "device physics" of DNA strand displacement (Chapter 2). We present a unified view of strand displacement biophysics and kinetics by studying the process at multiple levels of detail, using an intuitive model of a random walk on a 1-dimensional energy landscape, a secondary structure kinetics model with single base-pair steps, and a coarse-grained molecular model that incorporates three-dimensional geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Our findings are consistent with previously measured or inferred rates for hybridization, fraying, and branch migration, and provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

In Chapters 3 and 4, we identify and overcome the crucial experimental challenges involved in using our general DNA-based technology for engineering dynamical behaviors in the test tube. In this process, we identify important design rules that inform our choice of molecular motifs and our algorithms for designing and verifying DNA sequences for our molecular implementation. We also develop flexible molecular strategies for "tuning" our reaction rates and stoichiometries in order to compensate for unavoidable non-idealities in the molecular implementation, such as imperfectly synthesized molecules and spurious "leak" pathways that compete with desired pathways.

We successfully implement three distinct autocatalytic reactions, which we then combine into a de novo chemical oscillator. Unlike biological networks, which use sophisticated evolved molecules (like proteins) to realize such behavior, our test tube realization is the first to demonstrate that Watson-Crick base pairing interactions alone suffice for oscillatory dynamics. Since our design pipeline is general and applicable to any CRN, our experimental demonstration of a de novo chemical oscillator could enable the systematic construction of CRNs with other dynamic behaviors.

Estimation of catalyst activity profiles and deactivation parameters from reactor operating data with applications to naphtha reforming

Techniques are developed for estimating activity profiles in fixed bed reactors and catalyst deactivation parameters from operating reactor data. These techniques are applicable, in general, to most industrial catalytic processes. The catalytic reforming of naphthas is taken as a broad example to illustrate the estimation schemes and to signify the physical meaning of the kinetic parameters of the estimation equations. The work is described in two parts. Part I deals with the modeling of kinetic rate expressions and the derivation of the working equations for estimation. Part II concentrates on developing various estimation techniques.

Part I: The reactions used to describe naphtha reforming are dehydrogenation and dehydroisomerization of cycloparaffins; isomerization, dehydrocyclization and hydrocracking of paraffins; and the catalyst deactivation reactions, namely coking on alumina sites and sintering of platinum crystallites. The rate expressions for the above reactions are formulated, and the effects of transport limitations on the overall reaction rates are discussed in the appendices. Moreover, various types of interaction between the metallic and acidic active centers of reforming catalysts are discussed as characterizing the different types of reforming reactions.

Part II: In catalytic reactor operation, the activity distribution along the reactor determines the kinetics of the main reaction and is needed for predicting the effect of changes in the feed state and the operating conditions on the reactor output. In the case of a monofunctional catalyst and of bifunctional catalysts in limiting conditions, the cumulative activity is sufficient for predicting steady reactor output. The estimation of this cumulative activity can be carried out easily from measurements at the reactor exit. For a general bifunctional catalytic system, the detailed activity distribution is needed for describing the reactor operation, and some approximation must be made to obtain practicable estimation schemes. This is accomplished by parametrization techniques using measurements at a few points along the reactor. Such parametrization techniques are illustrated numerically with a simplified model of naphtha reforming.

To determine long term catalyst utilization and regeneration policies, it is necessary to estimate catalyst deactivation parameters from the the current operating data. For a first order deactivation model with a monofunctional catalyst or with a bifunctional catalyst in special limiting circumstances, analytical techniques are presented to transform the partial differential equations to ordinary differential equations which admit more feasible estimation schemes. Numerical examples include the catalytic oxidation of butene to butadiene and a simplified model of naphtha reforming. For a general bifunctional system or in the case of a monofunctional catalyst subject to general power law deactivation, the estimation can only be accomplished approximately. The basic feature of an appropriate estimation scheme involves approximating the activity profile by certain polynomials and then estimating the deactivation parameters from the integrated form of the deactivation equation by regression techniques. Different bifunctional systems must be treated by different estimation algorithms, which are illustrated by several cases of naphtha reforming with different feed or catalyst composition.

Tendência temporal do índice de massa corporal em alistados do Brasil (1980-2005)

Apesar da crescente prevalência da obesidade em países desenvolvidos e em desenvolvimento, há pouca evidência da associação com fatores ambientais. Objetivos: Investigar a evolução temporal do IMC em jovens alistados do sexo masculino de 18 anos no Brasil entre 1980 e 2005; identificar pontos específicos de maior variância na série temporal e comparar pontos específicos no tempo, a evolução temporal do IMC com as mudanças socioeconômicas no Brasil. Métodos: O presente estudo explorou uma série temporal de 26 anos em homens brasileiros que se alistaram no período de 1980 a 2005. A amostra compreendeu cerca de 35-40% de todos os jovens brasileiros de 18 anos de idade. O peso corporal e a estatura foram obtidos no momento do exame médico durante o alistamento militar. Todas as mensurações antropométricas foram realizadas por pessoal especializado e treinado. As prevalências do sobrepeso e da obesidade foram calculadas com intervalos de confiança de 95%. Com a finalidade de testar a presença de heterocedasticidade na série do IMC, realizou-se o teste de Multiplicador de Lagrange (LM). Para os pontos no tempo, com oscilações acima da média do IMC, variáveis dummies foram testadas utilizando-se o modelo ARCH (Autoregressivo de Heterocedasticidade Condicionada), com um nível de significância de p <0,05. Para aqueles pontos no tempo com oscilações acima da média do IMC (anos de 1985, 1994 e 2000), variáveis dummy foram incluídos sob a hipótese foi de que a taxa de crescimento do IMC não fosse a mesma ao longo da série temporal. Para as possíveis explicações para os aumentos bruscos na curva do IMC, foram consideradas as alterações nos principais indicadores econômicos do Brasil (Instituto Brasileiro de Geografia e Estatística e Instituto de Pesquisa Econômica Aplicada). Os fatores econômicos analisados foram: taxa de inflação anual, produção de alimentos, pobreza (%), o consumo de refrigerantes e o rendimento médio anual. Resultados: A prevalência de sobrepeso também passou de 4,5%, em 1980, para 12,5%, em 2005, um aumento de 2,6 vezes, enquanto a prevalência de obesidade aumentou de 0,5%, em 1980, para 1,9%, em 2005, um aumento de quase 300%, mas por comparação internacional estão abaixo da média. Particularmente em 1985-6 e 1994-5, houve um aumento acentuado e significativo do IMC. Em 1985-6, a média do IMC aumentou de 21,4 kg/m2 para 21,5 kg/m2 e, em 1994-5, a média do IMC médio aumentou de 21,7 kg/m2 para 21,9 kg/m2. Nesses dois pontos (1985-1986 e 1994-1995) ocorreram logo após duas grandes mudanças políticas econômicas que aumentaram o poder de compra da população. Em 1985-6, as mudanças foram principalmente relacionadas a fatores econômicos, tais como: a redução do nível de desigualdade social; aumento da renda familiar; redução da pobreza; o controle da inflação; aumento do tempo assistindo televisão e aumento do consumo de alimentos. Em 1994-5, além das mudanças no poder de compra, houve uma modificação na atividade física obrigatória nas escolas. Conclusão: O presente estudo mostrou um aumento abrupto da obesidade na população de homens jovens no Brasil em duas ocasiões durante esta série temporal (anos de 1985-6 e 1994-5), quando uma possível redução no gasto calórico e aumento do consumo de alimentos da população foram observados.