Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Oceans”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34Linked to original sources
USGS2006

Overview of the Opportunity Mars Exploration Rover Mission to Meridiani Planum: Eagle Crater to Purgatory Ripple

The Mars Exploration Rover Opportunity touched down at Meridiani Planum in January 2004 and since then has been conducting observations with the Athena science payload. The rover has traversed more than 5 km, carrying out the first outcrop‐scale investigation of sedimentary rocks on Mars. The rocks of Meridiani Planum are sandstones formed by eolian and aqueous reworking of sand grains that are composed of mixed fine‐grained siliciclastics and sulfates. The siliciclastic fraction was produced by chemical alteration of a precursor basalt. The sulfates are dominantly Mg‐sulfates and also include Ca‐sulfates and jarosite. The stratigraphic section observed to date is dominated by eolian bedforms, with subaqueous current ripples exposed near the top of the section. After deposition, interaction with groundwater produced a range of diagenetic features, notably the hematite‐rich concretions known as “blueberries.” The bedrock at Meridiani is highly friable and has undergone substantial erosion by wind‐transported basaltic sand. This sand, along with concretions and concretion fragments eroded from the rock, makes up a soil cover that thinly and discontinuously buries the bedrock. The soil surface exhibits both ancient and active wind ripples that record past and present wind directions. Loose rocks on the soil surface are rare and include both impact ejecta and meteorites. While Opportunity's results show that liquid water was once present at Meridiani Planum below and occasionally at the surface, the environmental conditions recorded were dominantly arid, acidic, and oxidizing and would have posed some significant challenges to the origin of life.

Journal of Geophysical Research E: Planets↗
USGS2007

Visible and near-infrared multispectral analysis of rocks at Meridiani Planum, Mars, by the Mars Exploration Rover Opportunity

Multispectral measurements in the visible and near infrared of rocks at Meridiani Planum by the Mars Exploration Rover Opportunity's Pancam are described. The Pancam multispectral data show that the outcrops of the Burns formation consist of two main spectral units which in stretched 673, 535, 432 nm color composites appear buff- and purple-colored. These units are referred to as the HFS and LFS spectral units based on higher and lower values of 482 to 535 nm slope. Spectral characteristics are consistent with the LFS outcrop consisting of less oxidized, and the HFS outcrop consisting of more oxidized, iron-bearing minerals. The LFS surfaces are not as common and appear, primarily, at the distal ends of outcrop layers and on steep, more massive surfaces, locations that are subject to greater eolian erosion. Consequently, the HFS surfaces are interpreted as a weathering rind. Further inherent spectral differences between layer's and between different outcrop map units, both untouched and patches abraded by the rover's Rock Abrasion Tool, are also described. Comparisons of the spectral parameters of the Meridiani outcrop with a set of laboratory reflectance measurements of Fe3+-bearing minerals show that the field of outcrop measurements plots near the fields of hematite, ferrihydrite, poorly crystalline goethite, and schwertmannite. Rind and fracture fill materials, observed intermittently at outcrop exposures, are intermediate in their spectral character between both the HFS and LFS spectral classes and other, less oxidized, surface materials (basaltic sands, spherules, and cobbles). Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2007

Dust deposition on the Mars Exploration Rover Panoramic Camera (Pancam) calibration targets

The Panoramic Camera (Pancam) on the Mars Exploration Rover mission has acquired in excess of 20,000 images of the Pancam calibration targets on the rovers. Analysis of this data set allows estimates of the rate of deposition and removal of aeolian dust on both rovers. During the first 150-170 sols there was gradual dust accumulation on the rovers but no evidence for dust removal. After that time there is ample evidence for both dust removal and dust deposition on both rover decks. We analyze data from early in both rover missions using a diffusive reflectance mixing model. Assuming a dust settling rate proportional to the atmospheric optical depth, we derive spectra of optically thick layers of airfall dust that are consistent with spectra from dusty regions on the Martian surface. Airfall dust reflectance at the Opportunity site appears greater than at the Spirit site, consistent with other observations. We estimate the optical depth of dust deposited on the Spirit calibration target by sol 150 to be 0.44 ?? 0.13. For Opportunity the value was 0.39 ?? 0.12. Assuming 80% pore space, we estimate that the dust layer grew at a rate of one grain diameter per ???100 sols on the Spirit calibration target. On Opportunity the rate was one grain diameter per ???125 sols. These numbers are consistent with dust deposition rates observed by Mars Pathfinder taking into account the lower atmospheric dust optical depth during the Mars Pathfinder mission. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2007

