Erupted magma volume estimates at Santiaguito and Pacaya Volcanoes, Guatemala using digital elevation models

High resolution digital elevation models (DEMs) of Santiaguito and Pacaya volcanoes, Guatemala, were used to estimate volume changes and eruption rates between 1954 and 2001. The DEMs were generated from contour maps and aerial photography, which were analyzed in ArcGIS 9.0®. Because both volcanoes were growing substantially over the five decade period, they provide a good data set for exploring effective methodology for estimating volume changes. The analysis shows that the Santiaguito dome complex grew by 0.78 ± 0.07 km3 (0.52 ± 0.05 m3 s-1) over the 1954-2001 period with nearly all the growth occurring on the El Brujo (1958-75) and Caliente domes (1971-2001). Adding information from field data prior to 1954, the total volume extruded from Santiaguito since 1922 is estimated at 1.48 ± 0.19 km3. Santiaguito’s growth rate is lower than most other volcanic domes, but it has been sustained over a much longer period and has undergone a change toward more exogenous and progressively slower extrusion with time. At Santiaguito some of the material being added at the dome is subsequently transported downstream by block and ash flows, mudflows and floods, creating channel shifting and areas of aggradation and erosion. At Pacaya volcano a total volume of 0.21 ± 0.05 km3 was erupted between 1961 and 2001 for an average extrusion rate of 0.17 ± 0.04 m3 s-1. Both the Santiaguito and Pacaya eruption rate estimates reported here are minima, because they do not include estimates of materials which are transported downslope after eruption and data on ashfall which may result in significant volumes of material spread over broad areas. Regular analysis of high resolution DEMs using the methods outlined here, would help quantify the effects of fluvial changes to downstream populated areas, as well as assist in tracking hazards related to dome collapse and eruption.

Areal Geology of the Post-Paleozoic Surface of the Northwest Plains

Probably most of the area included in this report has been ex­amined to some extent by oil geologists, and most, if not all, of the important domes have been discovered and surveyed thoroughly. In parts of the area, the bedrock is covered by glacial drift or alluvium materi­al, but it is reasonable to believe that no new domal structure will be found. This means that surface examination alone will be insufficient in locating new oil fields, so future prospecting will be dependent, to a great extent, on studies of sub-surface stratigraphy.

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.

Psychoendocrine validation of a short measure for assessment of perceived stress management skills in different non-clinical populations

