SEARCH · USGS Science

Results for “The R Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,783 records · Page 99Linked to original sources

USGS2011

The influence of irrigation water on the hydrology and lake water budgets of two small arid-climate lakes in Khorezm, Uzbekistan

Little is known regarding the origins and hydrology of hundreds of small lakes located in the western Uzbekistan province of Khorezm, Central Asia. Situated in the Aral Sea Basin, Khorezm is a productive agricultural region, growing mainly cotton, wheat, and rice. Irrigation is provided by an extensive canal network that conveys water from the Amu Darya River (AD) throughout the province. The region receives on average 10 cm/year of precipitation, yet potential evapotranspiration exceeds this amount by about 15 times. It was hypothesized that the perennial existence of the lakes of interest depends on periodic input of excess irrigation water. This hypothesis was investigated by studying two small lakes in the region, Tuyrek and Khodjababa. In June and July 2008, surface water and shallow groundwater samples were collected at these lake systems and surrounding communities and analyzed for δ 2 H, δ 18 O, and major ion hydrochemistry to determine water sources. Water table and lake surface elevations were monitored, and the local aquifer characteristics were determined through aquifer tests. These data and climate data from a Class A evaporation pan and meteorological stations were used to estimate water budgets for both lakes. Lake evaporation was found to be about 0.7 cm/day during the study period. Results confirm that the waters sampled at both lake systems and throughout central Khorezm were evaporated from AD water to varying degrees. Together, the water budgets and stable isotope and major ion hydrochemistry data suggest that without surface water input from some source (i.e. excess irrigation water), these and other Khorezm lakes with similar hydrology may decrease in volume dramatically, potentially to the point of complete desiccation.

Khorezm↗

USGS2014

Identifying dominant controls on hydrologic parameter transfer from gauged to ungauged catchments: a comparative hydrology approach

Daily streamflow information is critical for solving various hydrologic problems, though observations of continuous streamflow for model calibration are available at only a small fraction of the world’s rivers. One approach to estimate daily streamflow at an ungauged location is to transfer rainfall–runoff model parameters calibrated at a gauged (donor) catchment to an ungauged (receiver) catchment of interest. Central to this approach is the selection of a hydrologically similar donor. No single metric or set of metrics of hydrologic similarity have been demonstrated to consistently select a suitable donor catchment. We design an experiment to diagnose the dominant controls on successful hydrologic model parameter transfer. We calibrate a lumped rainfall–runoff model to 83 stream gauges across the United States. All locations are USGS reference gauges with minimal human influence. Parameter sets from the calibrated models are then transferred to each of the other catchments and the performance of the transferred parameters is assessed. This transfer experiment is carried out both at the scale of the entire US and then for six geographic regions. We use classification and regression tree (CART) analysis to determine the relationship between catchment similarity and performance of transferred parameters. Similarity is defined using physical/climatic catchment characteristics, as well as streamflow response characteristics (signatures such as baseflow index and runoff ratio). Across the entire US, successful parameter transfer is governed by similarity in elevation and climate, and high similarity in streamflow signatures. Controls vary for different geographic regions though. Geology followed by drainage, topography and climate constitute the dominant similarity metrics in forested eastern mountains and plateaus, whereas agricultural land use relates most strongly with successful parameter transfer in the humid plains.

Journal of Hydrology↗

USGS2024

Evaluating seawater intrusion forecast uncertainty under climate change in the Pajaro Valley, California

