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.

TECHNOLOGICAL ADVANCEMENT IN REMOTE SENSING AND INNOVATIVE APPLICATION TO VOLCANOLOGY

The integration of remote monitoring techniques at different scales is of crucial importance for monitoring of volcanoes and assessment of the associated hazard. In this optic, technological advancement and collaboration between research groups also play a key role. Vhub is a community cyberinfrastructure platform designed for collaboration in volcanology research. Within the Vhub framework, this dissertation focuses on two research themes, both representing novel applications of remotely sensed data in volcanology: advancement in the acquisition of topographic data via active techniques and application of passive multi-spectral satellite data to monitoring of vegetated volcanoes. Measuring surface deformation is a critical issue in analogue modelling of Earth science phenomena. I present a novel application of the Microsoft Kinect sensor to measurement of vertical and horizontal displacements in analogue models. Specifically, I quantified vertical displacement in a scaled analogue model of Nisyros volcano, Greece, simulating magmatic deflation and inflation and related surface deformation, and included the horizontal component to reconstruct 3D models of pit crater formation. The detection of active faults around volcanoes is of importance for seismic and volcanic hazard assessment, but not a simple task to be achieved using analogue models. I present new evidence of neotectonic deformation along a north-south trending fault from the Mt Shasta debris avalanche deposit (DAD), northern California. The fault was identified on an airborne LiDAR campaign of part of the region interested by the DAD and then confirmed in the field. High resolution LiDAR can be utilized also for geomorphological assessment of DADs, and I describe a size-distance analysis to document geomorphological aspects of hummock in the Shasta DAD. Relating the remote observations of volcanic passive degassing to conditions and impacts on the ground provides an increased understanding of volcanic degassing and how satellite-based monitoring can be used to inform hazard management strategies in nearreal time. Combining a variety of satellite-based spectral time series I aim to perform the first space-based assessment of the impacts of sulfur dioxide emissions from Turrialba volcano, Costa Rica, on vegetation in the surrounding environment, and establish whether vegetation indices could be used more broadly to detect volcanic unrest.

EXPLORATION OF THE MTSAT2 SATELLITE CAPABILITIES FOR REAL TIME DETECTION AND CHARACTERIZATION OF VOLCANIC EMISSIONS

In this report, we attempt to define the capabilities of the infrared satellite remote sensor, Multifunctional Transport Satellite-2 (MTSAT-2) (i.e. a geosynchronous instrument), in characterizing volcanic eruptive behavior in the highly active region of Indonesia. Sulfur dioxide data from NASA's Ozone Monitoring Instrument (OMI) (i.e. a polar orbiting instrument) are presented here for validation of the processes interpreted using the thermal infrared datasets. Data provided from two case studies are analyzed specifically for eruptive products producing large thermal anomalies (i.e. lava flows, lava domes, etc.), volcanic ash and SO2 clouds; three distinctly characteristic and abundant volcanic emissions. Two primary methods used for detection of heat signatures are used and compared in this report including, single-channel thermal radiance (4-µm) and the normalized thermal index (NTI) algorithm. For automated purposes, fixed thresholds must be determined for these methods. A base minimum detection limit (MDL) for single-channel thermal radiance of 2.30E+05 Wm- 2sr-1m-1 and -0.925 for NTI generate false alarm rates of 35.78% and 34.16%, respectively. A spatial comparison method, developed here specifically for use in Indonesia and used as a second parameter for detection, is implemented to address the high false alarm rate. For the single-channel thermal radiance method, the utilization of the spatial comparison method eliminated 100% of the false alarms while maintaining every true anomaly. The NTI algorithm showed similar results with only 2 false alarms remaining. No definitive difference is observed between the two thermal detection methods for automated use; however, the single-channel thermal radiance method coupled with the SO2 mass abundance data can be used to interpret volcanic processes including the identification of lava dome activity at Sinabung as well as the mechanism for the dome emplacement (i.e. endogenous or exogenous). Only one technique, the brightness temperature difference (BTD) method, is used for the detection of ash. Trends of ash area, water/ice area, and their respective concentrations yield interpretations of increased ice formation, aggregation, and sedimentation processes that only a high-temporal resolution instrument like the MTSAT-2 can analyze. A conceptual model of a secondary zone of aggregation occurring in the migrating Kelut ash cloud, which decreases the distal fine-ash component and hazards to flight paths, is presented in this report. Unfortunately, SO2 data was unable to definitively reinforce the concept of a secondary zone of aggregation due to the lack of a sufficient temporal resolution. However, a detailed study of the Kelut SO2 cloud is used to determine that there was no climatic impacts generated from this eruption due to the atmospheric residence times and e-folding rate of ~14 days for the SO2. This report applies the complementary assets offered by utilizing a high-temporal and a high-spatial resolution satellite, and it demonstrates that these two instruments can provide unparalleled observations of dynamic volcanic processes.

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.

