Converting NAD83 GPS heights into NAVD88 elevations with LVGEOID, a hybrid geoid height model for the Long Valley volcanic region, California

A GPS survey of leveling benchmarks done in Long Valley Caldera in 1999 showed that the application of the National Geodetic Survey (NGS) geoid model GEOID99 to tie GPS heights to historical leveling measurements would significantly underestimate the caldera ground deformation (known from other geodetic measurements). The NGS geoid model was able to correctly reproduce the shape of the deformation, but required a local adjustment to give a realistic estimate of the magnitude of the uplift. In summer 2006, the U.S. Geological Survey conducted a new leveling survey along two major routes crossing the Long Valley region from north to south (Hwy 395) and from east to west (Hwy 203 – Benton Crossing). At the same time, 25 leveling bench marks were occupied with dual frequency GPS receivers to provide a measurement of the ellipsoid heights. Using the heights from these two surveys, we were able to compute a precise geoid height model (LVGEOID) for the Long Valley volcanic region. Our results show that although the LVGEOID and the latest NGS GEOID03 model practically coincide in areas outside the caldera, there is a difference of up to 0.2 m between the two models within the caldera. Accounting for this difference is critical when using the geoid height model to estimate the ground deformation due to magmatic or tectonic activity in the calder

The Influence of Viscoelastic Crustal Rheologies on Volcanic Ground Deformation: Insights from Models of Pressure and Volume Change

