the december 2002 july 2003 effusive event at stromboli volcano italy insights into the shallow plumbing system by petrochemical studies abstract on 28 december 2002 the persistent strombolian activity at stromboli was interrupted by the sudden onset of lava emission onto the sciara del fuoco slope a horseshoe shaped depression on the nw flank of the volcano the effusive episode went on until 22 july 2003 and produced a cumulative volume of lavas of 11 106 m3 ranking the event as the largest occurred in the past 30 years the eruptive vents were mainly located in the ne sector of the sciara del fuoco depression at an elevation of 550 600 m a s l on 30 december the eastern portion of the sciara del fuoco collapsed producing a tsunamigenic landslide on 5 april 2003 a paroxysmal eruption occurred at the summit craters during which crystal poor pumiceous products were emitted the paroxysm did not interrupt the lava emission the strombolian activity at the summit craters gradually resumed starting from march 2003 and fully recovered by the end of july periodic sampling of the lava at the active vents was carried out during the entire effusive event all the analysed samples are shoshonitic basalts sio2 48 5 50 4 wt k2o 2 1 2 4 wt in the range of composition observed in the products erupted during the past 20 years they bear about 50 vol zoned crystals of plagioclase an90 60 diopside augite and olivine fo70 73 in a compositionally homogeneous shoshonitic groundmass bulk rock major and trace element contents measured with different methodologies at different laboratories show only minor variations sr and nd isotope ratios are close to those of the crystal rich scoriae erupted in the previous years despite of small compositional variations related to the emptying of the zoned topmost portion of the conduit the data show that the switch from strombolian explosive to effusive activity is not associated with changes of the textures and composition of the erupted products slight but somehow systematic variation of trace elements and isotope ratios between products erupted before and after the 5 april eruption are likely accounted by limited mixing between the fresh volatile rich crystal poor magma erupted as pumice during the paroxysm and the volatile poor crystal rich magma feeding the lava flow the uniform composition of the erupted lava indicates the presence of a large volume of well mixed crystal rich homogeneous magma residing in the shallow plumbing system of the volcano two possible trigger mechanism of the effusive event are here proposed i a discrete input of fresh magma into the lower part of the shallow magmatic system occurred some months before the eruption and was followed by crystallisation degassing and mixing with the crystal rich shallow magma which re homogenised the system these processes eventually led to the rise of the magma at a higher level and failure of the conduit walls ii the onset of the effusion may represent the consequence of a gradual rise of the magma level in the conduits occurred in the past two decades the process of progressive refilling initiated after the 1985 effusive eruption ultimately culminated on 28 december 2002 with the failure of the conduit wall and the opening of vents 100 150 m below the summit craters

on estimating crystal shape for crystal size distribution analysis abstract the accurate determination of three dimensional crystal habit long intermediate short axis is a requirement to enable correct stereological conversion of 2d crystal size measurements to true 3d crystal size distributions csds in this contribution we introduce a database and spreadsheet program which provides an objective estimate of true crystal habit from raw 2d measurements the database compares the sample s 2d measurements with 2d shape curves for random sections through 703 different habits the output gives the five best match curves and corresponding crystal habits based on a least squares fit between sample and database the minimum sample size required to give good shape estimates from 2d data is tested using random runs of sections through known shapes with different population densities at minimum 75 crystal sections are required to robustly determine crystal habit for csd measurements if crystals are tabular in shape with more acicular shapes requiring a minimum of 250 sections suggesting a sample size of 250 sections to be used for shape determination in natural examples where the true 3d shape is unknown the crystal habit of the best match shape can then be used to provide a more robust csd analysis which in turn can be used to investigate magmatic processes with more certainty

paleomagnetism and 40ar 39ar chronology of lavas from meseta del lago buenos aires patagonia 1 as part of a larger project investigating paleosecular variation in south america 36 paleomagnetic and geochronology sites were sampled on meseta del lago buenos aires 47 s 289 e 200 km east of the active plate boundary in southern patagonia basaltic lavas sampled range from late miocene to late pleistocene yet most of the individual lavas are younger than 3 3 ma thirty eight isochron ages determined via 40ar 39ar incremental heating experiments frame a detailed stratigraphy the isochron ages range from 66 9 4 1 ka for scoria cones and a youthful lava flow in the rio pinturas valley to 10 12 1 35 ma for thick plateau forming lavas exposed along the se edge of the meseta the last 3 3 myr is characterized by seven episodes of volcanism at ca 3 2 3 0 ma 2 4 ma 1 7 ma 1 35 ma 1 0 ma 750 ka 430 330 ka and 110 ka the bulk of lavas forming the surface of the meseta erupted in the last 1 7 my all sites have stable magnetization and after step demagnetization using either thermal or alternating field techniques yield characteristic directions held by magnetite and or titanomagnetite ten sites have distinct transitional directions defined by pole latitudes 55 and associated ages indicate possible connections to known reversals within the matuyama chron including the onset of the jaramillo subchron 1 016 0 01 ma the cobb mountain subchron 1 25 0 03 ma the ontong java 1 event 1 37 0 03 ma and the termination of the olduvai subchron 1 72 0 02 ma remaining sites are divided into normal 14 sites and reversed polarity 12 sites at 95 confidence mean normal and reversed directions overlap the mean direction for 26 sites is i 63 0 d 3 4 95 5 4 which is indistinguishable from the expected geocentric axial dipole paleosecular variation measured by the dispersion of virtual geomagnetic poles about the rotation axis is higher than expected at 20 0 model g using data for the past 5 myr predicts a dispersion of 17 this discrepancy may be due to true dispersion of the field in southern south america or it may be an artifact of inadequate sampling

crystal size distributions of plagioclase in lavas from the july august 2001 mount etna eruption during the 2001 eruption of mount etna two independent vent systems simultaneously erupted two different lavas the upper vents system uv opened between 3100 and 2650 m a s l emitted products that are markedly porphyritic and rich in plagioclase while the lower vents system lv opened at 2100 and 2550 m a s l emitted products that are sparsely porphyritic with scarce plagioclase in this study the crystal size distributions csds of plagioclase were measured for a series of 14 samples collected from all the main flows of the 2001 eruption the coefficient of r 2 determination was used to evaluate the goodness of fit of linear models to the csds and the results are represented as a grid of r 2 values by using a numerical code developed ad hoc r 2 diagrams suggest that the 2001 products can be separated into two main groups with slightly different characteristics plagioclase csds from the uvs can be modeled by three straight lines with different slopes while the plagioclase csds from the lvs are largely concave we have interpreted the csds of the uvs as representing three different populations of plagioclases i the large phenocrysts type i which started to crystallize at lower cooling rate in a deep reservoir from 13 to 8 months before eruption onset ii the phenocrysts type ii which crystallized largely during continuous degassing in a shallow reservoir and iii the microlites which crystallized during magma ascent immediately prior to the eruption the plagioclase csd curves for the lvs lava are interpreted to reflect strong and rapid changes in undercooling induced by strong and sudden degassing

