Flux Transport Dynamos: From Kinematics to Dynamics

Over the past several decades, Flux-Transport Dynamo (FTD) models have emerged as a popular paradigm for explaining the cyclic nature of solar magnetic activity. Their defining characteristic is the key role played by the mean meridional circulation in transporting magnetic flux and thereby regulating the cycle period. Most FTD models also incorporate the so-called Babcock-Leighton (BL) mechanism in which the mean poloidal field is produced by the emergence and subsequent dispersal of bipolar active regions. This feature is well grounded in solar observations and provides a means for assimilating observed surface flows and fields into the models in order to forecast future solar activity, to identify model biases, and to clarify the underlying physical processes. Furthermore, interpreting historical sunspot records within the context of FTD models can potentially provide insight into why cycle features such as amplitude and duration vary and what causes extreme events such as Grand Minima. Though they are generally robust in a modeling sense and make good contact with observed cycle features, FTD models rely on input physics that is only partially constrained by observation and that neglects the subtleties of convective transport, convective field generation, and nonlinear feedbacks. Here we review the formulation and application of FTD models and assess our current understanding of the input physics based largely on complementary 3D MHD simulations of solar convection, dynamo action, and flux emergence.

Tectono-geomorphic and environmental set-up deduced during deposition of Mio-Pleistocene sediments in NW Himalaya, India

