Paleomagnetic Study of the Portage Lake Volcanics exposed in the Quincy Mine

A detailed paleomagnetic and rock-magnetic investigation was conducted on thirty six basaltic flows of the ~1095 Ma Portage Lake Volcanics. The flows were sampled along the East Adit of the Quincy Mine (Hancock, MI). Thirty two flows yielded well-defined primary magnetization directions carried by magnetite. A secondary magnetization component carried by hematite was also found in twenty nine flows. After correction for serial correlation between the flows, nineteen independent mean directions were calculated. The corresponding paleomagnetic pole is located at 25.5 °N, 182.1 °W (A95 = 3.5°). The new pole overlaps with the pole from the ~1087 Ma Lake Shore Traps suggesting a standstill of the North American plate during that time period. The low angular dispersion of virtual geomagnetic poles (S = 7.9°) suggests that the flows were erupted within a short time period, or that the strength of geomagnetic secular variation was lower than that of the recent field.

STATIC AND DYNAMIC STRESS CHANGE AT 27 VOLCANOES OF THE CENTRAL AMERICAN VOLCANIC ARC AFTER THE MW. 7.6 COSTA RICA EARTHQUAKE OF 5 SEPTEMBER 2012

Large earthquakes may strongly influence the activity of volcanoes through static and dynamic processes. In this study, we quantify the static and dynamic stress change on 27 volcanoes in Central America, after the Mw 7.6 Costa Rica earthquake of 5 September 2012. Following this event, 8 volcanoes showed signs of activity. We calculated the static stress change due to the earthquake on hypothetical faults under these volcanoes with Coulomb 3.3. For the dynamic stress change, we computed synthetic seismograms to simulate the waveforms at these volcanoes. We then calculated the Peak Dynamic Stress (PDS) from the modeled peak ground velocities. The resulting values are from moderate to minor changes in stress (10-1-10-2 MPa) with the PDS values generally an order of magnitude larger than the static stress change. Although these values are small, they may be enough to trigger a response by the volcanoes, and are on the order of stress changes implicated in many other studies of volcano and earthquake triggering by large earthquakes. This study provides insight into the poorly-constrained mechanism for remote triggering.

COMPARISON OF NON-HEATING PALEOINTENSITY TECHNIQUES USING BASALTS FROM LEMPTÉGY VOLCANO, FRANCE AND SYNTHETIC MAGNETITE-BEARING SAMPLES

Data of the strength of Earth’s magnetic field (paleointensity) in the geological past are crucial for understanding the geodynamo. Conventional paleointensity determination methods require heating a sample to a high temperature in one or more steps. Consequently, many rocks are unsuitable for these methods due to a heating-induced experimental alteration. Alternative non-heating paleointensity methods are investigated to assess their effectiveness and reliability using both natural samples from Lemptégy Volcano, France, and synthetic samples. Paleointensity was measured from the natural and synthetic samples using the Pseudo-Thellier, ARM, REM, REMc, REM’, and Preisach methods. For the natural samples, only the Pseudo-Thellier method was able to produce a reasonable paleointensity estimate consistent with previous paleointensity data. The synthetic samples yielded more successful estimates using all the methods, with the Pseudo-Thellier and ARM methods producing the most accurate results. The Pseudo-Thellier method appears to be the best alternative to the heating-based paleointensity methods.

USING AUTO- AND CROSS-CORRELATIONS FROM SEISMIC NOISE TO MONITOR VELOCITY CHANGES AT VILLARRICA VOLCANO, CHILE