paleomagnetic secular variation of the honolulu volcanic series 33 700 ka o ahu hawaii abstract thirty years after the precursory study j geophys res 13 1972 2129 19 dated sites with ages between 33 and 700 ka from the honolulu volcanic series hvs southeastern part of the island of o ahu in hawaii have been revisited a minimum of 16 samples per site have been studied using alternating field and thermal demagnetization both techniques yielded similar results coeval flows could be identified by relying on their k ar ages similar rock magnetic parameters and identical directions reliable paleomagnetic directions have been obtained for all sites two site mean directions which deviate significantly from the geocentric axial dipole gad correspond to periods of low dipole intensity and may thus reflect a larger proportion of non dipole components the average value of the site mean declinations is 358 6 while the mean inclination of 34 5 is 3 5 lower than the value of the geocentric dipole at the site the north geographic pole is largely contained within the uncertainties surrounding the position of the mean virtual geomagnetic poles vgp latitude the deviation of 2 3 in latitude is further reduced by 1 if we consider the existence of a standing small quadrupole component that would represent 5 of the dipole the angular dispersion of the virtual geomagnetic poles 12 is only 1 7 smaller than expected at the site latitude and thus not a firm indication of low secular variation thus there is no significant difference with the global time averaged field for the same period other recent results from the same island for an older time period display similar characteristics but with reverse and normal polarities that are not perfectly antipodal it is suggested that this could result from changes in the mean dipole field intensity during these time intervals this points out to the importance of scrutinizing the entire field vector in any analysis of paleosecular changes these results do not bring any evidence for a persistent field anomaly under hawaii

crystal size distribution csd in rocks and the kinetics and dynamics of crystallization crystal size distributions csds measured in metamorphic rocks yield quantitative information about crystal nucleation and growth rates growth times and the degree of overstepping t of reactions during metamorphism csds are described through use of a population density function n dn dl where n is the cumulative number of crystals per unit volume and l is a linear crystal size plots of ln n vs l for olivine pyroxene and magnetite in high temperature 1000 c basalt hornfelses from the isle of skye define linear arrays indicating continuous nucleation and growth of crystals during metamorphism using the slope and intercept of these linear plots in conjunction with growth rate estimates we infer minimum mineral growth times of less than 100 years at t 10 c and nucleation rates between 10 4 and 10 1 cm3 s garnet and magnetite in regionally metamorphosed pelitic schists from south central maine have csds which are bell shaped we interpret this form to be the result of two processes 1 initial continuous nucleation and growth of crystals and 2 later loss of small crystals due to annealing the large crystals in regional metamorphic rocks retain the original size frequency distribution and may be used to obtain quantitative information on the original conditions of crystal nucleation and growth the extent of annealing increases with increasing metamorphic grade and could be used to estimate the duration of annealing conditions if the value of a rate constant were known finally the different forms of crystal size distributions directly reflect differences in the thermal histories of regional vs contact metamorphosed rocks because contact metamorphism involves high temperatures for short durations resulting csds are linear and unaffected by annealing similar to those produced by crystallization from a melt because regional metamorphism involves prolonged cooling from high temperatures primary linear csds are later modified by annealing to bell shapes

magma dynamics during the 2007 stromboli eruption aeolian islands italy mineralogical geochemical and isotopic data abstract after the 6 month long effusive event of 2002 2003 a new lava effusion occurred at stromboli between 27 february and 2 april 2007 despite the different durations approximately the same volume of magma was emitted in both eruptions in the order of 107 m3 a paroxysmal eruption occurred at the summit craters in both the 2002 2003 and 2007 episodes during which a significant amount of low porphyritic lp volatile rich magma was erupted in both cases the paroxysm did not interrupt the lava emission here we present compositional data including texture mineralogy chemistry and sr and nd isotope ratios of bulk rock groundmass and separated minerals of lavas erupted in 2007 together with chemistry and sr and nd isotope composition of the pumices emitted during the 15 march paroxysm as a whole the lavas have the same texture and chemistry that characterize the highly porphyritic hp products usually erupted at stromboli during normal strombolian activity and effusive events compared to the previous hp products the 2007 lavas show minor but systematic mineralogical and isotopic variations which are consistent with a modest increase of the magma supply rate of the volcano compositional variations during the entire duration of the event are very modest glass chemistry changes in lavas erupted in the second half of march can be explained by the minor mixing between the volatile rich lp magma rising through the shallow magmatic system during the 15 march paroxysm and the degassed residing hp magma a first conclusion of this study is that there is no compositional evidence supporting major changes in the magma dynamics of the volcano accompanying the effusive activity as also suggested for the 2002 2003 event the activity of stromboli is controlled by a steady state feeding system in which refilling mixing degassing and crystallization at shallow level continuously operate with modest oscillations in the magma supply rate switching between normal strombolian and effusive activity is related to periods of relatively more vigorous refilling of the shallow system leading to progressive pressure increase in the upper conduits associated with only minor compositional variations in the erupted products

crystal recycling in the steady state system of the active stromboli volcano a 2 5 ka story inferred from in situ sr isotope and trace element data in situ sr isotope data by microdrilling coupled with major and trace element analyses have been performed on plagioclase and clinopyroxene from seven samples collected during the 2002 2003 eruptive crisis at stromboli volcano aeolian islands italy on 28 december 2002 the persistent moderate explosive activity was broken by an effusive event lasting about 7 months a more violent explosion paroxysm occurred on 5 april 2003 two magma types were erupted namely a volatile poor and highly porphyritic magma hp magma poured out as scoria or lava and a volatile rich phenocryst poor magma lp magma found as pumice lp magma differs from the hp magma also for its slightly less evolved chemistry the groundmass composition and the lower sr isotope ratios micro sr isotope data show the presence of zoned minerals in strong isotope disequilibrium as previously found in products erupted in 1984 1985 and 1996 ad with 87sr 86sr values generally decreasing from cores to rims of minerals only some outer rims testify for equilibrium with the host groundmass the internal mineral zones with high sr isotope ratios 0 70665 0 70618 are interpreted as antecrysts crystallised during the previous activity and recycled in the present day system since the opening shoshonitic activity of the recent period which occurred at about 2 5 ka ago this result has implications for the dynamics of the present day plumbing system of stromboli at intermediate pressure about 2 3 km depth and allows us to propose a model whereby an hp magma reservoir is directly interconnected at the bottom with a cumulate crystal much reservoir efficient mixing between residing hp and input lp magmas can occur in this reservoir due to more similar rheological characteristics of the two magmas than in the conduit where crystallisation is enhanced by degassing antecrysts and possibly melts re enter in the hp magma reservoir both from the bottom recycled by ascending lp magmas crossing the crystal mush and from the top recycled by descending degassed and dense hp magma residual of the periodic strombolian explosions at the surface the isotope variation measured in the groundmasses allows calculating the proportion of the lp magma entering the shallow hp magma reservoir at 20 from this proportion we estimate that the total volume of lp magma input during 2002 2003 closely matches the magma volume erupted in the effusive event suggesting a steady state system at broadly constant volume the comparison with estimates of the lp magma volume ejected by the paroxysm indicates that the lp magma amount directly reaching the surface during the 5 april paroxysm is minimal with respect to that entering the system

