Dynamics and kinematics of eruptive activity at Fuego volcano, Guatemala 2005-2009

Volcanoes are the surficial expressions of complex pathways that vent magma and gasses generated deep in the Earth. Geophysical data record at least the partial history of magma and gas movement in the conduit and venting to the atmosphere. This work focuses on developing a more comprehensive understanding of explosive degassing at Fuego volcano, Guatemala through observations and analysis of geophysical data collected in 2005 – 2009. A pattern of eruptive activity was observed during 2005 – 2007 and quantified with seismic and infrasound, satellite thermal and gas measurements, and lava flow lengths. Eruptive styles are related to variable magma flux and accumulation of gas. Explosive degassing was recorded on broadband seismic and infrasound sensors in 2008 and 2009. Explosion energy partitioning between the ground and the atmosphere shows an increase in acoustic energy from 2008 to 2009, indicating a shift toward increased gas pressure in the conduit. Very-long-period (VLP) seismic signals are associated with the strongest explosions recorded in 2009 and waveform modeling in the 10 – 30 s band produces a best-fit source location 300 m west and 300 m below the summit crater. The calculated moment tensor indicates a volumetric source, which is modeled as a dike feeding a SW-dipping (35°) sill. The sill is the dominant component and its projection to the surface nearly intersects the summit crater. The deformation history of the sill is interpreted as: 1) an initial inflation due to pressurization, followed by 2) a rapid deflation as overpressure is explosively release, and finally 3) a reinflation as fresh magma flows into the sill and degasses. Tilt signals are derived from the horizontal components of the seismometer and show repetitive inflation deflation cycles with a 20 minute period coincident with strong explosions. These cycles represent the pressurization of the shallow conduit and explosive venting of overpressure that develops beneath a partially crystallized plug of magma. The energy released during the strong explosions has allowed for imaging of Fuego’s shallow conduit, which appears to have migrated west of the summit crater. In summary, Fuego is becoming more gas charged and its summit centered vent is shifting to the west - serious hazard consequences are likely.

Experimental results on gravity driven fully condensing flows in vertical tubes, their agreement with theory, and their differences with shear driven flows' boundary-condition sensitivities

This doctoral thesis presents the experimental results along with a suitable synthesis with computational/theoretical results towards development of a reliable heat transfer correlation for a specific annular condensation flow regime inside a vertical tube. For fully condensing flows of pure vapor (FC-72) inside a vertical cylindrical tube of 6.6 mm diameter and 0.7 m length, the experimental measurements are shown to yield values of average heat transfer co-efficient, and approximate length of full condensation. The experimental conditions cover: mass flux G over a range of 2.9 kg/m2-s ≤ G ≤ 87.7 kg/m2-s, temperature difference ∆T (saturation temperature at the inlet pressure minus the mean condensing surface temperature) of 5 ºC to 45 ºC, and cases for which the length of full condensation xFC is in the range of 0 < xFC < 0.7 m. The range of flow conditions over which there is good agreement (within 15%) with the theory and its modeling assumptions has been identified. Additionally, the ranges of flow conditions for which there are significant discrepancies (between 15 -30% and greater than 30%) with theory have also been identified. The paper also refers to a brief set of key experimental results with regard to sensitivity of the flow to time-varying or quasi-steady (i.e. steady in the mean) impositions of pressure at both the inlet and the outlet. The experimental results support the updated theoretical/computational results that gravity dominated condensing flows do not allow such elliptic impositions.

Modeling spontaneous three-dimensional polymerization

For human beings, the origin of life has always been an interesting and mysterious matter, particularly how life arose from inorganic matter through natural processes. Polymerization is always involved in such processes. In this paper we built what we refer to as ideal and physical models to simulate spontaneous polymerization based on certain physical principles. As the modeling confirms, without taking external energy, small and simple inorganic molecules formed bigger and more complicated molecules, which are necessary ingredients of all living organisms. In our simulations, we utilized actual ranges of parameters according to their experimentally observed values. The results from the simulations led to a good agreement with the nature of polymerization. After sorting out through all the models that were built, we arrived at a final model that, it is hoped, can be used to simply and efficiently describe spontaneous polymerization using only three parameters: the dipole moment, the distance between molecules, and the temperature.

DEVELOPMENT OF AN OPTIMIZATION MODEL FOR BIOFUEL FACILITY SIZE AND LOCATION AND A SIMULATION MODEL FOR DESIGN OF A BIOFUEL SUPPLY CHAIN

