Hydrodynamical simulations of Galactic fountains - II. Evolution of multiple fountains

The ejection of gas out of the disc in late-type galaxies is related to star formation and is mainly due to the explosion of Type II supernovae (SN II). In a previous paper, we considered the evolution of a single Galactic fountain, that is, a fountain powered by a single SN cluster. Using three-dimensional hydrodynamical simulations, we studied in detail the fountain flow and its dependence with several factors, such as the Galactic rotation, the distance to the Galactic centre and the presence of a hot gaseous halo. As a natural followup, this paper investigates the dynamical evolution of multiple generations of fountains generated by similar to 100 OB associations. We have considered the observed size-frequency distribution of young stellar clusters within the Galaxy in order to appropriately fuel the multiple fountains in our simulations. Most of the results of the previous paper have been confirmed, like for example the formation of intermediate velocity clouds above the disc by the multiple fountains. Also, this work confirms the localized nature of the fountain flows: the freshly ejected metals tend to fall back close to the same Galactocentric region where they are delivered. Therefore, the fountains do not change significantly the radial profile of the disc chemical abundance. The multiple fountain simulations also allowed us to consistently calculate the feedback of the star formation on the halo gas. We found that the hot gas gains about 10 per cent of all the SN II energy produced in the disc. Thus, the SN feedback more than compensate for the halo radiative losses and allow for a quasi steady-state disc-halo circulation to exist. Finally, we have also considered the possibility of mass infall from the intergalactic medium and its interaction with the clouds that are formed by the fountains. Though our simulations are not suitable to reproduce the slow rotational pattern that is typically observed in the haloes around the disc galaxies, they indicate that the presence of an external gas infall may help to slow down the rotation of the gas in the clouds and thus the amount of angular momentum that they transfer to the coronal gas, as previously suggested in the literature.

Integrated marine gravity field in the Brazilian coast from altimeter-derived sea surface gradient and shipborne gravity

Sea surface gradients derived from the Geosat and ERS-1 satellite altimetry geodetic missions were integrated with marine gravity data from the National Geophysical Data Center and Brazilian national surveys. Using the least squares collocation method, models of free-air gravity anomaly and geoid height were calculated for the coast of Brazil with a resolution of 2` x 2`. The integration of satellite and shipborne data showed better statistical results in regions near the coast than using satellite data only, suggesting an improvement when compared to the state-of-the-art global gravity models. Furthermore, these results were obtained with considerably less input information than was used by those reference models. The least squares collocation presented a very low content of high-frequency noise in the predicted gravity anomalies. This may be considered essential to improve the high resolution representation of the gravity field in regions of ocean-continent transition. (C) 2010 Elsevier Ltd. All rights reserved.

Skull Modularity in Neotropical Marsupials and Monkeys: Size Variation and Evolutionary Constraint and Flexibility

An organism is built through a series of contingent factors, yet it is determined by historical, physical, and developmental constraints. A constraint should not be understood as an absolute obstacle to evolution, as it may also generate new possibilities for evolutionary change. Modularity is, in this context, an important way of organizing biological information and has been recognized as a central concept in evolutionary biology bridging on developmental, genetics, morphological, biochemical, and physiological studies. In this article, we explore how modularity affects the evolution of a complex system in two mammalian lineages by analyzing correlation, variance/covariance, and residual matrices (without size variation). We use the multivariate response to selection equation to simulate the behavior of Eutheria and Metharia skulls in terms of their evolutionary flexibility and constraints. We relate these results to classical approaches based on morphological integration tests based on functional/developmental hypotheses. Eutherians (Neotropical primates) showed smaller magnitudes of integration compared with Metatheria (didelphids) and also skull modules more clearly delimited. Didelphids showed higher magnitudes of integration and their modularity is strongly influenced by within-groups size variation to a degree that evolutionary responses are basically aligned with size variation. Primates still have a good portion of the total variation based on size; however, their enhanced modularization allows a broader spectrum of responses, more similar to the selection gradients applied (enhanced flexibility). Without size variation, both groups become much more similar in terms of modularity patterns and magnitudes and, consequently, in their evolutionary flexibility. J. Exp. Zool. (Mol. Dev. Evol.) 314B:663-683, 2010. (C) 2010 Wiley-Liss, Inc.

