two dimensional boundary fitted circulation model in spherical coordinates a two dimensional vertically averaged unsteady circulation model using a nonorthogonal boundary fitted technique was developed in spherical coordinates for predicting sea level and currents in estuarine and shelf waters both the dependent and independent variables are transformed into a curvilinear coordinate system the governing equations are solved by a semiimplicit method in which the elevations are solved implicitly and the vertically averaged velocities are solved explicitly the model employs a space staggered grid system and a three level time discretization truncation errors are second order both in space and time the model was tested against analytic solutions for a standing wave in a closed basin tidal circulation in a simple rectangular channel with an irregular grid system and various degrees of rotation and tidal flow in an annular section channel with quadratic bottom topography the model was also tested against steady state wind induced setup in a closed irregular basin with cons

advances in fourier transform infrared spectroscopy of natural glasses from sample preparation to data analysis fourier transform infrared spectroscopy ftir is an analytical technique utilized to measure the concentrations of h and c species in volcanic glasses water and co2 are the most abundant volatile species in volcanic systems water is present in magmas in higher concentrations than co2 and is also more soluble at lower pressures and therefore it is the dominant volatile forming bubbles during volcanic eruptions dissolved water affects both phase equilibria and melt physical properties such as density and viscosity therefore water is important for understanding magmatic processes additionally quantitative measurements of different volatile species using ftir can be achieved at high spatial resolution recent developments in analytical equipment such as synchrotron light sources and the development of focal plane array fpa detectors allow higher resolution measurements and the acquisition of concentration maps these new capabilities are being used to characterize spatial gradients or lack thereof around bubbles and other textural features which in turn lead to new insights into the behavior of volcanic feeder systems here practical insights about sample preparation and analysis of the distribution and speciation of volatiles in volcanic glasses using ftir spectroscopy are discussed new advances in the field of ftir analysis produce reliable data at high spatial resolution that can be used to produce datasets on the distribution dissolution and diffusion of volatiles in volcanic materials

facies distribution of ejecta in analog tephra rings from experiments with single and multiple subsurface explosions the volume grain size and depositional facies of material deposited outside an explosion crater ejecta are sensitive to the depth of the explosion the explosion energy and the presence or absence of a crater before the explosion we detonate buried chemical explosives as an analog for discrete volcanic explosions in experiments to identify unique characteristics of proximal medial and distal ejecta facies and their distribution from a range of scaled depths in undisturbed and cratered ground ejecta are here discussed in terms of three facies 1 proximal ejecta which form a constructional landform around a crater 2 medial ejecta which form a continuous sheet deposit that thins much more gradually with distance and 3 distal ejecta that are deposited as isolated clasts the extent of proximal ejecta away from the crater relative to crater size is not sensitive to scaled depth but the volume proportion of proximal ejecta to the total ejecta deposit is sensitive to the presence of a crater and scaled depth medial ejecta distribution and volume contributions are both sensitive to the presence of a crater and to scaled depth distal ejecta distance is dependent on scaled depth and the presence of a crater while the volume proportion of distal ejecta is less sensitive to scaled depth or presence of a crater experimental facies and deposit structures inferred from observations of jet dynamics are used to suggest facies associations anticipated from eruptions dominated by explosions of different scaled depth configurations we emphasize that significant differences in tephra ring deposits can result from the effects of scaled depth and preexisting craters on ejecta dynamics and are not necessarily related to fundamental differences in explosion mechanisms or degree of magma fragmentation

modeling the tide and wind induced circulation in buzzards bay abstract hydrodynamic model application to buzzards bay is performed using a three dimensional boundary fitted hydrodynamic model in this study the model is forced with observed tidal harmonic constants along the open boundaries and winds on the surface the main focus of the present study is to model the detailed wind and tide induced circulation in buzzards bay the observed surface elevations and currents given in butman b signell r shoukimas p beardsley r c 1988 current observations in buzzards bay 1982 1986 open file report 88 5 united states geological survey and the tide and current harmonics given in signell r p 1987 tide and wind forced currents in buzzards bay massachusetts technical report wh 87 15 woods hole oceanographic institution woods hole massachusetts are used to validate the model predictions the calibrated model is then used to study the relative contributions of tidal and wind forcing on the instantaneous and residual circulation in buzzards bay the amplitudes and phases of the principal tidal constituents at 10 tidal stations in buzzards bay obtained from a harmonic analysis of a 60 day simulation compare well with the observed data the predicted amplitude and phase of the m 2 tidal constituent of surface elevations at these stations are respectively within 4 cm and 5 of the observed data the errors in the model predicted m 2 harmonic principal current speeds are less than 6 cm s and the principal current directions and phases are within 14 of the observations the observed surface elevations and currents given in butman b signell r shoukimas p beardsley r c 1988 current observations in buzzards bay 1982 1986 open file report 88 5 united states geological survey are used to validate the model predicted low frequency surface elevations and currents the model predictions in low frequency surface elevations at woods hole closely follow the trends seen in the observations with a correlation coefficient of 0 735 but fail to capture some of the peak surges seen in the observations the model predicted low frequency currents in the east west direction at stations in buzzards bay compare well with the observations with the correlation coefficient exceeding 0 811 and the model capturing the trends seen in the observations for the most part however the model predicted north south velocities does not compare well with the observations the model predictions agree with the observations that the tidal currents in vineyard sound lagged the currents in buzzards bay by more than 3 h the interaction of wind stress with large bathymetric gradients was shown to cause many vortices in buzzards bay as seen from the model predictions model simulations show that the winds play a more dominant role than the tides in the generation of the barotropic residual currents in buzzards bay while the model predicted tide induced residual current was seen to be small

transport capacity of pyroclastic density currents experiments and models of substrate flow interaction 1 one of the more distinctive features of many ignimbrites is the presence of large lithics some greater than meter scale and pumices that have been transported great distances 10 km from the eruptive vent sometimes over steep terrain and expanses of water in many cases these particles have been transported much further than can be explained by aerodynamic forces and ballistic trajectories we examine the forces responsible for transport of large clasts and examine in detail the momentum transfer occurring when particles interact with their boundaries we performed a suite of experiments and numerical simulations to quantify the mass and momentum transfer that occurs when particles interact with a pumice bed substrate and with water substrate two geologically motivated flow end members we find that clasts transported in dilute currents are particularly sensitive to the nature of the boundary and while large particles can skip several times on a water substrate they travel less far than particles that impact pumice bed substrates all else being equal large particles in dense pyroclastic density currents are themselves relatively insensitive to the details of their boundaries however one of the most important ways boundary conditions influence large particles is not through direct interaction but by changing the local concentration of fine particles momentum transfer from fine particles to large particles appears to be required to transport large clasts great distances if initially dense flows become dilute during transport then the transport capacity of large particles in the flow is substantially decreased

vesiculation in magmas from stromboli and implications for normal strombolian activity and paroxysmal explosions in basaltic systems 1 we performed a series of x ray tomographic experiments and lattice boltzmann permeability simulations on pyroclastic products from explosive activity at stromboli between december 2004 and may 2006 we reconstructed the 3 d textures of vesicles to investigate the relationship between the nature of vesiculation in the erupted products and the dynamics of gas transport in the shallow conduit in order to derive implications for the eruptive behavior of basaltic volcanoes scoriae from normal strombolian explosions display remarkably consistent vesicle volume distributions fit by power laws with an exponent of 1 0 2 we ascribe the origin of such distributions to the combined effect of coalescence and continuous nucleation events in the steady state shallow magma system that supplies normal strombolian activity volume distributions and textures of vesicles in pumice clasts from the 5 april 2003 and 15 march 2007 paroxysmal activity are markedly different from those in the scoriae besides a power law function with a higher exponent portions of these distributions can be also fit by an exponential function suggesting the attempt of the system to reach near equilibrium conditions the investigated pumice clasts also lack the large connecting vesicles responsible for the development of degassing pathways in the stromboli magma that erupts the scoriae this testifies to a decreased degassing efficiency of the magma associated with paroxysmal explosions and potential overpressure buildup at depth by comparison with degassing experiments on basaltic melts we derive a time constraint on the order of minutes to hours for the incubation of paroxysms at stromboli

externally triggered renewed bubble nucleation in basaltic magma the 12 october 2008 eruption at halema uma u overlook vent kīlauea hawai i usa 1 from october 2008 until present dozens of small impulsive explosive eruptions occurred from the overlook vent on the southeast side of halema uma u crater at k lauea volcano usa these eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column here we use microtextural observations and data from clasts erupted during the well characterized 12 october 2008 explosive eruption at halema uma u to extend existing models of eruption triggering we present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets we measure the size and number density of bubbles preserved in juvenile ejecta using 2d images and x ray microtomography our data suggest that accumulations of large bubbles with diameters of 50 m to at least millimeters existed at shallow levels within the conduit prior to the 12 october 2008 explosion furthermore a high number density of small bubbles 50 m is measured in the clasts implying very rapid nucleation of bubbles visual observations combined with preexisting geophysical data suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that 1 nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and 2 expanded the preexisting bubbles driving upward acceleration the trigger of eruption and bubble nucleation is thus external to the degassing system

dynamics of coupled diffusive and decompressive bubble growth in magmatic systems bubble growth in an ascending parcel of magma is controlled both by diffusion of oversaturated volatiles and decompression as the magma rises we have developed a numerical model which explores the processes involved in water exsolution from basaltic and rhyolitic melts rising at a constant rate from magma chamber depths of 4 and 1 km while the model does not attempt to simulate natural eruptions it sheds light on the processes which control eruptive behavior under various conditions ascent rates are defined such that a constant rate of decompression dp dt is maintained a variety of initial ascent rates are considered in the model from 1 m s to 100 m s for basalts and from a few centimeters per second to 10 m s for rhyolite at the base of the conduit the model results indicate that for any reasonable ascent rate basaltic melt degasses at a rate sufficient to keep the dissolved volatile concentration at equilibrium with the decreasing ambient pressure rhyolitic melt reaches the surface at equilibrium if its ascent rate is less than 1 m s but it can erupt with high oversaturation at greater ascent rates the latter may lead to explosive eruptions if the ascent rate of rhyolite is 10 m s or more then melt barely degasses at all in the conduit and erupts with the highest oversaturation possible for the case of slow magma rise bubble growth is limited by decompression for the case of rapid magma rise bubble growth is limited by diffusion the results of our simple model do not accurately simulate natural volcanic eruptions but suggest that subsequent more complex models may be able to simulate eruptions using the insights regarding diffusive and decompressive bubble growth processes explored in this study numerical modeling of volcanic degassing may eventually lead to better prediction of eruption timing energetics and hazards of active volcanoes

tidal computations for morecambe bay summary a numerical finite difference technique for computing tides in areas of shoaling water has been established and used to calculate the m2 tide in morecambe bay an essential feature of the method is the representation of drying banks which appear as the water level falls on the ebb and which are re submerged on the flood as the level rises again since the non linear advective acceleration should be important in shallow water schemes for incorporating these terms are examined and the influence of advection on the tide assessed in a series of tests the open sea tides applied to the outer boundary of the numerical bay model were taken first from existing tidal charts then from calculations with a larger scale encompassing model and finally from observations of current in each case the computed tidal elevations within the bay turned out to be in satisfactory agreement with those obtained from observations however it is significant that in the first case inaccuracy of the input tidal distribution deduced from charts gives rise to erroneous currents near the boundary