A range of societal issues have been caused by fossil fuel consumption in the transportation sector in the United States (U.S.), including health related air pollution, climate change, the dependence on imported oil, and other oil related national security concerns. Biofuels production from various lignocellulosic biomass types such as wood, forest residues, and agriculture residues have the potential to replace a substantial portion of the total fossil fuel consumption. This research focuses on locating biofuel facilities and designing the biofuel supply chain to minimize the overall cost. For this purpose an integrated methodology was proposed by combining the GIS technology with simulation and optimization modeling methods. The GIS based methodology was used as a precursor for selecting biofuel facility locations by employing a series of decision factors. The resulted candidate sites for biofuel production served as inputs for simulation and optimization modeling. As a precursor to simulation or optimization modeling, the GIS-based methodology was used to preselect potential biofuel facility locations for biofuel production from forest biomass. Candidate locations were selected based on a set of evaluation criteria, including: county boundaries, a railroad transportation network, a state/federal road transportation network, water body (rivers, lakes, etc.) dispersion, city and village dispersion, a population census, biomass production, and no co-location with co-fired power plants. The simulation and optimization models were built around key supply activities including biomass harvesting/forwarding, transportation and storage. The built onsite storage served for spring breakup period where road restrictions were in place and truck transportation on certain roads was limited. Both models were evaluated using multiple performance indicators, including cost (consisting of the delivered feedstock cost, and inventory holding cost), energy consumption, and GHG emissions. The impact of energy consumption and GHG emissions were expressed in monetary terms to keep consistent with cost. Compared with the optimization model, the simulation model represents a more dynamic look at a 20-year operation by considering the impacts associated with building inventory at the biorefinery to address the limited availability of biomass feedstock during the spring breakup period. The number of trucks required per day was estimated and the inventory level all year around was tracked. Through the exchange of information across different procedures (harvesting, transportation, and biomass feedstock processing procedures), a smooth flow of biomass from harvesting areas to a biofuel facility was implemented. The optimization model was developed to address issues related to locating multiple biofuel facilities simultaneously. The size of the potential biofuel facility is set up with an upper bound of 50 MGY and a lower bound of 30 MGY. The optimization model is a static, Mathematical Programming Language (MPL)-based application which allows for sensitivity analysis by changing inputs to evaluate different scenarios. It was found that annual biofuel demand and biomass availability impacts the optimal results of biofuel facility locations and sizes.

Meta-Environment and Executable Meta-Language using Smalltalk: an Experience Report

Object-oriented modelling languages such as EMOF are often used to specify domain specific meta-models. However, these modelling languages lack the ability to describe behavior or operational semantics. Several approaches have used a subset of Java mixed with OCL as executable meta-languages. In this experience report we show how we use Smalltalk as an executable meta-language in the context of the Moose reengineering environment. We present how we implemented EMOF and its behavioral aspects. Over the last decade we validated this approach through incrementally building a meta-described reengineering environment. Such an approach bridges the gap between a code-oriented view and a meta-model driven one. It avoids the creation of yet another language and reuses the infrastructure and run-time of the underlying implementation language. It offers an uniform way of letting developers focus on their tasks while at the same time allowing them to meta-describe their domain model. The advantage of our approach is that developers use the same tools and environment they use for their regular tasks. Still the approach is not Smalltalk specific but can be applied to language offering an introspective API such as Ruby, Python, CLOS, Java and C#.

Improving the dependability of research in personality and social psychology: Recommendations for research and educational practice

In this article, the Society for Personality and Social Psychology (SPSP) Task Force on Publication and Research Practices offers a brief statistical primer and recommendations for improving the dependability of research. Recommendations for research practice include (a) describing and addressing the choice of N (sample size) and consequent issues of statistical power, (b) reporting effect sizes and 95% confidence intervals (CIs), (c) avoiding “questionable research practices” that can inflate the probability of Type I error, (d) making available research materials necessary to replicate reported results, (e) adhering to SPSP’s data sharing policy, (f) encouraging publication of high-quality replication studies, and (g) maintaining flexibility and openness to alternative standards and methods. Recommendations for educational practice include (a) encouraging a culture of “getting it right,” (b) teaching and encouraging transparency of data reporting, (c) improving methodological instruction, and (d) modeling sound science and supporting junior researchers who seek to “get it right.”