Coordinated analyses of orbital and spirit rover data to characterize surface materials on the cratered plains of Gusev Crater, Mars

Comparison of the Mars Exploration Rover Spirit's Pancam (0.4 to 1.0 ??m) and Mars Express Observatoire pour la Mineralogie l'Eau, les Glaces et l'Activite?? (OMEGA) (0.4 to 2.5 ??m) spectral reflectance data over Spirit's traverses shows that Gusev cratered plains are dominated by nanophase ferric-oxide-rich dust covering weakly altered basaltic sands. This interpretation is also consistent with both observations from OMEGA data covering plains beyond the traverse region and interpretations of data from the other payload instruments on the Spirit Rover. OMEGA observations of relatively low albedo regions where dust has presumably been stripped by dust devils show negative spectral reflectance slopes from 1.5 to 2.5 ??m and moderately masked spectral features which are indicative of olivine or pyroxene. High-albedo regions north and south of the Spirit landing site have flat spectral reflectance slopes and few spectral features, although all spectra have a nanophase ferric-oxide absorption edge between 0.4 and 0.75 ??m. Comparison of THEMIS-derived thermal inertia values with OMEGA-derived spectral parameters shows that although the dust cover can be optically thick (0.4 to 2.5 ??m wavelength region) in some areas, it is not thick enough (???1 cm) to mask the thermal inertia of the underlying substrate for areas included in this study. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2006

Form drag in rivers due to small-scale natural topographic features: 1. Regular sequences

Small-scale topographic features are commonly found on the boundaries of natural rivers, streams, and floodplains. A simple method for determining the form drag on these features is presented, and the results of this model are compared to laboratory measurements. The roughness elements are modeled as Gaussian-shaped features defined in terms of three parameters: a protrusion height, H; a streamwise length scale, ??; and a spacing between crests, ??. This shape is shown to be a good approximation to a wide variety of natural topographic bank features. The form drag on an individual roughness element embedded in a series of identical elements is determined using the drag coefficient of the individual element and a reference velocity that includes the effects of roughness elements further upstream. In addition to calculating the drag on each element, the model determines the spatially averaged total stress, skin friction stress, and roughness height of the boundary. The effects of bank roughness on patterns of velocity and boundary shear stress are determined by combining the form drag model with a channel flow model. The combined model shows that drag on small-scale topographic features substantially alters the near-bank flow field. These methods can be used to improve predictions of flow resistance in rivers and to form the basis for fully predictive (no empirically adjusted parameters) channel flow models. They also provide a foundation for calculating the near-bank boundary shear stress fields necessary for determining rates of sediment transport and lateral erosion.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Form drag in rivers due to small-scale natural topographic features: 2. Irregular sequences

The size, shape, and spacing of small-scale topographic features found on the boundaries of natural streams, rivers, and floodplains can be quite variable. Consequently, a procedure for determining the form drag on irregular sequences of different-sized topographic features is essential for calculating near-boundary flows and sediment transport. A method for carrying out such calculations is developed in this paper. This method builds on the work of Kean and Smith (2006), which describes the flow field for the simpler case of a regular sequence of identical topographic features. Both approaches model topographic features as two-dimensional elements with Gaussian-shaped cross sections defined in terms of three parameters. Field measurements of bank topography are used to show that (1) the magnitude of these shape parameters can vary greatly between adjacent topographic features and (2) the variability of these shape parameters follows a lognormal distribution. Simulations using an irregular set of topographic roughness elements show that the drag on an individual element is primarily controlled by the size and shape of the feature immediately upstream and that the spatial average of the boundary shear stress over a large set of randomly ordered elements is relatively insensitive to the sequence of the elements. In addition, a method to transform the topography of irregular surfaces into an equivalently rough surface of regularly spaced, identical topographic elements also is given. The methods described in this paper can be used to improve predictions of flow resistance in rivers as well as quantify bank roughness.