dynamics and energetics of bubble growth in magmas analytical formulation and numerical modeling we have developed a model of diffusive and decompressive growth of a bubble in a finite region of melt which accounts for the energetics of volatile degassing and melt deformation as well as the interactions between magmatic system parameters such as viscosity volatile concentration and diffusivity on the basis of our formulation we constructed a numerical model of bubble growth in volcanic systems we conducted a parametric study in which a saturated magma is instantaneously decompressed to one bar and the sensitivity of the system to variations in various parameters is examined variations of each of seven parameters over practical ranges of magmatic conditions can change bubble growth rates by 2 4 orders of magnitude our numerical formulation allows determination of the relative importance of each parameter controlling bubble growth for a given or evolving set of magmatic conditions an analysis of the modeling results reveals that the commonly invoked parabolic law for bubble growth dynamics r t1 2 is not applicable to magma degassing at low pressures or high water oversaturation but that a logarithmic relationship r log t is more appropriate during active bubble growth under certain conditions a second aspect of our study involved a constant decompression bubble growth model in which an initially saturated magma was subjected to a constant rate of decompression model results for degassing of initially water saturated rhyolitic magma with a constant decompression rate show that oversaturation at the vent depends on the initial depth of magma ascent on the basis of decompression history explosive eruptions of silicic magmas are expected for magmas rising from chambers deeper than 2 km for ascent rates 1 5 m s 1

delayed disequilibrium degassing in rhyolite magma decompression experiments and implications for explosive volcanism recent numerical models and analog shock tube experiments show that disequilibrium degassing during magma ascent may lead to violent vesiculation very near the surface in this study a series of decompression experiments using crystal free rhyolite melt were conducted to examine the development of large supersaturations due to delayed homogenous spontaneous bubble nucleation melts were saturated at 900 c and 200 mpa with either 5 2 wt dissolved h2o or with 4 2 wt h2o and 640 ppm co2 and isothermally decompressed at linear rates of either 0 003 0 025 or 8 5 mpa s to final pressures between 25 and 175 mpa additional isobaric saturation experiments 900 c 200 25 mpa using pure h2o or mixed h2o co2 fluids establish reference equilibrium solubility curves values homogenous nucleation is triggered in both h2o only and h2o co2 experiments once the supersaturation pressure pss reaches 120 150 mpa and the melt contains two times its equilibrium water contents bubble number density and nucleation rate depend on the supersaturation pressure with values on the order of 102 cm3 and 1 cm3 s for pss 120 mpa 106 cm3 and 103 105 cm3 s for pss 130 150 mpa and 107 cm3 and 106 cm3 s for pss 160 175 mpa nucleation rates are consistent with classical nucleation theory and infer an activation energy for nucleation of 1 5 10 18 j nucleus a critical bubble radius of 2 10 9 m and an effective surface tension for rhyolite at 5 2 wt h2o and 900 c of 0 10 0 11 n m the long nucleation delay limits the time available for subsequent diffusion such that disequilibrium dissolved h2o and co2 contents persist to the end of our runs the disequilibrium degassing paths inferred from our experiments contrast markedly with the equilibrium or quasi equilibrium paths found in other studies where bubble nucleation occurs heterogenously on crystals or other discontinuities in the melt at low pss homogenous and heterogenous nucleation rates are comparable however as are bubble number densities so that at a given decompression rate it appears that nucleation mechanism rather than nucleation rate determines degassing efficiency by fixing the pressure depth at which vesiculation commences and hence the time available for equilibration prior to eruption although real systems are probably never truly crystal free our results show that rhyolitic magmas containing up to 104 crystals cm3 and perhaps as high as 106 crystals cm3 are controlled by homogenous rather than heterogenous nucleation during ascent

a 3d imaging textural characterization of pyroclastic products from the 1538 ad monte nuovo eruption campi flegrei italy abstract the explosive volcanic event of the 1538 ad monte nuovo eruption volcanic explosivity index vei 2 in the campi flegrei high risk caldera italy has a strategic significance in the framework of volcanology and volcanic hazard of caldera forming magmatic systems in fact it represents the last and unique historical eruption of the highly populated phlegraean restless caldera and its precursory and eruptive phenomena are well known because they were described in detail by contemporaneous eyewitnesses in this study a set of samples representative of the complete stratigraphic sequence of the monte nuovo eruption was characterized using phase contrast synchrotron radiation computed microtomography and quantitatively investigated through the development of a new protocol for 3d textural analysis of highly vesiculated volcanic rocks previous studies of products from this eruption are available in the literature mostly based on 2d imaging techniques and thus provide a useful data set for comparison the 3d textural measurements allow us to investigate the subvolcanic processes mechanisms and timing of magma degassing that occurred during magma ascent in the conduit for each stage of the eruption and their relationship with the variations in the eruptive style described in the contemporaneous accounts of the eruption this information is fundamental for the definition of a volcanic eruption scenario for such low vei events often recurrent in the history of the caldera and is useful both for hazard assessment and emergency planning

new geochemical insights into volcanic degassing magma degassing plays a fundamental role in controlling the style of volcanic eruptions whether a volcanic eruption is explosive or effusive is of crucial importance to approximately 500 million people living in the shadow of hazardous volcanoes worldwide studies of how gases exsolve and separate from magma prior to and during eruptions have been given new impetus by the emergence of more accurate and automated methods to measure volatile species both as volcanic gases and dissolved in the glasses of erupted products the composition of volcanic gases is dependent on a number of factors the most important being magma composition and the depth of gas melt segregation prior to eruption this latter parameter has proved difficult to constrain in the past yet is arguably the most critical for controlling eruptive style spectroscopic techniques operating in the infrared have proved to be of great value in measuring the composition of gases at high temporal resolution such methods when used in tandem with microanalytical geochemical investigations of erupted products are leading to better constraints on the depth at which gases are generated and separated from magma a number of recent studies have focused on transitions between explosive and effusive activity and have led to a better understanding of gas melt segregation at basaltic volcanoes other studies have focused on degassing during intermediate and silicic eruptions important new results include the recognition of fluxing by deep derived gases which buffer the amount of dissolved volatiles in the melt at shallow depths and the observation of gas flow up permeable conduit wall shear zones which may be the primary mechanism for gas loss at the cusp of the most explosive and unpredictable volcanic eruptions in this paper i review current and future directions in the field of geochemical studies of volcanic degassing processes and illustrate how the new insights are beginning to change the way in which we understand and classify volcanic eruptions

time and temperature in petroleum formation application of lopatin s method to petroleum exploration n v lopatin in the soviet union has developed a method for taking both time and temperature into account as factors in thermal maturation of kerogen lopatin s time temperature index of maturity tti values correlate with the thermal regimes corresponding to generation and preservation of hydrocarbons because such information is potentially of great interest for oil exploration a calibration and evaluation have been made of lopatin s method within the limitations of the data presently available the following statements can be made 1 the rate of the chemical reactions involved in thermal maturation of organic material appears to double with every 10 c 18 f rise in temperature 2 threshold values of lopatin s time temperature index of maturity tti are 15 onset of oil generation 75 peak oil generation 160 end oil generation 500 40 oil preservation deadline 1 000 50 oil preservation deadline 1 500 wet gas preservation deadline 65 000 dry gas preservation deadline 3 tti values calculated from lopatin reconstructions consistently agree with other maturation parameters commonly used by petroleum geochemists potential applications of lopatin s method for oil exploration include timing of oil generation calculation of volume of hydrocarbons generated within a basin and determination of economic deadlines

residual circulation in the gulf of suez egypt abstract a three dimensional non linear hydrodynamical numerical model has been used to study the water movement in the shallow gulf of suez at the northern end of the red sea the computational grid used to schematize the gulf has a horizontal resolution of 3 3 km and levels at 5 15 30 and 60 m in the vertical the simulations were carried out separately for tidal forcing and wind forcing further in each case stratified conditions were considered to represent the state of the gulf during late summer periods in the case of wind induced circulations a uniform wind speed of 5 m s 1was applied over the gulf blowing from nw n and ne directions the results suggest that the dynamics of the gulf of suez are governed by the interaction of wind and density forcing only occasionally for instance in the vicinity of headlands and islands tidal residuals are of similar importance for the overall residual currents it was found that the typical quasi baroclinic two layered inverse estuarine circulation observed in the gulf of suez can be produced under the combined action of m2 tide and winds blowing steadily from the north eastern direction north western and northern winds produced the reverse pattern

eruption dynamics of hawaiian style fountains the case study of episode 1 of the kīlauea iki 1959 eruption hawaiian eruptions are characterized by fountains of gas and ejecta sustained for hours to days that reach tens to hundreds of meters in height quantitative analysis of the pyroclastic products from the 1959 eruption of k lauea iki k lauea volcano hawai i provides insights into the processes occurring during typical hawaiian fountaining activity this short lived but powerful eruption contained 17 fountaining episodes and produced a cone and tephra blanket as well as a lava lake that interacted with the vent and fountain during all but the first episode of the eruption the focus of this paper microtextural analysis of hawaiian fountaining products from this opening episode is used to infer vesiculation processes within the fountain and shallow conduit vesicle number densities for all clasts are high 106 107 cm 3 post fragmentation expansion of bubbles within the thermally insulated fountain overprints the pre fragmentation bubble populations leading to a reduction in vesicle number density and increase in mean vesicle size however early quenched rims of some clasts with vesicle number densities approaching 107 cm 3 are probably a valid approximation to magma conditions near fragmentation the extent of clast evolution from low vesicle to melt ratio and corresponding high vesicle number density to higher vesicle to melt ratio and lower vesicle number density corresponds to the length of residence time within the fountain

application of three dimensional boundary fitted circulation model to providence river a three dimensional boundary fitted circulation model by muin and spaulding in 1996b was applied to simulate tidal and wind driven flow in upper narragansett bay and the providence and seekonk rivers the model employs a split mode solution technique the exterior mode vertically averaged was solved using a semiimplicit method the interior mode vertical structure was solved explicitly except for the vertical diffusion term that is implicit the vertical eddy viscosity was calculated from a one equation turbulence model model predicted tidal currents for the principal constituents were in good agreement with field observations the im d2 dominated the tidal currents accounting for 80 of the total tidal energy a combination of bottom friction coefficient of 0 0025 and constant 0 3 in the mixing length formulation of the turbulence model provided the best fit to the available observations the model predicted a cooscillating wave pattern with the im d2 water elevation about three hours out of phase with the current an approximate one hour phase difference between the currents in the deep dredged channel and the shallow water on either side was predicted wind driven flow was simulated under constant wind forcing of 0 01 n m ua 0 1 dyne cm along and across the river axis the steady state flow patterns are consistent with observations showing flows against the wind direction in the dredged channel the model predicted steady state currents in the dredged channel have maximum values of 3 4 cm s for winds along the channel axis and 1 cm s for winds across the channel axis