Climate change and climate variability impacts such as rising sea levels have the potential to exacerbate seawater intrusion and the strain on coastal freshwater resources in already stressed groundwater basins such as those in the Pajaro Valley groundwater basin, California. Regional hydrologic models are often coupled with climate projections to forecast future hydrologic conditions and inform adaptive resources management strategies. However, there is high uncertainty in the future projections of water resources due to uncertainties from downscaling global general circulation models (GCMs) to local scale climate change projections, future land use changes, and the inherent uncertainty of developed hydrologic models. Future climate projections and the magnitude of their influence on modeled hydrologic drivers are highly variable. Therefore, to develop a forecast model, an ensemble of different projections can be used to capture a wider range of basin responses and the associated uncertainties in the modeled forecasts. Understanding the reliability and uncertainty of forecasts is important for developing climate adaptation strategies such as developing protective thresholds, particularly at the basin scale where the impacts are felt, and adaptation is implemented. To demonstrate this, an uncertainty analysis of groundwater level and seawater intrusion forecasts for the Pajaro Valley groundwater basin was performed using an ensemble of three future climate projections with the Pajaro Valley Integrated Hydrologic Model (PVIHM) and the first-order second moment (FOSM) method. FOSM uncertainty analysis of hydrologic forecasts across a multi-GCM climate ensemble provides an upper and lower bound of potential impacts of climate change on sustainability targets related to mitigating seawater intrusion. The groundwater level forecasts’ narrow range of variability can help policymakers in adaptation planning by constraining possible outcomes to a focused range for risk-management decisions. However, less than one-third of groundwater level forecasts met the current protection thresholds to prevent chronic lowering of groundwater. Therefore, sustainability targets may need to be reassessed. Relative to groundwater level changes, the seawater intrusion forecasts had larger uncertainty due to the GCM climate projections and the simulated hydrologic response that were compounded by the propagation of scaling and bias from the GCMs and model simplifications in simulating the coastal boundary.

California↗

USGS2025

Characterization of DOC in riparian wetland pore water and its interference in stable water isotope analysis of δ2H and δ18O

Wetlands play an important role in carbon retention, which is greatly dependent on hydrological conditions. Therefore, the interest in wetland hydrology has increased over recent decades. In wetland hydrology, stable water isotopes ( δ 2 H, δ 18 O) have found increasing use, as laser-based instruments have become readily available and allow for faster analysis and lower costs. However, the use of laser-based determinations of stable hydrogen and oxygen isotopes in wetland pore waters might be challenging because of dissolved organic carbon (DOC), especially as wetland pore waters often contain elevated DOC concentrations. To study the potential interference of DOC on laser-based stable water isotopes determination, we determined concentration, fluorescence and absorbance characteristics of DOC in pore waters collected from six restored and six near-natural riparian wetlands. The water samples’ δ 2 H and δ 18 O values were determined by two laser absorption spectrometers (LAS) and, for reference, in parallel by dual-inlet isotope-ratio mass spectrometry (DI-IRMS). The two LAS methods showed significant deviations from the DI-IRMS-determinations. Variations in the specific UV absorbance index (SUVA 254 , the quotient of the absorbance at 254 nm and the DOC concentration), which indicates the aromaticity of the DOC, partially accounted for interference for δ 2 H measurements. Elevated SUVA 254 index values were especially linked to restored wetlands. Studies that investigate pore waters with high DOC aromaticity but where the DOC is not characterized, and which require a high accuracy of the δ 2 H-values for interpretation of the hydrological system, may be severely affected by DOC interference.

Journal of Hydrology↗

USGS2011

OSL dating of a Pleistocene maar: Birket Ram, the Golan heights

Direct dating of maars and their phreatomagmatic deposits is difficult due to the dominance of lithic (host rock) fragments and glassy particles of the juvenile magma. In this paper we demonstrate that optically stimulated luminescence (OSL) dating can be successfully used for age determination of phreatomagmatic deposits. We studied the tuff deposit of Birket Ram, a basanitic maar located at the northern edge of the Golan heights on the western Arabian plate. The maar is underlain by a thick section of Pleistocene basalts, and currently hosts a small lake. It is filled by approximately 90m of lacustrine sediments with radiocarbon ages extrapolated to 108ka at the base. OSL was applied to quartz grains extracted from tuffs and paleosols in order to set the time frame of the phreatomagmatism at the site. A maximum age constraint of 179??13ka was determined for a paleosol that underlies the maar ejecta. Quartz grains from two layers in the tuff section yielded a direct age of 129??6ka for the phreatomagmatic eruption. A younger age of 104??7ka, which was determined for a tuff layer underlying a basaltic flow, was attributed to thermal resetting during the lava emplacement. This was confirmed by an 40Ar/39Ar age of 101??3ka determined on the overlying basalt. The internal consistency of the OSL ages and the agreement with the 40Ar/39Ar age determination as well as with previous estimates demonstrates the potential of OSL for maar dating. ?? 2010 Elsevier B.V.

Journal of Volcanology and Geothermal Research↗

USGS2010