Journal of Geophysical Research F: Earth Surface↗
USGS2007

Glacier changes in southeast Alaska and northwest British Columbia and contribution to sea level rise

The digital elevation model (DEM) from the 2000 Shuttle Radar Topography Mission (SRTM) was differenced from a composite DEM based on air photos dating from 1948 to 1987 to detennine glacier volume changes in southeast Alaska and adjoining Canada. SRTM accuracy was assessed at ??5 in through comparison with airborne laser altimetry and control locations measured with GPS. Glacier surface elevations lowered over 95% of the 14,580 km2 glacier-covered area analyzed, with some glaciers thinning as much as 640 in. A combination of factors have contributed to this wastage, including calving retreats of tidewater and lacustrine glaciers and climate change. Many glaciers in this region are particularly sensitive to climate change, as they have large areas at low elevations. However, several tidewater glaciers that had historically undergone calving retreats are now expanding and appear to be in the advancing stage of the tidewater glacier cycle. The net average rate of ice loss is estimated at 16.7 ?? 4.4 km3/yr, equivalent to a global sea level rise contribution of 0.04 ?? 0.01 mm/yr. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2007

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗
USGS2007

Modeling englacial radar attenuation at Siple Dome, West Antarctica, using ice chemistry and temperature data

The radar reflectivity of an ice-sheet bed is a primary measurement for discriminating between thawed and frozen beds. Uncertainty in englacial radar attenuation and its spatial variation introduces corresponding uncertainty in estimates of basal reflectivity. Radar attenuation is proportional to ice conductivity, which depends on the concentrations of acid and sea-salt chloride and the temperature of the ice. We synthesize published conductivity measurements to specify an ice-conductivity model and find that some of the dielectric properties of ice at radar frequencies are not yet well constrained. Using depth profiles of ice-core chemistry and borehole temperature and an average of the experimental values for the dielectric properties, we calculate an attenuation rate profile for Siple Dome, West Antarctica. The depth-averaged modeled attenuation rate at Siple Dome (20.0 ?? 5.7 dB km-1) is somewhat lower than the value derived from radar profiles (25.3 ?? 1.1 dB km-1). Pending more experimental data on the dielectric properties of ice, we can match the modeled and radar-derived attenuation rates by an adjustment to the value for the pure ice conductivity that is within the range of reported values. Alternatively, using the pure ice dielectric properties derived from the most extensive single data set, the modeled depth-averaged attenuation rate is 24.0 ?? 2.2 dB km-1. This work shows how to calculate englacial radar attenuation using ice chemistry and temperature data and establishes a basis for mapping spatial variations in radar attenuation across an ice sheet. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Quantity-activity relationship of denitrifying bacteria and environmental scaling in streams of a forested watershed