estimating kinetic parameters for organic reactions from geological data an example from the gippsland basin australia abstract an examination has been made of strategies to estimate kinetic parameters activation energy and pre exponential factor for organic reactions from geological data least squares regression and a direct search algorithm were applied to a set of sterane isomerization data from volador 1 an oil exploration well in the gippsland basin australia the best fit kinetic parameters depend on the number and quality of data points used as well as the thermal history because both the data and thermal history have intrinsic uncertainties a range of kinetic parameters can explain the data to a specified level of misfit for reasonable variations in the thermal history the effective kinetic parameters for the data from volador 1 are e 75 100kj mol and a 10 1 10 5 s it is observed that the range of acceptable parameter pairs has a linear relation of the form ln a a be and the slope b is dependent on the thermal history small departures of either parameter from this relation lead to significant variations in the predicted values for the reaction progress

the 15 march 2007 explosive crisis at stromboli volcano italy assessing physical parameters through a multidisciplinary approach basaltic volcanoes are dominated by lava emission and mild explosive activity nevertheless many basaltic systems exhibit from time to time poorly documented and little understood violent explosions a short lived multiblast explosive crisis paroxysmal explosion occurred on 15 march 2007 during an effusive eruptive crisis at stromboli italy the explosive crisis which started at 20 38 14 ut had a total duration of 5 min the combined use of multiparametric data collected by the permanent instrumental networks seismic acoustic and thermal records and a field survey carried out immediately after the event enabled us to constrain the eruptive dynamics and quantify physical parameters the eruption consisted of three major pulses in the first lithic blocks and ash were ejected at speeds of 100 155 m s and 130 210 m s respectively the high solid load of the eruptive jet resulted in the partial collapse of the column with the formation of a small volume pyroclastic density current the second 12 s long pulse emitted 2 2 2 7 107 kg of tephra mass discharge rate 1 9 2 3 106 kg s forming a 3 km high convective plume dispersing tephra up to the west coast and a dilute density current with limited dispersal downslope of the craters a final 30 s long phase formed a scoria flow with a volume of 1 5 1 7 104 m3 mass discharge rate 5 9 6 7 105 kg s a total runout of 200 m and a velocity of 45 m s the total gas volume involved in the explosion was 1 3 1 9 104 m3 with an initial overpressure of 7 9 0 4 mpa we compared the 15 march 2007 event with historical paroxysms in particular with that of 5 april 2003 which was remarkably similar

vesiculation of basaltic magma during eruption vesicle size distributions in vent lavas from the pu9u9o9o kupaianaha eruption of kilauea volcano are used to estimate nucleation and growth rates of h 2 o rich gas bubbles in basaltic magma nearing the earth9s surface 120 m depth by using well constrained estimates for the depth of volatile exsolution and magma ascent rate nucleation rates of 35 9 events cm 3 s 1 and growth rates of 3 2 x 10 4 cm s are determined directly from size distribution data the results are consistent with diffusion controlled growth as predicted by a parabolic growth law this empirical approach is not subject to the limitations inherent in classical nucleation and growth theory and provides the first direct measurement of vesiculation kinetics in natural settings in addition perturbations in the measured size distributions are used to examine bubble escape accumulation and coalescence prior to the eruption of magma

application of biological markers in the recognition of palaeohypersaline environments the hydrocarbon fractions of a marl sample from a messinian basin in the northern appenines and four oils rozel point oil utah and three seep oils from sicily were studied

using so 2 camera imagery and seismicity to examine degassing and gas accumulation at kīlauea volcano may 2010 abstract so 2 camera measurements at k lauea volcano hawaii in may of 2010 captured two occurrences of lava lake rise and fall within the halema um au crater summit vent during high lava stands we observed diminished so 2 emission rates and decreased seismic tremor similar events at k lauea have been described as the result of sporadic degassing following gas accumulation beneath a mostly impermeable lava lake surface incorporation of so 2 camera data into a multi parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of these high stand events with accumulated gas release upon lake level drop compensating for the gas deficit reached during accumulation

shallowly driven fluctuations in lava lake outgassing gas pistoning kīlauea volcano abstract lava lakes provide ideal venues for directly observing and understanding the nature of outgassing in basaltic magmatic systems k lauea volcano s summit lava lake has persisted for several years during which seismic and infrasonic tremor amplitudes have exhibited episodic behavior associated with a rise and fall of the lava surface gas pistoning since 2010 the outgassing regime of the lake has been tied to the presence or absence of gas pistoning during normal behavior no gas pistoning the lake is in a spattering regime consisting of higher tremor amplitudes and gas emissions in comparison gas piston events are associated with an abrupt rise in lava level up to 20 m during which the lake enters a non spattering regime with greatly decreased tremor and gas emissions we study this episodic behavior using long term multidisciplinary monitoring of the lake including seismicity infrasound gas emission and geochemistry and time lapse camera observations the non spattering regime i e rise phase of a gas piston cycle reflects gas bubbles accumulating near the top of the lake perhaps as a shallow foam while spattering regimes represent more efficient decoupling of gas from the lake we speculate that the gas pistoning might be controlled by time varying porosity and or permeability in the upper portions of the lava lake which may modulate foam formation and collapse competing models for gas pistoning such as deeply sourced gas slugs or dynamic pressure balances are not consistent with our observations unlike other lava lakes which have cyclic behavior that is thought to be controlled by deeply sourced processes external to the lake itself we show an example of lava lake fluctuations driven by cycles of activity at shallow depth and close to the lake s surface these observations highlight the complex and unsteady nature of outgassing from basaltic magmatic systems

dynamics of degassing at kilauea volcano hawaii at kilauea volcano hawaii the recent long lived eruptions of mauna ulu and pu u o o have occurred in two major stages defining a characteristic eruptive pattern the first stage consists of cyclic changes of activity between episodes of fire fountaining and periods of quiescence or effusion of vesicular lava the second stage consists only of continuous effusion of lava we suggest that these features reflect the dynamics of magma degassing in a chamber which empties into a narrow conduit in the volcano chamber gas bubbles rise through magma and accumulate at the roof in a foam layer the foam flows toward the conduit and its shape is determined by a dynamic balance between the input of bubbles from below and the output into the conduit the foam thickness is proportional to lq e2 l g 1 4 where l and l are the viscosity and density of magma e is the gas volume fraction in the foam g is the acceleration of gravity and q is the gas flux the bubbles in the foam deform under the action of buoyancy and the maximum permissible foam thickness is hc 2 e lgr where is the coefficient of surface tension and r is the original bubble radius if this critical thickness is reached the foam collapses into a large gas pocket which erupts into the conduit foam accumulation then resumes and a new cycle begins the attainment of the foam collapse threshold requires a gas flux in excess of a critical value which depends on viscosity surface tension and bubble size hence two different eruption regimes are predicted 1 alternating regimes of foam buildup and collapse leading to the periodic eruption of large gas volumes and 2 steady foam flow at the roof leading to continuous bubbly flow in the conduit the essential result is that the continuous process of degassing can lead to discontinuous eruptive behavior data on eruption rates and repose times between fountaining phases from the 1969 mauna uiu and the 1983 1986 pu u o o eruptions yield constraints on three key variables the area of the chamber roof must be a few tens of square kilometers with a minimum value of about 8 km2 magma reservoirs of similar dimensions are imaged by seismic attenuation tomography below the east rift zone close to the roof the gas volume fraction is a few percent and the gas bubbles have diameters lying between 0 1 and 0 6 mm these estimates are close to the predictions of models for bubble nucleation and growth in basaltic melts as well as to the observations on deep submarine basalts the transition between cyclic and continuous activity occurs when the mass flux of gas becomes lower than a critical value of the order of 103 kg s in this model changes of eruptive regime reflect changes in the amount and size of bubbles which reach the chamber roof

organic matter maturation under the influence of a deep intrusive heat source a natural experiment for quantitation of hydrocarbon generation and expulsion from a petroleum source rock toarcian shale northern germany four shallow boreholes were drilled in the hils syncline northern germany in order to determine quantitatively the amount of hydrocarbons generated and expelled during maturation of a typical kerogen type ii bearing source rock the holes penetrated the carbonceous lias shales posidonia shale lower toarcian and part of the adjacent dogger and lias mudstones the maturity of the organic matter in the cores recovered ranges from immature 0 48 r0 to overmature 1 45 r0 due to location of the hils syncline in the vicinity of the vlotho massif which is deep seated intrusive body facies variations of the lias within the short geographical distances in the study area are negligible organic matter mass balance calculations were based on detailed organic geochemical analyses of residual material in the lias shales kerogen bitumen etc and on the evidence of a uniform initial composition of these sediments in the study area dead carbon determinations supported this latter criterion but were not used as a parameter in the calculations about 50 of the initial kerogen was transformed into oil gas and inorganic compounds during the vitrinite reflectance increase from 0 48 to 0 88 ro and only marginally more during the maturity increase from 0 88 to 1 45 ro only a small portion of the generated material remained in the source rock even at a relatively early stage of generation 0 68 ro expulsion efficiency of oil plus gas reached a value of 86 at the end of the main generation stage 0 88 ro

tectonic history from vitrinite reflectance summary the degree of organic metamorphism of carbonaceous matter is expressed in terms of the reflectance of vitrinite vitrinite reflectance is a function of temperature and time a model relating vitrinite reflectance to temperature and time is proposed from published vitrinite reflectance temperature time nomograms a plot of the logarithm of vitrinite reflectance versus depth gives a linear curve for vitrinite reflectance greater than 0 3 per cent in a sedimentary basin which has had a constant geothermal gradient throughout its tectonic history the geothermal gradient c km 1 is equal to the slope of the curve multiplied by 194 8 the model is also applied to sedimentary basins with decreasing geothermal gradient with time i e cooling continental margin basins the vitrinite reflectance model is demonstrated to be consistent with observed reflectance data and tectonic information from the sydney bowen surat and gippsland basins in eastern australia

very long period conduit oscillations induced by rockfalls at kilauea volcano hawaii 1 eruptive activity at the summit of kilauea volcano hawaii beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period vlp seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the halemaumau pit crater we use raw data recorded with an 11 station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 august 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal to noise ratio to determine the source centroid location and source mechanism we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account we apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of green s functions representing the seismometer response to translation and tilt ground motions this method enables a robust description of the source mechanism over the period range 1 1000 s the vlp signals associated with the rockfalls originate in a source region 1 km below the eastern perimeter of the halemaumau pit crater the observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack dike dipping 80 to the north intersecting a north striking crack another dike dipping 65 to the east each rockfall is marked by a similar step like inflation trailed by decaying oscillations of the volumetric source attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column assuming a simple lumped parameter representation of the shallow magmatic system the observed pressure and volume variations can be modeled with the following attributes rockfall volume 200 4500 m3 length of magma column 120 210 m diameter of pipe connecting the halemaumau pit crater to the subjacent dike system 6 m average thickness of the two underlying dikes 3 6 m and effective magma viscosity 30 210 pa s most rockfalls occur during episodes of sustained deflation of the kilauea summit the mass loss rate in the shallow magmatic system is estimated to be 1400 15 000 kg s 1 based on measurements of the temporal variation of vlp period in the two large rockfalls that occurred on 29 august 2012

relationships between thermal maturity indices calculated using arrhenius equation and lopatin method implications for petroleum exploration thermal maturity can be calculated with time temperature indices tti based on the arrhenius equation using kinetics applicable to a range of types ii and iii kerogens these ttis are compared with tti calculations based on the lopatin method and are related theoretically and empirically via vitrinite reflectance to the petroleum generation window the ttis for both methods are expressed mathematically as integrals of temperature combined with variable linear heating rates for selected temperature intervals heating rates control the thermal maturation trends of buried sediments relative to arrhenius ttis lopatin ttis tend to underestimate thermal maturity at high heating rates and overestimate it as low heating rates complex burial histories applicable to a range of tectonic environments illustrate the different exploration decisions that might be made on the basis of independent results of these two thermal maturation models