Satellite observations of a surtseyan eruption: Hunga Ha'apai, Tonga

On March 17, 2009, a surtseyan eruption occurred around Hunga Ha'apai Island, Tonga. A series of observations from the high-spatial resolution Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), and the high-temporal resolution Moderate Resolution Imaging Spectroradiometer (MODIS), was used to estimate the magnitude, location, start time, and duration of the eruption, and measure the evolving characteristics of the new Hunga Ha'apai Island. The eruption start time was estimated to be between 01:50 and 11:10 local time, on March 17, 2009 (i.e., between 12:50 and 22:10 UTC, March 16). The initial explosive phase lasted 3–5 days and consisted of multiple steam and tephra explosions from two distinct vent sources, one on the northwest side, and another about 100 m off the south shore of the pre-existing island. The eruption plume reached 4.0 to 7.6 km altitude above sea level, and tephra added new land around each of the vents, initially tripling the area of the pre-existing island. The next phase of steaming from newly formed crater lakes around the vents lasted a few days. Three warm crater lakes formed initially, but disappeared with time as the shoreline eroded. After ~ 2 months, vegetation that was initially buried by tephra was recovering; after ~ 10 months, the size of the island had eroded down to ~ twice that of the pre-existing island, and the one remaining crater lake had a temperature of ~ 68 ± 17 °C. The volume of erupted material was estimated to be ~ 0.0176 km 3 and the volcanic explosivity index (VEI) was estimated to be VEI = 2.

Hunga Ha'apai Island↗

USGS2016

Seismo-acoustic evidence for an avalanche driven phreatic eruption through a beheaded hydrothermal system: An example from the 2012 Tongariro eruption

The 6 August 2012 Te Maari eruption comprises a complex eruption sequence including multiple eruption pulses, a debris avalanche that propagated ~ 2 km from the vent, and the formation of a 500 m long, arcuate chasm, located ~ 300 m from the main eruption vent. The eruption included 6 distinct impulses that were coherent across a local infrasound network marking the eruption onset at 11:52:18 (all times UTC). An eruption energy release of ~ 3 × 1012 J was calculated using a body wave equation for radiated seismic energy. A similar calculation based on the infrasound record, shows that ~ 90% of the acoustic energy was released from three impulses at onset times 11:52:20 (~ 20% of total eruption energy), 11:52:27 (~ 50%), and 11:52:31 (~ 20%). These energy impulses may coincide with eyewitness accounts describing an initial eastward directed blast, followed by a westward directed blast, and a final vertical blast. Pre-eruption seismic activity includes numerous small unlocatable micro-earthquakes that began at 11:46:50. Two larger high frequency earthquakes were recorded at 11:49:06 and 11:49:21 followed directly by a third earthquake at 11:50:17. The first event was located within the scarp based on an arrival time location from good first P arrival times and probably represents the onset of the debris avalanche. The third event was a tornillo, characterised by a 0.8 Hz single frequency resonance, and has a resonator attenuation factor of Q ~ 40, consistent with a bubbly fluid filled resonator. This contrasts with a similar tornillo event occurring 2.5 weeks earlier having Q ~ 250–1000, consistent with a dusty gas charged resonator. We surmise from pre-eruption seismicity, and the observed attenuation change, that the debris avalanche resulted from the influx of fluids into the hydrothermal system, causing destabilisation and failure. The beheaded hydrothermal system may have then caused depressurisation frothing of the remaining gas charged system leading to the onset of explosive activity.

Journal of Volcanology and Geothermal Research↗

USGS2016

Shallow and deep controls on lava lake surface motion at Kīlauea Volcano