Detailed pedofacies characterization along-with lithofacies investigations of the Mio-Pleistocene Siwalik sediments exposed in the Ramnagar sub-basin have been studied so as to elucidate variability in time and space of fluvial processes and the role of intra- and extra-basinal controls on fluvial sedimentation during the evolution of the Himalayan foreland basin (HFB). Dominance of multiple, moderately to strongly developed palaeosol assemblages during deposition of Lower Siwalik (similar to 12-10.8 Ma) sediments suggest that the HFB was marked by Upland set-up of Thomas et al. (2002). Activity of intra-basinal faults on the uplands and deposition of terminal fans at different times caused the development of multiple soils. Further, detailed pedofacies along-with lithofacies studies indicate prevalence of stable tectonic conditions and development of meandering streams with broad floodplains. However, the Middle Siwalik (similar to 10.8-4.92 Ma) sub-group is marked by multistoried sandstones and minor mudstone and mainly weakly developed palaeosols, indicating deposition by large braided rivers in the form of megafans in a Lowland set-up of Thomas et al. (2002). Significant change in nature and size of rivers from the Lower to Middle Siwalik at similar to 10 Ma is found almost throughout of the basin from Kohat Plateau (Pakistan) to Nepal because the Himalayan orogeny witnessed its greatest tectonic upheaval at this time leading to attainment of great heights by the Himalaya, intensification of the monsoon, development of large rivers systems and a high rate of sedimentation, hereby a major change from the Upland set-up to the Lowland set-up over major parts of the HFB. An interesting geomorphic environmental set-up prevailed in the Ramnagar sub-basin during deposition of the studied Upper Siwalik (similar to 4.92 to <1.68 Ma) sediments as observed from the degree of pedogenesis and the type of palaeosols. In general, the Upper Siwalik sub-group in the Ramnagar sub-basin is subdivided from bottom to top into the Purmandal sandstone (4.92-4.49 Ma), Nagrota (4.49-1.68 Ma) and Boulder Conglomerate (<1.68 Ma) formations on the basis of sedimentological characters and change in dominant lithology. Presence of mudstone, a few thin gravel beds and dominant sandstone lithology with weakly to moderately developed palaeosols in the Purmandal sandstone Fm. indicates deposition by shallow braided fluvial streams. The deposition of mudstone dominant Nagrota Fm. with moderately to some well developed palaeosols and a zone of gleyed palaeosols with laminated mudstones and thin sandstones took place in an environment marked by numerous small lakes, water-logged regions and small streams in an environment just south of the Piedmont zone, perhaps similar to what is happening presently in the Upland region/the Upper Gangetic plain. This area is locally called the `Trai region' (Pascoe, 1964). Deposition of Boulder Conglomerate Fm. took place by gravelly braided river system close to the Himalayan Ranges. Activity along the Main Boundary Fault led to progradation of these environments distal-ward and led to development of on the whole a coarsening upward sequence. (C) 2014 Elsevier B.V. All rights reserved.

On failure modes and strength characterization of brittle disordered materials under uniaxial compression and tension

With a newly developed Material Failure Process Analysis code (MFPA(2D)), influence of hetero geneity on fracture processes and strength characterization of brittle disorder materials such as rock or concrete is numerically studied under uniaxial compression and tension conditions. It is found th at, due to the heterogeneity of the disordered material, relatively more diffused micro-fractures appear in the early stage of loading. Different from homogeneous materials such as glass, macro-crack nucleation starts well before the peak stress is reached and the crack propagation and coalescence can be traced, which can be taken as a precursory to predict the macro-fracture of the material. The presence of residual strength in the post-peak region and the resemblance in the stress-strain curves between tension and compression are significant results and are found to be dependent on the heterogeneity of the specimens. Examples showing the tentative applications of MFPA(2D) in modeling failure of composite materials and rock or civil engineering problem are also given in this paper.

Energy conversions and storage caused by an unsteady poloidal flow in active solar regions

In this paper we discuss coupling processes between a magnetic field and an unsteady plasma motion, and analyze the features of energy storage and conversions in active region. It is pointed out that the static force-free field is insufficient for a discussion of storage processes, and also the pure unsteady plasma rotation is not a perfect approach. In order to analyze the energy storage, we must consider the addition of poloidal plasma motion. The paper shows that because the unsteady poloidal flow is added and coupling occurs between the magnetic field and both the toroidal and the poloidal plasma flows, an unsteady process is maintained which changes the force-free factor with time. Hence, the energy in the lower levels can be transferred to the upper levels, and a considerable energy can be stored in the active region. Finally, another storage process is given which is due to the pure poloidal flow. The article shows that even if there is no twisted magnetic line of force, the energy in the lower levels may still be transferred to the upper levels and stored there.

Low streamflow in the Myakka River Basin Area in Florida

The Sarasota-Manatee area is a water-short area and the study was undertaken in 1963 in order to determine the storage capability and discharge rates of the Myakka water shed. It was found that many of the streams of the water shed were virtually dry during part of every year. However, the basins of the Myakka lakes, through which the river flows offer some storage potential, that if properly developed would provide a continuance drift of about seven million gallons of water per day of good quality water that would be high in color and temperature upon occasion. With reasonable treatment some of this water could be used to meet the present needs of the rapidly expanding coastal areas. (PDF contains 40 pages.)

Invasion of Egeria into the Hawkesbury-Nepean River, Australia

The Hawkesbury-Nepean River in New South Wales (NSW), Australia, is the largest river system in the Sydney metropolitan area, and it drains most of the developing areas to the west. This catchment is under increasing pressure from urban expansion and the river frequently experiences extended periods of low flows due to a combination of extensive river regulation and the Australian temperate climate. Added to this, the river and several of its tributaries receive treated sewage and stormwater from various sources. Habitats and biota within the Hawkesbury-Nepean River catchment have been altered since European settlement and many introduced species have spread throughout the terrestrial and aquatic environment (Recher et al. 1993). Submersed macrophyte assemblages within the river have undergone significant changes in their distribution and abundance due to eutrophication, habitat alteration and changes to river flows (Recher et al 1993). Anecdotal evidence and some early unpublished studies suggest that egeria (Egeria densa Planchon), introduced from South America as an aquarium plant, was present in the Hawkesbury-Nepean River prior to 1980. Sainty (1973) reported a persistent and troublesome infestation over a number of years at Wallacia in the upper Nepean River. Here, as part of a larger study on the ecology of macrophyte and invertebrate assemblages associated with anthropogenic disturbance in the Hawkesbury-Nepean River, we document the rapid spread of egeria since 1994. Significant increases in egeria biomass were also found, and we present preliminary evidence which suggest that the native ribbonweed, vallisneria (Vallisneria americana Michx.) is being displaced.

Seasonal water quality monitoring in the Klamath River estuary, 1991-1994

The California Department of Fish and Game's Natural Stocks Assessment Project (NSAP) collected water quality data at high tides on a monthly basis from February 1991 to October 1994, and during low tides from March 1992 to June 1994 in the Klamath River estuary to describe water quality conditions. NSAP collected data on water temperature, dissolved oxygen, salinity, depth of saltwedge, and Klamath River flow. Klamath River flows ranged from 44.5 cubic meters per second (1570 cfs) in August 1994 to 3832.2 cubic meters per second (135,315 cfs) in March 1993. Saltwater was present in the estuary primarily in the summer and early fall and generally extended 2 to 3 miles upstream. Surface water temperatures ranged from 6-8° C in the winter to 20-24° C in the summer. Summer water temperatures within the saltwedge were generally 5 to 8° C cooler than the surface water temperature. Dissolved oxygen in the estuary was generally greater than 6 to 7 ppm year-round. A sand berm formed at the mouth of the river each year in the late summer or early fall which raised the water level in the estuary and reduced tidal fluctuation so that the Klamath estuary became essentially a lagoon. I hypothesize the formation of the sand berm may increase the production of the estuary and help provide favorable conditions for rearing juvenile chinook salmon.

San Francisco Bay regional sediment management strategy development

The San Francisco Bay Conservation and Development Commission (BCDC), in continued partnership with the San Francisco Bay Long Term Management Strategies (LTMS) Agencies, is undertaking the development of a Regional Sediment Management Plan for the San Francisco Bay estuary and its watershed (estuary). Regional sediment management (RSM) is the integrated management of littoral, estuarine, and riverine sediments to achieve balanced and sustainable solutions to sediment related needs. Regional sediment management recognizes sediment as a resource. Sediment processes are important components of coastal and riverine systems that are integral to environmental and economic vitality. It relies on the context of the sediment system and forecasting the long-range effects of management actions when making local project decisions. In the San Francisco Bay estuary, the sediment system includes the Sacramento and San Joaquin delta, the bay, its local tributaries and the near shore coastal littoral cell. Sediment flows from the top of the watershed, much like water, to the coast, passing through rivers, marshes, and embayments on its way to the ocean. Like water, sediment is vital to these habitats and their inhabitants, providing nutrients and the building material for the habitat itself. When sediment erodes excessively or is impounded behind structures, the sediment system becomes imbalanced, and rivers become clogged or conversely, shorelines, wetlands and subtidal habitats erode. The sediment system continues to change in response both to natural processes and human activities such as climate change and shoreline development. Human activities that influence the sediment system include flood protection programs, watershed management, navigational dredging, aggregate mining, shoreline development, terrestrial, riverine, wetland, and subtidal habitat restoration, and beach nourishment. As observed by recent scientific analysis, the San Francisco Bay estuary system is changing from one that was sediment rich to one that is erosional. Such changes, in conjunction with increasing sea level rise due to climate change, require that the estuary sediment and sediment transport system be managed as a single unit. To better manage the system, its components, and human uses of the system, additional research and knowledge of the system is needed. Fortunately, new sediment science and modeling tools provide opportunities for a vastly improved understanding of the sediment system, predictive capabilities and analysis of potential individual and cumulative impacts of projects. As science informs management decisions, human activities and management strategies may need to be modified to protect and provide for existing and future infrastructure and ecosystem needs. (PDF contains 3 pages)

Investigation of Fundamental Processes Governing Secondary Organic Aerosol Formation in Laboratory Chambers

Our understanding of the processes and mechanisms by which secondary organic aerosol (SOA) is formed is derived from laboratory chamber studies. In the atmosphere, SOA formation is primarily driven by progressive photooxidation of SOA precursors, coupled with their gas-particle partitioning. In the chamber environment, SOA-forming vapors undergo multiple chemical and physical processes that involve production and removal via gas-phase reactions; partitioning onto suspended particles vs. particles deposited on the chamber wall; and direct deposition on the chamber wall. The main focus of this dissertation is to characterize the interactions of organic vapors with suspended particles and the chamber wall and explore how these intertwined processes in laboratory chambers govern SOA formation and evolution.

A Functional Group Oxidation Model (FGOM) that represents SOA formation and evolution in terms of the competition between functionalization and fragmentation, the extent of oxygen atom addition, and the change of volatility, is developed. The FGOM contains a set of parameters that are to be determined by fitting of the model to laboratory chamber data. The sensitivity of the model prediction to variation of the adjustable parameters allows one to assess the relative importance of various pathways involved in SOA formation.

A critical aspect of the environmental chamber is the presence of the wall, which can induce deposition of SOA-forming vapors and promote heterogeneous reactions. An experimental protocol and model framework are first developed to constrain the vapor-wall interactions. By optimal fitting the model predictions to the observed wall-induced decay profiles of 25 oxidized organic compounds, the dominant parameter governing the extent of wall deposition of a compound is identified, i.e., wall accommodation coefficient. By correlating this parameter with the molecular properties of a compound via its volatility, the wall-induced deposition rate of an organic compound can be predicted based on its carbon and oxygen numbers in the molecule.

Heterogeneous transformation of δ-hydroxycarbonyl, a major first-generation product from long-chain alkane photochemistry, is observed on the surface of particles and walls. The uniqueness of this reaction scheme is the production of substituted dihydrofuran, which is highly reactive towards ozone, OH, and NO3, thereby opening a reaction pathway that is not usually accessible to alkanes. A spectrum of highly-oxygenated products with carboxylic acid, ester, and ether functional groups is produced from the substituted dihydrofuran chemistry, thereby affecting the average oxidation state of the alkane-derived SOA.

The vapor wall loss correction is applied to several chamber-derived SOA systems generated from both anthropogenic and biogenic sources. Experimental and modeling approaches are employed to constrain the partitioning behavior of SOA-forming vapors onto suspended particles vs. chamber walls. It is demonstrated that deposition of SOA-forming vapors to the chamber wall during photooxidation experiments can lead to substantial and systematic underestimation of SOA. Therefore, it is likely that a lack of proper accounting for vapor wall losses that suppress chamber-derived SOA yields contribute substantially to the underprediction of ambient SOA concentrations in atmospheric models.

Force chains, friction, and flow: Behavior of granular media across length scales

We study the behavior of granular materials at three length scales. At the smallest length scale, the grain-scale, we study inter-particle forces and "force chains". Inter-particle forces are the natural building blocks of constitutive laws for granular materials. Force chains are a key signature of the heterogeneity of granular systems. Despite their fundamental importance for calibrating grain-scale numerical models and elucidating constitutive laws, inter-particle forces have not been fully quantified in natural granular materials. We present a numerical force inference technique for determining inter-particle forces from experimental data and apply the technique to two-dimensional and three-dimensional systems under quasi-static and dynamic load. These experiments validate the technique and provide insight into the quasi-static and dynamic behavior of granular materials.

At a larger length scale, the mesoscale, we study the emergent frictional behavior of a collection of grains. Properties of granular materials at this intermediate scale are crucial inputs for macro-scale continuum models. We derive friction laws for granular materials at the mesoscale by applying averaging techniques to grain-scale quantities. These laws portray the nature of steady-state frictional strength as a competition between steady-state dilation and grain-scale dissipation rates. The laws also directly link the rate of dilation to the non-steady-state frictional strength.

At the macro-scale, we investigate continuum modeling techniques capable of simulating the distinct solid-like, liquid-like, and gas-like behaviors exhibited by granular materials in a single computational domain. We propose a Smoothed Particle Hydrodynamics (SPH) approach for granular materials with a viscoplastic constitutive law. The constitutive law uses a rate-dependent and dilation-dependent friction law. We provide a theoretical basis for a dilation-dependent friction law using similar analysis to that performed at the mesoscale. We provide several qualitative and quantitative validations of the technique and discuss ongoing work aiming to couple the granular flow with gas and fluid flows.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

Modelagem hidrológica da bacia hidrográfica do rio Bengalas, Nova Friburgo, RJ, utilizando o potencial de geotecnologias na definição de áreas de risco de inundação.

A presente pesquisa concentra-se no estudo hidrológico utilizando o potencial das geotecnologias na modelagem do escoamento na bacia do rio Bengalas, cujo rio principal de mesmo nome, corta o município de Nova Friburgo, RJ, no sentido Sul-Norte. Esse município, um dos mais importantes da região serrana, sofre frequentemente com inundações e deslizamentos, onde, dados históricos e acontecimentos recentes mostram que a ocupação inadequada de encostas e calhas dos rios são as áreas destacada e negativamente afetadas. A metodologia tem suporte no uso de um SIG, extraindo informações, que por sua vez, serão entrada de dados na fase de modelagem, e reforçando a apresentação dos resultados das simulações através de mapas. Ela está divida basicamente em três etapas: "SIG", "Modelagem" e "Suporte à Tomada de Decisão/Simulação". Esse primeiro estudo permitiu compor um banco de dados geográfico com as características fisiográficas da bacia; a seleção criteriosa de uma modelagem matemática e encadeamento de seus parâmetros com os componentes do ciclo hidrológico; realizar a calibração do modelo de transformação chuva-vazão, Soil Conservation Service (CN); e simular a passagem dos volumes gerados pela precipitação efetiva na calha do rio Bengalas, com o objetivo de identificar e analisar as áreas suscetíveis a inundações na porção central da cidade de Nova Friburgo. Modelagem dessa natureza vem sendo empregada, principalmente, no gerenciamento de recursos hídricos, onde a tomada de decisões embasada nos resultados de simulações computacionais, contribuem para evitar prejuízos materiais e financeiros, e ainda, perdas de vidas humanas em áreas de risco, neste caso, aquelas suscetíveis a inundações. Analisando os resultados encontrados temos que a área suscetível à inundação para uma chuva com tempo de recorrência de 50 anos, o mais crítico estudado, seja de aproximadamente 1,0 km, distribuídos nos seus 8,5km na região central de Nova Friburgo-RJ, sendo está, ora delimitada, prioritariamente edificada.