crystalmom a new model for the evolution of crystal size distributions in magmas with the quadrature based method of moments nucleation and growth of crystals and the resulting crystal size distribution play a fundamental role in controlling the physical properties of magmas and consequently the dynamics of the eruptions in the past decades laboratory experiments demonstrated that size and shape of crystals strongly control the physical properties of magma and lava additionally natural and experimental samples are usually characterized in terms of their crystal size distribution to link it with physical processes that are not directly observable such as cooling or decompression mechanisms in this paper we present crystalmom a new predictive model based on the quadrature based method of moments developed for studying the kinetic of crystallization in volcanic systems the quadrature based method of moments well established in the field of chemical engineering represents a mesoscale modelling approach that rigorously simulates the space time evolution of a distribution of particles by considering its moments the method is applied here for the first time for studying the equilibrium disequilibrium crystallization in magma modelling the temporal evolution of the moments of a crystal size distribution the model verified against numerical and experimental data represents a valuable tool to infer the cooling and decompression rates from the crystal size distribution observed in natural samples

dynamics of the december 2002 flank failure and tsunami at stromboli volcano inferred by volcanological and geophysical observations 1 the 2002 effusive flank eruption at stromboli volcano started on december 28 after several months of strong explosive activity at the summit craters on december 30 the seismic network recorded two large flank failures and associated tsunami waves this is the first time that a flank collapse and tsunami and their associated phenomena have been recorded by a multi disciplinary monitoring system volcanological and geophysical monitoring as well as thermal surveys performed immediately before and after the failure allowed us to define and interpret the sequence of events the still on going eruption has provided for the first time the opportunity to look into the dynamics of stromboli s effusive eruptions flank failure and landslide formation and their potential hazard

onset of the persistent activity at stromboli volcano italy stromboli known worldwide as the lighthouse of the mediterranean is commonly believed to have been in a state of persistent activity for the past 2000 2500 years however historical sources older than 1000 a d are not accurate enough to assess if the activity of the volcano was exactly the same as we see at present in order to attempt to identify the onset of the present eruptive regime and assess if it has been maintained with the same characteristics through time stratigraphic and radiometric studies of the recent tephra deposits were undertaken up to 4 m deep stratigraphic trenches dug at a height of approximately 500 m on the ne flank of the volcano exposed a conformable tephra pile containing charcoal fragments one of the most interesting finds was the discovery of a 7 cm thick weathered bed rich in organic matter thin palaeosol approximately 3 m below the surface the sequence underneath the palaeosol consists of decimetre thick lapilli fallout beds alternating with ash deposits bearing small charcoals with calibrated ages of between the fourth century b c and the first century a d the sequence above the palaeosol is charcoal free and consists of coarse ash deposits with discrete centimetre thick lapilli fallout beds composed of crystal poor golden pumice and subordinate crystal rich black scoriae similar to scoria pumice pairs emitted during the more energetic explosions of the present day activity the data collected indicate that between the third and seventh centuries a d after a period of quiescence the activity resumed with an eruptive style identical to the present one we conclude that the lighthouse of the mediterranean actually began its activity in a period much later than previously thought

combining csd and isotopic microanalysis magma supply and mixing processes at stromboli volcano aeolian islands italy integrating isotopic microanalysis with other analytical techniques creates powerful new methodologies for understanding the evolution of rock samples at the sub grain scale here we present crystal size distribution csd data for a 26 000 year old sample from stromboli volcano and accompanying isotopic microanalysis of the phenocrysts a technique called the icsd plot is introduced which given stated assumptions allows the integration of both sets of data to generate timelines of isotopic evolution through the volcanic system the combined approach is powerful allowing investigation of the magma supply mixing crystallisation and contamination processes prior to eruption of a volcanic sample for stromboli volcano the combined analysis suggests that the change in magma type following a cone collapse took roughly five years to complete similar to the timescale of changes seen in recent decades

dynamics of magma mixing and degassing recorded in plagioclase at stromboli aeolian archipelago italy crystal rich materials scoriae and lava flows emitted during the 1985 2000 activity of stromboli were taken into consideration for systematic study of bulk rock matrix glass chemistry and in particular for the study of chemical and textural zoning of plagioclase the most abundant mineral phase over the considered time period bulk rock composition remained fairly constant in both major sio2 49 2 50 9 wt and k2o 1 96 2 18 wt and trace elements the quite constant chemistry of matrix glasses also indicates that the degree of crystallization of magma was maintained at around 50 vol plagioclase ranges in composition between an62 and an88 and is characterized by alternating 10 100 m thick bytownitic and labradoritic concentric layers although the dominant and representative plagioclase of scoriae is an68 the labradoritic layers an62 70 show small scale 1 5 m oscillatory zoning are free of inclusions and appear to record episodes of slow crystal growth in equilibrium with a degassed liquid having the composition of the matrix glass in contrast the bytownitic layers an70 an88 are patchy zoned show sieve structure with abundant micrometric glass inclusions and voids and are attributed to rapid crystal growth

using quantitative textural analysis to understand the emplacement of shallow level rhyolitic laccoliths a case study from the halle volcanic complex germany received february 1 2002 accepted november 8 2002in qualitatively homogeneous magmatic bodies quantitativetextural analysis such as crystal size distribution modalabundance and spatial distribution pattern analyses allowstheir internal heterogeneity to be measured and interpreted inthis study these methods are applied to samples from a 300mdrill core through one of the porphyritic rhyolitic laccoliths petersberg unit of the 300ma halle volcanic complex germany qualitatively thegeochemicallyhomogeneouspetersbergunit does not show much textural variation quantitatively however the crystal size distributions of the three most commonphenocryst phases orthoclase plagioclase and quartz suggestcontinuouscrystalgrowthduringmagmaascentandemplacement but different growth histories of the phenocryst phases throughoutthe genesis of the laccolith in situ cooling did not affect thephenocryst population size distributions of the phenocrysts varyon a centimetre to decimetre scale but are similar on the scale ofthe laccolith the modal abundance of the phenocryst phases isvery similar throughout the drill core quantification of thespatial distribution of phenocrysts however reveals a trend forclusteringtowardstheinteriororupperpartofthelaccolith whichis attributed to flow and shear processes during emplacement anddiscontinuities in the interior relating to the intrusion of differentmagmapulses circularstatisticsoftheorientationoflongaxesofcrystals reveala weakalignment oftheorthoclaseand plagioclasephenocrystsonthesamplescaleasaresultofflowinthemagmainspite of little acicularity in general laccoliths can be fedby several pulses of magma without major cooling betweenbatches