SIZE AS A LINE OF LEAST RESISTANCE II: DIRECT SELECTION ON SIZE OR CORRELATED RESPONSE DUE TO CONSTRAINTS?

Evolutionary change in New World Monkey (NWM) skulls occurred primarily along the line of least resistance defined by size (including allometric) variation (g(max)). Although the direction of evolution was aligned with this axis, it was not clear whether this macroevolutionary pattern results from the conservation of within population genetic covariance patterns (long-term constraint) or long-term selection along a size dimension, or whether both, constraints and selection, were inextricably involved. Furthermore, G-matrix stability can also be a consequence of selection, which implies that both, constraints embodied in g(max) and evolutionary changes observed on the trait averages, would be influenced by selection Here, we describe a combination of approaches that allows one to test whether any particular instance of size evolution is a correlated by-product due to constraints (g(max)) or is due to direct selection on size and apply it to NWM lineages as a case study. The approach is based on comparing the direction and amount of evolutionary change produced by two different simulated sets of net-selection gradients (beta), a size (isometric and allometric size) and a nonsize set. Using this approach it is possible to distinguish between the two hypotheses (indirect size evolution due to constraints or direct selection on size), because although both may produce an evolutionary response aligned with g(max), the amount of change produced by random selection operating through the variance/covariance patterns (constraints hypothesis) will be much smaller than that produced by selection on size (selection hypothesis). Furthermore, the alignment of simulated evolutionary changes with g(max) when selection is not on size is not as tight as when selection is actually on size, allowing a statistical test of whether a particular observed case of evolution along the line of least resistance is the result of selection along it or not. Also, with matrix diagonalization (principal components [PC]) it is possible to calculate directly the net-selection gradient on size alone (first PC [PC1]) by dividing the amount of phenotypic difference between any two populations by the amount of variation in PC1, which allows one to benchmark whether selection was on size or not

Scale architecture of the palmar and plantar epidermis of Polychrus acutirostris Spix, 1825 (Iguania, Polychrotidae) and its relationship to arboreal locomotion

In tetrapod squamates, the diversity of micro-ornamentations of the epidermis of the contact areas of hands and feet is generally associated with constraints and modalities related to locomotion. Polychrus acutirostris is a medium-sized lizard that occurs in open heterogeneous habitats in South America, such as the cerrados, caatingas, and fallow lands. It progresses slowly on branches of various diameters in its arboreal environment. It can also move more rapidly on the ground. The hands and feet are prehensile and may be considered an adaptation for grasping and climbing. Epidermal surfaces from the palmar and plantar areas of the hands and feet of P. acutirostris were prepared for SEM examination, and studied at various magnifications. They show three major levels of complexity: (1) scale types, organized in gradients of size and imbrication, (2) scalar ornamentations, organized by increasing complexity and polarity, and (3) presence of Oberhautchen showing typically iguanian honeycomb micro-ornamentations. The shape and surface structure of the scales with their pattern of micro-ornamental peaks, which improve grip, and the grasping hands and feet indicate that P. acutirostris is morpho-functionally specialized for arboreality. (C) 2009 Elsevier GmbH. All rights reserved.

Untangling associations between chironomid taxa in Neotropical streams using local and landscape filters