the structure of basaltic scoria and reticulite and inferences for vesiculation foam formation and fragmentation in lava fountains in this investigation pyroclast structures are used to constrain degassing in basaltic lava fountains vesicle size shape number density interconnectedness and packing character are quantified and related to 1 the kinetics of bubble nucleation and growth 2 the structural evolution of magmatic foams and 3 the influence of vesiculation rate on magma fragmentation measurements made on a diverse suite of pyroclasts from kilauea volcano indicate that basaltic foams evolve through an initially disordered closed celled spherical state to a well ordered open celled polyhedral state as the vesicularity rises from 75 to 98 the structural changes occur rapidly 10 s in the conduit and fountain in response to an intense vesiculation burst vesicle size distribution systematics indicate bubble nucleation rates 2 104 events cm 3s 1 that are approximately three orders of magnitude greater than those found for effusive eruptive activity bubble growth rates 9 10 4 cm s exceed effusive estimates by a factor of 3 the observed runaway rate of bubble production indicates strong supersaturations at the onset of nucleation we speculate that the rise speed of the magma as it reaches the level where significant volatile exsolution begins determines the intensity of the vesiculation burst and hence the vigor of the eruption rapid expansion and acceleration of the magma under these conditions may provide the impetus for fragmentation

biological marker characteristics of oils and asphalts from carbonate source rocks in a rapidly subsiding graben dead sea israel abstract a detailed gc ms study of biological marker compounds in the saturated and aromatic hydrocarbon fractions of oils and asphalts from the dead sea area israel provided decisive information to the solution of a long lasting controversy by showing that the asphalts are products of early generation in an immature stage from the same type of carbonate source rock which generated more mature oils the asphalts are not biodegraded residues of the oils oils from six different wells and asphalts from wells outcrops and a floating block from the dead sea all have very similar sterane and triterpane patterns they all lack rearranged steranes diasteranes indicating a carbonate source matrix and compare reasonably well with a sample of upper cretaceous bituminous chalk from nebi musa the main difference between the oils and the asphalts is a significantly higher triaromatic to mono plus triaromatic steroid hydrocarbon ratio in the former this is explained as a result of rapid subsidence and heating of their source rock close to the deep parts of the dead sea graben the oils thus were generated in the more deeply buried source rock blocks under the graben fill whereas the asphalts either originate from an immature source rock section closer to the graben rims or represent an earlier phase of generation and expulsion this study also provides general information on the evolution of biological markers in carbonate source rocks low activation energy processes like isomerisation of steranes appear to occur much faster at low temperatures than in shales the high sulfur content and less cross linking of the biogenic organic matter into a complex kerogen structure are suggested to be responsible for this care should be taken when using only sterane isomerisation to assess the maturity of hydrocarbons from carbonate rocks and of carbonate derived oils

abrupt transitions during sustained explosive eruptions examples from the 1912 eruption of novarupta alaska plinian ignimbrite activity stopped briefly and abruptly 16 and 45 h after commencement of the 1912 novarupta eruption defining three episodes of explosive volcanism before finally giving way after 60 h to effusion of lava domes we focus here on the processes leading to the termination of the second and third of these three episodes early erupted pumice from both episodes show a very similar range in bulk vesicularity but the modal values markedly decrease and the vesicularity range widens toward the end of episode iii clasts erupted at the end of each episode represent textural extremes at the end of episode ii clasts have very thin glass walls and a predominance of large bubbles whereas at the end of episode iii clasts have thick interstices and more small bubbles quantitatively all clasts have very similar vesicle size distributions which show a division in the bubble population at 30 m vesicle diameter and cumulative number densities ranging from 107 109 cm 3 patterns seen in histograms of volume fraction and the trends in the vesicle size data can be explained by coalescence signatures superimposed on an interval of prolonged nucleation and free growth of bubbles compared to experimental data for bubble growth in silicic melts the high 1912 number densities suggest homogeneous nucleation was a significant if not dominant mechanism of bubble nucleation in the dacitic magma the most distinct clast populations occurred toward the end of plinian activity preceding effusive dome growth distributions skewed toward small sizes thick walls and teardrop vesicle shapes are indicative of bubble wall collapse marking maturation of the melt and onset of processes of outgassing the data suggest that the superficially similar pauses in the 1912 eruption which marked the ends of episodes ii and iii had very different causes through episode iii the trend in vesicle size data reflects a progressive shift in the degassing process from rapid magma ascent and coupled gas exsolution to slower ascent with partial open system outgassing as a precursor to effusive dome growth no such trend is visible in the episode ii clast assemblages we suggest that external changes involving failure of the conduit vent walls are more likely to have effected the break in explosive activity at 45 h

time scales of foam stability in shallow conduits insights from analogue experiments volcanic systems can exhibit periodical trends in degassing activity characterized by a wide range of time scales understanding the dynamics that control such periodic behavior can provide a picture of the processes occurring in the feeding system towards this end we analyzed the periodicity of outgassing in a series of decompression experiments performed on analogue material argon saturated silicone oil plus glass beads fibers scaled to serve as models of basaltic magma to define the effects of liquid viscosity and crystal content on the time scale of outgassing we investigated both 1 pure liquid systems at differing viscosities 100 and 1000 pa s 2 particle bearing suspensions diluted and semi diluted the results indicate that under dynamic conditions e g decompressive bubble growth fluid ascent within the conduit the periodicity of foam disruption may be up to several orders of magnitude less than estimates based on the analysis of static conditions this difference in foam disruption time scale is inferred to result from the contribution of bubble shear and bubble growth to inter bubble film thinning the presence of particles in the semi diluted regime is further linked to shorter bubble bursting times likely resulting from contributions of the presence of a solid network and coalescence processes to the relative increase in bubble break up rates finally it is argued that these experiments represent a good analogue of gas piston activity i e the periodical rise and fall of a basaltic lava lake surface implying a dominant role for shallow foam accumulation as a source process for these phenomena this article is protected by copyright all rights reserved

a global synthesis of lava lake dynamics abstract active lava lakes represent a variety of open vent volcanism in which a sizable body of lava accumulates at the top of the magma column constrained by the vent and or crater geometry the longevity of lava lakes reflects a balancing of cooling and outgassing occurring at the surface and input of hot and gas rich magma from below due to their longevity and relative accessibility lava lakes provide a natural laboratory for studying fundamental volcanic processes such as degassing convection and cooling this article examines all seven lakes that existed at the time of writing in 2018 located in the pacific antarctica africa and south and central america these lakes span all tectonic environments and a range of magma compositions we focus on analysis of the lake surface motion using image velocimetry which reveals both similarities and contrasts in outgassing and lake dynamics when comparing the different lakes we identify two categories of lake behavior organized erta ale nyiragongo k lauea after 2011 and erebus and chaotic villarrica masaya marum this division does not map directly to lake size viscosity gas emission rate or temperature instead when examined together we find that the lakes follow a linear relationship between average surface speed and the ratio of total gas flux to lake surface area this relationship points to the combined importance of both flux and lake size in addition to the total volume of gas emission and suggests that a shared deep mechanism controls the supply of heat and gas to all lakes on the other hand the differences between chaotic and organized lakes highlight the important role of the geometry of the conduit lake transition which superimposes a shallow signal on that of the deep circulation the spatial patterns of surface motion we document suggest that the release of gas bubbles at chaotic lakes is more efficient i e bubbles are less likely to be retained and recycled compared with organized lakes in addition the data presented here indicate that the solidified crust of organized lakes plays a role in regulating convection and outgassing in lava lakes

thermal effects of magma degassing abstract the thermodynamics of diffusive bubble growth is dominated by latent heat of water exsolution vaporization and the work of gas expansion p dv work a numerical assessment of the effects of these on cooling of bubbly magma water in albite indicates that a magma can exsolve volatiles at equilibrium or with varying degrees of oversaturation depending on decompression history and degassing kinetics heat of water exsolution vaporization from magmatic melts has not previously been experimentally or analytically determined and we have evaluated it from others thermodynamic functions known for the albite water system heat of exsolution is small for pressures over 100 mpa but can reach 20 kj mole at 1 2 mpa 10 20 bar oversaturation degassing at 10 20 bar can cause cooling of albite by 8 k wt of exsolved water the thermal effect of equilibrium degassing depends on the starting pressure of decompression because it follows the solubility curve for a saturated melt at 100 mpa 4 wt water it can cool a magma at least by 35 k before it disrupts into spray as gas volume fraction exceeds 0 8 although we do not specifically consider the dynamics of bubble interactions in this study our results suggest that the dynamic cooling rate could be significant at the vent of an erupting magma column and under extreme eruption conditions can lead to glass formation around fast growing bubbles of magmatic foam this could cause fragmentation into fine ash by brittle failure conversely bubble wall cooling even before glass formation can serve to reduce diffusive volatile flux into bubbles decreasing the overall cooling of the system and lead to solidification of oversaturated magmas as occasionally observed in eruption products

vesiculation and fragmentation history in a submarine scoria cone forming eruption an example from nishiizu izu peninsula japan an uplifted 50 m thick half dissected submarine emplaced below wave base scoria cone occurs as dipping beds in coastal outcrops at nishiizu on the izu peninsula in japan concentrically outward dipping weakly stratified ungraded framework supported thin to very thick beds consist of brown coarse tuff to scoria lapilli tuff with outsized fluidal bombs throughout accessory lithic clasts chiefly occur in the lowermost visible beds scoria bombs have quenched margins weak bread crust textures and their vesicle number densities decrease inward which is indicative of fast surface cooling composite textures in the scoria bombs indicate recycling and agglutination of quenched and semi molten pyroclasts at the submarine vent in contrast to weak concentric gradations in vesicle size distribution in the bombs lapilli have asymmetrical gradients in vesicle size distribution indicating that they are fragments of coarser quenched lumps three grain size modes characterise the nishiizu brown scoria with coarse magma lumps ejected during magmatic fragmentation and quench jointed upon contact with seawater to be subsequently fragmented into lapilli and coarse ash by various styles of fragmentation where seawater plays a critical role the cone was constructed by slow moving fallout fed granular flow creep fed directly by suspension settling focused at the crater rim but extending onto the cone flanks with only minor resedimentation by granular flows nishiizu deposits yield an exceptional record of eruption and sedimentation dynamics during submarine cone building activity and in this study we compare their vesiculation and fragmentation mechanisms with those of potential subaerial analogues