Redistribution and attitudes toward immigrants: A European perspective

We investigate whether the dependence of immigrants on welfare benefits leads to opposition to further immigration by natives and immigrants in a pooled cross-section of 21 European countries for the 2004-2010 period. Explicitly controlling for the dependence of immigrants and natives on benefits we find that higher benefit take-up rates among immigrants than among natives lead to less favourable attitudes of natives towards immigration. Interestingly, we do not find similar stylised facts for immigrants' attitudes towards immigration.

Studies of pseudoachondroplasia chondrocytes and knockdown of the cartilage oligomeric matrix protein

Cartilage oligomeric matrix protein (COMP) is a large, homopentameric, extracellular matrix glycoprotein. Mutations in COMP cause two skeletal dysplasias: pseudoachondroplasia (PSACH) and multiple epiphyseal dysplasia (EMD1). These dwarfing conditions are caused by retention of misfolded mutant COMP with type IX collagen and matrilin-3 (MATN3) in the rough endoplasmic reticulum (rER) of the chondrocyte. These proteins form a matrix in the rER that continues to expand until it fills the entire cell, eventually causing cell death. Interestingly, loss of COMP in COMP null mice does not affect normal bone development or growth, suggesting that elimination of COMP (wildtype and mutant) expression may prevent PSACH. The hypothesis of these studies was that a hammerhead ribozyme could eliminate or knockdown COMP mRNA expression in PSACH chondrocytes . To test this hypothesis, a human chondrocyte model system that recapitulates the PSACH chondrocyte phenotype was developed by over-expressing mutant (mt-) COMP in normal chondrocytes using a recombinant adenovirus. Chondrocytes over-expressing mt-COMP developed giant rER cisternae containing COMP, type IX collagen and MATN3. Deconvolution microscopy and computer modeling showed that these proteins formed an ordered matrix surrounding a type II pro-collagen core. Additionally, the results show that a hammerhead ribozyme, ribozyme 56 (Ribo56) reduced over-expressed mt-COMP in COS cells and endogenous COMP in normal chondrocytes and mt-COMP in three PSACH chondrocytes cell line (with different mutations) by 40-70%. Altogether, these studies show that the PSACH cellular phenotype can be created in vitro and that the mt-COMP protein burden can be reduced by the presence of a COMP-specific ribozyme. Future studies will focus on designing ribozymes or short interfering RNA (siRNA) technologies that will result in better knockdown of COMP expression as well as the temporal constraints imposed by the PSACH phenotype. ^

Age determination and stable isotope record of foraminifera of sediment core SO82_5-2

Surface and deepwater paleoclimate records in Irminger Sea core SO82-5 (59°N, 31°W) and Icelandic Sea core PS2644 (68°N, 22°W) exhibit large fluctuations in thermohaline circulation (THC) from 60 to 18 calendar kyr B.P., with a dominant periodicity of 1460 years from 46 to 22 calendar kyr B.P., matching the Dansgaard-Oeschger (D-O) cycles in the Greenland Ice Sheet Project 2 (GISP2) temperature record [Grootes and Stuiver, 1997, doi:10.1029/97JC00880]. During interstadials, summer sea surface temperatures (SSTsu) in the Irminger Sea averaged to 8°C, and sea surface salinities (SSS) averaged to ~36.5, recording a strong Irminger Current and Atlantic THC. During stadials, SSTsu dropped to 2°-4°C, in phase with SSS drops by ~1-2. They reveal major meltwater injections along with the East Greenland Current, which turned off the North Atlantic deepwater convection and hence the heat advection to the north, in harmony with various ocean circulation and ice models. On the basis of the IRD composition, icebergs came from Iceland, east Greenland, and perhaps Svalbard and other northern ice sheets. However, the southward drifting icebergs were initially jammed in the Denmark Strait, reaching the Irminger Sea only with a lag of 155-195 years. We also conclude that the abrupt stadial terminations, the D-O warming events, were tied to iceberg melt via abundant seasonal sea ice and brine water formation in the meltwater-covered northwestern North Atlantic. In the 1/1460-year frequency band, benthic ?18O brine water spikes led the temperature maxima above Greenland and in the Irminger Sea by as little as 95 years. Thus abundant brine formation, which was induced by seasonal freezing of large parts of the northwestern Atlantic, may have finally entrained a current of warm surface water from the subtropics and thereby triggered the sudden reactivation of the THC. In summary, the internal dynamics of the east Greenland ice sheet may have formed the ultimate pacemaker of D-O cycles.