We used the Green's functions from auto-correlations and cross-correlations of seismic ambient noise to monitor temporal velocity changes in the subsurface at Villarrica volcano in the Southern Andes of Chile. Campaigns were conducted from March to October 2010 and February to April 2011 with 8 broadband and 6 short-period stations, respectively. We prepared the data by removing the instrument response, normalizing with a root-mean-square method, whitening the spectra, and filtering from 1 to 10 Hz. This frequency band was chosen based on the relatively high background noise level in that range. Hour-long auto- and cross-correlations were computed and the Green's functions stacked by day and total time. To track the temporal velocity changes we stretched a 24 hour moving window of correlation functions from 90% to 110% of the original and cross correlated them with the total stack. All of the stations' auto-correlations detected what is interpreted as an increase in velocity in 2010, with an average increase of 0.13%. Cross-correlations from station V01, near the summit, to the other stations show comparable changes that are also interpreted as increases in velocity. We attribute this change to the closing of cracks in the subsurface due either to seasonal snow loading or regional tectonics. In addition to the common increase in velocity across the stations, there are excursions in velocity on the same order lasting several days. Amplitude decreases as the station's distance from the vent increases suggesting these excursions may be attributed to changes within the volcanic edifice. In at least two occurrences the amplitudes at stations V06 and V07, the stations farthest from the vent, are smaller. Similar short temporal excursions were seen in the auto-correlations from 2011, however, there was little to no increase in the overall velocity.

A Surface Displacement Analysis for Volcan Pacaya from October 2001 through March 2013 by Means of 3-D Modeling of Precise Position GPS Data

Volcán Pacaya is one of three currently active volcanoes in Guatemala. Volcanic activity originates from the local tectonic subduction of the Cocos plate beneath the Caribbean plate along the Pacific Guatemalan coast. Pacaya is characterized by generally strombolian type activity with occasional larger vulcanian type eruptions approximately every ten years. One particularly large eruption occurred on May 27, 2010. Using GPS data collected for approximately 8 years before this eruption and data from an additional three years of collection afterwards, surface movement covering the period of the eruption can be measured and used as a tool to help understand activity at the volcano. Initial positions were obtained from raw data using the Automatic Precise Positioning Service provided by the NASA Jet Propulsion Laboratory. Forward modeling of observed 3-D displacements for three time periods (before, covering and after the May 2010 eruption) revealed that a plausible source for deformation is related to a vertical dike or planar surface trending NNW-SSE through the cone. For three distinct time periods the best fitting models describe deformation of the volcano: 0.45 right lateral movement and 0.55 m tensile opening along the dike mentioned above from October 2001 through January 2009 (pre-eruption); 0.55 m left lateral slip along the dike mentioned above for the period from January 2009 and January 2011 (covering the eruption); -0.025 m dip slip along the dike for the period from January 2011 through March 2013 (post-eruption). In all bestfit models the dike is oriented with a 75° westward dip. These data have respective RMS misfit values of 5.49 cm, 12.38 cm and 6.90 cm for each modeled period. During the time period that includes the eruption the volcano most likely experienced a combination of slip and inflation below the edifice which created a large scar at the surface down the northern flank of the volcano. All models that a dipping dike may be experiencing a combination of inflation and oblique slip below the edifice which augments the possibility of a westward collapse in the future.

DOING MORE WITH SHORT PERIOD DATA: DETERMINING MAGNITUDES FROM CLIPPED AND OVER-RUN SEISMIC DATA AT MOUNT ST. HELENS

How can we calculate earthquake magnitudes when the signal is clipped and over-run? When a volcano is very active, the seismic record may saturate (i.e., the full amplitude of the signal is not recorded) or be over-run (i.e., the end of one event is covered by the start of a new event). The duration, and sometimes the amplitude, of an earthquake signal are necessary for determining event magnitudes; thus, it may be impossible to calculate earthquake magnitudes when a volcano is very active. This problem is most likely to occur at volcanoes with limited networks of short period seismometers. This study outlines two methods for calculating earthquake magnitudes when events are clipped and over-run. The first method entails modeling the shape of earthquake codas as a power law function and extrapolating duration from the decay of the function. The second method draws relations between clipped duration (i.e., the length of time a signal is clipped) and the full duration. These methods allow for magnitudes to be determined within 0.2 to 0.4 units of magnitude. This error is within the range of analyst hand-picks and is within the acceptable limits of uncertainty when quickly quantifying volcanic energy release during volcanic crises. Most importantly, these estimates can be made when data are clipped or over-run. These methods were developed with data from the initial stages of the 2004-2008 eruption at Mount St. Helens. Mount St. Helens is a well-studied volcano with many instruments placed at varying distances from the vent. This fact makes the 2004-2008 eruption a good place to calibrate and refine methodologies that can be applied to volcanoes with limited networks.