geochemical effects of petroleum migration and expulsion from toarcian source rocks in the hils syncline area nw germany for mature oil prone source rocks of toarcian age in the hils syncline area nw germany the geochemical effects of primary migration of petroleum are shown expulsion has occurred at high efficiencies primarily via a network of macro and micro fractures the effects of this fracture controlled expulsion with respect to gross and molecular composition of c15 soluble organic matter are documented samples from fractured intervals have expelled total c15 soluble organic matter more efficiently than adjacent unfractured intervals a systematic comparison of samples taken at cm intervals away from macro fractures reveals compositional gradients which are interpreted to reflect preferential depletion of saturated and aromatic hydrocarbons from the source rock shale into the fractures extracts of shale samples directly at fracture walls are enriched in nso compounds and highly polar asphaltene like fractions the composition of residual petroleum fluids extracted from the fracture fill materials is variable in some cases there is a close match with that of the shale at some distance away from the fractures in other cases there is geochemical evidence to suggest that petroleum fluids were chemically altered during lateral flow along the fractures based on available geological and geochemical evidence a conceptual model is proposed to explain fracture development and fluid flow processes rapid heating by a nearby intrusive body at shallow depths when the toarcian sediments still contained about 15 20 pore water has played an important role overpressuring of aqueous pore fluids with abundant dissolved carbonate caused intense fracturing of the rich intervals in the mature toarcian sections finally the closely interdependent relationships between petroleum generation and petroleum migration processes are documented especially the evolution of these relationships with maturity progress from immature 0 48 rm to mature 0 68 and 0 88 rm and overmature 1 45 rm stages is shown for example in the maturity interval from 0 68 to 0 88 rm most of the measured shale compaction and porosity decrease is achieved by expulsion of petroleum for some samples there is an inverse relationship between kerogen quality and extractable hydrocarbon yields suggesting that expulsion efficiencies were higher from the better quality intervals because the petroleum achieved higher saturations and hence higher relative permeabilities

biomarker transformations as constraints for the depositional environment and for maximum temperatures during burial of opalinus clay and posidonia shale in northern switzerland abstract samples from two argillaceous formations opalinus clay and posidonia shale of near identical maturity from northern switzerland were subjected to a geochemical characterisation of organic matter and to confined system pyrolysis experiments throughout the study area the characteristics of organic matter are similar indicating a spatially homogeneous sedimentary facies posidonia shale contains marine organic matter deposited in a reducing environment while a predominantly terrigenous source and a more oxidising environment of deposition was identified for opalinus clay in the western and central parts of the study area organic maturity is close to the onset of oil generation in the easternmost part a higher maturity has been reached due to a deeper burial below thick tertiary molasse deposits isothermal pyrolysis experiments were conducted at temperatures between 250 and 390 c over 24 h bitumen yields increase along similar pathways for both opalinus clay and posidonia shale but the maximum values are displaced by 10 20 c data pertaining to maturity were determined from gc ms analyses of saturated hydrocarbons and specific attention was given to c 29 sterane and c 32 hopane isomerisation ratios the evolution of these parameters with rising temperature is slightly different in the two formations which is attributed to the contrasting organic facies the pyrolysis data together with literature data from natural basins were used to calculate kinetic parameters for c 29 sterane and c 32 hopane assuming a single step isomerisation scheme according to the arrhenius law the resulting values based on pyrolysis data alone are very similar to those based on the combination of pyrolysis and natural data activation energies are similar in both formations while the frequency factors are up to one order of magnitude higher for posidonia shale when compared to opalinus clay for the benken site maximum temperature during cretaceous burial was calculated on the basis of the kinetic data using the tti approach the resulting temperatures of 75 80 c are 5 10 c below those derived in the literature from apatite fission track analysis vitrinite reflectance and basin modelling

biomarker geochemistry of crude oils and lower paleozoic source rocks in the tarim basin western china an oil source rock correlation study abstract in order to select biomarker parameters that have potential to quantitatively distinguish between the cambrian lower ordovician and upper ordovician source rocks in the tarim basin with a total organic carbon content 0 5 42 lower paleozoic shale samples were collected from a number of wells and outcrops in the tarim basin northwestern china the biomarkers in the extracts of the source rocks were analyzed by gas chromatography mass spectrometry the results show that c23 and c30 were dominant in the tricyclic terpanes and pentacyclic triterpanes respectively and c27 and c29 were both dominant in c27 29 steranes of the source rock samples the distribution of n alkanes isoprenoids and terpanes indicate that the source rock samples have been slightly to moderately biodegraded and they originated from a reducing depositional environment with considerable marine algal input in comparatively upper ordovician source rock samples presented lower values of c28 c29 regular steranes and higher values of c29 20s 20s 20r and c29 abnormally lower values of c29 20s 20s 20r and c29 were observed in the cambrian lower ordovician strata the set of data in this study show that homohopane index including c31homohopane c30hopane c31 35homohopanes c30hopane and the ratios related to the triaromatic steroids including c26 20s c28 20s tas c27 20r c28 20r c26 20r c27 20s c28 20s and c28 c26 c27 tas could serve as potential biomarker parameters to distinguish between the two groups of source rocks in the tarim basin the values of those parameters in above order in cambrian lower ordovician source rocks were 0 3 0 5 0 85 and 0 8 1 9 1 2 respectively the correlation between the source rock samples and 19 oils obtained from the tadong tabei and tazhong uplifts in the tarim basin show that most of the oils in tabei uplift are probably from upper ordovician source rocks the oils in tadong uplift are more likely from cambrian lower ordovician source rocks and the tazhong uplift is more like a mixed oil region with oils from both of the two groups of source rocks

linking variable explosion style and magma textures during 2002 at stromboli volcano italy microtextural characteristics of fresh ejecta from stromboli volcano were examined from three periods of differing eruption style and intensity in 2002 activity shifted from relatively weak and infrequent ash charged explosions during january through may into two broad cycles of waxing activity in june through late september and late september through december followed by the onset on 28 december of the 2002 2003 effusive eruption analyzed sets of lapilli from may september october and 28 december show contrasts in the physical properties of magma resident in the shallow conduit during this range of activity three distinct textures are observed among the analyzed pyroclasts low density ld with an abundance of subspherical bubbles the presence of large irregularly shaped bubbles and a light to transparent glass matrix transitional texture tt with an intermediate number of subspherical bubbles a high frequency of large irregularly shaped bubbles and a honey colored glass matrix and high density hd with sparse relatively small bubbles conspicuous large irregular bubbles and a dark glass matrix observational and quantitative data density vesicle size indicate that these textures are linked through variable residence time in stromboli s shallow conduit with an ongoing evolution from ld to hd magma calculations suggest that residual ld magma will evolve to hd texture in a period of hours to days contrasting amounts of the ld tt and hd magmas are present in each sample with the most tt in may the most ld in september october and the most hd in december this implies that the shallow magma had a different rheology at each collection period the viscosity of ld and hd magmas are calculated to be in the range of 2 000 to 2 600 and 3 000 to 5 000 pa s respectively which with their changing proportions must have implications for rates of bubble slug ascent and processes of fragmentation this study suggests that an increasing maturity of magma in stromboli s shallow conduit with resultant increase in viscosity feeds back to reduce the intensity of explosions whereas a steady flux of ld magma favors more powerful explosions

vesiculation of high fountaining hawaiian eruptions episodes 15 and 16 of 1959 kīlauea iki the 1959 summit eruption of k lauea volcano produced the highest recorded hawaiian fountain in hawai i quantitative analysis of closely spaced samples from the final two high fountaining episodes of the eruption result in a fine scale textural study of pyroclasts and provide a record of postfragmentation processes as clast vesicularity increases the vesicle number density decreases and vesicle morphology shifts from small and round to larger and more irregular the shift in microtexture corresponds to greater degrees of postfragmentation expansion of clasts with higher vesicularity we suggest the range of clast morphologies in the deposit is related to thermal zonation within a hawaiian fountain where the highest vesicularity clasts traveled in the center and lowest traveled along the margins vesicle number densities are greatest in the highest fountaining episode and therefore scale with intensity of activity major element chemical analyses and fasciculate crystal textures indicate microlite rich zones within individual clasts are portions of recycled lava lake material that were incorporated into newly vesiculating primary melt

pyroclast characteristics of a subaqueous to emergent surtseyan eruption black point volcano california abstract we present the first ever quantitative analysis of pyroclast textures from subaqueous deposits of a surtseyan emergent volcano the volcano black point formed in mono lake s pleistocene precursor lake russell at 13 300 500 bp from an explosive basaltic eruption involving the interaction of shallow lake water 105 m with vesiculating magma the bulk of the edifice comprises unconsolidated tephra and represents submerged stages of eruption thin summit deposits formed during a final phase of emergent eruption and produced a group of low tuff rings lapillus clasts 16 mm to 32 mm and fine ash grains 4 0 0625 mm from the subaqueous deposits have textures that reveal a conduit history of bubble nucleation growth collapse intense microlite crystallisation and magma water interaction lapilli vesicularity ranges broadly to high values up to 89 while vesicle number densities 5 4 102 4 8 104 vesicles per mm3 are comparable to those published for strombolian and plinian basaltic eruptions 1 104 vesicles per mm3 vesicle size and volume distributions suggest continuous bubble nucleation but the suite of data as a whole indicates that a late stage nucleation event played a significant role in shaping the eruption ash grain morphology and surface features indicative of magma water interactions attest to water s role in fragmentation overall the pyroclast textures are consistent with water modification of a rather intense basaltic eruption which might have been of a violent strombolian or plinian style if the eruption had been dry and occurred subaerially

terrasar x interferometry reveals small scale deformation associated with the summit eruption of kīlauea volcano hawai i 1 on 19 march 2008 a small explosive eruption at the summit of k lauea volcano hawai i heralded the formation of a new vent along the east wall of halema uma u crater in the ensuing years the vent widened due to collapses of the unstable rim and conduit wall some collapses impacted an actively circulating lava pond and resulted in small explosive events we used synthetic aperture radar data collected by the terrasar x satellite a joint venture between the german aerospace center dlr and eads astrium to identify and analyze small scale surface deformation around the new vent during 2008 2012 lidar data were used to construct a digital elevation model to correct for topographic phase allowing us to generate differential interferograms with a spatial resolution of about 3 m in k lauea s summit area these interferograms reveal subsidence within about 100 m of the rim of the vent small baseline subset time series analysis suggests that the subsidence rate is not constant and over time may provide an indication of vent stability and potential for rim and wall collapse information with obvious hazard implications the deformation is not currently detectable by other space or ground based techniques

dynamics of an open basaltic magma system the 2008 activity of the halema uma u overlook vent kīlauea caldera on march 19 2008 a small explosive event accompanied the opening of a 35 m wide vent overlook vent on the southeast wall of halema uma u crater in k lauea caldera initiating an eruptive period that extends to the time of writing the peak of activity in 2008 consisted of alternating background open system outgassing and spattering punctuated by sudden short lived weak explosions triggered by collapses of the walls of the vent and conduit near daily sampling of the tephra from this open system along with exceptionally detailed observations allow us to study the dynamics of the activity during two eruptive sequences in late 2008 each sequence includes background activity preceding and following one or more explosions in september and october 2008 respectively componentry analyses were performed for daily samples to characterise the diversity of the ejecta nine categories of pyroclasts were identified in all the samples including wall rock fragments the six categories of juvenile clasts can be grouped in three classes based on vesicularity 1 poorly 2 uniformly highly to extremely and 3 heterogeneously highly vesicular the wall rock and juvenile clasts show dissimilar grainsize distributions reflecting different fragmentation mechanisms the wall rock particles formed by failure of the vent and conduit walls above the magma free surface and were then passively entrained in the eruptive plume the juvenile componentry reveals consistent contrasts in degassing and fragmentation processes before during and after the explosive events we infer a crude layering developed in the shallow melt in terms of both rheology and bubble and volatile contents beneath a convecting free surface during background activity a tens of centimetres thick viscoelastic surface layer was effectively outgassed and relatively cool while at depths of less than 100 m the melt remained slightly supersaturated in volatiles and actively vesiculating decoupled metre sized bubbles rising through the column burst through the free surface frequently ejecting fragments of the outgassed upper layer when the surface was abruptly perturbed by the rock falls existing mm sized bubbles expanded leading to the acceleration of adjacent melt upward and consecutive explosions while renewed nucleation created a minor population of 10 micron sized bubbles after each explosive event in september october 2008 this layering was re established but with decreasing vigour suggesting that the magma batch as a whole was becoming progressively depleted in dissolved volatiles