BACKGROUND: We investigated the psychometric properties of a short questionnaire for combined assessment of different perceived stress management skills in the general population and tested whether scores relate to physiological stress reactivity. METHODS: For psychometric evaluation, we determined the factor structure of the questionnaire and investigated its measurement invariance in the participant groups and over time in three different independent samples representing the general population (total N=332). Reliability was tested by estimating test-retest reliability, internal consistency, and item reliabilities. We examined convergent and criterion validity using selected criterion variables. For endocrine validation, 35 healthy non-smoking and medication-free men in a laboratory study and 35 male and female employees in a workplace study underwent an acute standardized psychosocial stress task. We assessed stress management skills and measured salivary cortisol before and several times up to 60 min (workplace study) and 120 min (laboratory study) after stress. Potential confounders were controlled. RESULTS: The factor structure of the questionnaire consists of five scales reflecting acceptably distinct stress management skills such as cognitive strategies, use of social support, relaxation strategies, anger regulation, and perception of bodily tension. This factor structure was stable across participant groups and over time. Internal consistencies, item reliabilities, and test-retest reliabilities met established statistical requirements. Convergent and criterion validity were also established. In both endocrine validation studies, higher stress management skills were independently associated with lower cortisol stress reactivity (p's<.029). CONCLUSIONS: Our findings suggest that the questionnaire has good psychometric properties and that it relates to subjective psychological and objective physiological stress indicators. Therefore, the instrument seems a suitable measure for differential assessment of stress management skills in the general population.

Changes in emotion recognition relate to hypertension status

Background: Emotional processing in essential hypertension beyond self-report questionnaire has hardly been investigated. The aim of this study is to examine associations between hypertension status and recognition of facial affect. Methods: 25 healthy, non-smoking, medication-free men including 13 hypertensive subjects aged between 20 and 65 years completed a computer-based task in order to examine sensitivity of recognition of facial affect. Neutral faces gradually changed to a specific emotion in a pseudo-continuous manner. Slides of the six basic emotions (fear, sadness, disgust, happiness, anger, surprise) were chosen from the „NimStim Set“. Pictures of three female and three male faces were electronically morphed in 1% steps of intensity from 0% to 100% (36 sets of faces with 100 pictures each). Each picture of a set was presented for one second, ranging from 0% to 100% of intensity. Participants were instructed to press a stop button as soon as they recognized the expression of the face. After stopping a forced choice between the six basic emotions was required. As dependent variables, we recorded the emotion intensity at which the presentation was stopped and the number of errors (error rate). Recognition sensitivity was calculated as emotion intensity of correctly identified emotions. Results: Mean arterial pressure was associated with a significantly increased recognition sensitivity of facial affect for the emotion anger (ß = - .43, p = 0.03*, Δ R2= .110). There was no association with the emotions fear, sadness, disgust, happiness, and surprise (p’s > .0.41). Mean arterial pressure did not relate to the mean number of errors for any of the facial emotions. Conclusions: Our findings suggest that an increased blood pressure is associated with increased recognition sensitivity of facial affect for the emotion anger, if a face shows anger. Hypertensives perceive facial anger expression faster than normotensives, if anger is shown.

Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

Ar-40/Ar-39 Evidence for Middle Proterozoic (1300-1500 Ma) Slow Cooling of the Southern Black Hills, South Dakota, Midcontinent, North America: Implications for Early Proterozoic P-T Evolution and Posttectonic Magmatism

Ar-40/Ar-39 total gas and plateau dates from muscovite and biotite in the southern Black Hills, South Dakota, provide evidence for a period of Middle Proterozoic slow cooling. Early Proterozoic (1600-1650 Ma) mica dates were obtained from metasedimentary rocks located in a synformal structure between the Harney Peak and Bear Mountain domes and also south of Bear Mountain. Metamorphic rocks from the dome areas and undeformed samples of the similar to 1710 Ma Harney Peak Granite (HPG) yield Middle Proterozoic mica dates (similar to 1270-1500 Ma). Two samples collected between the synform and Bear Mountain dome yield intermediate total gas mica dates of similar to 1550 Ma. We suggest two end-member interpretations to explain the map pattern of cooling ages: (1) subhorizontal slow cooling of an area which exhibits variation in mica Ar retention intervals or (2) mild folding of a Middle Proterozoic (similar to 1500 Ma) similar to 300 degrees C isotherm. According to the second interpretation, the preservation of older dates between the domes may reflect reactivation of a preexisting synformal structure (and downwarping of relatively cold rocks) during a period of approximately east-west contraction and slow uplift during the Middle Proterozoic. The mica data, together with hornblende data from the Black Hills published elsewhere, indicate that the ambient country-rock temperature at the 3-4 kbar depth of emplacement of the HPG was between 350 degrees C and 500 degrees C, suggesting that the average upper crustal geothermal gradient was 25 degrees-40 degrees C/km prior to intrusion. The thermochronologic data suggest HPG emplacement was followed by a similar to 200 m.y. period of stability and tectonic quiescence with little uplift. We propose that crust thickened during the Early Proterozoic was uplifted and erosionally(?) thinned prior to similar to 1710 Ma and that the HPG magma was emplaced into isostatically stable crust of relatively normal thickness. We speculate that uplift and crustal thinning prior to HPG intrusion was the result of differential thinning of the subcrustal lithosphere beneath the Black Hills. If so, this process would have also caused an increase in mantle heat flux across the Moho and triggered vapor-absent melting of biotite to produce the HPG magma. This scenario for posttectonic granite generation is supported, in part, by the fact that in the whole of the Black Hills, the HPG is spatially associated with the deepest exposed Early Proterozoic country rock.

The deglaciation history of the Simplon region (southern Swiss Alps) constrained by 10Be exposure dating of ice-molded bedrock surfaces

The deglaciation history of the Swiss Alps after the Last Glacial Maximum involved the decay of several ice domes and the subsequent disintegration of valley glaciers at high altitude. Here we use bedrock exposure dating to reconstruct the temporal and spatial pattern of ice retreat at the Simplon Pass (altitude: ∼2000 m) located 40 km southwest of the ‘Rhône ice dome’. Eleven 10Be exposure ages from glacially polished quartz veins and ice-molded bedrock surfaces cluster tightly between 13.5 ± 0.6 ka and 15.4 ± 0.6 ka (internal errors) indicating that the Simplon Pass depression became ice-free at 14.1 ± 0.4 ka (external error of mean age). This age constraint is interpreted to record the melting of the high valley glaciers in the Simplon Pass region during the warm Bølling–Allerød interstadial shortly after the Oldest Dryas stadial. Two bedrock samples collected a few hundred meters above the pass depression yield older 10Be ages of 17.8 ± 0.6 ka and 18.0 ± 0.6 ka. These ages likely reflect the initial downwasting of the Rhône ice dome and the termination of the ice transfluence from the ice dome across the Simplon Pass toward the southern foreland. There, the retreat of the piedmont glacier in Val d’Ossola was roughly synchronous with the decay of the Rhône ice dome in the interior of the mountain belt, as shown by 10Be ages of 17.7 ± 0.9 ka and 16.1 ± 0.6 ka for a whaleback at ∼500 m elevation near Montecrestese in northern Italy. In combination with well-dated paleoclimate records derived from lake sediments, our new age data suggest that during the deglaciation of the European Alps the decay of ice domes was approximately synchronous with the retreat of piedmont glaciers in the foreland and was followed by the melting of high-altitude valley glaciers after the transition from the Oldest Dryas to the Bølling–Allerød, when mean annual temperatures rose rapidly by ∼3 °C.

Deep drilling of glaciers: Russian projects in the Arctic (1975-1995)

The data collection "Deep Drilling of Glaciers: Soviet-Russian projects in Arctic, 1975-1995" was collected by the following basic considerations: - compilation of deep (>100 m) drilling projects on Arctic glaciers, using data of (a) publications; (b) archives of IGRAN; (c) personal communication of project participants; - documentation of parameters, references. Accuracy of data and techniques applied to determine different parameters are not evaluated. The accuracy of some geochemical parameters (up to 1984 and heavy metalls) is uncertain. Most reconstructions of ice core age and of annual layer thickness are discussed; - digitizing of published diagrams (in case, when original numerical data were lost) and subsequent data conversion to equal range series and adjustment to the common units. Therefore, the equal-range series were calculated from original data or converted from digitized chart values as indicated in the metadata. For the methodological purpose, the equal-range series obtained from original and reconstructed data were compared repeatedly; the systematic difference was less then 5-7%. Special attention should be given to the fact, that the data for individual ice core parameters varies, because some parameters were originally measured or registered. Parameters were converted in equal-range series using 2 m steps; - two or more parameter values were determined, then the mean-weighted (i.e. accounting the sample length) value is assigned to the entire interval; - one parameter value was determined, measured or registered independently from the parameter values in depth intervals which over- and underlie it, then the value is assigned to the entire interval; - one parameter value was determined, measured or registered for two adjoining depth intervals, then the specific value is assigned to the depth interval, which represents >75% of sample length ; if each of adjoining depth intervals represents <75% of sample length, then the correspondent parameter value is assigned to both intervals of depth. This collection of ice core data (version 2000) was made available through the EU funded QUEEN project by S.M. Arkhipov, Moscow.

The 2008 DGFI Realization of the ITRS: DTRF2008 (data)

The DTRF2008 is a realization of the International Terrestrial Reference System ITRS. The DTRF2008 consists of station positions and velocities of global distributed observing stations of the space geodetic observation techniques VLBI, SLR, GPS and DORIS. The DTRF2008 was released in May 2010 and includes the observation data of the techniques up to and including 2008. The observation data are processed and submitted by the corresponding international services: IGS (International GNSS Service, http://igscb.jpl.nasa.gov) IVS (International VLBI Service, http://ivscc.gsfc.nasa.gov) ILRS (International Laser Ranging Service, http://ilrs.gsfc.nasa.gov) IDS (International DORIS Service, http://ids-doris.org). The DTRF2008 is an independent ITRS realization, which is computed on the basis of the same input data as the ITRF2008 (IGN, Paris). Both realizations differ with respect to their computation strategies: while the ITRF2008 is based on the combination of solutions, the DTRF2008 is computed by the combination of normal equations. The DTRF2008 comprises the coordinates of 559 GPS-, 106 VLBI-, 122 SLR- and 132 DORIS-stations. The reference epoch is 1.1.2005, 0h UTC. The Earth Orientation Parameters (EOP) - that means the coordinates of the terrestrial and the celestial pole, UT1-UTC and the Length of Day (LOD) - were simultaneously estimated with the station coordinates. The EOP time series cover the period of 1983 to 2008. The station names are the official IERS indications: cdp numbers or 4-character IDs and DOMES numbers (http://itrf.ensg.ign.fr/doc_ITRF/iers_sta_list.txt). The solution is available in different file formats (SINEX and SSC), see below. A detailed description of the solution is given by Seitz M. et al. (2012). The results of a comparison of DTRF2008 and ITRF2008 is given by Seitz M. et al. (2013). More information as well as residual time series of the station positions can be made available by request.

Geochemistry and isotopes at DSDP Hole 79-546

At Site 546, below the Mazagan Escarpment at a water depth of 4 km, 36 m of salt rock was cored from the top of one of a field of salt domes. The core was studied by thin section and a variety of geochemical procedures. The salt rock contains 0.1 to 3% carnallite and lesser amounts of sylvite and polyhalite, which with the corresponding high level of bromide place it within the potash evaporite facies. The bromide profile is of a dominantly marine evaporite deposited in moderately shallow brine which, however, was not repeatedly desiccated. A mineralogical argument suggests that the brine surface was not below sea level. An average of about 5% elastics, with dispersed anhydrite, darken the salt rock to deep shades of red, brown, and gray green. Most of the included materials are in highly deformed boudins or dispersions in the salt rock that has also undergone cataclasis in a subsequent, probably tectonic, deformation. The salt rock is slightly deficient in anhydrite, and the usual separate beds and laminae of anhydrite are virtually absent. Stable isotope ratios of sulfur and oxygen in the sulfate are clearly derived from sea water of Permian to Scythian age, in contrast to the late Triassic or Early Jurassic age of evaporites onshore in Morocco and Portugal and the corresponding evaporites offshore Maritime Canada. In contrast to those evaporites off the axis of Atlantic rifting, the salt at Site 546 may have been deposited in a very early central rift fed by marine waters from Tethys through the Gibraltar or South Atlas fracture zones.

(Supplemet S1) Silicon oxide and Potassium + Sodium oxide from the PACManus hydrothermal area

This study presents a systematic analysis and interpretation of autonomous underwater vehicle-based microbathymetry combined with remotely operated vehicle (ROV) video recordings, rock analyses and temperaturemeasurements within the PACManus hydrothermal area located on Pual Ridge in the Bismarck Sea of eastern Manus Basin. The data obtained during research cruise Magellan-06 and So-216 provides a framework for understanding the relationship between the volcanism, tectonismand hydrothermal activity. PACManus is a submarine felsic vocanically-hosted hydrothermal area that hosts multiple vent fields locatedwithin several hundredmeters of one another but with different fluid chemistries, vent temperatures and morphologies. The total area of hydrothermal activity is estimated to be 20,279m**2. Themicrobathymetrymaps combinedwith the ROV video observations allow for precise high-resolution mapping estimates of the areal extents of hydrothermal activity.We find the distribution of hydrothermal fields in the PACManus area is primarily controlled by volcanic features that include lava domes, thick andmassive blocky lava flows, breccias and feeder dykes. Spatial variation in the permeability of local volcanic facies appears to control the distribution of venting within a field.We define a three-stage chronological sequence for the volcanic evolution of the PACManus based on lava flow morphology, sediment cover and lava SiO2 concentration. In Stage-1, sparsely to moderately porphyritic dacite lavas (68-69.8 wt.% SiO2) erupted to form domes or cryptodomes. In Stage-2, aphyric lava with slightly lower SiO2 concentrations (67.2-67.9 wt.% SiO2) formed jumbled and pillowed lava flows. In the most recent phase Stage-3, massive blocky lavaswith 69 to 72.5wt.% SiO2were erupted throughmultiple vents constructing a volcanic ridge identified as the PACManus neovolcanic zone. The transition between these stages may be gradual and related to progressive heating of a silicic magma following a recharge event of hot, mantle-derived melts.