1. Analyses of species association have major implications for selecting indicators for freshwater biomonitoring and conservation, because they allow for the elimination of redundant information and focus on taxa that can be easily handled and identified. These analyses are particularly relevant in the debate about using speciose groups (such as the Chironomidae) as indicators in the tropics, because they require difficult and time-consuming analysis, and their responses to environmental gradients, including anthropogenic stressors, are poorly known. 2. Our objective was to show whether chironomid assemblages in Neotropical streams include clear associations of taxa and, if so, how well these associations could be explained by a set of models containing information from different spatial scales. For this, we formulated a priori models that allowed for the influence of local, landscape and spatial factors on chironomid taxon associations (CTA). These models represented biological hypotheses capable of explaining associations between chironomid taxa. For instance, CTA could be best explained by local variables (e.g. pH, conductivity and water temperature) or by processes acting at wider landscape scales (e.g. percentage of forest cover). 3. Biological data were taken from 61 streams in Southeastern Brazil, 47 of which were in well-preserved regions, and 14 of which drained areas severely affected by anthropogenic activities. We adopted a model selection procedure using Akaike`s information criterion to determine the most parsimonious models for explaining CTA. 4. Applying Kendall`s coefficient of concordance, seven genera (Tanytarsus/Caladomyia, Ablabesmyia, Parametriocnemus, Pentaneura, Nanocladius, Polypedilum and Rheotanytarsus) were identified as associated taxa. The best-supported model explained 42.6% of the total variance in the abundance of associated taxa. This model combined local and landscape environmental filters and spatial variables (which were derived from eigenfunction analysis). However, the model with local filters and spatial variables also had a good chance of being selected as the best model. 5. Standardised partial regression coefficients of local and landscape filters, including spatial variables, derived from model averaging allowed an estimation of which variables were best correlated with the abundance of associated taxa. In general, the abundance of the associated genera tended to be lower in streams characterised by a high percentage of forest cover (landscape scale), lower proportion of muddy substrata and high values of pH and conductivity (local scale). 6. Overall, our main result adds to the increasing number of studies that have indicated the importance of local and landscape variables, as well as the spatial relationships among sampling sites, for explaining aquatic insect community patterns in streams. Furthermore, our findings open new possibilities for the elimination of redundant data in the assessment of anthropogenic impacts on tropical streams.

Specific leaf areas of the tank bromeliad Guzmania monostachia perform distinct functions in response to water shortage

Leaves comprise most of the vegetative body of tank bromeliads and are usually subjected to strong longitudinal gradients. For instance, while the leaf base is in contact with the water accumulated in the tank, the more light-exposed middle and upper leaf sections have no direct access to this water reservoir. Therefore, the present study attempted to investigate whether different leaf portions of Guzmania monostachia, a tank-forming C(3)-CAM bromeliad, play distinct physiological roles in response to water shortage, which is a major abiotic constraint in the epiphytic habitat. Internal and external morphological features, relative water content, pigment composition and the degree of CAM expression were evaluated in basal, middle and apical leaf portions in order to allow the establishment of correlations between the structure and the functional importance of each leaf region. Results indicated that besides marked structural differences, a high level of functional specialization is also present along the leaves of this bromeliad. When the tank water was depleted, the abundant hydrenchyma of basal leaf portions was the main reservoir for maintaining a stable water status in the photosynthetic tissues of the apical region. In contrast, the CAM pathway was intensified specifically in the upper leaf section, which is in agreement with the presence of features more suitable for the occurrence of photosynthesis at this portion. Gas exchange data indicated that internal recycling of respiratory CO(2) accounted for virtually all nighttime acid accumulation, characterizing a typical CAM-idling pathway in the drought-exposed plants. Altogether, these data reveal a remarkable physiological complexity along the leaves of G. monostachia, which might be a key adaptation to the intermittent water supply of the epiphytic niche. (C) 2009 Elsevier GmbH. All rights reserved.

Exact Search Directions for Optimization of Linear and Nonlinear Models Based on Generalized Orthonormal Functions

A novel technique for selecting the poles of orthonormal basis functions (OBF) in Volterra models of any order is presented. It is well-known that the usual large number of parameters required to describe the Volterra kernels can be significantly reduced by representing each kernel using an appropriate basis of orthonormal functions. Such a representation results in the so-called OBF Volterra model, which has a Wiener structure consisting of a linear dynamic generated by the orthonormal basis followed by a nonlinear static mapping given by the Volterra polynomial series. Aiming at optimizing the poles that fully parameterize the orthonormal bases, the exact gradients of the outputs of the orthonormal filters with respect to their poles are computed analytically by using a back-propagation-through-time technique. The expressions relative to the Kautz basis and to generalized orthonormal bases of functions (GOBF) are addressed; the ones related to the Laguerre basis follow straightforwardly as a particular case. The main innovation here is that the dynamic nature of the OBF filters is fully considered in the gradient computations. These gradients provide exact search directions for optimizing the poles of a given orthonormal basis. Such search directions can, in turn, be used as part of an optimization procedure to locate the minimum of a cost-function that takes into account the error of estimation of the system output. The Levenberg-Marquardt algorithm is adopted here as the optimization procedure. Unlike previous related work, the proposed approach relies solely on input-output data measured from the system to be modeled, i.e., no information about the Volterra kernels is required. Examples are presented to illustrate the application of this approach to the modeling of dynamic systems, including a real magnetic levitation system with nonlinear oscillatory behavior.

