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

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

Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Geophysical and geological analysis of the Collingwood Member of the Trenton Formation

The Collingwood Member is a mid to late Ordovician self-sourced reservoir deposited across the northern Michigan Basin and parts of Ontario, Canada. Although it had been previously studied in Canada, there has been relatively little data available from the Michigan subsurface. Recent commercial interest in the Collingwood has resulted in the drilling and production of several wells in the state of Michigan. An analysis of core samples, measured laboratory data, and petrophysical logs has yielded both a quantitative and qualitative understanding of the formation in the Michigan Basin. The Collingwood is a low permeability and low porosity carbonate package that is very high in organic content. It is composed primarily of a uniformly fine grained carbonate matrix with lesser amounts of kerogen, silica, and clays. The kerogen content of the Collingwood is finely dispersed in the clay and carbonate mineral phases. Geochemical and production data show that both oil and gas phases are present based on regional thermal maturity. The deposit is richest in the north-central part of the basin with thickest deposition and highest organic content. The Collingwood is a fairly thin deposit and vertical fractures may very easily extend into the surrounding formations. Completion and treatment techniques should be designed around these parameters to enhance production.

Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

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.

Seismic observations of regional gas blowdown at Teal South

A re-examination of seismic time-lapse data from the Teal South field provides support for a previously proposed model of regional pressure decline and the associated liberation of gas from nearby reservoirs due to the production from the only reservoir among them that is under production. The use of a specific attribute, instantaneous amplitude, and a series of time slices, however, provides increased detail in understanding fluid migration into or out of the reservoirs, and the path taken by pressure changes across faults. The regional decrease of pressure due to production in one reservoir has dramatic effects in nearby untapped reservoirs, one of which appears to exhibit evidence for the escape, and possible re-trapping nearby, of hydrocarbons from a spill point. The influx of water into the producing reservoir is also evidenced by a decrease in amplitude at one end of the oil-water contact.

Morphologic changes in patients with geographic atrophy assessed with a novel spectral OCT-SLO combination

To investigate the appearance of geographic atrophy in high-resolution optical coherence tomography (OCT) images, the fundus autofluorescence (FAF) pattern, and infrared images simultaneously recorded with a novel combined OCT-scanning laser ophthalmology (SLO) system.

DIFFERENTIATING PATH AND SOURCE PROCESSES IN COMPLEX GEOLOGIC STRUCTURE THROUGH PASSIVE SEISMIC STUDIES: BERING GLACIER, ALASKA AND VILLARRICA VOLCANO, CHILE

Measuring shallow seismic sources provides a way to reveal processes that cannot be directly observed, but the correct interpretation and value of these signals depend on the ability to distinguish source from propagation effects. Furthermore, seismic signals produced by a resonating source can look almost identical to those produced by impulsive sources, but modified along the path. Distinguishing these two phenomena can be accomplished by examining the wavefield with small aperture arrays or by recording seismicity near to the source when possible. We examine source and path effects in two different environments: Bering Glacier, Alaska and Villarrica Volcano, Chile. Using three 3-element seismic arrays near the terminus of the Bering Glacier, we have identified and located both terminus calving and iceberg breakup events. We show that automated array analysis provided a robust way to locate icequake events using P waves. This analysis also showed that arrivals within the long-period codas were incoherent within the small aperture arrays, demonstrating that these codas previously attributed to crack resonance were in fact a result of a complicated path rather than a source effect. At Villarrica Volcano, seismometers deployed from near the vent to ~10 km revealed that a several cycle long-period source signal recorded at the vent appeared elongated in the far-field. We used data collected from the stations nearest to the vent to invert for the repetitive seismic source, and found it corresponded to a shallow force within the lava lake oriented N75°E and dipping 7° from horizontal. We also used this repetitive signal to search the data for additional seismic and infrasonic properties which included calculating seismic-acoustic delay times, volcano acoustic-seismic ratios and energies, event frequency, and real-time seismic amplitude measurements. These calculations revealed lava lake level and activity fluctuations consistent with lava lake level changes inferred from the persistent infrasonic tremor.

Differential optical densities of intraretinal spaces

PURPOSE: To test the hypothesis that hyporeflective spaces in the neuroretina found on optical coherence tomography (OCT) examination have different optical reflectivities according to whether they are associated with exudation or degeneration. METHODS: Retrospective analysis of eyes with idiopathic perifoveal telangiectasia (IPT), diabetic macular edema (DME), idiopathic central serous chorioretinopathy (CSC), retinitis pigmentosa (RP), or cone dystrophy (CD) and eyes of healthy control subjects. OCT scans were performed. Raw scan data were exported and used to calculate light reflectivity profiles. Reflectivity data were acquired by projecting three rectangular boxes, each 50 pixels long and 5 pixels wide, into the intraretinal cystoid spaces, centrally onto unaffected peripheral RPE, and onto the prefoveolar vitreous. Light reflectivity in the retinal pigment epithelium (RPE), vitreous, and intraretinal spaces for the different retinal conditions and control subjects were compared. RESULTS: Reflectivities of the vitreous and the RPE were similar among the groups. Hyporeflective spaces in eyes with exudation (DME, RP, and CSC) had higher reflectivity compared with the mean reflectivity of the vitreous, whereas the cystoid spaces in the maculae of the eyes without exudation (CD and IPT) had a lower reflectivity than did the normal vitreous. CONCLUSIONS: Analysis of the light reflectivity profiles may be a tool to determine whether the density of hyporeflective spaces in the macula is greater or less than that of the vitreous, and may be a way to differentiate degenerative from exudative macular disease.