ar ar ages from transitionally magnetized lavas on la palma canary islands and the geomagnetic instability timescale 1 a detailed study of 43 lava flows comprising two stratigraphic sequences exposed along the north and south walls of barranco de los tilos on the island of la palma canary islands reveals a complex temporally segmented record of geodynamo behavior that contains no less than three distinct geomagnetic events the matuyama brunhes m b reversal is recorded in five transitionally magnetized lava flows from the north tn section the isochrons obtained from three of the lower four m b lavas are defined by 14 incremental heating experiments that together with a previous age determination yielded a weighted mean of 798 4 6 2 ka all uncertainties 2s in addition a 780 3 10 3 ka isochron was determined for the overlying transitionally magnetized flow indicating that it was erupted during a distinctly younger portion of the transition near the base of the south ts section one finds a sequence of weakly magnetized flows associated with virtual geomagnetic pole vgp positions in the southwest indian ocean between latitudes 56 s and 65 s suggesting instability of the geomagnetic field beyond that of typical secular variation 40 ar 39 ar isochrons from three of these flows defined by 11 separate incremental heating experiments gave a weighted mean of 822 2 8 7 ka this anomalous field behavior recorded 24 11 kyr prior to the m b reversal may coincide with an event featured in several marine sediment records directly above two normal polarity flows 40 ar 39 ar isochrons of 751 9 18 1 ka and 675 0 15 7 ka are nine transitionally magnetized lavas having magnetization directions associated with low to midlatitude vgps spanning 23 60 n these flows are then capped by a single flow possessing normal polarity based on 12 incremental heating experiments 40 ar 39 ar isochrons of five of these nine lavas along with the uppermost flow gave a weighted mean age of 580 2 7 8 ka for this period of transitional to normal field behavior from these same transitional lavas quidelleur et al 1999 reported three unspiked k ar ages with a weighted mean of 602 24 ka and proposed a new event called the la palma excursion however the 40 ar 39 ar age presented here is three times more precise than the k ar age and is indistinguishable at the 95 confidence level from the 40 ar 39 ar age of a lava from the snake river plain idaho that originally defined the big lost event transitional field behavior of similar age observed in astronomically dated marine cores further establishes that the big lost event recorded at la palma was indeed global in extent rigorous temporal and geomagnetic constraints for several additional periods of geomagnetic field instability during the last several million years will comprise a geomagnetic instability timescale that can be factored confidently into models of the dynamo process index terms 1560 geomagnetism and paleomagnetism time variations secular and long term 1520 geomagnetism and paleomagnetism magnetostratigraphy 1035 geochemistry geochronology 1513 geomagnetism and paleomagnetism geomagnetic excursions keywords argonargon dating lavas excursion paleomagnetism reversal

experimental constraints on bubble interactions in rhyolite melts implications for vesicle size distributions we have studied interactions between bubbles of two distinct size classes in rhyolite melts experimentally decompressed between 200 and 80 mpa the first set of decompression bubbles has a size range rdec of 1 11 m and is formed from nucleation and growth upon isothermal decompression of the melt the larger populations of hydration bubbles are on average 30 40 m in radius rhyd and are formed from pore spaces present that were filled with water vapor during the saturation runs prior to the decompression experiments the first type of interaction results in the elongation of decompression bubbles oriented radially around the larger hydration bubbles the degree of elongation increases both as a function of distance and with increasing ratio of hydration to decompression bubble size rhyd rdec the second type of interaction studied results in a reduction of the size of decompression bubbles located within a range of distances from 10 to 65 m from a hydration bubble surface relative to the modal size of the unaffected bubbles in the same sample in addition within an average distance of 10 m melt next to the hydration bubble surface is depleted in decompression bubbles our results indicate that concentration gradients in the melt are probably responsible for bubble size reduction and the depleted zones because the predicted time scales for ostwald ripening are much longer than those of the experiments these effects persist even to the lowest ending pressures studied 80 mpa which indicates that size distributions of small bubbles may be affected by concentration gradients in the depleted melt shell surrounding large bubbles large bubbles present in an ascending magma prior to a subsequent nucleation event could therefore affect the growth of the smaller bubble population occurring within the depleted melt shell of the larger bubbles and produce a bimodal vesicle size distribution elongated decompression bubbles may be strained as a result of melt flowing away from the much larger hydration bubbles as they grow estimates of capillary number ca plotted against deformation df indicate that bubbles in water rich rhyolite melts are deformable even at small sizes 1 m and small values of ca our results show a different trend of df with ca than previous studies in non geological systems predict indicating that viscosity effects may be important the preservation of deformation textures depends strongly on relaxation time explaining the lack of deformation textures in less viscous natural lavas

three magmatic components in the 1973 eruption of eldfell volcano iceland evidence from plagioclase crystal size distribution csd and geochemistry abstract the 1973 eruption of eldfell volcano iceland appears to have been a short simple event but textural and geochemical evidence suggest that it may have had three different magmatic components the first erupted fissure magmas were chemically evolved rich in plagioclase 18 and had shallow straight crystal size distribution csd curves the early lavas were less evolved chemically had lower plagioclase contents 13 and steeper slightly concave up csds the late lavas were chemically similar to the early lavas but even richer in plagioclase than the initial magmas 24 and had the steepest csds there was no chemical evidence for plagioclase fractionation but compositional diversity could be produced by clinopyroxene fractionation which must have occurred at depth we propose that the eruption started with old coarsened ostwald ripened magma left over from a previous eruption possibly that which produced surtsey island ten years earlier the early flows may be mixtures of small amounts of this old magma with a new low crystallinity uncoarsened magma or a completely new magma the late flows are another new magma from depth chemically similar to the early flows but which has grown plagioclase under increasing saturation undercooling perhaps during its ascent all three magmatic components may have originated from the same parent but had varying degrees of clinopyroxene fractionation plagioclase nucleation and growth and coarsening

textural analysis of magmatic enclaves from the kameni islands santorini greece abstract the kameni islands have been the focus of historic volcanic activity on santorini greece the islands comprise a series of dacite lavas the erupted compositions of which have remained approximately constant for the last 2200 years the dacite lava flows of nea kameni contain a variety of mafic enclaves distinguished by differences in texture composition and mechanical behaviour the mafic enclaves represent chilled fragments of mafic magma that crystallised rapidly on contact with the cooler more silicic host magma the absence of chilled margins on the enclaves suggests that the replenishing magma initially formed a dense layer at the base of the chamber groundmass textures are interpreted as an indication of the degree of undercooling between the enclaves and the host groundmass textural information obtained from the enclaves has been used to determine the relative volumes of replenishing magmas injected into the host magma before eruption which indicates that the volume of erupted magma is directly proportional to the volume of replenishing magma emplaced in the chamber prior to eruption which thus may act as a potential eruption trigger the variety of replenishing magmas implies the existence of a complex conduit and chamber system beneath santorini volcano andesitic enclaves with linear crystal size distributions are thought to have originated as an aphyric andesitic melt expelled from a crystal mush by filter pressing disequilibrium phenocryst assemblages in the lavas and phenocryst bearing enclaves provide evidence for the entrainment and recycling of older phenocryst populations