The spatial variability of subreach denitrification rates in streams was evaluated with respect to controlling environmental conditions, molecular examination of denitrifying bacteria, and dimensional analysis. Denitrification activities ranged from 0 and 800 ng-N gsed-1 d-1 with large variations observed within short distances (<50 m) along stream reaches. A log-normal probability distribution described the range in denitrification activities and was used to define low (16% of the probability distributibn), medium (68%), and high (16%) denitrification potential groups. Denitrifying bacteria were quantified using a competitive polymerase chain reaction (cPCR) technique that amplified the nirK gene that encodes for nitrite reductase. Results showed a range of nirK quantities from 103 to 107 gene-copy-number gsed.-1 A nonparametric statistical test showed no significant difference in nirK quantifies among stream reaches, but revealed that samples with a high denitrification potential had significantly higher nirK quantities. Denitrification activity was positively correlated with nirK quantities with scatter in the data that can be attributed to varying environmental conditions along stream reaches. Dimensional analysis was used to evaluate denitrification activities according to environmental variables that describe fluid-flow properties, nitrate and organic material quantities, and dissolved oxygen flux. Buckingham's pi theorem was used to generate dimensionless groupings and field data were used to determine scaling parameters. The resulting expressions between dimensionless NO3- flux and dimensionless groupings of environmental variables showed consistent scaling, which indicates that the subreach variability in denitrification rates can be predicted by the controlling physical, chemical, and microbiological conditions. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2007

State factor relationships of dissolved organic carbon and nitrogen losses from unpolluted temperate forest watersheds

We sampled 100 unpolluted, old-growth forested watersheds, divided among 13 separate study areas over 5 years in temperate southern Chile and Argentina, to evaluate relationships among dominant soil-forming state factors and dissolved carbon and nitrogen concentrations in watershed streams. These watersheds provide a unique opportunity to examine broad-scale controls over carbon (C) and nitrogen (N) biogeochemistry in the absence of significant human disturbance from chronic N deposition and land use change. Variations in the ratio dissolved organic carbon (DOC) to nitrogen (DON) in watershed streams differed by underlying soil parent material, with average C:N = 29 for watersheds underlain by volcanic ash and basalt versus C:N = 73 for sedimentary and metamorphic parent materials, consistent with stronger adsorption of low C:N hydrophobic materials by amorphous clays commonly associated with volcanic ash and basalt weathering. Mean annual precipitation was related positively to variations in both DOC (range: 0.2-9.7 mg C/L) and DON (range: 0.008-0.135 mg N/L) across study areas, suggesting that variations in water volume and concentration may act synergistically to influence C and N losses across dry to wet gradients in these forest ecosystems. Dominance of vegetation by broadleaf versus coniferous trees had negligible effects on organic C and N concentrations in comparison to abiotic factors. We conclude that precipitation volume and soil parent material are important controls over chemical losses of dissolved organic C and N from unpolluted temperate forest watersheds. Our results raise the possibility that biotic imprints on watershed C and N losses may be less pronounced in naturally N-poor forests than in areas impacted by land use change and chronic N deposition. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2007

Nutrient (N, P) loads and yields at multiple scales and subbasin types in the Yukon River basin, Alaska

Loads and yields of dissolved and particulate nitrogen (N) and phosphorus (P) were measured and modeled at three locations on the Yukon River (YR) and on the Tanana and Porcupine Rivers in Alaska during 2001-2005. Total export of N and P upstream of Yukon Delta averaged 120 Gg N a-1 and 56 Gg P a-1, respectively, with 43.5% of total N (TN) as dissolved organic N, and 98% of total P (TP) as particulate phosphorus. Approximately half of the annual export of TN and TP occurred during spring. Hydrologic yields ofTN (5.6-13.3 mmol N m-2 a-1) and TP (0.8-9.0 mmol P m-2 a-1) were least in the Porcupine basin and greatest in the Tanana basin and were proportional to water yield. Comparison of current and historical dissolved organic matter (DOM) export from the basin indicates decreased DON export with respect to total water discharge during summer and autumn in recent decades. Any possible climate-related change in annual water discharge will result in proportional changes in N and P export.

Journal of Geophysical Research G: Biogeosciences↗
USGS2007

The role of historical fire disturbance in the carbon dynamics of the pan-boreal region: A process-based analysis