Lava lakes provide a rare window into magmatic behavior, and lake surface motion has been used to infer deeper properties of the magmatic system. At Halema'uma'u Crater, at the summit of Kīlauea Volcano, multidisciplinary observations for the past several years indicate that lava lake surface motion can be broadly divided into two regimes: 1) stable and 2) unstable. Stable behavior is driven by lava upwelling from deeper in the lake (presumably directly from the conduit) and is an intrinsic process that drives lava lake surface motion most of the time. This stable behavior can be interrupted by periods of unstable flow (often reversals) driven by spattering – a shallowly-rooted process often extrinsically triggered by small rockfalls from the crater wall. The bursting bubbles at spatter sources create void spaces and a localized surface depression which draws and consumes surrounding surface crust. Spattering is therefore a location of lava downwelling, not upwelling. Stable (i.e. deep, upwelling-driven) and unstable (i.e. shallow, spattering-driven) behavior often alternate through time, have characteristic surface velocities, flow directions and surface temperature regimes, and also correspond to changes in spattering intensity, outgassing rates, lava level and seismic tremor. These results highlight that several processes, originating at different depths, can control the motion of the lava lake surface, and long-term interdisciplinary monitoring is required to separate these influences. These observations indicate that lake surface motion is not always a reliable proxy for deeper lake or magmatic processes. From these observations, we suggest that shallow outgassing (spattering), not lake convection, drives the variations in lake motion reported at Erta 'Ale lava lake.

Hawaii↗

USGS2017

Monitoring gas and heat emissions at Norris Geyser Basin, Yellowstone National Park, USA based on a combined eddy covariance and Multi-GAS approach

We quantified gas and heat emissions in an acid-sulfate, vapor-dominated area (0.04-km 2 ) of Norris Geyser Basin, located just north of the 0.63 Ma Yellowstone Caldera and near an area of anomalous uplift. From 14 May to 3 October 2016, an eddy covariance system measured half-hourly CO 2 , H 2 O and sensible ( H ) and latent ( LE ) heat fluxes and a Multi-GAS instrument measured (1 Hz frequency) atmospheric H 2 O, CO 2 and H 2 S volumetric mixing ratios. We also measured soil CO 2 fluxes using the accumulation chamber method and temperature profiles on a grid and collected fumarole gas samples for geochemical analysis. Eddy covariance CO 2 fluxes ranged from − 56 to 885 g m − 2 d − 1 . Using wavelet analysis, average daily eddy covariance CO 2 fluxes were locally correlated with average daily environmental parameters on several-day to monthly time scales. Estimates of CO 2 emission rate from the study area ranged from 8.6 t d − 1 based on eddy covariance measurements to 9.8 t d − 1 based on accumulation chamber measurements. Eddy covariance water vapor fluxes ranged from 1178 to 24,600 g m − 2 d − 1 . Nighttime H and LE were considered representative of hydrothermal heat fluxes and ranged from 4 to 183 and 38 to 504 W m − 2 , respectively. The total hydrothermal heat emission rate ( H + LE + radiant) estimated for the study area was 11.6 MW and LE contributed 69% of the output. The mean ± standard deviation of H 2 O, CO 2 and H 2 S mixing ratios measured by the Multi-GAS system were 9.3 ± 3.1 parts per thousand, 467 ± 61 ppmv, and 0.5 ± 0.6 ppmv, respectively, and variations in the gas compositions were strongly correlated with diurnal variations in environmental parameters (wind speed and direction, atmospheric temperature). After removing ambient H 2 O and CO 2 , the observed variations in the Multi-GAS data could be explained by the mixing of relatively H 2 O-CO 2 -H 2 S-rich fumarole gases with CO 2 -rich and H 2 O-H 2 S-poor soil gases. The fumarole H 2 O/CO 2 and CO 2 /H 2 S end member ratios (101.7 and 27.1, respectively, on average) were invariant during the measurement period and fell within the range of values measured in direct fumarole gas samples. The soil gas H 2 O/CO 2 end member ratios (~ 15–30) were variable and low relative to the fumarole end member, likely resulting from water vapor loss during cooling and condensation in the shallow subsurface, whereas the CO 2 /H 2 S end member ratio was high (~ 160), presumably related to transport of CO 2 -dominated soil gas emissions mixed with trace fumarolic emissions to the Multi-GAS station. Nighttime eddy covariance ratios of H 2 O to CO 2 flux were typically between the soil gas and fumarole end member H 2 O/CO 2 ratios defined by Multi-GAS measurements. Overall, the combined eddy covariance and Multi-GAS approach provides a powerful tool for quasi-continuous measurements of gas and heat emissions for improved volcano-hydrothermal monitoring.

Wyoming↗

USGS2019

Eruption and fountaining dynamics of selected 1985–1986 high fountaining episodes at Kīlauea volcano, Hawai'i, from quantitative vesicle microtexture analysis

