Barium-rich fluids and melts in a subduction environment (La Corea and Sierra del Convento mélanges,eastern Cuba)

Whole-rock compositions of muscovite-bearing amphibolite, trondhjemite, pegmatite and quartz-muscovite rocks from Sierra del Convento and La Corea mélanges (eastern Cuba), as well as mineralogy, record complex circulation of Ba-rich fluids and melts in the subduction environment. Partial melting of fluid-fluxed, MORB-derived amphibolite produced trondhjemite magmas that crystallized at depth, in some cases evolving into pegmatites. Qtz-Ms rocks probably crystallized from primary fluids derived from subducted sediments. All these rocks have elevated concentrations of large-ion lithophile elements, especially Ba (up to 11,810 ppm), presumably released from slab sediments by fluids and/or melts. Fluid–rock interaction produced crystallization of phengite in parental amphibolites. The phengite crystallized in all types of rocks is rich in Ba, with concentric zoning, characterized by Ba-rich cores and Ba-poor rims, indicating a compatible behaviour of Ba in the studied systems. Zoning in phengite is governed primarily by the celadonite (tschermak) exchange vector ((Mg,Fe)Si-(Al)_-2), with more moderate contributions of the celsian (BaAl-(KSi)_-1) and paragonite (NaK_-1) exchange vectors. Late remobilization of Ba at relatively low temperature formed retrograde celsian. The compatible behaviour of Ba in the studied rocks strengthens the importance of the stability of phengite for the transfer of LIL elements from the subduction to the volcanic arc environments.     

Characterizing the severe turbulence environments associated with commercial aviation accidents. A real-time turbulence model (RTTM) designed for the operational prediction of hazardous aviation turbulence environments

Summary In this paper, we will focus on the real-time prediction of environments that are predisposed to producing moderate-severe (hazardous) aviation turbulence. We will describe the numerical model and its postprocessing system that is designed for said prediction of environments predisposed to severe aviation turbulence as well as presenting numerous examples of its utility. The purpose of this paper is to demonstrate that simple hydrostatic precursor circulations organize regions of preferred wave breaking and turbulence at the nonhydrostatic scales of motion. This will be demonstrated with a hydrostatic numerical modeling system, which can be run in real time on a very inexpensive university computer workstation employing simple forecast indices. The forecast system is designed to efficiently support forecasters who are directing research aircraft to measure the environment immediately surrounding turbulence. The numerical model is MASS version 5.13, which is integrated over three different grid matrices in real-time on a university workstation in support of NASA-Langley’s B-757 turbulence research flight missions. The model horizontal resolutions are 60, 30, and 15 km and the grids are centered over the region of operational NASA-Langley B-757 turbulence flight missions. The postprocessing system includes several turbulence-related products including four turbulence forecasting indices, winds, streamlines, turbulence kinetic energy, and Richardson numbers. Additionally there are convective products including precipitation, cloud height, cloud mass fluxes, lifted index, and K-index. Furthermore, soundings, sounding parameters, and Froude number plots are also provided. The horizontal cross section plot products are provided from 16,000–46,000 feet in 2,000 feet intervals. Products are available every three hours at the 60 and 30 km grid interval and every 1.5 hours at the 15 km grid interval. The model is initialized from the NWS ETA analyses and integrated two times a day.     

The rms–flux relations in different branches in Cyg X-2

In this paper, the rms-flux (root mean square-flux) relation along the Z-track of the bright Z-Source Cyg X-2 is analyzed using the observational data of Rossi X-ray Timing Explorer (RXTE). Three types of rms-flux relations, i.e. positive, negative, and ‘arch’-like correlations are found in different branches. The rms is positively correlated with flux in normal branch (NB), but anti-correlated in the vertical horizontal branch (VHB). The rms–flux relation shows an ‘arch’-like shape in the horizontal branch (HB). We also try to explain this phenomenon using existing models.     

Functional responses for zooplankton feeding on multiple resources: a review of assumptions and biological dynamics

Modelers often need to quantify the rates at which zooplankton consume a variety of species, size classes and trophic types. Implicit in the equations used to describe the multiple resource functional response (i.e. how nutritional intake varies with resource densities) are assumptions that are not often stated, let alone tested. This is problematic because models are sensitive to the details of these formulations. Here, we enable modelers to make more informed decisions by providing them with a new framework for considering zooplankton feeding on multiple resources. We define a new classification of multiple resource responses that is based on preference, selection and switching, and we develop a set of mathematical diagnostics that elucidate model assumptions. We use these tools to evaluate the assumptions and biological dynamics inherent in published multiple resource responses. These models are shown to simulate different resource preferences, implied single resource responses, changes in intake with changing resource densities, nutritional benefits of generalism, and nutritional costs of selection. Certain formulations are further shown to exhibit anomalous dynamics such as negative switching and sub-optimal feeding. Such varied responses can have vastly different ecological consequences for both zooplankton and their resources; inappropriate choices may incorrectly quantify biologically-mediated fluxes and predict spurious dynamics. We discuss how our classes and diagnostics can help constrain parameters, interpret behaviors, and identify limitations to a formulation's applicability for both regional (e.g. High-Nitrate-Low-Chlorophyll regions comprising large areas of the Pacific) and large-scale applications (e.g. global biogeochemical or climate change models). Strategies for assessing uncertainty and for using the mathematics to guide future experimental investigations are also discussed.     