[1] Wildfire is a common occurrence in ecosystems of northern high latitudes, and changes in the fire regime of this region have consequences for carbon feedbacks to the climate system. To improve our understanding of how wildfire influences carbon dynamics of this region, we used the process-based Terrestrial Ecosystem Model to simulate fire emissions and changes in carbon storage north of 45°N from the start of spatially explicit historically recorded fire records in the twentieth century through 2002, and evaluated the role of fire in the carbon dynamics of the region within the context of ecosystem responses to changes in atmospheric CO 2 concentration and climate. Our analysis indicates that fire plays an important role in interannual and decadal scale variation of source/sink relationships of northern terrestrial ecosystems and also suggests that atmospheric CO 2 may be important to consider in addition to changes in climate and fire disturbance. There are substantial uncertainties in the effects of fire on carbon storage in our simulations. These uncertainties are associated with sparse fire data for northern Eurasia, uncertainty in estimating carbon consumption, and difficulty in verifying assumptions about the representation of fires that occurred prior to the start of the historical fire record. To improve the ability to better predict how fire will influence carbon storage of this region in the future, new analyses of the retrospective role of fire in the carbon dynamics of northern high latitudes should address these uncertainties.

Journal of Geophysical Research G: Biogeosciences↗
USGS2008

Biogeochemical legacy of prescribed fire in a giant sequoia - Mixed conifer forest: A 16-year record of watershed balances

The effects of prescription burning on watershed balances of major ions in mixed conifer forest were examined in a 16-year paired catchment study in Sequoia National Park, California. The objective was to determine whether fire-related changes in watershed balances persist as long as estimated low-end natural fire-return intervals (???10 years), and whether cumulative net export caused by fire could deplete nutrient stocks between successive fires. Inputs (wet + dry deposition) and outputs (stream export) of N, S, Cl-, HCO3-, Ca2+, Mg2+, Na+, K+, H+, and SiO2 were measured for 7 years preceding, and 9 years following, a prescribed burn of one of the catchments. After fire, runoff coefficients increased by 7% (in dry years) to 35% (in wet years). Inorganic N was elevated in stream water for 3 years after fire. Increased export of water, SO42-, Cl-,SiO2, and base cations continued through the end of the study. Pools and processes attributed to fire led to the cumulative loss, per hectare, of 1.2 kg N, 16 kg S, 25 kg Cl-, 130 kg Ca2+, 19 kg Mg2+, 71 kg Na+, 29 kg K+ and 192 kg Si, above that predicted by prefire regression equations relating export in the paired catchments. This additional export equaled <1% of the N, up to one-third of the Ca and Mg, and up to three-fourths of the K, contained in the forest floor prior to combustion. Changes in watershed balances indicated that low-end natural fire-return intervals may prevent complete reaccumulation of several elements between fires. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2008

Constitutive relationships and physical basis of fault strength due to flash heating

We develop a model of fault strength loss resulting from phase change at asperity contacts due to flash heating that considers a distribution of contact sizes and nonsteady state evolution of fault strength with displacement. Laboratory faulting experiments conducted at high sliding velocities, which show dramatic strength reduction below the threshold for bulk melting, are well fit by the model. The predicted slip speed for the onset of weakening is in the range of 0.05 to 2 m/s, qualitatively consistent with the limited published observations. For this model, earthquake stress drops and effective shear fracture energy should be linearly pressure-dependent, whereas the onset speed may be pressure-independent or weakly pressure-dependent. On the basis of the theory, flash weakening is expected to produce large dynamic stress drops, small effective shear fracture energy, and undershoot. Estimates of the threshold slip speed, stress drop, and fracture energy are uncertain due to poor knowledge of the average ontact dimension, shear zone thickness and gouge particle size at seismogenic depths. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Monte Carlo method for determining earthquake recurrence parameters from short paleoseismic catalogs: Example calculations for California