Tephra from the early Hawaiian fountaining episodes of the ongoing eruption of Pu'u 'Ō'ō in the East Rift Zone (ERZ) of Kīlauea provides an opportunity to study the vesicle microtextures of pyroclasts erupted from a single vent over a prolonged period of time. We report the results of microtextural analysis of pyroclasts from five of Pu'u 'Ō'ō's high (>200 m) Hawaiian fountaining episodes (episodes 32, 37, 40, 44 and 45) erupted during 1985–1986. This analysis was carried out to constrain the parameters that led to large variations in fountain height at Pu'u 'Ō'o, and the extent to which pyroclast residence times in the fountain modified microtextures. Our results confirm the finding of Stovall et al., 2011 , Stovall et al., 2012 that pyroclasts from a single Hawaiian fountain can vary greatly in texture (from bubbly to foamy), and have vesicle volume densities (N m v ) and vesicle to melt ratios (V G /V L ) that vary by an order of magnitude. This range in vesicle texture and population is due to extensive growth and coalescence of vesicles within the fountain after fragmentation. Only one pyroclast from four of five episodes was found to have textures interpreted as indicative of the vesicle population near the moment of fragmentation: bubbly texture, high density (typically >500 kg m −3 ), high N m v (2.2 × 10 6 to 4.4 × 10 6 ), and low V G /V L of 2.06 to 4.65. We demonstrate a linear correlation between Δ(V G /V L ) and peak fountain height across a range of Hawaiian fountains from Kilauea. This correlation could be used to infer peak heights of unobserved Hawaiian fountaining eruptions after further testing using well-recorded events.

Hawai'i↗

USGS1999

Liquid chromatographic determination of benzocaine and N-acetylbenzocaine in the edible fillet tissue from rainbow trout

A method was developed for determining benzocaine and N-acetylbenzocaine concentrations in fillet tissue of rainbow trout. The method involves extracting the analytes with acetonitrile, removing lipids or hydrophobic compounds from the extract with hexane, and providing additional clean-up with solid-phase extraction techniques. Analyte concentrations are determined using reversed-phase high-performance liquid chromatographic techniques with an isocratic mobile phase and UV detection. The accuracy (range, 92 to 121%), precision (R.S.D., <14%), and sensitivity (method quantitation limit, <24 ng/g) for each analyte indicate the usefulness of this method for studies characterizing the depletion of benzocaine residues from fish exposed to benzocaine. Copyright (C) 1999.

Journal of Chromatography A↗

USGS1998

Infiltration and hydraulic connections from the Niagara River to a fractured-dolomite aquifer in Niagara Falls, New York

The spatial distribution of hydrogen and oxygen stable-isotope values in groundwater can be used to distinguish different sources of recharge and to trace groundwater flow directions from recharge boundaries. This method can be particularly useful in fractured-rock settings where multiple lines of evidence are required to delineate preferential flow paths that result from heterogeneity within fracture zones. Flow paths delineated with stable isotopes can be combined with hydraulic data to form a more complete picture of the groundwater flow system. In this study values of ??D and ??18O were used to delineate paths of river-water infiltration into the Lockport Group, a fractured dolomite aquifer, and to compute the percentage of fiver water in groundwater samples from shallow bedrock wells. Flow paths were correlated with areas of high hydraulic diffusivity in the shallow bedrock that were delineated from water-level fluctuations induced by diurnal stage fluctuations in man-made hydraulic structures. Flow paths delineated with the stable-isotope and hydraulic data suggest that fiver infiltration reaches an unlined storm sewer in the bedrock through a drainage system that surrounds aqueducts carrying river water to hydroelectric power plants. This finding is significant because the storm sewer is the discharge point for contaminated groundwater from several chemical waste-disposal sites and the cost of treating the storm sewer's discharge could be reduced if the volume of infiltration from the river were decreased.The spatial distribution of hydrogen and oxygen stable-isotope values in groundwater can be used to distinguish different sources of recharge and to trace groundwater flow directions from recharge boundaries. This method can be particularly useful in fractured-rock settings where multiple lines of evidence are required to delineate preferential flow paths that result from heterogeneity within fracture zones. Flow paths delineated with stable isotopes can be combined with hydraulic data to form a more complete picture of the groundwater flow system. In this study values of ??D and ??18O were used to delineate paths of river-water infiltration into the Lockport Group, a fractured dolomite aquifer, and to compute the percentage of river water in groundwater samples from shallow bedrock wells. Flow paths were correlated with areas of high hydraulic diffusivity in the shallow bedrock that were delineated from water-level fluctuations induced by diurnal stage fluctuations in man-made hydraulic structures. Flow paths delineated with the stable-isotope and hydraulic data suggest that river infiltration reaches an unlined storm sewer in the bedrock through a drainage system that surrounds aqueducts carrying river water to hydroelectric power plants. This finding is significant because the storm sewer is the discharge point for contaminated groundwater from several chemical waste-disposal sites and the cost of treating the storm sewer's discharge could be reduced if the volume of infiltration from the river were decreased.