Geochemical characteristics and LA-ICP-MS zircon U-Pb age of the Lanhualing granite in Ningguo,Anhui ProvinceEI北大核心CSCD

The Lanhualing tungsten-molybdenum deposit is a skarn-type deposit located in Ningguo county, Anhui province. This deposit is mainly hosted in the Yinzhubu Formation and the Yanwashan Formation of Ordovician, and genetically related to the Lanhualing granite. The Lanhualing granite belongs to high-K calc-alkaline series with high alkali (Na2O+K2O=4.00%-7.03%), SiO2 (67.87%-74.92%) and MgO (0.62%-1.23%) contents. The granitic rocks show right-dipping chondrite normalized REE patterns with weak δEu anomalies. The granitic rocks are relatively enriched in large ion lithophile elements (LILE) and depleted in high field strength elements (HFSE). The ore-bearing granite was dated at 148.17±0.94 Ma by LA-ICP-MS zircon U-Pb method. The late Yanshanian is the main tungsten mineralization epoch in the South Anhui-north Jiangxi area; and indeed, the Dongyuan, Zhuxi, Yangchuling, Dahutang and other large and super-large tungsten deposits were formed in this period. Geochemical Characteristics of the Lanhualing granite indicate a crustal source but with mantle input under tectonic regime of compression thickening. ©, 2015, Science Press. All right reserved.     

Developments in Jovian Radio Emissions Tomography and Observations Techniques

Jupiter radio emission is known to be the most powerful nonthermal planetary radiation. In recent years specifically space-based observations allow us to permanently cover a large frequency band(from 100 kHz up to 40 MHz combined with ground-based telescopes)of the Jovian spectrum. The Plasma and Wave Science experiment onboard Galileo enables the observation of Jovian kilometric and hectometric emissions; Wind/WAVES and ground-based telescopes (mainly Decametric Array in Nancay, France, and UTR-2 in Kharkov, Ukraine) cover also hectometric and mainly decametric emissions. Specific geometrical configurations between Cassini approaching Jupiter and Wind spacecraft orbiting Earth, with Galileo orbiting Jupiter and Wind, in combination with ground-based observations provide a new approach to perform Jovian radio tomography. The tomography technique is used to analyze ray paths of Jovian radio emission observed in different directions (e.g. solar and anti-solar direction) and for different declination of Earth. The developments of Jovian radio emission tomography in recent years treated refraction effects and its connection to the local magnetic field in the radio source as well as the radio wave propagation through the Io torus and the terrestrial ionosphere. Most recently ground-based multi-site and simultaneous Jupiter decametric radio observations by means of digital spectropolarimeter and waveform receiver provide the basis of a new data analysis treatment. The above addressed topics are without exemption deeply connected to the plasma structures the radio waves are generated in and propagating through. This revised version was published online in July 2006 with corrections to the Cover Date.     