Paleoearthquake observations often lack enough events at a given site to directly define a probability density function (PDF) for earthquake recurrence. Sites with fewer than 10-15 intervals do not provide enough information to reliably determine the shape of the PDF using standard maximum-likelihood techniques (e.g., Ellsworth et al., 1999). In this paper I present a method that attempts to fit wide ranges of distribution parameters to short paleoseismic series. From repeated Monte Carlo draws, it becomes possible to quantitatively estimate most likely recurrence PDF parameters, and a ranked distribution of parameters is returned that can be used to assess uncertainties in hazard calculations. In tests on short synthetic earthquake series, the method gives results that cluster around the mean of the input distribution, whereas maximum likelihood methods return the sample means (e.g., NIST/SEMATECH, 2006). For short series (fewer than 10 intervals), sample means tend to reflect the median of an asymmetric recurrence distribution, possibly leading to an overestimate of the hazard should they be used in probability calculations. Therefore a Monte Carlo approach may be useful for assessing recurrence from limited paleoearthquake records. Further, the degree of functional dependence among parameters like mean recurrence interval and coefficient of variation can be established. The method is described for use with time-independent and time-dependent PDFs, and results from 19 paleoseismic sequences on strike-slip faults throughout the state of California are given.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Finding concealed active faults: Extending the southern Whidbey Island fault across the Puget Lowland, Washington

The southern Whidbey Island fault zone (SWIF), as previously mapped using borehole data, potential field anomalies, and marine seismic reflection surveys, consists of three subparallel, northwest trending strands extending ∼100 km from near Vancouver Island to the northern Puget Lowland. East of Puget Sound, the SWIF makes landfall between the cities of Seattle and Everett but is concealed beneath a thick mantle of young glacial deposits and vegetation. A ∼20‐km‐wide, northwest trending swath of subparallel, low‐amplitude aeromagnetic anomalies crosses this region of the Puget Lowland and is on strike with the SWIF. The most prominent aeromagnetic anomaly, the Cottage Lake lineament, extends at least 18 km and lies approximately on strike with the SWIF on Whidbey Island. Subtle scarps and topographic lineaments on Pleistocene surfaces, visible on high‐resolution lidar topography at a number of locations along the SWIF, lie on or near these magnetic anomalies. In the field, scarps exhibit northeast‐side‐up and vertical relief of 1 to 5 m. Excavations across several lidar scarps lying on or near magnetic anomalies show evidence for multiple folding and faulting events since deglaciation, most likely above buried reverse/oblique faults. Excavations in areas away from magnetic anomalies do not show evidence of tectonic deformation. In total, paleoseismological evidence suggests that the SWIF produced at least four earthquakes since deglaciation about 16,400 years ago, the most recent less than 2700 years ago.

Washington↗
USGS2008

Variability of passive gas emissions, seismicity, and deformation during crater lake growth at White Island Volcano, New Zealand, 2002-2006

We report on 4 years of airborne measurements of CO2, SO2, and H2S emission rates during a quiescent period at White Island volcano, New Zealand, beginning in 2003. During this time a significant crater lake emerged, allowing scrubbig processes to be investigated. CO2 emissions varied from a baseline of 250 to >2000 t d-1 and demonstrated clear annual cycling that was consistent with numbers of earthquake detections and annual changes in sea level. The annual variability was found to be most likely related to increases in the strain on the volcano during sea level highs, temporarily causing fractures to reduce in size in the upper conduit. SO2 emissions varied from 0 to >400 t d-1 and were clearly affected by scrubbing processes within the first year of take development. Scrubbing caused increases of SO42- and Cl- in lake waters, and the ratio of carbon to total sulphur suggested that elemental sulphur deposition was also significant in the lake during the first year. Careful measurements of the lake level and chemistry allowed estimates of the rate of H2O(g) and HCl(g) input into the lake and suggested that the molar abundances of major gas species (H2O, CO2, SO2, and HCl) during this quiescent phase were similar to fumarolic ratios observed between earlier eruptive periods. The volume of magma estimated from CO2 emissions (0.0 15-0.04 km3) was validated by Cl- increases in the lake, suggesting that the gas and magma are transported from deep to shallow depths as a closed system and likely become open in the upper conduit region. The absence of surface deformation further leads to a necessity of magma convection to supply and remove magma from the degassing depths. Two models of convection configurations are discussed. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