Journal of Hydrology↗

USGS1998

Analysis and simulation of reactive transport of metal contaminants in ground water in Pinal Creek Basin, Arizona

Large-scale mining activities have generated a plume of acidic ground water more than 15 km long in the regional aquifer of the Pinal Creek Basin. A one-dimensional reactive-transport model was developed using PHREEQC to aid in the analysis of transport and chemical processes in the plume and to determine the uses and limitations of this type of modeling approach. In 1984, the acidic part of the plume had a pH as low as 3.4 and contained milligram-per-liter concentrations of iron, copper, aluminum and other metals. From 1984 to 1994, concentrations of contaminants in the alluvial aquifer in Pinal Creek Basin, Arizona, decreased as a result of mixing, recharge, remedial pumping and chemical reactions. For reactions involving gypsum and rhodochrosite, the equilibrium modeling assumption of a local geochemical equilibrium was generally valid. From 1984 to 1990, water along the simulated flow path was at equilibrium or slightly supersaturated with gypsum, and gypsum equilibria controlled dissolved concentrations of calcium and sulfate. Beginning in 1991, water in the acidic part of the plume became increasingly undersaturated with respect to gypsum, indicating that the gypsum available for dissolution in the aquifer may have been completely consumed by about 1991. Rhodochrosite precipitation was thought responsible for the measured attenuation in dissolved manganese in the neutralized zone. For reactions involving calcite, the assumption of a local geochemical equilibrium was generally not valid. Dissolution of calcite in the transition zone was not sufficient to establish equilibrium although, following neutralization, the calcite saturation index decreased to −1.2 in 1986. Calcite undersaturation decreased along the flow path in the neutralized zone, and equilibrium was attained about 7 km downgradient of the transition zone. The assumption of a local geochemical equilibrium was not valid for oxidation–reduction reactions that involved iron oxides and manganese oxides. Kinetically controlled oxidation–reduction reactions continued in the acidic part of the flow path for years following the passage of the transition zone. Although the equilibrium approach helped to provide an increased understanding of contaminant transport at Pinal Creek, future work will require a kinetic modeling approach to more accurately simulate selected reactions between the plume and aquifer materials.

Arizona↗

USGS1999

The types of data needed for assessing the environmental and human health impacts of coal

Coal is one of the most important sources of energy. Its worldwide use will continue to expand during the next several decades, particularly in rapidly developing countries such as China and India. Unfortunately, coal use may bring with it environmental and human health costs. Many of the environmental and health problems attributed to coal combustion are due to mobilization of potentially toxic elements. Some of these problems could be minimized or even avoided if comprehensive databases containing appropriate coal quality information were available to decision makers so that informed decisions could be made regarding coal use. Among the coal quality parameters that should be included in these databases are: C, H, N, O, pyritic sulfur, organic sulfur, major, minor, and trace element concentrations, modes of occurrence of environmentally sensitive elements, cleanability, mineralogy, organic chemistry, petrography, and leachability.Coal is one of the most important sources of energy. Its worldwide use will continue to expand during the next several decades, particularly in rapidly developing countries such as China and India. Unfortunately, coal use may bring with it environmental and human health costs. Many of the environmental and health problems attributed to coal combustion are due to mobilization of potentially toxic elements. Some of these problems could be minimized or even avoided if comprehensive databases containing appropriate coal quality information were available to decision makers so that informed decisions could be made regarding coal use. Among the coal quality parameters that should be included in these databases are: C, H, N, O, pyritic sulfur, organic sulfur, major, minor, and trace element concentrations, modes of occurrence of environmentally sensitive elements, cleanability, mineralogy, organic chemistry, petrography, and leachability.

