Dust masses and abundances in planetary nebulae: likely strong elemental depletion in low-abundance sources

An analysis is undertaken of the relation between dust/gas mass ratios and elemental abundances within planetary nebulae (PNe). It is found that M _DUST/ M _GAS is broadly invariant with abundance, and similar to the values observed in asymptotic giant branch (AGB)-type stars. However, it is noted that the masses of dust observed in low-abundance PNe are similar to the masses of heavy elements observed in the gas phase. This is taken to imply that levels of elemental depletion must be particularly severe, and extend to many more species than have been identified so far. In particular, given that levels of C and O depletion are likely to be large, then this probably implies that species such as Fe, S, Si and Mg are depleted as well. There is already evidence for depletion of Fe, Si and Mg in individual PNe. It follows that whilst quoted abundances may accurately reflect gas-phase conditions, they are likely to be at variance with intrinsic abundances in low Z _N nebulae.
Finally, we note that there appears to be a variation in dust/gas mass ratios with galactocentric distance, with gradient similar to that observed for several elemental abundances. This may represent direct evidence for a correlation between dust/gas mass ratios and nebular abundances.     

A comparative study of geothermal and meteorological records of climate change in Kamchatka

To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski. Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing 20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature (SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response to the changes occurred on the surface.     

Compositional Zoning of the Bishop Tuff

Compositional data for >400 pumice clasts, organized accordingto eruptive sequence, crystal content, and texture, providenew perspectives on eruption and pre-eruptive evolution of the>600 km³ of zoned rhyolitic magma ejected as the Bishop Tuffduring formation of Long Valley caldera. Proportions and compositionsof different pumice types are given for each ignimbrite packageand for the intercalated plinian pumice-fall layers that eruptedsynchronously. Although withdrawal of the zoned magma was lesssystematic than previously realized, the overall sequence displaystrends toward greater proportions of less evolved pumice, morecrystals (0·5–24 wt %), and higher FeTi-oxide temperatures(714–818°C). No significant hiatus took place duringthe 6 day eruption of the Bishop Tuff, nearly all of which issuedfrom an integrated, zoned, unitary reservoir. Shortly beforeeruption, however, the zoned melt-dominant portion of the chamberwas invaded by batches of disparate lower-silica rhyolite magma,poorer in crystals than most of the resident magma but slightlyhotter and richer in Ba, Sr, and Ti. Interaction with residentmagma at the deepest levels tapped promoted growth of Ti-richrims on quartz, Ba-rich rims on sanidine, and entrapment ofnear-rim melt inclusions relatively enriched in Ba and CO₂.Varied amounts of mingling, even in higher parts of the chamber,led to the dark gray and swirly crystal-poor pumices sparselypresent in all ash-flow packages. As shown by FeTi-oxide geothermometry,the zoned rhyolitic chamber was hottest where crystal-richest,rendering any model of solidification fronts at the walls orroof unlikely. The main compositional gradient (75–195ppm Rb; 0·8–2·2 ppm Ta; 71–154 ppmZr; 0·40–1·73% FeO*) existed in the melt,prior to crystallization of the phenocryst suite observed, whichincluded zircon as much as 100 kyr older than the eruption.The compositions of crystals, though themselves largely unzoned,generally reflect magma temperature and the bulk compositionalgradient, implying both that few crystals settled or were transportedfar and that the observed crystals contributed little to establishingthat gradient. Upward increases in aqueous gas and dissolvedwater, combined with the adiabatic gradient (for the 5 km depthrange tapped) and the roofward decline in liquidus temperatureof the zoned melt, prevented significant crystallization againstthe roof, consistent with dominance of crystal-poor magma earlyin the eruption and lack of any roof-rind fragments among theBishop ejecta, before or after onset of caldera collapse. Amodel of secular incremental zoning is advanced wherein numerousbatches of crystal-poor melt were released from a mush zone(many kilometers thick) that floored the accumulating rhyoliticmelt-rich body. Each batch rose to its own appropriate levelin the melt-buoyancy gradient, which was self-sustaining againstwholesale convective re-homogenization, while the thick mushzone below buffered it against disruption by the deeper (non-rhyolitic)recharge that augmented the mush zone and thermally sustainedthe whole magma chamber. Crystal–melt fractionation wasthe dominant zoning process, but it took place not principallyin the shallow melt-rich body but mostly in the pluton-scalemush zone before and during batchwise melt extraction. KEY WORDS: Bishop Tuff; ignimbrite; magma zonation; mush model; rhyolite     

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.     

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.     

The geometry and petrogenesis of hydrothermal dolomites at Navan, Ireland

The dolomites at Navan, Ireland, formed in Courceyan peritidal and shallow-shelf limestones. The dolomite body has a plume-like geometry, cross-cutting both lithological boundaries and diagenetic barriers generated by sea-floor cementation and emergence. The dolomitizing fluids rose parallel to major faults to diffuse laterally through the succession, controlled by variations in permeability that reflect both facies variation on deposition and pre-dolomitization diagenesis. Cathodoluminescent zones reveal three principal stages of dolomite emplacement, separated by dissolution surfaces, with each stage reflecting several changes in the character of dolomitizing solutions. The predominance of dull zones indicates burial rather than surface conditions. The dolomites formed some time after burial in response to an areally limited hydrothermally-driven flow. Isotope values (σ¹⁸O of — 6σ6 to — 10.4%_δ and σ¹³C of — 0σ2 to +2σ5%_δ PDB), and fluid inclusion data, suggest that these fluids had compositions similar to those of Carboniferous seawater. However, these became hotter with time, with temperature increasing from 60 to 160δC. The Navan dolomites are closely associated with Europe's largest zinc-lead deposit. The distribution of the ores follows the same trend as that of the dolomites and paragenetic relationships indicate that dolomitization and mineralization were temporally and genetically related.     

Pillow talk

Three distinct types of pillows and pillow lava sequences with different modes of origin have been recognized in the extrusive sequences comprising the upper parts of ophiolite complexes that represent the mafic portion of the floor of an Early Cretaceous back-arc basin in southern Chile. One type of pillow formed by non-explosive submarine effusion. A second type formed by magmatic intrusion into pre-existing aquagene tuff formed by explosive eruption. The third type of pillow occurs within dikes, forming pillowed dikes, possibly as a result of vapor streaming within a cooling dike. Where studied in southern Chile, aquagene tuffs and intrusive pillows decrease and water-lain pillows increase in relative abundance from north to south. This variation corresponds with a north-to-south decrease in both the relative volume of extrusives to extensional dikes and the range and volume of differentiated rocks, suggesting a southward increase in rate of extension relative to rate of magma supply within the spreading ridges at which the ophiolites formed. In the northern part of the original basin where the rate of extension was small relative to the rate of magma supply, magma remained in magma chambers longer, resulting in a greater range and volume of differentiated rocks. The larger volume of more differentiated, cooler and more viscous magmas, in conjunction with the likelihood of more violent eruption of volatile-rich differentiates, may have been responsible for the large volume of aquagene tuff in the northern part of the original basin. These observations in southern Chile suggest that ophiolites which contain a great abundance of aquagene tuffs and intrusive pillow lavas formed in tectonic environments in which the rate of extension was small relative to the rate of magma supply (island arcs, embryonic marginal basins). Ophiolites with predominantly water-lain pillowed and massive lavas formed in tectonic environments in which the rate of extension was large relative to the rate of magma supply (mid-ocean ridges, mature back-arc basins). Thus geologic field data may supplement geochemical data as a tool in distinguishing the original igneous-tectonic environments in which ophiolites originate.     

Inverse Analysis of Deep Excavation Using Differential Evolution Algorithm

This paper presents the applications of the differential evolution (DE) algorithm in back analysis of soil parameters for deep excavation problems. A computer code, named Python‐based DE, is developed and incorporated into the commercial finite element software ABAQUS, with a parallel computing technique to run an FE analysis for all trail vectors of one generation in DE in multiple cores of a cluster, which dramatically reduces the computational time. A synthetic case and a well‐instrumented real case, that is, the Taipei National Enterprise Center (TNEC) project, are used to demonstrate the capability of the proposed back‐analysis procedure. Results show that multiple soil parameters are well identified by back analysis using a DE optimization algorithm for highly nonlinear problems. For the synthetic excavation case, the back‐analyzed parameters are basically identical to the input parameters that are used to generate synthetic response of wall deflection. For the TNEC case with a total of nine parameters to be back analyzed, the relative errors of wall deflection for the last three stages are 2.2, 1.1, and 1.0%, respectively. Robustness of the back‐estimated parameters is further illustrated by a forward prediction. The wall deflection in the subsequent stages can be satisfactorily predicted using the back‐analyzed soil parameters at early stages. Copyright © 2014 John Wiley & Sons, Ltd.