Anisotropy of thermal stresses in confined dusty plasmas

Anisotropy of thermal stresses in confined dusty plasmas is considered. It is shown that in a multi-component low-temperature plasma containing electrons, ions and dust, the complicated dependence of the ion viscosity on ion temperature gradients leads to a plasma equilibrium state with anisotropic pressure. This pressure anisotropy can be of the order of the ion pressure in some limiting cases, in which the ion Larmor radius or the ion mean free path are of the order of the characteristic length of the plasma nonuniformity. For a sufficiently large dust number density, they contribute to the plasma pressure anisotropy and to its spatial dependence. Currently, it is not yet clear whether this equilibrium state is stable or not. Under these conditions, some convective plasma flows can arise in confinement devices. Therefore, this question needs special consideration.

Hyperfine interaction studies with (181)Ta and (111)Cd probes in the compound Ti(2)Ag

By using the time-differential perturbed angular correlation technique, the electric field gradients (EFG) at (181)Hf/(181)Ta and (111)In/(111)Cd probe sites in the MoSi(2)-type compound Ti(2)Ag have been measured as a function of temperature in the range from 24 to 1073 K. Ab initio EFG calculations have been performed within the framework of density functional theory using the full-potential augmented plane wave + local orbitals method as implemented in the WIEN2k package. These calculations allowed assignments of the probe lattice sites. For Ta, a single well-defined EFG with very weak temperature dependence was established and attributed to the [4(e)4mm] Ti site. For (111)Cd probes, two of the three measured EFGs are well defined and correlated with substitutional lattice sites, i.e. both the [4(e)4mm] Ti site and the [2(a)4/mmm] Ag site.

Rock magnetism in soils of the Pampean plain. Buenos Aires province, Argentina. Linking climate and magnetic behaviour

Results of the study of Argiudolls in two localities (Zarate and Veronica) of the Pampean plain, Argentina, are presented in this contribution. This is a typical area covered by loess. The magnetic studies carried out allowed to determine the presence of detrital magnetite and titanomagnetite, as well as maghemite, pedogenic goethite and superparamagnetic particles (SP). In Veronica soils, a depletion of ferromagnetic minerals is recorded The dominant process in these soils has been the reductive loss of detrital magnetite and titanomagnetite. This is associated with a greater degree of evolution of the soil, which is determined by the concentration and type of detected clays. The higher clay concentration in these soils facilitated reducing conditions and a greater loss of detrital magnetic particles. The loss is reinforced under poor drainage conditions. In the poorly drained soil of Zarate, a concentration of magnetic particles is observed in the Bt horizon, which is associated with an illuviation process. The well drained soil of the same locality shows neoformation of SP particles. These particles would have an ephemeral life until a new wet period in the annual cycle occurs. Although some characteristics of the magnetic signal appear reinforced by the conditions of drainage, this aspect does not seem to be too significant, at least in the Pampean region with low topographic gradients.

High T/P evolution and metamorphic ages of the migmatitic basement of northern Sierras Pampeanas, Argentina: Characterization of a mid-crustal segment of the Famatinian belt

New petrologic, thermobarometric and U-Pb monazite geochronologic information allowed to resolve the metamorphic evolution of a high temperature mid-crustal segment of an ancient subduction-related orogen. The EI Portezuelo Metamorphic-Igneous Complex, in the northern Sierras Pampeanas, is mainly composed of migmatites that evolved from amphibolite to granulite metamorphic facies, reaching thermal peak conditions of 670-820 degrees C and 4.5-5.3 kbar. The petrographic study combined with conventional and pseudosection thermobarometry led to deducing a short prograde metamorphic evolution within migmatite blocks. The garnet-absent migmatites represent amphibolite-facies rocks, whereas the cordierite-garnet-K-feldspar-sillimanite migmatites represent higher metamorphic grade rocks. U-Pb geochronology on monazite grains within leucosome record the time of migmatization between approximate to 477 and 470 Ma. Thus, the El Portezuelo Metamorphic-Igneous Complex is an example of exhumed Early Ordovician anatectic middle crust of the Famatinian mobile belt. Homogeneous exposure of similar paleo-depths throughout the Famatinian back-arc and isobaric cooling paths suggest slow exhumation and consequent longstanding crustal residence at high temperatures. High thermal gradients uniformly distributed in the Famatinian back-arc can be explained by shallow convection of a low-viscosity asthenosphere promoted by subducting-slab dehydration. (C) 2011 Elsevier Ltd. All rights reserved.