the problem of heat flow density determination from inaccurate data abstract we present an inverse method for processing rough data from oil wells in order to evaluate the heat flow density assuming that the evaluations of temperatures and thermal resistivities present gaussian uncertainties we give formulas for computing the optimal value of the heat flow density and its standard deviation a synthetic example using likely assumptions on various uncertainties leads to the prediction that a standard deviation of about 10 mw m 2 could apply to heat flow estimates from oil well data for typical situations

pigeonholing pyroclasts insights from the 19 march 2008 explosive eruption of kīlauea volcano we think conventionally of volcanic explosive eruptions as being triggered in one of two ways by release and expansion of volatiles dissolved in the ejected magma magmatic explosions or by transfer of heat from magma into an external source of water phreatic or phreatomagmatic explosions we document here an event where neither magma nor an external water source was involved in explosive activity at k lauea instead the eruption was powered by the expansion of decoupled magmatic volatiles released from deeper magma which was not ejected by the eruption and the trigger was a collapse of near surface wall rocks that then momentarily blocked that volatile flux mapping of the advected fall deposit a day after this eruption has highlighted the difficulty of constraining deposit edges from unobserved or prehistoric eruptions of all magnitudes our results suggest that the dispersal area of advected fall deposits could be miscalculated by up to 30 of the total raising issues for accurate hazard zoning and assessment eruptions of this type challenge existing classification schemes for pyroclastic deposits and explosive eruptions and in the past have probably been interpreted as phreatic explosions where the eruptive mechanism has been assumed to involve flashing of groundwater to steam

bubble nucleation in rhyolitic melts experiments at high pressure temperature and water content abstract we report the first measurement of bubble nucleation in hydrated rhyolitic melts in response to pressure release two rhyolitic obsidians one containing less than 1 of microlites of fe ti oxides and the other about 20 of various crystals were hydrated at 150 mpa and 780 850 c after saturation was reached 5 3 5 5 wt water pressure was lowered and the samples were allowed to nucleate and grow bubbles for various amounts of time before the final rapid quenching of the experiments the results demonstrate the importance of heterogeneous nucleation microlites of fe ti oxides are very efficient as sites for bubble nucleation in their presence modest nucleation was observed even after decompression by 6 10 8 bubbles cm 3 in the absence of microlites no nucleation occurred at p mpa at p 10 mpa bubbles also nucleated on crystals of biotite zircon and apatite modest nucleation 10 3 10 5 cm 3 took place even in crystal free samples but it was still heterogeneous when p exceeded 80 mpa nucleation in crystal free samples became extensive 10 5 10 7 cm 3 the lack of correlation of bubble density with either time or decompression suggests that nucleation was still heterogeneous nucleation rates were controlled mainly by the availability of sites rates were faster than 10 6 cm 3 s 1 when microlites were present and faster than 10 5 cm 3 s 1 in the absence of microlites at p 70 mpa narrow size distributions in most samples suggests that nucleation took place immediately after the pressure drop the experimental data we present here indicate that the presence or absence of efficient nucleation sites can lead to two distinct modes of bubble formation when a large number of efficient sites e g fe ti oxide are present bubble nucleation requires very little supersaturation and to a good approximation gas and magma are in equilibrium in magmas that are crystal free or contain crystals that are inefficient at nucleating bubbles very high degrees of supersaturation are required in order to initiate nucleation these two modes of exsolution may lead to contrasting styles of convection pressure build up and eruption

fragmentation mechanisms associated with explosive lava water interactions in a lacustrine environment rootless cones form when partially outgassed lava interacts explosively with external water the explosions represent an end member system that can elucidate mechanisms of explosive magma water interactions in the absence of magmatic fragmentation induced by outgassing the proportion of finely fragmented ejecta i e ash generated in rootless explosions may contribute significantly to the energy of the explosion even if the ash volume is small relative to coarser ejecta laboratory experiments indicate that the degree of melt water mixing and energy release are proportional to the abundance of blocky grains fragmented by brittle disintegration which effectively contribute thermal energy to the system to constrain the mechanisms and dynamics of rootless explosive activity we assess the nature and modes of fragmentation and ejecta characteristics through morphological textural and density analysis of rootless tephra associated with a p hoehoe lava flow in a lacustrine lake basin environment we observe strong correlations between the mean grain size and the mass percentage of both blocky negative power law trend and fluidal positive logarithmic tephra clasts of all sizes we interpret these trends as scale dependent fragmentation behavior due to the decreasing efficacy of hydrodynamic fragmentation as it occurs over finer scales especially over the ash size range additionally all analyzed beds contain fine ash sized blocky and mossy clasts which are thought to be diagnostic of a high transfer rate of thermal to mechanical energy characteristic of molten fuel coolant interactions these results agree with a recent model of rootless cone formation prior fragmentation theory and scaled laboratory experiments and therefore provide a field based analog for future experimental and modeling efforts

kinetic method for modeling vitrinite reflectance a kinetic method for the prediction of vitrinite reflectance is outlined its application allows the effect of spatially variable syndepositional and postdepositional processes to be considered through geologic time the model can be used for deterministic and probablilistic analyses and allows the determination of the temporal and spatial relationship between hydrocarbon generation and vitrinite reflectance in a basin

sensitivity of biomarker properties to depositional environment and or source input in the lower toarcian of sw germany abstract a sequence of lower toarcian sediments which are highly contrasting in their depositional environment over a 5 m depth interval has been investigated in detail for the variability of geochemical properties which are used for source characterization and oil to source correlation lower liassic marls are representative of a well aerated shallow water environment in contrast to the immediately overlying bituminous shales which are deposited under extremely reducing conditions various properties have been found to vary considerably within these two units amongst the most important are carbon isotopic composition of the kerogen the pristane phytane ratios nickel as opposed to vanadyl porphyrins and the c27 dia regular steranes although maturation within the profile does not change some of the maturation dependent biomarker properties such as the monoaromatic steroid side chain cracking and the tm ts ratio exhibit large changes which can be assigned to diagenetic processes another maturation dependent property the 20s 20r epimerization of c29 steranes exhibits smaller changes which could also be due to early diagenetic processes the study suggests that reducing and oxidizing conditions i e eh and ph in the sediment exert an influence on several biomarker precursor product pathways in maturation studies initial variations due to depositional conditions have therefore to be taken into account

role of temperature and time on organic metamorphism geologic notes organic metamorphism kerogen to hydrocarbon transformation is a function of both temperature and exposure time the effect of temperature is exponential of that of time linear at very low temperatures 266 f 130 c the reaction proceeds at such a high rate that time also plays no important role the effect of time is noticeable in the range of 158 to 212 f 70 to 100 c where intermediate conversion rates prevail stopping the heating process termination of burial freezes the rate constant of the reaction and allows organic metamorphism to continue at a linear rather than exponential rate this effect is demonstrated with the help of a simple geologic model for the first example a source rock was subjected to continuous linear heating whereas for the second the heating was terminated at a certain temperature level the model was analyzed by using the oleum scale a modified lopatin approach

oil generation in the michigan basin a biological marker and carbon isotope approach abstract palaeozoic strata in the michigan basin produced crude oils which despite their rather small volume and economic value are interesting because of the considerable age of the basin and because of the opportunity to study generation and migration of oils within a relatively simple but ancient geological setting based on n alkane profiles biological marker distributions and carbon isotope ratios the oils belong to three main families of different genetic origin and a few less important mixed types silurian oils from salina and niagara limestones have broad n alkane distributions and abundant isoprenoid hydrocarbons a strong phytane over pristane predominance and the lack of diasteranes indicate a carbonate source for these oils chemical maturity parameters show that they are more mature than the oils from the other main families from which they are also clearly distinguished by carbon isotope ratios of hydrocarbon fractions and single n alkanes oils found in the calcareous ordovician trenton formation contain n alkane cyclohexylalkane and alkyl phenanthrene distributions typical of immature oils many oils from the devonian dundee reservoirs are very similar in overall composition to the trenton oils devonian traverse oils are considered to be mainly from a devonian source with some contribution of ordovician type oil and are of intermediate maturity based on calculations using kinetic parameters of biological marker reactions considerably deeper subsidence of the devonian source rock in the past is implied

natural and artificial maturation of biological markers in a toarcian shale from northern germany steroid and terpenoid biological marker evolution was studied in a maturity sequence of lias toarcian shales and adjacent strata from three shallow boreholes in the hils syncline northern germany excellent facies control within the study area ensures the comparability of the organic matter type among the stratigraphically equivalent samples of the natural maturation series in the early hydrocarbon generation stage 0 68 r0 a strong facies dependence of some sterane and hopane isomerisation values was found within the same borehole at only a few metres depth difference the extent of sterane isomerisation at c 14 and c 17 was close to 80 in the toarcian shale and between about 55 and 65 in the adjacent dogger and lias sediments similar differences were observed for the 17 h hopane moretane and 22 29 30 trisnor 18 h neohopane 22 29 30 trisnor 17 h hopane ratios whereas the 20s 20s 20r sterane isomerisation ratio did not show this difference the observed differences mostly disappear by the end of the main oil generation phase 0 88 ro where several reactions appear to have reached completion hydrous pyrolysis experiments by variations of time and temperature were used to simulate the natural evolution of the biological markers results for sterane isomerisation in the high temperature range were compared to those published by mackenzie and mckenzie geol mag 120 417 470 1983 for natural sediments kinetic data derived for this reaction from the hydrous pyrolysis results alone or from the combination of hydrous pyrolysis and natural data sets differ markedly from those determined by mackenzie and mckenzie 1983 the hydrous pyrolysis data alone yield an activation energy for sterane isomerisation at c 20 of 131 kj mol the combined data of 170 kj mol the frequency factors are in the 105 108 s 1 range the significance of the data is discussed in detail

occurrence of 20r and 20s δ8 14 and δ14 5α h sterenes and the origin of 5α h 14β h 17β h steranes in an immature sediment abstract series of 20r and 20s 8 14 and 145 h sterenes have been identified in an immature marl from a hypersaline deposit northern apennines italy by co injection studies with authentic compounds these alkenes co occur with 13 17 5 h spirosterenes and transient 20r and 20s 14 5 h 17 h sterenes the source of these isomeric steroid alkenes is probably 7 and 8 14 5 h sterols via 2 7 2 8 14 and 2 14 5 h steradienes reduction processes in the sediment via the 14 5 h 17 h sterenes can give rise to the c 20 isomeric 5 h 14 h 17 h steranes present in anomalously high relative abundances from a thermal maturity viewpoint this diagenetic pathway is supported by laboratory rearrangement hydrogenation experiments of 5 h cholest 7 ene and from considerations of carbon number distributions of various steroid hydrocarbons in the sediment extract