a case history of paroxysmal explosion at stromboli timing and dynamics of the april 5 2003 event abstract on april 5 2003 stromboli volcano italy produced the most violent explosion of the past 50 years the event was exceptionally well documented thanks to the presence on the island of several scientists and a large number of instruments deployed over the preceding months to monitor the effusive eruption that began in december 2002 integration of visual documentation deposit features and geophysical data allowed an accurate reconstruction of the explosive event and its dynamics the eruption consisted of a 8 min long explosive event which evolved through four phases whose timing was precisely recorded by an infrared thermometer located about 450 m from the summit crater phases 2 and 3 lasted 39 and 42 s respectively both had an impulsive character were responsible for ejecting almost the entire mass of the pyroclastic products phases 1 and 4 represented respectively a short lived precursory event and a waning tale during phase 2 meter sized ballistic blocks were launched with velocities of 170 m s to altitudes of up to 1400 m above the craters these fell on the volcano flanks and on the village of ginostra about 2 km distant from the vent a vertical jet rose above the craters which developed to feed a convective plume that reached a height of up to 4 km the calculated mass of the phase 2 fallout deposit and mass discharge rate were 1 1 1 4 108 kg and 2 8 3 6 106 kg s respectively during phase 3 a scoria flow deposit with an estimated volume of 0 9 1 1 104 m3 was erupted from the same vent that fed the ongoing sustained lava flow the average mass discharge rate for this phase was 2 5 3 1 105 kg s products emitted during phases 2 and 3 consisted of lithic and fresh magmatic material in similar proportions the juvenile fraction consisted of a deep originated almost aphyric highly vesicular pumice mingled with a shallow derived crystal rich moderately vesicular scoria similarities with the eruption dynamics of other historical paroxysms at stromboli makes the april 5 2003 explosion representative of these highly energetic events that constitute the most hazardous volcanic phenomena at stromboli volcano

paleomagnetism of the southwestern u s a recorded by 0 5 ma igneous rocks 1 the issue of permanent nondipole contributions to the time averaged field lies at the very heart of paleomagnetism and the study of the ancient geomagnetic field in this paper we focus on paleomagnetic directional results from igneous rocks of the southwestern u s a in the age range 0 5 ma and investigate both the time averaged field and its variability about the mean value several decades of work in the southwestern united states have resulted in the publication of paleomagnetic data from over 800 individual paleomagnetic sites as part of a new investigation of the san francisco volcanics we collected paleomagnetic samples from 47 lava flows many of which have been previously dated the new data combined with published data are highly scattered contributions to the scatter were considered and we find that removal of data sets from tectonically active areas and judicious selection according to fisher s 1953 precision parameter results in an axially symmetric data distribution with normal and reverse modes that are indistinguishable from antipodal monte carlo simulations suggest that a minimum of 5 samples per site are needed to estimate the precision parameter sufficiently accurately to allow its use as a determinant of data quality numerical simulations from statistical paleosecular variation models indicate the need for several hundred paleomagnetic sites to get an accurate determination of the average field direction and are also used to investigate the directional bias that results from averaging unit vectors rather than using the full field vector average directions for the southwestern u s a show small deviations from a geocentric axial dipole field but these cannot be considered statistically significant virtual geomagnetic pole vgp dispersions are consistent with those from globally distributed observations analyzed by mcelhinny and mcfadden 1997 however a systematic investigation of the effect of imposing a cutoff on vgps with large deviations from the geographic axis indicates that while it may reduce bias in calculating the average direction such a procedure can result in severe underestimates of the variance in the geomagnetic field a more satisfactory solution would be to use an unbiased technique for joint estimation of the mean direction and variance of the field distribution

multiple sector collapses at stromboli volcano italy how they work this paper demonstrates that four large sector collapses have affected the nw flank of the stromboli volcano in the past 13 ka alternating with growth phases in order to contribute to the evaluation of the critical conditions which trigger lateral collapses a reconstruction of the geometry of each collapse of the volcano edifice in the four stages that preceded the relative collapse events is also presented and a computation of the landslide volume this reconstruction is based on new field data plotted in three dimensions prior to the initial 13 ka collapse the volcano was 1125 100 m high above sea level the collapse had a volume of 2 23 0 87 km3 whereas the pre collapse volcano volume was 218 8 7 7 km3 the next edifice that failed was 900 70 m high a s l the collapse volume was 1 0 54 km3 with a precollapse volcano volume of 201 4 5 4 km3 the edifice then grew to 1000 60 m a s l the third collapse had a volume of 1 08 0 39 km3 and occurred within a volcano with a volume of 209 1 4 6 km3 this was followed by a new growth phase followed by the last collapse with a volume of 0 73 0 22 km3 the volcano volume was about the same as the present one the present active crater zone is at 780 m a s l in the first three collapses sliding surfaces cut the main magma conduit in the last collapse the upper scarp coincided with the conduit location dyking along a main ne trending weakness zone across the volcano summit exerted a lateral force for collapse inception the decrease of the landslide volumes with the age and the concentric scarps of the four collapses suggest that the younger sliding planes tended to become more superficial and to decrease the areal extent this is interpreted as due to a successively weaker eruptive products from dominantly lavas to dominantly pyroclastics and b the feedback effects between collapses and dykes that injected along the lateral segments of the first collapse slide plane

effects of magma storage and ascent on the kinetics of crystal growth the size distributions of crystals of olivine plagioclase and oxides of the 1991 93 eruption at mt etna italy are analyzed the simultaneous collection of this information for different minerals gives precious insight into the cooling history of lavas three distinct episodes are detectable a storage of the magma in a deep reservoir characterized by nearly constant and low nucleation and growth rates near to equilibrium an ascent phase with an ever increasing nucleation rate related to volatile exsolution and finally a quenching phase in addition to geochemical and geophysical evidence the similarity of the crystal size distributions of the present eruption with those of previous ones of this century makes it possible to exclude that crystal size distributions of etnean lavas are due to mixing of different populations this strongly suggests that the main features of the volcano feeding system have not changed despite observed variations in the magma output rates

dissolution crystallization kinetics recorded in the 2002 2003 lavas of stromboli italy on 28 december 2002 new vents opened on the flanks of stromboli just below the summit craters interrupting the persistent activity of the volcano with a 7 month long effusive eruption we here report on the plagioclase size distribution plgsd in lava samples collected following the chronology of the 2002 2003 eruption data reveal a linear plgsd similar to that found in samples of normal stromboli activity indicating that the switch from strombolian explosive to effusive activity is not associated with changes in texture nevertheless the crystal size distribution slopes and intercepts exhibit slight sinusoidal temporal variations that are here ascribed to a magma supply mechanism able to induce resonance in the crystal size distribution with an amplitude that depends on the supply rate