Comparison of magnetic properties and petrography between dykes and lava flows from La Cienega, New Mexico and Thunder Bay area, Canada

Data on the evolution of geomagnetic paleointensity are crucial for understanding the geodynamo and Earth’s thermal history. Although basaltic flows are preferred for paleointensity experiments, quickly cooled mafic dykes have also been used. However, the paleointensity values obtained from the dykes are systematically lower than those from lava flows. This bias may originate from the difference in cooling histories and resultant magnetic mineralogies of extrusive and intrusive rocks. To explore this hypothesis, the magnetic mineralogy of two feeder dyke-lave flow systems, from Thunder Bay (Canada) and La Cienega (New-Mexico), has been studied using magnetic and microscopy methods. Within each system, the flow and dyke show different stages of deuteric oxidation of titanomagnetite, but the oxidation stages also differ between the two systems. It is concluded that the tested hypothesis is viable, but the relationships between the magnetic and mineralogical properties of flows and dykes are complex and need a further investigation.

Sonic log prediction in carbonates

This work is conducted to study the complications associated with the sonic log prediction in carbonate logs and to investigate the possible solutions to accurately predict the sonic logs in Traverse Limestone. Well logs from fifty different wells were analyzed to define the mineralogy of the Traverse Limestone by using conventional 4-mineral and 3-mineral identification approaches. We modified the conventional 3-mineral identification approach (that completely neglects the gamma ray response) to correct the shale effects on the basis of gamma ray log before employing the 3-mineral identification. This modification helped to get the meaningful insight of the data when a plot was made between DGA (dry grain density) and UMA (Photoelectric Volumetric Cross-section) with the characteristic ternary diagram of the quartz, calcite and dolomite. The results were then compared with the 4-mineral identification approach. Contour maps of the average mineral fractions present in the Traverse Limestone were prepared to see the basin wide mineralogy of Traverse Limestone. In the second part, sonic response of Traverse Limestone was predicted in fifty randomly distributed wells. We used the modified time average equation that accounts for the shale effects on the basis of gamma ray log, and used it to predict the sonic behavior from density porosity and average porosity. To account for the secondary porosity of dolomite, we subtracted the dolomitic fraction of clean porosity from the total porosity. The pseudo-sonic logs were then compared with the measured sonic logs on the root mean square (RMS) basis. Addition of dolomite correction in modified time average equation improved the results of sonic prediction from neutron porosity and average porosity. The results demonstrated that sonic logs could be predicted in carbonate rocks with a root mean square error of about 4μsec/ft. We also attempted the use of individual mineral components for sonic log prediction but the ambiguities in mineral fractions and in the sonic properties of the minerals limited the accuracy of the results.

A PETROPHYSICAL EVALUATION OF THE TRENTON-BLACK RIVER FORMATION OF THE MICHIGAN BASIN