Geochemistry of dikes and lavas from the north wall of the Hess Deep Rift

An investigation of ~1-m.y.-old dikes and lavas from the north wall of the Hess Deep Rift (2°15'N, 101°30'W) collected during Alvin expeditions provides a detailed view of the evolution of fast spreading oceanic crust. The study area encompasses 25 km of an east-west flow line, representing ~370,000 years of crustal accretion at the East Pacific Rise. Samples analyzed exhibit depleted incompatible trace element abundances and ratios [(La/Sm)N < 1]. Indices of fractionation (MgO), and incompatible element ratios (La/Sm, Nb/Ti) show no systematic trends along flow line. Rather, over short (<4 m) and long (~25 km) distances, significant variations are observed in major and trace element concentrations and ratios. Modeling of these variations attests to the juxtaposition of dikes of distinct parental magma compositions. These findings, combined with studies of segmentation of the subaxial magma chamber and lateral magma transport in dikes along rift-dominated systems, suggest a more realistic model of the magmatic system underlying the East Pacific Rise relative to the commonly assumed twodimensional model. In this model, melts from a heterogeneous mantle feed distinct portions of a segmented axial magma reservoir. Dikes emanating from these distinct reservoirs transport magma along axis, resulting in interleaved dikes and host lavas with different evolutionary histories. This model suggests the use of axial or flow line lava compositions to infer the evolution of axial magma chambers should be approached with caution because dikes may never erupt lava or may transport magma significant distances along axis and erupt lavas far from their axial magma chamber of origin.

Chemical and Sr isotopic compositions of 3 step leachates and residues of bulk samples, zeolites, and smectites from Lower Miocene sediments at ODP Site 166-1007

Massive clinoptilolite authigenesis was observed at about 1105 meters below sea floor (mbsf) in lower Miocene wellcompacted carbonate periplatform sediments from the Great Bahama Bank [Ocean Drilling Program, ODP Leg 166, Site 1007]. The diagenetic assemblage comprises abundant zeolite crystallized within foraminifer tests and sedimentary matrix, as well as Mg smectites. In carbonate-rich deposits, the formation of the zeolite requires a supply of silica. Thus, the objective of the study is to determine the origin of the silica supply, its diagenetic evolution, and consequently the related implications on interpretation of the sedimentary record, in terms of local or global paleoceanographic change. For lack of evidence for any volcaniclastic input or traces of Si-enriched deep fluids circulation, an in situ biogenic source of silica is validated by isotopic data and chemical modeling for the formation of such secondary minerals in shallow-water carbonate sequences. Geochemical and strontium isotopic data clearly establish the marine signature of the diagenetic zeolite, as well as its contemporaneous formation with the carbonate deposition (Sr model ages of 19.6-23.2 Ma). The test of saturation for the pore fluids specifies the equilibrium state of the present mineralogical assemblage. Seawater-rock modeling specifies that clinoptilolite precipitates from the dissolution of biogenic silica, which reacts with clay minerals. The amount of silica (opal-A) involved in the reaction has to be significant enough, at least 10 wt.%, to account for the observed content of clinoptilolite occurring at the most zeolite-rich level. Modeling also shows that the observed amount of clinoptilolite (~19%) reflects an in situ and short-term reaction due to the high reactivity of primary biogenic silica (opal-A) until its complete depletion. The episodic occurrence of these well-lithified zeolite-rich levels is consistent with the occurrence of seismic reflectors, particularly the P2 seismic sequence boundary located at 1115 mbsf depth and dated as 23.2 Ma. The age range of most zeolitic sedimentary levels (biostratigraphic ages of 21.5-22 Ma) correlates well with that of the early Miocene glaciation Mi-1 and Mi-1a global events. Thus, the clinoptilolite occurrence in the shallow carbonate platform environment far from volcanogenic supply, or in other sensitive marine areas, is potentially a significant new proxy for paleoproductivity and oceanic global events, such as the Miocene events, which are usually recognized in deep-sea pelagic sediments and high latitude deposits.

Sulfur contents and S, C and O isotopic compositions in serpentinized peridotites at ODP Site 153-920