Volcanic forcing for climate modeling: a new microphysics-based data set covering years 1600–present

As the understanding and representation of the impacts of volcanic eruptions on climate have improved in the last decades, uncertainties in the stratospheric aerosol forcing from large eruptions are now linked not only to visible optical depth estimates on a global scale but also to details on the size, latitude and altitude distributions of the stratospheric aerosols. Based on our understanding of these uncertainties, we propose a new model-based approach to generating a volcanic forcing for general circulation model (GCM) and chemistry–climate model (CCM) simulations. This new volcanic forcing, covering the 1600–present period, uses an aerosol microphysical model to provide a realistic, physically consistent treatment of the stratospheric sulfate aerosols. Twenty-six eruptions were modeled individually using the latest available ice cores aerosol mass estimates and historical data on the latitude and date of eruptions. The evolution of aerosol spatial and size distribution after the sulfur dioxide discharge are hence characterized for each volcanic eruption. Large variations are seen in hemispheric partitioning and size distributions in relation to location/date of eruptions and injected SO2 masses. Results for recent eruptions show reasonable agreement with observations. By providing these new estimates of spatial distributions of shortwave and long-wave radiative perturbations, this volcanic forcing may help to better constrain the climate model responses to volcanic eruptions in the 1600–present period. The final data set consists of 3-D values (with constant longitude) of spectrally resolved extinction coefficients, single scattering albedos and asymmetry factors calculated for different wavelength bands upon request. Surface area densities for heterogeneous chemistry are also provided.

Effects of air pollution on asthma hospitalization rates in different age groups in metropolitan cities of Korea

Many studies have shown relationships between air pollution and the rate of hospital admissions for asthma. A few studies have controlled for age-specific effects by adding separate smoothing functions for each age group. However, it has not yet been reported whether air pollution effects are significantly different for different age groups. This lack of information is the motivation for this study, which tests the hypothesis that air pollution effects on asthmatic hospital admissions are significantly different by age groups. Each air pollutant's effect on asthmatic hospital admissions by age groups was estimated separately. In this study, daily time-series data for hospital admission rates from seven cities in Korea from June 1999 through 2003 were analyzed. The outcome variable, daily hospital admission rates for asthma, was related to five air pollutants which were used as the independent variables, namely particulate matter <10 micrometers (μm) in aerodynamic diameter (PM10), carbon monoxide (CO), ozone (O3), nitrogen dioxide (NO2), and sulfur dioxide (SO2). Meteorological variables were considered as confounders. Admission data were divided into three age groups: children (<15 years of age), adults (ages 15-64), and elderly (≥ 65 years of age). The adult age group was considered to be the reference group for each city. In order to estimate age-specific air pollution effects, the analysis was separated into two stages. In the first stage, Generalized Additive Models (GAMs) with cubic spline for smoothing were applied to estimate the age-city-specific air pollution effects on asthmatic hospital admission rates by city and age group. In the second stage, the Bayesian Hierarchical Model with non-informative prior which has large variance was used to combine city-specific effects by age groups. The hypothesis test showed that the effects of PM10, CO and NO2 were significantly different by age groups. Assuming that the air pollution effect for adults is zero as a reference, age-specific air pollution effects were: -0.00154 (95% confidence interval(CI)= (-0.0030,-0.0001)) for children and 0.00126 (95% CI = (0.0006, 0.0019)) for the elderly for PM 10; -0.0195 (95% CI = (-0.0386,-0.0004)) for children for CO; and 0.00494 (95% CI = (0.0028, 0.0071)) for the elderly for NO2. Relative rates (RRs) were 1.008 (95% CI = (1.000-1.017)) in adults and 1.021 (95% CI = (1.012-1.030)) in the elderly for every 10 μg/m3 increase of PM10 , 1.019 (95% CI = (1.005-1.033)) in adults and 1.022 (95% CI = (1.012-1.033)) in the elderly for every 0.1 part per million (ppm) increase of CO; 1.006 (95%CI = (1.002-1.009)) and 1.019 (95%CI = (1.007-1.032)) in the elderly for every 1 part per billion (ppb) increase of NO2 and SO2, respectively. Asthma hospital admissions were significantly increased for PM10 and CO in adults, and for PM10, CO, NO2 and SO2 in the elderly.^

Actividad antioxidante de los compuestos fenólicos contenidos en las semillas de la vid (vitis vinifera l.)