Impact of optic media opacities and image compression on quantitative analysis of optical coherence tomography

To analyze the impact of opacities in the optical pathway and image compression of 32-bit raw data to 8-bit jpg images on quantified optical coherence tomography (OCT) image analysis.

Macular thickness measurements in healthy eyes using six different optical coherence tomography instruments

To compare central retinal thickness (CRT) measurements in healthy eyes by different commercially available OCT instruments and to compare the intersession reproducibility of such measurements.

A new method to monitor visual field defects caused by photoreceptor degeneration by quantitative optical coherence tomography

PURPOSE: To correlate the dimension of the visual field (VF) tested by Goldman kinetic perimetry with the extent of visibility of the highly reflective layer between inner and outer segments of photoreceptors (IOS) seen in optical coherence tomography (OCT) images in patients with retinitis pigmentosa (RP). METHODS: In a retrospectively designed cross-sectional study, 18 eyes of 18 patients with RP were examined with OCT and Goldmann perimetry using test target I4e and compared with 18 eyes of 18 control subjects. A-scans of raw scan data of Stratus OCT images (Carl Zeiss Meditec, AG, Oberkochen, Germany) were quantitatively analyzed for the presence of the signal generated by the highly reflective layer between the IOS in OCT images. Starting in the fovea, the distance to which this signal was detectable was measured. Visual fields were analyzed by measuring the distance from the center point to isopter I4e. OCT and visual field data were analyzed in a clockwise fashion every 30 degrees , and corresponding measures were correlated. RESULTS: In corresponding alignments, the distance from the center point to isopter I4e and the distance to which the highly reflective signal from the IOS can be detected correlate significantly (r = 0.75, P < 0.0001). The greater the distance in VF, the greater the distance measured in OCT. CONCLUSIONS: The authors hypothesize that the retinal structure from which the highly reflective layer between the IOS emanates is of critical importance for visual and photoreceptor function. Further research is warranted to determine whether this may be useful as an objective marker of progression of retinal degeneration in patients with RP.

Quantitative OCT analysis of idiopathic perifoveal telangiectasia

PURPOSE: To identify and quantitate specific changes in optical coherence tomography (OCT) images of patients with type 2 idiopathic perifoveal telangiectasia (IPT). METHODS: In a prospectively designed, observational, case-control study, 28 eyes of 14 consecutive patients with IPT were examined with OCT and compared with eyes of 14 unaffected control subjects. Light reflectivity profiles of raw scan data of OCT images were quantitatively analyzed for differences in distance between different retinal reflectivity layers and their respective reflectivities. Maculae were examined in four separate regions: (1) central fovea, (2) nasal perifovea, (3) temporal perifovea, and (4) outside the fovea. RESULTS: Retinal thinning, shortening of the photoreceptor outer segments and loss of reflectivity of the photoreceptor ellipsoid region were found in the central foveal region as well as the nasal and temporal perifoveal regions in eyes with IPT. In addition, increased reflectivity of the outer nuclear layer was found in a sharply demarcated area of the inferotemporal perifoveal region in all affected eyes. Retinal tissue located more than 2000 mum away from the foveola was indistinguishable from that in normal eyes. CONCLUSIONS: Quantitative OCT analysis shows unique and specific changes in the photoreceptors of the central macula in IPT which can be detected from first clinical presentation. These changes may be of use as an additional diagnostic tool. Correlation of the findings in the outer nuclear layer with histologic studies may help identify the nature of the reflectivity increase and define more clearly the type of damage sustained by the photoreceptors in this condition.

Decreased visual function after patchy loss of retinal pigment epithelium induced by low-dose sodium iodate

PURPOSE: To correlate damage to the retinal pigment epithelium (RPE) with decreased visual function after the systemic administration of sodium iodate (NaIO(3)). METHODS: Damage was produced in mice by injection of 15, 25, or 35 mg/kg NaIO(3). Visual function was assessed with the cued water maze (WM) behavioral test and the optokinetic reflex (OKR) measurement at different times after injection. Autofluorescence in whole eye flatmounts was quantified, and hematoxylin and eosin staining of paraffin sections was performed to assess changes in the outer retina. RESULTS: After 15 mg/kg NaIO(3), cued WM test results were normal, whereas OKR measurements were significantly decreased at all times. Focal RPE loss began on day 21, but no significant damage to the outer nuclear layer was observed. After 25 mg/kg NaIO(3), the cued WM test was transitionally reduced and the OKR measurement again decreased at all times. Large areas of RPE loss occurred on day 14 with a reduced outer nuclear layer on the same day. With 35 mg/kg NaIO(3), the cued WM test was reduced beginning on day 14 with complete obliteration of the OKR beginning on day 3, large areas of RPE loss on the same day, and a reduced outer nuclear layer on day 7. CONCLUSIONS: Stable, patchy RPE loss was observed with a low concentration of NaIO(3). The OKR measurement showed changes in visual function earlier than the cued WM test and before histologic findings were observed.