The opaque mineralogy and the contents and isotope compositions of sulfur in serpentinized peridotites from the MARK (Mid-Atlantic Ridge, Kane Fracture Zone) area were examined to understand the conditions of serpentinization and evaluate this process as a sink for seawater sulfur. The serpentinites contain a sulfur-rich secondary mineral assemblage and have high sulfur contents (up to 1 wt.%) and elevated d34S_sulfide (3.7 to 12.7?). Geochemical reaction modeling indicates that seawater-peridotite interaction at 300 to 400°C alone cannot account for both the high sulfur contents and high d34S_sulfide. These require a multistage reaction with leaching of sulfide from subjacent gabbro during higher temperature (~400°C) reactions with seawater and subsequent deposition of sulfide during serpentinization of peridotite at ~300°C. Serpentinization produces highly reducing conditions and significant amounts of H2 and results in the partial reduction of seawater carbonate to methane. The latter is documented by formation of carbonate veins enriched in 13C (up to 4.5?) at temperatures above 250°C. Although different processes produce variable sulfur isotope effects in other oceanic serpentinites, sulfur is consistently added to abyssal peridotites during serpentinization. Data for serpentinites drilled and dredged from oceanic crust and from ophiolites indicate that oceanic peridotites are a sink for up to 0.4 to 6.0 mln ton seawater S per year. This is comparable to sulfur exchange that occurs in hydrothermal systems in mafic oceanic crust at midocean ridges and on ridge flanks and amounts to 2 to 30% of the riverine sulfate source and sedimentary sulfide sink in the oceans. The high concentrations and modified isotope compositions of sulfur in serpentinites could be important for mantle metasomatism during subduction of crust generated at slow spreading rates.

Cenozoic oxygen isotopic values for benthic foraminifera from DSDP and ODP low and high latitude marine sediment cores

The climate during the Cenozoic era changed in several steps from ice-free poles and warm conditions to ice-covered poles and cold conditions. Since the 1950s, a body of information on ice volume and temperature changes has been built up predominantly on the basis of measurements of the oxygen isotopic composition of shells of benthic foraminifera collected from marine sediment cores. The statistical methodology of time series analysis has also evolved, allowing more information to be extracted from these records. Here we provide a comprehensive view of Cenozoic climate evolution by means of a coherent and systematic application of time series analytical tools to each record from a compilation spanning the interval from 4 to 61 Myr ago. We quantitatively describe several prominent features of the oxygen isotope record, taking into account the various sources of uncertainty (including measurement, proxy noise, and dating errors). The estimated transition times and amplitudes allow us to assess causal climatological-tectonic influences on the following known features of the Cenozoic oxygen isotopic record: Paleocene-Eocene Thermal Maximum, Eocene-Oligocene Transition, Oligocene-Miocene Boundary, and the Middle Miocene Climate Optimum. We further describe and causally interpret the following features: Paleocene-Eocene warming trend, the two-step, long-term Eocene cooling, and the changes within the most recent interval (Miocene-Pliocene). We review the scope and methods of constructing Cenozoic stacks of benthic oxygen isotope records and present two new latitudinal stacks, which capture besides global ice volume also bottom water temperatures at low (less than 30°) and high latitudes. This review concludes with an identification of future directions for data collection, statistical method development, and climate modeling.

(Figures S2a-c) Phase diagrams of the upper crust, the lower crust and the upper mantle based on averaged compositions of the oceanic crust and mantle

The age of the subducting Nazca Plate off Chile increases northwards from 0 Ma at the Chile Triple Junction (46°S) to 37 Ma at the latitude of Valparaíso (32°S). Age-related variations in the thermal state of the subducting plate impact on (a) the water influx to the subduction zone, as well as on (b) the volumes of water that are released under the continental forearc or, alternatively, carried beyond the arc. Southern Central Chile is an ideal setting to study this effect, because other factors for the subduction zone water budget appear constant. We determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the Chile Trench. The water release under forearc and arc is determined by coupling FEM thermal models of the subducting plate with stability fields of water-releasing mineral reactions for upper and lower crust and hydrated mantle. Results show that both the influx of water stored in, and the outflux of water released from upper crust, lower crust and mantle vary drastically over segment boundaries. In particular, the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments, but their release flux to the forearc is only about one fourth of the latter. This high variability over a subduction zone of < 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes. This article is protected by copyright. All rights reserved.