El propósito de este trabajo es obtener un antioxidante natural a partir de las semillas de uva (Vitis vinifera L.), para emplear en alimentos. Para ello se compararon distintos solventes para la extracción de fenoles de las semillas de la uva, de modo de obtener el extracto más concentrado en compuestos activos con la mínima degradación de su poder antioxidante durante el proceso de obtención. La concentración de fenoles totales de los extractos se determinó por el método Folin Ciocalteu. El poder reductor de los extractos se midió empleando el método de Oyaizu. Una vez seleccionado el solvente más adecuado para la extracción, se analizó la cinética de extracción, optimizando el tiempo de tratamiento. El extracto fue concentrado al vacío, y se veríficó la conservación del poder reductor en el extracto concentrado, por el método de Oyaizu. El extracto de semillas concentrado y sin concentrar se empleó en un sistema real sujeto a oxidación, tal como el jugo de manzanas. El grado de oxidación del jugo se midió por el método de Özoglu. El extracto concentrado fue deshidratado por secado en lecho de espuma y por liofilizado. En ambos casos se verificó el efecto del tratamiento de secado sobre el poder reductor. Finalmente, se evaluó la actividad antioxidante del extracto líquido concentrado de semillas de vid, respecto de otros antioxidantes comerciales como ácido ascórbico y dióxido de azufre. El sustrato oxidable fue el jugo de manzanas, y el grado de oxidación se midió por el método de Özoglu. El análisis estadístico de los datos se realizó mediante el análisis de la varianza; cuando no fue posible emplear el mencionado análisis, debido a que no se verificaban los supuestos básicos para su aplicación, se empleó la prueba de Kruskal –Wallis. En todos los casos, se utilizó el programa Statgraphics plus ®4.0. Para obtener un extracto antioxidante a partir de semillas de vid se utilizó una ex-tracción con agua a 90ºC, durante 4 horas. La relación sólido- líquido empleada fue de 1g de semillas enteras por 10 ml de solvente. El extracto obtenido presentaba una concentración de 12,587 mg de fenoles totales por gramo de semillas de uva extractadas y un poder reductor de 1,290 unidades. Como consecuencia del análisis de la cinética de extracción, el tiempo de tratamiento se redujo de 4 horas a 3 horas. La concentración del extracto se realizó al vacío a 60ºC, verificándose un aumen-to del poder reductor en el extracto concentrado, comprobado sobre jugo de manzanas. Comparando el extracto concentrado y el extracto sin concentrar se observa que la concentración de fenoles totales aumentó 29,57 veces, mientras que el poder reductor aumentó 37,39 veces. El deshidratado del extracto por medio del lecho de espuma permitió conservar el poder reductor del mismo, no ocurrió lo mismo en el deshidratado por liofilizado, donde se produjo un deterioro del poder reductor. Para un mismo contenido de fenoles totales agregado al jugo de manzanas, el ex-tracto líquido sin concentrar produjo un 28,4% de inhibición de la oxidación, mientras que el de extracto líquido concentrado produjo un 51,5 % de inhibición de la oxidación del jugo de manzanas. El extracto de semillas de vid, aplicado como antioxidante en jugo de manzanas, inhibió el desarrollo de la oxidación en un 31,51%, considerando 24 horas el tiempo de tratamiento. Este desempeño supera al ácido ascórbico, que en iguales condiciones, inhibió el desarrollo de la oxidación en un 2,6%. Pero en las condiciones de tra-bajo, el dióxido de azufre resulta mejor antioxidante que ambos, ya que logró inhibir el desarrollo de la oxidación en un 97,40 %.

Extracto de semillas de vid (Vitis vinifera L.) con actividad antioxidante : concentración, deshidratación y comparación con antioxidantes de uso comercial

Los objetivos de este trabajo fueron: 1) Evaluar un proceso de concentración y deshidratación de un extracto polifenólico de semillas de vid (Vitis vinifera L.). El proceso debía conservar el poder reductor del extracto. 2) Valorar la eficiencia del antioxidante obtenido por comparación con antioxidantes de uso habitual en un producto vegetal susceptible de ser oxidado, tal como el jugo de manzanas. El extracto de semillas de vid fue concentrado a 60°C en un concentrador rotativo al vacío de laboratorio. La concentración de los compuestos fenólicos fue determinada por el método de Folin-Ciocalteu modificado. La actividad antioxidante se determinó midiendo el poder reductor, por el método de Oyaizu. La oxidación del jugo de manzanas se determinó por el método de Özoglu. El extracto concentrado obtenido fue deshidratado por dos métodos: liofilización y secado en lecho de espuma. A igual concentración fenólica, el extracto concentrado mostró mayor capacidad antioxidante que el extracto simple. El extracto secado en lecho de espuma conservó la misma capacidad antioxidante que el extracto concentrado. El extracto liofilizado experimentó una pérdida significativa de la actividad antioxidante. El extracto concentrado de semillas de vid inhibió la oxidación del jugo de manzanas en un 31,51%, el ácido ascórbico en un 2,60% y el dióxido de azufre en un 97,40%.

Geochemistry of oceanic crust of the Cocos Plate