excitation and resonance of acoustic gravity waves in a column of stratified bubbly magma oscillations of magma in volcanic conduits are thought to be the source of certain seismic and infrasonic signals observed near active volcanoes however the multiphase and stratified nature of magma within the conduit complicates the calculation of resonant modes that is required to interpret observations here we present a linearized mathematical framework to describe small amplitude oscillations and waves in a stably stratified column of two phase magma liquid melt and gas bubbles with a traction free upper surface a lava lake we explore the role of time dependent mass exchange between the phases depth varying fluid properties and gravity on the modes of oscillation of inviscid magma within an axisymmetric vertical conduit non equilibrium phase exchange which we refer to as bubble growth and resorption bgr is parameterized by introduction of a kinetic time scale quantifying mass exchange between the liquid and gas phases that evolves the mixture towards a state of thermodynamic equilibrium using a provably stable finite difference method we solve the eigenvalue problem for the resonance frequencies decay rates and spatial structure of the conduit eigenmodes the numerical method is then extended to time domain simulations of waves excited by internal volumetric sources in the conduit or forces applied to the surface of the lava lake we connect time dependent wave propagation simulations to the modal analysis by identifying the primary modes that are excited by representative excitation processes waves propagating through bubbly magma are dispersive and their behaviour is determined by three dimensionless parameters one quantifies the importance of buoyancy and gravitational restoring forces relative to compressibility the second quantifies differences between fluid properties e g mixture compressibility under equilibrium and non equilibrium conditions and the third compares the wave period to the bgr time scale pronounced depth variations in background fluid properties such as the transition from liquid melt with dissolved volatiles at the high pressures at depth to bubbly magma above the gas exsolution depth segment the conduit into distinct regions the longest period modes which are expressed with the largest amplitudes for typical excitation processes are most sensitive to the length of the bubbly region and properties of the bubbly magma within it while the boundary condition at the bottom of the conduit determines whether the fundamental mode is affected by the total conduit length modes localized above the exsolution depth are remarkably insensitive to the overall conduit length our analysis suggests that parameters affecting eruption style such as total volatile content and kinetic time scales of bgr along with excitation source characteristics are imprinted on long period seismic and infrasonic signals at active volcanoes

organic geochemistry of the lower toarcian posidonia shale in nw europe abstract a suite of core and outcrop samples of the lower toarcian posidonia shale and time equivalent deposits was collected from multiple locations including southern germany swabian alb and franconian alb northern germany the east coast of great britain the netherlands and luxembourg in order to study major differences and similarities with respect to depositional environment and redox potential based on molecular geochemical data all samples are immature or in the early oil window with the highest thermal maturity occurring in the profile from great britain about 0 7 vitrinite reflectance molecular geochemical maturity parameters were used to quantify maturation in detail during shale deposition in the west cleveland basin great britain h 2 s was immediately consumed in the water column by terrigenous input supplying abundant reactive iron to form pyrite these samples are characterized by low dibenzothiophene phenanthrene ratios 95 of 2 and h 2 s in contrast sulfidic euxinic bottom waters possibly with enhanced salinity and stratified water column prevailed due to a more restricted iron supply in luxembourg northern germany and the franconian alb of southern germany a more oxygenated dysoxic bottom water column associated with anoxic sediments i e anoxic dysoxic boundary close to the sediment water interface prevailed in the swabian alb of southern germany despite observed changes within each vertical profile at each site large scale regional trends in depositional conditions can be clearly identified for the early toarcian

microscopic and sedimentologic evidence for the generation and migration of hydrocarbons in toarcian source rocks of different maturities abstract this study is part of a comprehensive research program on the generation and migration of petroleum from the organic rich posidonia shale lower toarcian jurassic of the hils syncline northern germany four cores of the complete posidonia shale with varying amounts of underlying pliensbachian and overlying aalenian mudstones provided the sample base depending upon their distance from a presumed subjacent heat source the vlotho massif the cores display progressive levels of increasing thermal maturation with r 0 values 0 48 0 68 0 88 and 1 45 which span the limits of the oil generation window primary stratigraphic lithologic and petrographic attributes of the four sections are all very similar each includes a lower marlstone member about 5 m thick overlain by 11 31 m of calcareous clay shale clay minerals calcite pyrite kerogen and quartz are the principal constituents thus there is no evidence that depositional conditions of the posidonia shale varied significantly over the sampled area the omnipresence of framboidal pyrite clearly formed from the activity of sulphate reducing bacteria at or near the depositional interface is strong evidence that uniform reducing anoxic conditions were persistent at all localities this conclusion supports the assumption that initially each section was identical in overall properties and that significant differences observed in geochemistry organic macerals and other properties are directly related to the variable levels of thermal maturation kerogen from the least mature section r o 0 48 is type ii and composed of alginite a alginite b liptodetrinite bituminite and some terrestrial sporinite vitrinite and inertinite macerals with increasing maturation the alginite macerals exhibit a progressive change in morphology and fluorescence properties alginite of the least mature section displays a strong green fluorescence at 0 68 r o alginite macerals exhibit yellow orange fluorescence and are structurally well preserved although geochemical data indicate generation and loss of large amounts of oil at 0 88 r o hydrocarbon generation is almost complete and few primary alginites remain they have been transformed into secondary liptinites designated as metaalginite which have similar fluorescence colours at the highest maturation level 1 45 r o metaalginites are only weakly fluorescent with a shift toward red colours calcite and or bitumen filled bedding plane micro and macrofractures occur exclusively in the organic rich posidonia shale at maturity levels of 0 68 r o and above such features are totally absent in the immature posidonia section and in the adjacent organic poor mudstones this observation indicates that the fractures developed in direct response to overpressuring i e pore pressure exceeding lithologic load related to oil and gas generation in the source rock the fractures undoubtedly provided important expulsion conduits for generated products

continuous gravity measurements reveal a low density lava lake at kīlauea volcano hawai i abstract on 5 march 2011 the lava lake within the summit eruptive vent at k lauea volcano hawai i began to drain as magma withdrew to feed a dike intrusion and fissure eruption on the volcano s east rift zone the draining was monitored by a variety of continuous geological and geophysical measurements including deformation thermal and visual imagery and gravity over the first 14 hours of the draining the ground near the eruptive vent subsided by about 0 15 m gravity dropped by more than 100 gal and the lava lake retreated by over 120 m we used gps data to correct the gravity signal for the effects of subsurface mass loss and vertical deformation in order to isolate the change in gravity due to draining of the lava lake alone using a model of the eruptive vent geometry based on visual observations and the lava level over time determined from thermal camera data we calculated the best fit lava density to the observed gravity decrease to our knowledge the first geophysical determination of the density of a lava lake anywhere in the world our result 950 300 kg m 3 suggests a lava density less than that of water and indicates that k lauea s lava lake is gas rich which can explain why rockfalls that impact the lake trigger small explosions knowledge of such a fundamental material property as density is also critical to investigations of lava lake convection and degassing and can inform calculations of pressure change in the subsurface magma plumbing system

continuous gravity and tilt reveal anomalous pressure and density changes associated with gas pistoning within the summit lava lake of kīlauea volcano hawai i

cyclic spattering seismic tremor and surface fluctuation within a perched lava channel kīlauea volcano in late 2007 a perched lava channel built up to 45 m above the preexisting surface developed during the ongoing eruption near pu u cone on k lauea volcano s east rift zone the lava channel was segmented into four pools extending over a total of 1 4 km from late october to mid december a cyclic behavior consisting of steady lava level rise terminated by vigorous spattering and an abrupt drop in lava level was commonly observed in pool 1 we use geologic observations video time lapse camera images and seismicity to characterize and understand this cyclic behavior spattering episodes occurred at intervals of 40 100 min during peak activity and involved small 5 10 m high fountains limited to the margins of the pool most spattering episodes had fountains which migrated downchannel each spattering episode was associated with a rapid lava level drop of about 1 m which was concurrent with a conspicuous cigar shaped tremor burst with peak frequencies of 4 5 hz we interpret this cyclic behavior to be gas pistoning and this is the first documented instance of gas pistoning in lava well away from the deeper conduit our observations and data indicate that the gas pistoning was driven by gas accumulation beneath the visco elastic component of the surface crust contrary to other studies which attribute similar behavior to the periodic rise of gas slugs the gas piston events typically had a gas mass of about 2 500 kg similar to the explosions at stromboli with gas accumulation and release rates of about 1 1 and 5 7 kg s 1 respectively the time averaged gas output rate of the gas pistoning events accounted for about 1 2 of the total gas output rate of the east rift zone eruption

insights into shallow magmatic processes at kīlauea volcano hawaiʻi from a multiyear continuous gravity time series continuous gravity data collected near the summit eruptive vent at k lauea volcano hawai i during 2011 2015 show a strong correlation with summit area surface deformation and the level of the lava lake within the vent over periods of days to weeks suggesting that changes in gravity reflect variations in volcanic activity joint analysis of gravity and lava level time series data indicates that over the entire time period studied the average density of the lava within the upper tens to hundreds of meters of the summit eruptive vent remained low approximately 1000 1500 kg m3 the ratio of gravity change adjusted for earth tides and instrumental drift to lava level change measured over 15 day windows rose gradually over the course of 2011 2015 probably reflecting either 1 a small increase in the density of lava within the eruptive vent or 2 an increase in the volume of lava within the vent due to gradual vent enlargement superimposed on the overall time series were transient spikes of mass change associated with inflation and deflation of k lauea s summit and coincident changes in lava level the unexpectedly strong mass variations during these episodes suggest magma flux to and from the shallow magmatic system without commensurate deformation perhaps indicating magma accumulation within and withdrawal from void space a process that might not otherwise be apparent from lava level and deformation data alone continuous gravity data thus provide unique insights into magmatic processes arguing for continued application of the method at other frequently active volcanoes

shallow degassing events as a trigger for very long period seismicity at kīlauea volcano hawai i the first eruptive activity at k lauea volcano s summit in 25 years began in march 2008 with the opening of a 35 m wide vent in halema uma u crater the new activity has produced prominent very long period vlp signals corresponding with two new behaviors episodic tremor bursts and small explosive events both of which represent degassing events from the top of the lava column previous work has shown that vlp seismicity has long been present at k lauea s summit and is sourced approximately 1 km below halema uma u by integrating video observations infrasound and seismic data we show that the onset of the large vlp signals occurs within several seconds of the onset of the degassing events this timing indicates that the vlp is caused by forces sourced at or very near the lava free surface due to degassing transmitted down the magma column and coupling to the surrounding rock at 1 km depth

magma oscillations in a conduit reservoir system application to very long period vlp seismicity at basaltic volcanoes 2 data inversion and interpretation at kīlauea volcano