International Journal of Coal Geology↗

USGS2001

Progression of natural attenuation processes at a crude-oil spill site . I. Geochemical evolution of the plume

A 16-year study of a hydrocarbon plume shows that the extent of contaminant migration and compound-specific behavior have changed as redox reactions, most notably iron reduction, have progressed over time. Concentration changes at a small scale, determined from analysis of pore-water samples drained from aquifer cores, are compared with concentration changes at the plume scale, determined from analysis of water samples from an observation well network. The small-scale data show clearly that the hydrocarbon plume is growing slowly as sediment iron oxides are depleted. Contaminants, such as ortho -xylene that appeared not to be moving downgradient from the oil on the basis of observation well data, are migrating in thin layers as the aquifer evolves to methanogenic conditions. However, the plume-scale observation well data show that the downgradient extent of the Fe 2+ and BTEX plume did not change between 1992 and 1995. Instead, depletion of the unstable Fe (III) oxides near the subsurface crude-oil source has caused the maximum dissolved iron concentration zone within the plume to spread at a rate of approximately 3 m/year. The zone of maximum concentrations of benzene, toluene, ethylbenzene and xylene (BTEX) has also spread within the anoxic plume. In monitoring the remediation of hydrocarbon-contaminated ground water by natural attenuation, subtle concentration changes in observation well data from the anoxic zone may be diagnostic of depletion of the intrinsic electron-accepting capacity of the aquifer. Recognition of these subtle patterns may allow early prediction of growth of the hydrocarbon plume.

Journal of Contaminant Hydrology↗

USGS2003

The effect of entrapped nonaqueous phase liquids on tracer transport in heterogeneous porous media: Laboratory experiments at the intermediate scale

This work considers the applicability of conservative tracers for detecting high-saturation nonaqueous-phase liquid (NAPL) entrapment in heterogeneous systems. For this purpose, a series of experiments and simulations was performed using a two-dimensional heterogeneous system (10??1.2 m), which represents an intermediate scale between laboratory and field scales. Tracer tests performed prior to injecting the NAPL provide the baseline response of the heterogeneous porous medium. Two NAPL spill experiments were performed and the entrapped-NAPL saturation distribution measured in detail using a gamma-ray attenuation system. Tracer tests following each of the NAPL spills produced breakthrough curves (BTCs) reflecting the impact of entrapped NAPL on conservative transport. To evaluate significance, the impact of NAPL entrapment on the conservative-tracer breakthrough curves was compared to simulated breakthrough curve variability for different realizations of the heterogeneous distribution. Analysis of the results reveals that the NAPL entrapment has a significant impact on the temporal moments of conservative-tracer breakthrough curves. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

USGS1998

Recovery and validation of historical sediment quality data from coastal and estuarine areas: An integrated approach