This is the final version. Avaialble on open access from Wiley via the DOI in this recordInelastic rheological behaviour, such as viscoelasticity, is increasingly utilised in the modelling of volcanic ground deformation, as elevated thermal regimes induced by magmatic systems may necessitate the use of a mechanical model containing a component of time-dependent viscous behaviour. For the modelling of a given amplitude and footprint of ground deformation, incorporating a viscoelastic regime has been shown to reduce the magma reservoir overpressure requirements suggested by elastic models. This phenomenon, however, is restricted to pressure-based analyses and the associated creep behaviour. Viscoelastic materials exhibit additional constitutive time-dependent behaviours, determined by the stress and strain states, that are yet to be analysed in the context of volcanic ground deformation. By utilising a mechanically homogeneous model-space and distinct reservoir evolutions, we provide a comparison of three viscoelastic rheological models, including the commonly implemented Maxwell and Standard Linear Solid configurations, and their time-dependent behaviours from a fundamental perspective. We also investigate the differences between deformation timeseries resulting from a pressurisation or volume change; two contrasting approaches that are assumed to be equivalent through elastic modelling. Our results illustrate that the perceived influence of viscoelasticity is dependent on the mode of deformation, with stress-based pressurisation models imparting enhanced deformation relative to the elastic models, thus reducing pressure requirements. Strain-based volumetric models, however, exhibit reduced levels of deformation and may produce episodes of apparent ground subsidence induced by source inflation or vice versa, due to the relaxation of crustal stresses, dependent on whether the reservoir is modelled to be expanding or contracting, respectively.Natural Environment Research Council (NERC

Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon

<p>Abstract</p> <p>Background</p> <p>Several recent studies have demonstrated the use of Roche 454 sequencing technology for <it>de novo </it>transcriptome analysis. Low error rates and high coverage also allow for effective SNP discovery and genetic diversity estimates. However, genetically diverse datasets, such as those sourced from natural populations, pose challenges for assembly programs and subsequent analysis. Further, estimating the effectiveness of transcript discovery using Roche 454 transcriptome data is still a difficult task.</p> <p>Results</p> <p>Using the Roche 454 FLX Titanium platform, we sequenced and assembled larval transcriptomes for two butterfly species: the Propertius duskywing, <it>Erynnis propertius </it>(Lepidoptera: Hesperiidae) and the Anise swallowtail, <it>Papilio zelicaon </it>(Lepidoptera: Papilionidae). The Expressed Sequence Tags (ESTs) generated represent a diverse sample drawn from multiple populations, developmental stages, and stress treatments.</p> <p>Despite this diversity, > 95% of the ESTs assembled into long (> 714 bp on average) and highly covered (> 9.6× on average) contigs. To estimate the effectiveness of transcript discovery, we compared the number of bases in the hit region of unigenes (contigs and singletons) to the length of the best match silkworm (<it>Bombyx mori</it>) protein--this "ortholog hit ratio" gives a close estimate on the amount of the transcript discovered relative to a model lepidopteran genome. For each species, we tested two assembly programs and two parameter sets; although CAP3 is commonly used for such data, the assemblies produced by Celera Assembler with modified parameters were chosen over those produced by CAP3 based on contig and singleton counts as well as ortholog hit ratio analysis. In the final assemblies, 1,413 <it>E. propertius </it>and 1,940 <it>P. zelicaon </it>unigenes had a ratio > 0.8; 2,866 <it>E. propertius </it>and 4,015 <it>P. zelicaon </it>unigenes had a ratio > 0.5.</p> <p>Conclusions</p> <p>Ultimately, these assemblies and SNP data will be used to generate microarrays for ecoinformatics examining climate change tolerance of different natural populations. These studies will benefit from high quality assemblies with few singletons (less than 26% of bases for each assembled transcriptome are present in unassembled singleton ESTs) and effective transcript discovery (over 6,500 of our putative orthologs cover at least 50% of the corresponding model silkworm gene).</p

Anthropogenic and natural ground deformation in the Hengill geothermal area, Iceland

We investigate crustal deformation due to the extraction of water and steam from a high-enthalpy geothermal reservoir; a common occurrence, yet not well understood. The cause of this deformation can be a change in pressure or in temperature in the reservoir, both of which can be caused by extraction or injection of geothermal fluids. Our study area, the Hengill mountains in SW Iceland, is an active volcanic center and a plate triple junction that hosts two power plants producing geothermal energy. This combination of natural and anthropogenic processes causes a complex displacement field at the surface. We analyze geodetic data—Global Navigation Satellite System and Interferometric Synthetic Aperture Radar—to obtain the surface velocity field, which we then simulate using an inverse modeling approach. We focus on the deformation around the geothermal power plants but need to model the regional tectonic and volcanic deformation as well, because the signals are overlapping. We find that plate motion and a deep contracting body can explain the broad scale signal in the area. Local deformation near the two power plants, Hellisheidi and Nesjavellir, can be explained by extraction of geothermal fluids. We estimate reservoirs extending from 0.6 to 3.0 km depth at Hellisheidi, and 1.0 to 3.0 km depth at Nesjavellir for observed pressure decrease rates of 0.25 MPa/yr and 0.1 MPa/yr, respectively. We find that the main cause for the subsidence in the geothermal area is the observed pressure drawdown

Thermal structure of a gas-permeable lava dome and timescale separation in its response to perturbation

The thermal boundary layer at the surface of a volcanic lava dome is investigated through a continuum model of the thermodynamic advection diffusion processes resulting from magmatic gas flow through the dome matrix. The magmatic gas mass flux, porosity and permeability of the rock are identified as key parameters. New, theoretical, nonlinear steady-state thermal profiles are reported which give a realistic surface temperature of 210 degC for a region of lava dome surface through which a gas flux of 3.5 x 10-3 kg s-1 m-2 passes. This contrasts favourably with earlier purely diffusive thermal models, which cool too quickly. Results are presented for time-dependent perturbations of the steady states as a response to: changes in surface pressure, a sudden rockfall from the lava dome surface, and a change in the magmatic gas mass flux at depth. Together with a generalized analysis using the method of multiple scales, this identifies two characteristic time scales associated with the thermal evolution of a dome carapace: a short time scale of several minutes, over which the magmatic gas mass flux, density, and pressure change to a new quasi-steady-state, and a longer time scale of several days, over which the thermal profile changes to a new equilibrium distribution. Over the longer time scale the dynamic properties of the dome continue to evolve, but only in slavish response to the ongoing temperature evolution. In the light of this time scale separation, the use of surface temperature measurements to infer changes in the magmatic gas flux for use in volcanic hazard prediction is discussed

Integration of micro-gravity and geodetic data to constrain shallow system mass changes at Krafla Volcano, N Iceland

New and previously published micro-gravity data are combined with InSAR data, precise levelling and GPS measurements to produce a model for the processes operating at Krafla volcano, 20 years after its most recent eruption. The data have been divided into two periods: from 1990 to 1995 and from 1996 to 2003 and show that the rate of deflation at Krafla is decaying exponentially. The net micro-gravity change at the centre of the caldera is shown, using the measured Free Air Gradient, to be -85 μGal for the first and -100 μGal for the second period. After consideration of the effects of water extraction by the geothermal power station within the caldera, the net gravity decreases are -73 ± 17 μGal for the first and -65 ± 17 μGal for the second period. These decreases are interpreted in terms of magma drainage. Following a Mogi point source model we calculate the mass decrease to be ~2 x 1010 kg/yr reflecting a drainage rate of ~0.23 m3/s, similar to the ~0.13 m3/s drainage rate previously found at Askja volcano, N-Iceland. Based on the evidence for deeper magma reservoirs and the similarity between the two volcanic systems, we suggest a pressure-link between Askja and Krafla at deeper levels (at the lower crust or the crust-mantle boundary). After the Krafla fires, co-rifting pressure decrease of a deep source at Krafla stimulated the subsequent inflow of magma, eventually affecting conditions along the plate boundary in N-Iceland, as far away as Askja. We anticipate that the pressure of the deeper reservoir at Krafla will reach a critical value and eventually magma will rise from there to the shallow magma chamber, possibly initiating a new rifting episode. We have demonstrated that by examining micro-gravity and geodetic data, our knowledge of active volcanic systems can be significantly improved

Decadal to monthly timescales of magma transfer and reservoir growth at a caldera volcano

International audienceCaldera-forming volcanic eruptions are low-frequency, highimpact events capable of discharging tens to thousands of cubic kilometres of magma explosively on timescales of hours to days, with devastating effects on local and global scales1. Because no such eruption has been monitored during its long build-up phase, the precursor phenomena are not well understood. Geophysical signals obtained during recent episodes of unrest at calderas such as Yellowstone, USA, and Campi Flegrei, Italy, are difficult to interpret, and the conditions necessary for large eruptions are poorly constrained2,3. Here we present a study of pre-eruptive magmatic processes and their timescales using chemically zoned crystals from the 'Minoan' caldera-formingeruption of Santorini volcano,Greece4, which occurred in the late 1600s BC. The results provide insights into how rapidly large silicic systems may pass from a quiescent state to one on the edge of eruption5,6. Despite the large volume of erupted magma4 (40-60 cubic kilometres), and the 18,000-year gestation period between the Minoan eruption and the previous major eruption, most crystals in the Minoan magma record processes that occurred less than about 100 years before the eruption. Recharge of the magma reservoir by large volumes of silicic magma (and some mafic magma) occurred during the century before eruption, and mixing between different silicicmagmabatches was still taking place during the final months. Final assembly of large silicic magma reservoirs may occur on timescales that are geologically very short by comparison with the preceding repose period, with major growth phases immediately before eruption. These observations have implications for the monitoring of long-dormant, but potentially active, caldera systems

Seismic tremor at the 9°50′N East Pacific Rise eruption site

Ocean bottom seismic observations within the 9°50′N region of the East Pacific Rise indicate persistent, low-amplitude tremor activity throughout the October 2003 through February 2007 period of monitoring. These signals exhibit either monochromatic or polychromatic spectral characteristics, with a ∼6 Hz fundamental frequency and up to two harmonics. Individual events cannot be correlated between nearby (<1 km) stations, implying the presence of multiple, small-amplitude sources positioned within the shallow crust. Tremor exhibits a semidiurnal periodicity, with some stations recording activity during times of increasing tidal extension and others detecting tremor signals during times of increasing compression. The amplitude, duration, and rate of activity also correlate positively with fortnightly changes in the amplitude of the tides. These spatiotemporal patterns are consistent with tremor generation in response to tidally modulated fluid flow within a network of shallow cracks. Tremor energy flux is spatially and temporally heterogeneous; however, there are extended periods of greater and lesser activity that can be tracked across portions of the array. Despite their shallow crustal origin, changes in tremor amplitude and spectral character occur in the months prior to a major microearthquake swarm and inferred seafloor spreading event on 22 January 2006, with an increase in the degree of correlation between tremor activity and tidal strain in the weeks leading up to this event. After the spreading event, two eruption-surviving stations near the axis continue to show high rates of tremor activity, whereas these signals are suppressed at the single station recovered from the near-axis flanks. This off-axis quiescence may result from the dike-induced closing of cracks or perhaps from the emplacement of impermeable flows near the station

Exceptional river gorge formation from unexceptional floods

An understanding of rates and mechanisms of incision and knickpoint retreat in bedrock rivers is fundamental to perceptions of landscape response to external drivers, yet only sparse field data are available. Here we present eye witness accounts and quantitative surveys of rapid, amphitheatre-headed gorge formation in unweathered granite from the overtopping of a rock-cut dam spillway by small-moderate floods (~100–1,500 m3 s−1). The amount of erosion demonstrates no relationship with flood magnitude or bedload availability. Instead, structural pattern of the bedrock through faults and joints appears to be the primary control on landscape change. These discontinuities facilitate rapid erosion (>270 m headward retreat; ~100 m incision; and ~160 m widening over 6 years) principally through fluvial plucking and block topple. The example demonstrates the potential for extremely rapid transient bedrock erosion even when rocks are mechanically strong and flood discharges are moderate. These observations are relevant to perceived models of gorge formation and knickpoint retreat