Thermochronology of central Ribeira Fold Belt, SE Brazil: Petrological and geochronological evidence for long-term high temperature maintenance during Western Gondwana amalgamation

The studied sector of the central Ribeira Fold Belt (SE Brazil) comprises metatexites, diatexites, charnockites and blastomylonites. This study integrates petrological and thermochronological data in order to constrain the thermotectonic and geodynamic evolution of this Neoproterozoic-Ordovician mobile belt during Western Gondwana amalgamation. New data indicate that after an earlier collision stage at similar to 610 Ma (zircon, U-Pb age), peak metamorphism and lower crust partial melting, coeval with the main regional high grade D(1) thrust deformation, occurred at 572-562 Ma (zircon, U-Pb ages). The overall average cooling rate was low (<5 degrees C/Ma) from 750 to 250 degrees C (at similar to 455 Ma; biotite-WR Rb-Sr age), but disparate cooling paths indicate differential uplift between distinct lithotypes: (a) metatexites and blastomylonites show a overall stable 3-5 degrees C/Ma cooling rate; (b) charnockites and associated rocks remained at T>650 degrees C during sub-horizontal D(2) shearing until similar to 510-470 Ma (garnet-WR Sm-Nd ages) (1-2 degrees C/Ma), being then rapidly exhumed/cooled (8-30 degrees C/Ma) during post-orogenic D(3) deformation with late granite emplacement at similar to 490 Ma (zircon, U-Pb age). Cooling rates based on garnet-biotite Fe-Mg diffusion are broadly consistent with the geochronological cooling rates: (a) metatexites were cooled faster at high temperatures (6 degrees C/Ma) and slowly at low temperatures (0.1 degrees C/Ma), decreasing cooling rates with time; (b) charnockites show low cooling rates (2 degrees C/Ma) near metamorphic peak conditions and high cooling rates (120 degrees C/Ma) at lower temperatures, increasing cooling rates during retrogression. The charnockite thermal evolution and the extensive production of granitoid melts in the area imply that high geothermal gradients were sustained fora long period of time (50-90 Ma). This thermal anomaly most likely reflects upwelling of asthenospheric mantle and magma underplating coupled with long-term generation of high HPE (heat producing elements) granitoids. These factors must have sustained elevated crustal geotherms for similar to 100 Ma, promoting widespread charnockite generation at middle to lower crustal levels. (C) 2010 Elsevier B.V. All rights reserved.

Second-order negative-curvature methods for box-constrained and general constrained optimization

A Nonlinear Programming algorithm that converges to second-order stationary points is introduced in this paper. The main tool is a second-order negative-curvature method for box-constrained minimization of a certain class of functions that do not possess continuous second derivatives. This method is used to define an Augmented Lagrangian algorithm of PHR (Powell-Hestenes-Rockafellar) type. Convergence proofs under weak constraint qualifications are given. Numerical examples showing that the new method converges to second-order stationary points in situations in which first-order methods fail are exhibited.

Boundary and internal conditions for adjoint fluid-flow problems

The ever-increasing robustness and reliability of flow-simulation methods have consolidated CFD as a major tool in virtually all branches of fluid mechanics. Traditionally, those methods have played a crucial role in the analysis of flow physics. In more recent years, though, the subject has broadened considerably, with the development of optimization and inverse design applications. Since then, the search for efficient ways to evaluate flow-sensitivity gradients has received the attention of numerous researchers. In this scenario, the adjoint method has emerged as, quite possibly, the most powerful tool for the job, which heightens the need for a clear understanding of its conceptual basis. Yet, some of its underlying aspects are still subject to debate in the literature, despite all the research that has been carried out on the method. Such is the case with the adjoint boundary and internal conditions, in particular. The present work aims to shed more light on that topic, with emphasis on the need for an internal shock condition. By following the path of previous authors, the quasi-1D Euler problem is used as a vehicle to explore those concepts. The results clearly indicate that the behavior of the adjoint solution through a shock wave ultimately depends upon the nature of the objective functional.