heterogeneous bubble nucleation on fe ti oxide crystals in high silica rhyolitic melts the nucleation of h2o bubbles in magmas has been proposed as a trigger for volcanic eruptions to determine how bubbles nucleate heterogeneously in silicate melts experiments were carried out in which high silica rhyolitic melts were hydrated at 740 800 c and 50 175 mpa decompressed by 20 70 mpa and held at the lower pressures for 10 s before being quenched the hydration conditions were subliquidus and all samples contain blocky magnetite needle shaped hematite plagioclase magnetite is abundant at 800 c and high pressures whereas hematite becomes more abundant at lower temperatures and pressures bubbles nucleated in a single event in all samples with the number density nt of bubbles varying between 2 107 and 1 109 cm 3 at low degrees of supersaturation one to a few bubbles nucleate on faces of magnetite but at medium to high degrees of supersaturation multiple bubbles nucleate on single magnetite grains on hematite one to a few bubbles nucleated at the ends of the needle shaped crystals at medium supersaturations but formed along their entire lengths at high supersaturations nt increases as water diffusivity decreases indicating that the number of bubbles nucleated is influenced by their growth which depletes the melt with respect to h2o and lowers supersaturation if volcanic eruptions are triggered by bubble formation in magmas stored in shallow level magma chambers then the supersaturations needed for heterogeneous nucleation suggest that only small amounts of crystallization are needed whereas homogeneous nucleation is unlikely to trigger eruptions

stability of foams in silicate melts bubble coalescence and the spontaneous disruption of high porosity foams in silicate melts are the result of physical expulsion of interpore melt syneresis leading to bubble coalescence and diffusive gas exchange between bubbles melt expulsion can be achieved either along films between pairs of bubbles or along plateau borders which represent the contacts between 3 or more bubbles theoretical evaluation of these mechanisms is confirmed by experimental results enabling us to quantify the relevant parameters and determine stable bubble size and critical film thickness in a foam as a function of melt viscosity surface tension and time foam stability is controlled primarily by melt viscosity and time melt transport leading to coalescence of bubbles proceeds along inter bubble films for smaller bubbles and along plateau borders for larger bubbles thus the average bubble size accelerates with time in silicate melts the diffusive gas expulsion out of a region of foam is effective only for water and even then only at small length scales as the diffusion of co2 is negligible the results of our analyses are applicable to studies of vesicularity of lavas melt degassing and eruption mechanisms

biological marker compounds as indicators of the depositions history of the maoming oil shale abstract the eocene maoming oil shale from guangdong province occurs as a laterally uniform stratigraphic section typically 20 25 m thick from which the aliphatic hydrocarbon constituents of six representative samples were investigated using gc and c gc ms the sediments evaluated included the basal lignite a vitrinite lens from the overlying claystone and four intervals from the massive oil shale bed as expected the lignite and vitrinite differ markedly from the oil shales the lignite is dominated by bacterial hopanoids and components of higher plant origin including c 29 steroids and triterpenoids such as oleanenes visually the oil shale samples show corroded and degraded phytoclasts spores wispy particles of fluorescent organic material attributable to dinoflagellates and especially in the uppermost sample colonial algal bodies the distributions of biological markers in the oil shales show many features in common notably a dominance of dinoflagellate derived 4 methylsteroids and a significant proportion of higher plant derived n alkanes with marked odd over even carbon number predominance overall they exhibit several features that resemble characteristics of the messel shale the hydrocarbons of the lowest shale horizon suggest that there may have been a gradual transition between deposition of the original peat and the subsequent oil shales the aliphatic hydrocarbons of the uppermost shale are dominated by a number of c 31 and c 33 botryococcane homologues and other unusual branched alkanes possibly derived from green algae all of the samples are immature overall molecular and microscopic examination of the stratigraphic succession of the maoming oil shale suggests a shallow lacustrine environment within which peats were deposited this lake subsequently deepened to support abundant algal populations especially dinoflagellates culminating in a dominance of botryococcoid algae

triggering mechanism at the origin of paroxysms at stromboli aeolian archipelago italy the 5 april 2003 eruption 1 the 5 april 2003 paroxysm of stromboli occurred during an ongoing effusive episode that initiated on 28 december 2002 it mainly consisted in two a few seconds apart cannon like explosions followed by a vertical gas pyroclast jet the hk basaltic pumices that were generated mainly contain crystals that were drained back from the shallow parts of the plumbing system few euhedral to anhedral olivines and their melt inclusions testify to crystal fractionation and possible mixing between hk basaltic magmas variable in their volatile content extent of evolution and volume the most primitive volatile rich hk basalt term is sizable at micrometer scale only we propose that the triggering mechanism of the paroxysm was initiated at moderate pressures 240 mpa and was related to bubble driven ascent of magma blobs the latter rose through and interacted with overlying slightly more evolved melts and finally mingled with the shallow crystal rich magma just before the eruption

constraints upon the use of 24 ethylcholestane diastereomer ratios as indicators of the maturity of petroleum abstract we have measured the ratio 20 s 20 s 20 r the relative abundances of the 24 ethylcholestane diastereomers in nine pairs of closely adjacent samples of coals and shales which represent a range of geological ages and different sedimentary basins within some sample pairs the values were similar in others they were very different the differences between the values for coal and shale members of pairs are attributed to effects of sediment matrix type upon the processes which cause the changes in the ratio the kinetic parameters e and s for the changes in the 20 s 20 s 20 r ratio have been estimated using previously published measurements for the hod 1 well shale and jupiter 1 well coal and these values were found to be very different the process in the coal sequence has a much lower activation energy and is therefore less responsive to differences in heating rate regimes than is the case for the shale sequence values of the 20 s 20 s 20 r ratio in two other sequences are shown to be consistent with this explanation and the implications of these observations for the practice of using sterane isomer ratios for assessment of maturity of sediments and crude oils are discussed

determining flow type shear rate and shear stress in magmas from bubble shapes and orientations we interpret the shear environments that produced bubble textures in obsidian samples using the results of theoretical numerical and experimental studies on the deformation of bubbles in shear flows in particular we use the shapes and orientations of bubbles vesicles in obsidian to estimate shear rates and shear stresses and assess flow type simple vs pure shear this technique can be used to determine shear rates in volcanic conduits the origin of pyroclastic obsidian and the emplacement history and dynamics of obsidian flows a spherical bubble in a viscous fluid subjected to a low reynolds number steady flow field deforms until it reaches a steady shape and orientation bubble deformation is governed by the competing stresses from shearing that deforms and surface tension that rerounds the ratio of these stresses is the capillary number ca because the relationships among ca bubble orientation and shape differ for pure and simple shear we can distinguish between these flow types using preserved bubble geometries furthermore because ca is a function of shear rate we can use relationships between ca and the magnitude of deformation to calculate shear rates when melt viscosity and surface tension are known to demonstrate the potential of the technique we examine three obsidian samples chosen for diversity of origin and texture two of the samples have low crystallinities and banding defined by layers of different vesicularity bubble geometries in a spatter fed obsidian flow sample record deformation by pure shear whereas a juvenile obsidian clast from a pyroclastic fall deposit records predominantly simple shear a third sample from an obsidian flow has banding marked by variable concentrations of microlites and shows bubble deformation most consistent with dominantly simple shear and some bubble relaxation the highest shear rate and shear stress are recorded by the pyroclastic obsidian shear rate 0 01 s 1 shear stress 60 kpa the spatter fed obsidian flow sample deformed by pure shear records the lowest shear stress 4 9 kpa the obsidian flow sample deformed by simple shear records the lowest shear rates 10 6 6 10 6 9 s 1 despite high shear stresses 22 36 kpa because of the high viscosity of the degassed rhyolite 0 13 wt water

integrated petrochemical and geophysical data reveals thermal distribution of the feeding conduits at stromboli volcano italy 1 samples of scoriae erupted at stromboli volcano during its persistent strombolian activity were collected between 2005 and 2008 chemical and mineralogical compositions were obtained on products erupted from the three main crater sectors sw central and ne small chemical variations indicate a different degree of evolution coupled with small difference of magma temperature 10 c analysis of the acoustic data for the same time period as the scoria sampling indicates that puffing a persistent overpressurized bubble degassing was on average mainly observed at the central craters and at times moved to the ne sector the cross check of the two independent data sets allowed us to assess correlation between composition of products and puffing activity at vents the hotter products are always erupted from the vents where puffing occurs indicating that slightly higher temperature can be the expression of an enhanced two phase bubble flow dynamics

organic matter provenance and depositional environment of marine to continental mudstones and coals in eastern ordos basin china evidence from molecular geochemistry and petrology abstract cyclothems composed of interbedded mudstone coal and sandstone layers make up the taiyuan and shanxi formations in the late carboniferous to early permian in north china under a marine to continental depositional environment the cyclothems act as important fossil energy hosts such as coalbeds hydrocarbon source rocks and unconventional natural gas reservoirs organic geochemistry and petrology of mudstones and coals in the taiyuan and shanxi formations in the eastern ordos basin were studied to reveal the organic matter sources and paleoenvironments total organic carbon toc contents vary from 1 1 wt mudstone to 72 6 wt coal the samples are mainly within the oil window with the tmax values ranging from 433 to 469 c organic petrology and source biomarkers indicate that the mudstones were sourced from a mixed organic matter input and terrigenous organic matter predominates over aquatic organic matter the coals are mostly sourced by terrigenous organic matter inputs high concentrations of hopanes argue for a strong bacterial input some m z 217 mass chromatograms have peaks at the hopanes retention times as a result of high hopane to sterane ratios these hopane derived peaks do not interfere the identification of the steranes because the hopanes and the steranes have different retention times maturity dependent biomarkers demonstrate that the samples have been thermally mature which agree with the tmax values anomalously low c29 20s 20s 20r and c29 sterane ratios are present in all the samples and are interpreted as due to the terrigenous organic matter input or the coal related depositional environment in addition biomarkers and iron sulfide morphology indicate that the organic matter of the mudstones deposited in a proximal setting with shallow brackish fresh water bodies with consideration of preservation of organic matter the redox conditions are dysoxic redox oscillations resulted in the records of oxic conditions in some samples finally the coals and the mudstones mainly generate gas and have poor oil generative potential

isomerization and aromatization of hydrocarbons and the paleothermometry and burial history of alberta foreland basin the hydrocarbon fractions extracted from 34 cretaceous shales from a cross section of the alberta basin have been analyzed using gas chromatography mass spectrometry particular attention has been paid to the extents to which three temperature dependent geochemical reactions have occurred isomerization of steranes at c 20 isomerization of hopanes at c 22 and aromatization of c ring monoaromatic steroid hydrocarbons estimates of the kinetic constants of these reactions have been published previously therefore constraints may be put on the maximum temperature experienced by the sediments and the time for which this temperature existed the subsidence history of the alberta basin is simulated using a lithospheric flexural model in which the lithosphere is downwarped by the loading of overthrusts in the adjacent cordilleran orogene when overthrusting ceases erosion of the overthrusted region lightens the load and the lithosphere rebounds causing uplift and erosion of the foreland basin major overthrusting and therefore deep burial occurred during the laramide orogeny of late cretaceous to mid tertiary time the extent of chemical reactions suggests that the present maturity of the cretaceous shales is almost entirely a reflection of burial and heating during the laramide orogeny in a manner that is consistent with the lithospheric flexural model for the basin the chemical reactions suggest that the paleogeothermal gradient increased systematically away from the foothills of the rocky mountains along the profile studied that a similar trend exists today suggests this situation has persisted at least since the early stages of the laramide orogeny this variation is best explained by advective heat transfer water entering the basin in the faulted and topographically higher foothills is heated at depth thereby depressing the geothermal gradient this hot water then migrates updip and away from the foothills to enhance the geothermal gradients in the edmonton area and beyond the route of the migrating hot water is similar to that taken by hydrocarbons from devonian and cretaceous source rocks in the deep basin updip to the oil and tar sand accumulations of alberta