This work is conducted to study the geological and petrophysical features of the Trenton- Black River limestone formation. Log curves, crossplots and mineral identification methods using well-log data are used to determine the components and analyze changes in lithology. Thirty-five wells from the Michigan Basin are used to define the mineralogy of Trenton-Black River limestone. Using the different responses of a few log curves, especially gamma-ray, resistivity and neutron porosity, the formation tops for the Utica shale, the Trenton limestone, the Black River limestone and the Prairie du Chien sandstone are identified to confirm earlier authors’ work and provide a basis for my further work. From these, an isopach map showing the thickness of Trenton-Black River formation is created, indicating that its maximum thickness lies in the eastern basin and decreases gradually to the west. In order to obtain more detailed lithological information about the limestone formations at the thirty-five wells, (a) neutron-density and neutron-sonic crossplots, (b) mineral identification methods, including the M-N plot, MID plot, ϱmaa vs. Umaa MID plot, and the PEF plot, and (c) a modified mineral identification technique are applied to these wells. From this, compositions of the Trenton-Black River formation can be divided into three different rock types: pure limestone, partially dolomitized limestone, and shaly limestone. Maps showing the fraction of dolomite and shale indicate their geographic distribution, with dolomite present more in the western and southwestern basin, and shale more common in the north-central basin. Mineral identification is an independent check on the distribution found from other authors, who found similar distributions based on core descriptions. The Thomas Stieber method of analysis is best suited to sand-shale sequences, interpreting hree different distributions of shale within sand, including dispersed, laminated and structural. Since this method is commonly applied in clastic rocks, my work using the Thomas Stieber method is new, as an attempt to apply this technique, developed for clastics, to carbonate rocks. Based on the original assumption and equations with a corresponding change to the Trenton-Black River formation, feasibility of using the Thomas Stieber method in carbonates is tested. A graphical display of gamma-ray versus density porosity, using the properties of clean carbonate and pure shale, suggests the presence of laminated shale in fourteen wells in this study. Combined with Wilson’s study (2001), it is safe to conclude that when shale occurs in the Trenton-Black River formation, it tends to be laminated shale.

GROUND-BASED AND SATELLITE REMOTE SENSING OF PAROXYSMAL ERUPTIONS AT ETNA VOLCANO, 2011-2012.

Mt Etna's activity has increased during the last decade with a tendency towards more explosive eruptions that produce paroxysmal lava fountains. From January 2011 to April 2012, 25 lava fountaining episodes took place at Etna's New South-East Crater (NSEC). Improved understanding of the mechanism driving these explosive basaltic eruptions is needed to reduce volcanic hazards. This type of activity produces high sulfur dioxide (SO2) emissions, associated with lava flows and ash fall-out, but to date the SO2 emissions associated with Etna's lava fountains have been poorly constrained. The Ultraviolet (UV) Ozone Monitoring Instrument (OMI) on NASA's Aura satellite and the Atmospheric Infrared Sounder (AIRS) on Aqua were used to measure the SO2 loadings. Ground-based data from the Observatoire de Physique du Globe de Clermont-Ferrand (OPGC) L-band Doppler radar, VOLDORAD 2B, used in collaboration with the Italian National Institute of Geophysics and Volcanology in Catania (INGV-CT), also detected the associated ash plumes, giving precise timing and duration for the lava fountains. This study resulted in the first detailed analysis of the OMI and AIRS SO2 data for Etna's lava fountains during the 2011-2012 eruptive cycle. The HYSPLIT trajectory model is used to constrain the altitude of the observed SO2 clouds, and results show that the SO2 emission usually coincided with the lava fountain peak intensity as detected by VOLDORAD. The UV OMI and IR AIRS SO2 retrievals permit quantification of the SO2 loss rate in the volcanic SO2 clouds, many of which were tracked for several days after emission. A first attempt to quantitatively validate AIRS SO2 retrievals with OMI data revealed a good correlation for high altitude SO2 clouds. Using estimates of the emitted SO2 at the time each paroxysm, we observe a correlation with the inter-paroxysm repose time. We therefore suggest that our data set supports the collapsing foam (CF) model [1] as driving mechanism for the paroxysmal events at the NSEC. Using VOLDORAD-based estimates of the erupted magma mass, we observe a large excess of SO2 in the eruption clouds. Satellite measurements indicate that SO2 emissions from Etnean lava fountains can reach the lower stratosphere and hence could pose a hazard to aviation. [1] Parfitt E.A (2004). A discussion of the mechanisms of explosive basaltic eruptions. J. Volcanol. Geotherm. Res. 134, 77-107.