experimental study on the effects of temperature cycling on coarsening of plagioclase and olivine in an alkali basalt we conducted experiments on an alkali basalt at 1 atmosphere in order to investigate the effects of temperature cycling on crystal coarsening experiments at 1 150 c near the ni nio buffer indicate that coarsening of plagioclase and olivine crystals is greatly accentuated by temperature cycling for a given experiment duration crystal number density decreases with temperature cycle amplitude and average crystal size increases with increasing amplitude we observed little correlation between cycle period and crystal number density or average crystal size we suggest that dissolution and size proportional crystal growth during repeated heating and cooling decrease crystal number density and increase average crystal size these experiments indicate that the texture of silicate minerals can be modified by temperature cycling and that phenocrysts may develop quicker in silicate magmas when the magma temperature is cycled

transient explosions at open vent volcanoes the case of stromboli italy stromboli is a persistently active open vent basaltic volcano whose activity is controlled by the balance between magma supply outgassing and eruptive rates and is characterized by low intensity regular strombolian explosions however two types of large transient violent explosive eruptions suddenly occur with no clear precursory activity these explosions called major and paroxysmal depending on size cover a large variability in intensity and mag nitude but are all marked by short duration paroxysms have significantly larger intensities 106 kg s than major explosions 104 kg s and fundamental differences in the characteris tics composition crystallinity vesicularity of the erupted tephra suggesting that different sources feed these two eruption types paroxysms are generated by the explosive fragmen tation of low porphyricity lp magma mingled with high porphyricity resident magma in the shallow reservoir whereas major eruptions are likely associated with destabilization of the lower portion of the shallow magmatic system continuously hybridized by the arrival of lp magma in general the intensity of these explosions is related to the amount of the lp magma erupted 107 kg in paroxysms and 104 105 kg in major explosions suggesting that the magma plays a major role in the fragmentation mechanism despite its primary impor tance in the hazards of stromboli the total amount of magma erupted in these events in the past 10 years is less than 1 of the total mass erupted by the volcano

experimental quantification of plagioclase crystal size distribution during cooling of a basaltic liquid an experimental study has been conducted to constrain how thermal history controls crystal size distribution csd of plagioclase in cooling basalts data from all experiments are dominated by a log linear segment of decreasing number density with increasing crystal size consistent with observations in many natural rocks the slope of the csd is found to be a function of cooling rate faster cooling leading to greater slopes at constant cooling rate the csds flatten as temperature decreases and are sometimes kinked characteristics consistent with independent textural observations that crystal agglomeration contributes significantly to crystal growth a downturn is observed toward small size which suggests ostwald ripening furthermore we find that thermal history above the liquidus has a major influence on csds and on the temperatures of phase appearance we conclude that near the liquidus heterogeneous nucleation typically dominates

coupled degassing and crystallization experimental study at continuous pressure drop with application to volcanic bombs experiments on degassing of water saturated granite melts with a pressure drop from 100 and 450 mpa to 40 and 120 mpa respectively at temperatures close to feldspar liquidus 750 700 c were carried out to determine the modality of water exsolution and vesicle formation at the liquidus temperature pressure drop rates as small as approximately 100 bar day were used uniform space distributions of bubbles of exsolved water were obtained with starting glass containing a small fraction 0 5 vol of trapped air bubbles volume crystallization of feldspar was observed in degassed melts supplied with seeds bubble size distributions bsd measured in granite glasses after degassing are presented data on vesicle characteristics number radius area elongation were acquired on images digitized with standard software while the reconstruction of size distributions was performed with the schwartz saltikov unfolding procedure bubble size distributions of size classes in the range 5 1000 m were acquired with proper magnification and satisfactory statistical reliability of determined number densities the bsds of the experimental samples are compared with the results of measurements of rapidly degassed products of mt etna and vulcano island many particular features of the bubble nucleation and growth can be distinguished in an individual bsd however the general bsd of the whole data set including natural ones can be relatively well described with linear regression in bilogarithmic coordinates the slope of this regression is approximately 2 8 0 1 this dependence is in striking contrast with distributions theoretically predicted with classical nucleation models based on homogeneous nucleation of vesicles the theoretical distribution requires the occurrence of strong maxima that are not observed in our experimental and natural samples thus arguing for heterogeneous nucleation mechanisms

the rheology of crystal bearing basaltic magmas from stromboli and etna abstract our ability to forecast volcanic hazards at active volcanoes stems from our knowledge of parameters affecting eruption dynamics persistent activity displayed by basaltic volcanic systems such as etna and stromboli is governed mainly by the textural evolution and rheological properties of eruptive products here we investigate the high temperature 1131 c

an aeolian grainflow model for martian recurring slope lineae abstract recurring slope lineae rsl on mars have been enigmatic since their discovery their behavior resembles a seeping liquid but sources of water remain puzzling this work demonstrates that the properties of rsl are consistent with observed behaviors of martian and terrestrial aeolian processes specifically rsl are well explained as flows of sand that remove a thin coating of dust observed rsl properties are supportive of or consistent with this model which requires no liquid water or other exotic processes but rather indicates seasonal aeolian behavior these settings and behaviors resemble features observed by rovers and also explain the occurrence of many slope lineae on mars that do not meet the strict definition of rsl this indicates that rsl can be explained simply as aeolian features other processes may add complexities just as they could modify the behavior of any sand dune

crystal size distribution csd in rocks and the kinetics and dynamics of crystallization ii makaopuhi lava lake crystal size distribution csd theory has been applied to drill core samples from makaopuhi lava lake kilauea volcano hawaii plagioclase and fe ti oxide size distribution spectra were measured and population densities n were calculated and analyzed using a steady state crystal population balance equation n n0 exp l g slopes on ln n versus crystal size l plots determine the parameter g a product of average crystal growth rate g and average crystal growth time the intercept is j g where j is nucleation rate known temperature depth distributions for the lava lake provide an estimate of effective growth time allowing nucleation and growth rates to be determined that are independent of any kinetic model plagioclase growth rates decrease with increasing crystallinity 9 9 5 4 10 11 cm s as do plagioclase nucleation rates 33 9 1 6 10 3 cm3 s ilmenite growth and nucleation rates also decrease with increasing crystallinity 4 9 3 4 10 10 cm s and 15 2 2 10 3 cm3 s respectively magnetite growth and nucleation rates are also estimated from the one sample collected below the magnetite liquidus g 2 9 10 10 cm s j 7 6 10 2 cm3 s moments of the population density function were used to examine the change in crystallization rates with time preliminary results suggest that total crystal volume increases approximately linearly with time after 50 crystallization a more complete set of samples is needed for material with 50 crystals to define the entire crystallization history comparisons of calculated crystallization rates with experimental data suggests that crystallization in the lava lake occurred at very small values of undercooling this interpretation is consistent with proposed thermal models of magmatic cooling where heat loss is balanced by latent heat production to maintain equilibrium cooling