H2O activity in concentrated NaCl solutions at high pressures and temperatures measured by the brucite-periclase equilibrium

 H₂O activities in concentrated NaCl solutions were measured in the ranges 600°–900° C and 2–15 kbar and at NaCl concentrations up to halite saturation by depression of the brucite (Mg(OH)₂) – periclase (MgO) dehydration equilibrium. Experiments were made in internally heated Ar pressure apparatus at 2 and 4.2 kbar and in 1.91-cm-diameter piston-cylinder apparatus with NaCl pressure medium at 4.2, 7, 10 and 15 kbar. Fluid compositions in equilibrium with brucite and periclase were reversed to closures of less than 2 mol% by measuring weight changes after drying of punctured Pt capsules. Brucite-periclase equilibrium in the binary system was redetermined using coarsely crystalline synthetic brucite and periclase to inhibit back-reaction in quenching. These data lead to a linear expression for the standard Gibbs free energy of the brucite dehydration reaction in the experimental temperature range: ΔG° (±120J)=73418–134.95T(K). Using this function as a baseline, the experimental dehydration points in the system MgO−H₂O−NaCl lead to a simple systematic relationship of high-temperature H₂O activity in NaCl solution. At low pressure and low fluid densities near 2 kbar the H₂O activity is closely approximated by its mole fraction. At pressures of 10 kbar and greater, with fluid densities approaching those of condensed H₂O, the H₂O activity becomes nearly equal to the square of its mole fraction. Isobaric halite saturation points terminating the univariant brucite-periclase curves were determined at each experimental pressure. The five temperature-composition points in the system NaCl−H₂O are in close agreement with the halite saturation curves (liquidus curves) given by existing data from differential thermal analysis to 6 kbar. Solubility of MgO in the vapor phase near halite saturation is much less than one mole percent and could not have influenced our determinations. Activity concentration relations in the experimental P-T range may be retrieved for the binary system H₂O-NaCl from our brucite-periclase data and from halite liquidus data with minor extrapolation. At two kbar, solutions closely approach an ideal gas mixture, whereas at 10 kbar and above the solutions closely approximate an ideal fused salt mixture, where the activities of H₂O and NaCl correspond to an ideal activity formulation. This profound pressure-induced change of state may be characterized by the activity (a) – concentration (X) expression: a _{H 2}O=X _{H 2}O/(1+αX _NaCl), and a _NaCl=(1+α)^(1+α)[X _NaCl/(1+αX _NaCl)]^(1+α). The parameter α is determined by regression of the brucite-periclase H₂O activity data: α=exp[A–B/ϱ_{H 2}O ]-CP/T, where A=4.226, B=2.9605, C=164.984, and P is in kbar, T is in Kelvins, and ϱ_{H 2}O is the density of H₂O at given P and T in g/cm³. These formulas reproduce both the H₂O activity data and the NaCl activity data with a standard deviation of ±0.010. The thermodynamic behavior of concentrated NaCl solutions at high temperature and pressure is thus much simpler than portrayed by extended Debye-Hückel theory. The low H₂O activity at high pressures in concentrated supercritical NaCl solutions (or hydrosaline melts) indicates that such solutions should be feasible as chemically active fluids capable of coexisting with solid rocks and silicate liquids (and a CO₂-rich vapor) in many processes of deep crustal and upper mantle metamorphism and metasomatism. Received: 1 September 1995 / Accepted: 24 March 1996     

Radiation pressure effects in the oscillations of compressible rotating homogeneous spheroids

Earlier models of compressible, rotating, and homogeneous ellipsoids with gas pressure are generalized to include the presence of radiation pressure. Under the assumptions of a linear velocity field of the fluid and a bounded ellipsoidal surface, the dynamical behaviour of these models can be described by ordinary differential equations. These equations are used to study the finite oscillations of massive radiative models with masses 10M and 30M in which the effects of radiation pressure are expected to be important.Models with two different degrees of equilibrium are chosen: an equilibrium (i.e., dynamically stable) model with an initial asymmetric inward velocity, and a nonequilibrium model with a nonequilibrium central temperature and which falls inwards from rest. For each of these two degrees of equilibrium, two initial configurations are considered: rotating spheroidal and nonrotating spherical models.From the numerical integration of the differential equations for these models, we obtain the time evolution of their principal semi-diametersa ₁ anda ₃, and of their central temperatures, which are graphically displayed by making plots of the trajectories in the (a ₁,a ₃) phase space, and of botha ₁ and the total central pressureP _c against time.It is found that in all the equilibrium radiative models (in which radiation pressure is taken into account), the periods of the oscillations of botha ₁ andP _c are longer than those of the corresponding nonradiative models, while the reverse is true for the nonequilibrium radiative models. The envelopes of thea ₁ oscillations of the equilibrium radiative models also have much longer periods; this result also holds for the nonequilibrium models whenever the envelope is well defined. Further, as compared to the nonradiative models, almost all the radiative models collapse to smaller volumes before rebouncing, with the more massive model undergoing a larger collapse and attaining a correspondingly larger peakP _c.When the mass is increased, the dynamical behavior of the radiative model generally becomes more nonperiodic. The ratio of the central radiation pressure to the central gas pressure, which is small for low mass models, increases with mass, and at the center of the more massive model, the radiation pressure can be comparable in magnitude to the gas pressure. In all the radiative models, the average periods as well as the average amplitudes of both thea ₁ andP _c oscillations also increase with mass.When either rotation or radiation pressure effects or both are included in the equilibrium nonradiative model, the period of the envelope of thea ₁ oscillations is increased. The presence of rotation in the equilibrium radiative model, however, decreases this period.Some astrophysical implications of this work are briefly discussed.     

Spectral Variability of FSRQs

The optical variability of 29 flat spectrum radio quasars in SDSS Stripe 82 region are investigated by using DR7 released multi-epoch data. All FSRQs show variations with overall amplitude ranging from 0.24 mag to 3.46 mag in different sources. About half of FSRQs show a bluer-when-brighter trend, which is commonly observed for blazars. However, only one source shows a redder-when-brighter trend, which implies it is rare in FSRQs. In this source, the thermal emission may be responsible for the spectral behaviour.     

Lunar gravity field determination using SELENE same-beam differential VLBI tracking data

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from 1-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080±0.0009 km³/s² (10 sigma). The lunar degree 2 potential Love number k ₂ was also estimated, and has a value of 0.0255 ± 0.0016 (10 sigma as well).