New geochemical data from the Cocos Plate constrain the composition of the input into the Central American subduction zone and demonstrate the extent of influence of the Galápagos Hotspot on the Cocos Plate. Samples include sediments and basalts from Ocean Drilling Program (ODP) Site 1256 outboard of Nicaragua, gabbroic sills from ODP Sites 1039 and 1040, tholeiitic glasses from the Fisher Ridge off northwest Costa Rica, and basalts from the Galápagos Hotspot Track outboard of Central Costa Rica. Site 1256 basalts range from normal to enriched MORB in incompatible elements and have Pb and Nd isotopic compositions within the East Pacific Rise MORB field. The sediments have similar 206Pb/204Pb and only slightly more radiogenic 207Pb/204Pb and 208Pb/204Pb isotope ratios than the basalts. Altered samples from the subducting Galápagos Hotspot Track have similar Nd and Pb isotopic compositions to fresh Galápagos samples but have significantly higher Sr isotopic composition, indicating that the subduction input will have a distinct geochemical signature from Galápagos-type mantle material that may be present in the wedge beneath Costa Rica. Gabbroic sills from Sites 1039 and 1040 in East Pacific Rise (EPR) crust show evidence for influence of the Galápagos Hotspot ?100 km beyond the morphological hotspot track.

Geochemistry and ages at DSDP Holes 81-552A, 81-553A and 81-555

Glauconite-rich sediments have been encountered at two horizons during drilling in the southwest Rockall Plateau. The younger of these horizons lies at the base of a deep-sea ooze sequence and is of early or middle Miocene age. Glauconite formed in situ during periods of nondeposition related to strong bottom-water currents, in water depths of as much as 2500 m - five times greater than previously accepted limits for glauconite formation. The older horizon, of early Eocene age, is a record of the major transgression coincident with the separation of Rockall and Greenland. Isotopic age dating of the Miocene glauconites gives results in relatively close accord with their biostratigraphic age. However, an Eocene (NP12) glauconite gives a highly discrepant date (36.5 m.y. ago). One possible explanation is that the Eocene glauconites have continued to evolve after burial by the diagenetic uptake of potassium from the surrounding mud matrix, a possibility denied to the Miocene glauconites by the relative scarcity of available potassium in the nannofossil-foraminiferal ooze matrix.

Ash layer characteristics and major element oxide content of glasses from ODP Sites 162-982 and 162-985

The Cenozoic volcanic activity on Iceland has been recorded in North Atlantic sediments drilled during several Ocean Drilling Program (ODP)/Deep Sea Drilling Project legs (Legs 104, 151, 152, 162, and 163). Leg 162 (North Atlantic-Arctic Gateways II) recovered ash layers at Sites 982, 985, and 907 (Jansen, Raymo, Blum, et al., 1996, doi:10.2973/odp.proc.ir.162.1996). The revisited Site 907 was first drilled during Leg 151, and the ash from this site has been described in detail by Lacasse et al. (1996, doi:10.2973/odp.proc.sr.151.122.1996) and Werner et al. (1996, doi:10.2973/odp.proc.sr.151.123.1996). Site 982 is located within the Hatton-Rockall Basin on the Rockall Plateau, which is situated west of the British Isles. Site 985 is located northeast of Iceland at the foot of the eastern slope of the Iceland Plateau, adjacent to the Norwegian Basin. Here we report chemical analyses of Neogene tephra layers from Holes 982A, 983B, 982C, 985A, and 985B. The sedimentary sequence at Site 982 spans the lower Miocene-Holocene; Site 985 recovered sediments spanning the upper Oligocene-Holocene. Twenty-two distinct ash layers and ash-bearing sediments were sampled in Holes 982A-982C (Cores 162-982A-16H through 24H, 162-982B-14H through 56X, and 162-982C-15H through 27H), and 59 ash layers were sampled in Holes 985A and 985B (Cores 162-985A-11H through 59X, and 162-985B-11H through 14H). Almost 50% of the sampled ash is strongly altered (predominantly from Site 985). A cluster of altered thin layers in the lower Pliocene of Site 985 (top of Unit III) is remarkable.

Trade, economic growth and environment

The literature on trade openness, economic development, and the environment is largely inconclusive about the environmental consequences of trade. This study review previous studies focusing on treating trade and income as endogenous and estimating the overall impact of trade openness on environmental quality using the instrumental variables technique. The results show that whether or not trade has a beneficial effect on the environment varies depending on the pollutant and the country. Trade is found to benefit the environment in OECD countries. It has detrimental effects, however, on sulfur dioxide (SO2) and carbon dioxide (CO2) emissions in non-OECD countries, although it does lower biochemical oxygen demand (BOD) emissions in these countries. The results also find the impact is large in the long term, after the dynamic adjustment process, although it is small in the short term.