effusive to explosive transition during the 2003 eruption of stromboli volcano the persistent explosive activity of stromboli volcano italy ceased in december 2002 and correlated with the onset of a seven month long effusive eruption on the volcano flank from new vents that opened just below the summit craters we intensively monitored this effusive event collecting and interpreting in real time an extensive multiparametric geophysical data set the resulting data synergy allowed detailed insights into the conduit dynamics that drove the eruption and the transition back to the typical strombolian activity we present a direct link between gas flux magma volume flux and seismicity supporting a gas driven model whereby the balance between gas flux and gas overpressure determines whether the system will support effusive or explosive activity this insight enabled us to monitor the migration of the magma column up the conduit and to explain the onset of explosive activity

textural effects of steady state behaviour of the stromboli feeding system abstract the exceptionally persistent activity of stromboli volcano has lasted for at least 1400 years the normal strombolian activity is periodically interrupted by more energetic explosions 1 2 per year and by sporadic effusive episodes every 10 20 years normal activity and effusive episodes are characterized by crystal rich high k to shoshonitic basalts issuing from a volatile poor shallow system crystal poor pumice are emitted only during more violent explosions and are thought to derive from deep pulses of volatile rich magma shallow level degassing induces massive crystallization of deep pulses of feeding magma that continuously mixing with the resident one produces the crystal rich shoshonite of the persistent activity we examined the crystallization history of the crystal rich shallow reservoir using plagioclase crystal size distribution csd analysis of scoriae and lavas emitted in the past twenty years csds show a linear dependence from crystal size in the size interval 0 06 1 2 mm number density of larger crystals is biased by right hand truncation effects csds slopes and intercepts are quite constant during the whole considered time span revealing a system that is close to the equilibrium also from a kinetic point of view the linear crystal size distribution are reached by the system through episodes of growth and resorption respectively occurring in the degassed and undegassed magma during the continuous mixing in the feeding system plagioclase net growth rate 2 10 11 cm s results from a balance of growth 10 10 cm s and resorption episodes which induce spectacular zoning and resorption textures in crystals larger than 200 m csds of mafic phases cannot be accurately acquired on each single sample due to poor counting statistics the evaluation of pyroxene and olivine csd on the whole data set however confirms the conclusions acquired from plagioclase csds

magma dynamics beneath kameni volcano thera greece as revealed by crystal size and shape measurements size distributions of plagioclase crystals in series of recent porphyritic dacite lavas from kameni volcano greece can be modelled by mixing two populations of crystals each with overlapping linear crystal size distributions csd termed microlites and megacrysts the magmas bearing the microlites and megacrysts started to crystallise 6 13 and 24 96 years respectively before each eruption the dates of initiation of crystallisation of the megacrysts indicate that they are left overs of earlier injections of new magma into a shallow chamber some magma remains after each eruption and continues to crystallise new magma with few or no crystals is then introduced and the microlites crystallise from the mixed magma eruption followed 6 13 years after mixing such a model would suggest that some porphyritic magmas are products of a shallow magma chamber that is never completely emptied just topped up from time to time

geomagnetic paleosecular variation recorded in plio pleistocene volcanic rocks from possession island crozet archipelago southern indian ocean possession island in the crozet archipelago consists of volcanic units erupted mainly between 5 and 0 5 ma a paleomagnetic sampling was carried out along several sections distributed near the northern eastern and southeastern coasts a total of 45 independent flows were sampled 320 samples for each flow a precisely defined characteristic remanence direction was usually isolated after a careful progressive cleaning in alternating fields however particularly complex remanence behavior is often observed the magnetostratigraphy of the lava pile is quite simple with reversed rocks in the lower part and normal units in the upper part of two sections a third section is of normal polarity throughout its whole thickness including three excursional directions we did not find any intermediate directions between the normal and reverse magnetozones thus we have no evidence for the recording of the matuyama brunhes transition expected from a previous study watkins 1972 the amplitude of paleosecular variation estimated from between flow dispersion from the field of an axial dipole is 11 8 with 95 confidence limits between 9 3 deg and 14 0 deg this value is consistent with the general anisotropic statistical model for paleosecular variation of constable and johnson 1999

quantifying the building blocks of igneous rocks are clustered crystal frameworks the foundation most phenocryst populations in volcanic rocks and those preserved in shallow level igneous intrusions are clustered variously referred to as clots clumps or glomerocrysts these clusters of crystals are the building blocks that accumulate to form the high porosity touching crystal frameworks from which igneous cumulates form examination of touching crystal frameworks in olivine komatiite cumulates and experimental charges and plagioclase dominant crystal populations holyoke flood basalt connecticut usa reveal complex high porosity clustered crystal arrangements olivine touching frameworks in komatiite flows are interpreted to form in hundreds of days plagioclase frameworks are calculated to have formed in less than 17 years for a crystal growth rate of 1 x 10 10 mm s to less than 3 years for a growth rate of 5 x 10 10 mm s based on crystal size distributions the origin of crystal clusters is likely to involve either or a combination of heterogeneous nucleation remobilization of cumulate mushes or crystals sticking together during settling and or flow the spatial distribution pattern of clustered crystal frameworks from both natural and experimental examples constrains fields on spatial packing diagrams that allow the identification of touching and non touching crystal populations and further improve our understanding of crystal packing arrangements and cluster size distributions

recycling and re hydration of degassed magma inducing transient dissolution crystallization events at stromboli italy abstract intrusive degassing and recycling of degassed and dense magma at depth have been proposed for a long time at stromboli the brief explosive event that occurred at the summit craters on 9 january 2005 threw out bombs and lapilli that could be good candidates to illustrate recycling of shallow degassed magma at depth we present an extensive data set on both the textures and the mineral bulk rock and glassy matrix chemistry of the 9 jan products the latter have the common shoshonitic basaltic bulk composition of lavas and scoriae issued from typical strombolian activity in contrast they differ by the heterogeneous chemistry of their matrix glasses and their crystal textures that testify to crystal dissolution event s just prior magma crystallization upon ascent and eruption comparison between mineral paragenesis of the natural products and experimental phase equilibria suggest water induced magma re equilibration we propose that mineral dissolution is related to water enrichment of the recycled degassed magma via differential gas bubble transfer and to some extents its physical mixing with volatile rich magma blobs however all these features illustrate transient processes even though evidence of mineral dissolution is ubiquitous at stromboli its effect on the bulk magma chemistry is minor because of the subtle interplay between mineral dissolution and crystallization in magmas having comparable bulk chemistry