A comprehensive database of sediment chemistry and environmental parameters has been compiled for Boston Harbor and Massachusetts Bay. This work illustrates methodologies for rescuing and validating sediment data from heterogeneous historical sources. It greatly expands spatial and temporal data coverage of estuarine and coastal sediments. The database contains about 3500 samples containing inorganic chemical, organic, texture and other environmental data dating from 1955 to 1994. Cooperation with local and federal agencies as well as universities was essential in locating and screening documents for the database. More than 80% of references utilized came from sources with limited distribution (gray literature). Task sharing was facilitated by a comprehensive and clearly defined data dictionary for sediments. It also served as a data entry template and flat file format for data processing and as a basis for interpretation and graphical illustration. Standard QA/QC protocols are usually inapplicable to historical sediment data. In this work outliers and data quality problems were identified by batch screening techniques that also provide visualizations of data relationships and geochemical affinities. No data were excluded, but qualifying comments warn users of problem data. For Boston Harbor, the proportion of irreparable or seriously questioned data was remarkably small (<5%), although concentration values for metals and organic contaminants spanned 3 orders of magnitude for many elements or compounds. Data from the historical database provide alternatives to dated cores for measuring changes in surficial sediment contamination level with time. The data indicate that spatial inhomogeneity in harbor environments can be large with respect to sediment-hosted contaminants. Boston Inner Harbor surficial sediments showed decreases in concentrations of Cu, Hg, and Zn of 40 to 60% over a 17-year period.A comprehensive database of sediment chemistry and environmental parameters has been compiled for Boston Harbor and Massachusetts Bay. This work illustrates methodologies for rescuing and validating sediment data from heterogeneous historical sources. It greatly expands spatial and temporal data coverage of estuarine and coastal sediments. The database contains about 3500 samples containing inorganic chemical, organic, texture and other environmental data dating from 1995 to 1994. Cooperation with local and federal agencies as well as universities was essential in locating and screening documents for the database. More than 80% of references utilized came from sources with limited distribution (gray Task sharing was facilitated by a comprehensive and clearly defined data dictionary for sediments. It also served as a data entry template and flat file format for data processing and as a basis for interpretation and graphical illustration. Standard QA/QC protocols are usually inapplicable to historical sediment data. In this work outliers and data quality problems were identified by batch screening techniques that also provide visualizations of data relationships and geochemical affinities. No data were excluded, but qualifying comments warn users of problem data. For Boston Harbor, the proportion of irreparable or seriously questioned data was remarkably small (<5%), although concentration values for metals and organic contaminants spanned 3 orders of magnitude for many elements or compounds. Data from the historical database provide alternatives to dated cores for measuring changes in surficial sediment contamination level with time. The data indicate that spatial inhomogeneity in harbor environments can be large with respect to sediment-hosted contaminants. Boston Inner Harbor surficial sediments showed decreases in concentrations Cu, Hg, and Zn of 40 to 60% over a 17-year period.

Massachusetts↗

USGS2000

10,000 Years of explosive eruptions of Merapi Volcano, Central Java: archaeological and modern implications

Stratigraphy and radiocarbon dating of pyroclastic deposits at Merapi Volcano, Central Java, reveals ~10,000 years of explosive eruptions. Highlights include: (1) Construction of an Old Merapi stratovolcano to the height of the present cone or slightly higher. Our oldest age for an explosive eruption is 9630±60 14 C y B.P.; construction of Old Merapi certainly began earlier. (2) Collapse(s) of Old Merapi that left a somma rim high on its eastern slope and sent one or more debris avalanche(s) down its southern and western flanks. Impoundment of Kali Progo to form an early Lake Borobudur at ~3400 14 C y B.P. hints at a possible early collapse of Merapi. The latest somma-forming collapse occurred ~1900 14 C y B.P. The current cone, New Merapi, began to grow soon thereafter. (3) Several large and many small Buddhist and Hindu temples were constructed in Central Java between 732 and ~900 A.D. (roughly, 1400-1000 14 C y B.P.). Explosive Merapi eruptions occurred before, during and after temple construction. Some temples were destroyed and (or) buried soon after their construction, and we suspect that this destruction contributed to an abrupt shift of power and organized society to East Java in 928 A.D. Other temples sites, though, were occupied by "caretakers" for several centuries longer. (4) A partial collapse of New Merapi occurred <1130±50 14 C y B.P. Eruptions ~700-800 14 C y B.P. (12-14th century A.D.) deposited ash on the floors of (still-occupied?) Candi Sambisari and Candi Kedulan. We speculate but cannot prove that these eruptions were triggered by (the same?) partial collapse of New Merapi, and that the eruptions, in turn, ended "caretaker" occupation at Candi Sambisari and Candi Kedulan. A new or raised Lake Borobudur also existed during part or all of the 12-14th centuries, probably impounded by deposits from Merapi. (5) Relatively benign lava-dome extrusion and dome-collapse pyroclastic flows have dominated activity of the 20th century, but explosive eruptions much larger than any of this century have occurred many times during Merapi's history, most recently during the 19th century. Are the relatively small eruptions of the 20th century a new style of open-vent, less hazardous activity that will persist for the foreseeable future? Or, alternatively, are they merely low-level "background" activity that could be interrupted upon relatively short notice by much larger explosive eruptions? The geologic record suggests the latter, which would place several hundred thousand people at risk. We know of no reliable method to forecast when an explosive eruption will interrupt the present interval of low-level activity. This conclusion has important implications for hazard evaluation.

Journal of Volcanology and Geothermal Research↗