kinetics of crystal evolution as a probe to magmatism at stromboli aeolian archipelago italy abstract the kinetics of hydrated basaltic melts erupted during the present activity at stromboli have been studied to estimate the growth and dissolution rate of plagioclase specifically a high k basalt composition pst 9 has been studied to investigate magma and eruption dynamics at stromboli volcano by combining crystallization kinetics of plagioclase and csd measurements on natural samples from literature armienti et al 2007 fornaciai et al 2009 a series of water saturated decompression experiments over a range of final water pressure pf 75 5 mpa at constant temperature 1075 c show that plagioclase is systematically present from 50 to 5 mpa at water saturated conditions moreover these experiments show that anorthite an content decreases with decreasing p h 2 o reaching the same composition as the natural plagioclase in stromboli scoria at pressure below 20 mpa and that the plagioclase crystal fraction increases as the experimental conditions tend to lower final pressure plagioclase growth rate gl is observed to increase with undercooling for the pf investigated during decompression experiments except for the 75 mpa pf serie that only has two samples with the presence of plagioclase crystals the values of gl vary from 10 7 to 10 8 cm s for pf from 75 to 25 mpa while at pf from 10 to 5 mpa growth rates are approximately of 10 6 cm s a series of dissolution experiments at atmospheric pressure and over a range of temperature has been done for plagioclase t range of 1220 1240 c dissolution rate g for plagioclase 10 7 cm s tends to be slightly higher at higher temperature in the range of 1220 1240 c and appears to be time independent for the experimental durations investigated 10 30 h these trends could be related to development of a diffusion limited boundary layer adjacent to the dissolving crystal by comparison of the experimental data on plagioclase composition growth rates and dissolution in stromboli basalt it is possible to place the reservoir of the crystal rich magma in the upper part from 400 m to the surface of the volcanic conduit kinetic data of the plagioclase the most important phase of the shallow magmatic system of stromboli show that the magmatic processes are quite dynamic and in a relatively short time hours or several days the system can change considerably furthermore the results from this work combined with observations on natural samples help to improve our knowledge of the magma plumbing system of interactions between resident magma and new magmas the dynamics of volcanic activity of stromboli and the time scales of magmatic processes that change in a few hours to 1 month

paleosecular variation of brunhes age lava flows from british columbia canada 1 brunhes age lava flows have been sampled for paleosecular variation studies from several volcanic fields of southern british columbia silverthrone garibaldi lake mt meager clearwater wells gray park and kelowna area a total of 52 lava flows were sampled and 7 to 10 samples were drilled at each site previous radiometric studies indicate that the ages of these lava flows range from 2 3 to 760 ka stepwise thermal demagnetization 14 to 21 steps was carried out for all the samples in each site and af demagnetization was performed on one sample per site forty five sites were selected based on rigorous criteria 95 5 to calculate a mean direction d 356 9 i 70 2 95 2 8 that is statistically indistinguishable from the direction of the dipole field at the area i 68 3 virtual geomagnetic poles vgps do not show the far sided effect and the angular standard deviation is 17 5 a value in agreement with the paleomagnetic field dispersion for that latitude these high quality results are expected to improve the time averaged field taf and secular variation models

steady state volcanism paleoeffusion rates and magma system volume inferred from plagioclase crystal size distributions in mafic lavas dome mountain nevada abstract the dome mountain dm lavas located in the moat of the timber mountain caldera of the sw nevada volcanic field cap a thick sequence of welded tuffs massive and bedded tuffs and mafic and rhyolitic lava flows these 10 9 m y old flows comprise a pile of some twenty flows with a total thickness of about 300 m individual flows range in thickness from 3 to 20 m and are characterized by a basal vesicular zone a massive central region and an upper vesicular region containing in most cases cooling fractures the flows dip gently away from the summit of dm a prominent topographic feature of the region with an elevation of 2065 m silica contents range from 48 8 wt at the base of the pile to 59 4 wt at the top although most of the flows are andesitic crystal size distributions csds for plagioclase in these comagmatic lavas have been determined to provide information on magma storage times these times when coupled with phase equilibria can be used to infer cooling rates which are indicative of the relative dimensions of magmatic bodies and the vigor of various magmatic cooling processes most of the csds show no sign of crystal fractionation or accumulation assuming a growth rate of 10 10 cm s calculated mean crystal residence times for most lavas cluster in the range of 1 5 to 4 years two of the lower lavas have residence times in the range of 4 to 9 years the residence times correlate inversely with nucleation density and broadly correlate inversely with silica content high residence time low silica lavas occur at the base of the pile whereas low residence time high silica lavas occur at the top of the pile one interpretation of the residence time data is that the lavas are from an open system of more or less constant residence time but varying spatially in system locations e g a magma chamber in which steady state recharge and eruption may have been reached alternatively these lavas may represent small samples closely spaced in time of a slowly cooling large batch or reservoir of magma the greater residence times associated with the lower flows may indicate the increased amount of time necessary to form the conduit of the magma volcano system by the oldest magma or simply that part of the system with the highest liquidus temperature a minimum volume of 1 4 to 2 5 km 3 for the magma chamber is estimated for the dm system using the estimated residence times in a simplified cooling model of a sheet of aspect ratio r the actual system volume depends on r which is strictly unknown and thus the actual system volume may be much larger e g 200 km 3 a mean paleoeffusion rate of 0 05 km 3 yr is obtained for the lavas by assuming that the residence time of each lava corresponds to the amount of time elapsed since the eruption of the underlying flow

crystal size distributions and scaling laws in the quantification of igneous textures this study demonstrates the value of the combined use of scaling of the chemical evolution of cooling magmas have investigated the detailed relationships between temanalysis and crystal size distribution csd measurements in the interpretation of igneous textures and outlines applications of these perature phase abundance and phase composition the reaction series bowen 1922 and the cipw norm cross methods to other fundamentally important problems theoretical calculations and measurements of natural samples are used to et al 1903 have been supplemented by modern thermodynamic models that predict the compositions and abundcharacterize the relationship between igneous texture and cooling history the total number nt and mean length l of crystals ances of minerals in equilibrium with melt of a given composition and temperature e g ariskin et al 1993 in a sample are correlated through a scaling relationship of the form n t l this proportionality depends on the mineralogy ghiorso sack 1995 although there are a few areas of uncertainty the chemistry of igneous petrogenesis is of the rock and so a modal normalization factor c is introduced the csd slope s and intercept ln n are uniquely related to fairly well understood but chemical models cannot address the central probeach other through the equation ln n 4ln s ln c where s l this overall relationship allows the texture of a rock lem of igneous crystallization which is the physical process of generating a set of crystals from magma two to be related to the local duration of cooling through an arbitrary crystal growth model quantification of the link between the texture magmas with the same composition will normally produce the same minerals in the same proportions however of a rock and its cooling history makes it possible to predict spatial variations in texture we show an example of this method using depending on their thermal histories these two magmas can solidify into rocks with drastically different physical the sudbury igneous complex canada additional applications include relating the textures of volcanic rocks to spatial and temporal appearances e g basalt and gabbro the set of minerals is the same but what distinguishes them is the abundance variations in the magma chamber and extracting kinetic parameters from suites of comagmatic rocks and size of the constituent crystals the fine grained rocks in the chilled margins of most basic intrusions cooled rapidly and the coarse grained rocks in the interiors of plutons cooled slowly this has

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