SEARCH · USGS Science

Results for “Geothermics”

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 919 records · Page 51Linked to original sources

USGS2007

Hydrothermal systems and volcano geochemistry

The upward intrusion of magma from deeper to shallower levels beneath volcanoes obviously plays an important role in their surface deformation. This chapter will examine less obvious roles that hydrothermal processes might play in volcanic deformation. Emphasis will be placed on the effect that the transition from brittle to plastic behavior of rocks is likely to have on magma degassing and hydrothermal processes, and on the likely chemical variations in brine and gas compositions that occur as a result of movement of aqueous-rich fluids from plastic into brittle rock at different depths. To a great extent, the model of hydrothermal processes in sub-volcanic systems that is presented here is inferential, based in part on information obtained from deep drilling for geothermal resources, and in part on the study of ore deposits that are thought to have formed in volcanic and shallow plutonic environments.

Book chapter↗

USGS1988

The distribution, structure, and composition of freshwater ice deposits in Bolivian salt lakes

Freshwater ice deposits are described from seven, high elevation (4117-4730 m), shallow (mean depth <30 cm), saline (10-103 g l-1) lakes in the southwestern corner of Bolivia. The ice deposits range to several hundred meters in length and to 7 m in height above the lake or playa surface. They are located near the lake or salar margins; some are completely surrounded by water, others by playa deposits or salt crusts. Upper surfaces and sides of the ice deposits usually are covered by 20-40 cm of white to light brown, dry sedimentary materials. Calcite is the dominant crystalline mineral in these, and amorphous materials such as diatom frustules and volcanic glass are also often abundant. Beneath the dry overburden the ice occurs primarily as horizontal lenses 1-1000 mm thick, irregularly alternating with strata of frozen sedimentary materials. Ice represents from 10 to 87% of the volume of the deposits and yields freshwater (TFR <3 g l-1) when melted. Oxygen isotope ratios for ice are similar to those for regional precipitation and shoreline seeps but much lower than those for the lakewaters. Geothermal flux is high in the region as evidenced by numerous hot springs and deep (3.0-3.5 m) sediment temperatures of 5-10??C. This flux is one cause of the present gradual wasting away of these deposits. Mean annual air temperatures for the different lakes probably are all in the range of -2 to 4??C, and mean midwinter temperatures about 5??C lower. These deposits apparently formed during colder climatic conditions by the freezing of low salinity porewaters and the building up of segregation ice lenses. ?? 1988 Dr W. Junk Publishers.

Hydrobiologia↗

USGS1980

Stable isotope variations in the Quaternary epithermal calcite-fluorite deposit at Monte delle Fate near Cerveteri (Latium, central Italy)

Carbon, oxygen and hydrogen isotope variations have been measured in samples from the epithermal fluorite vein deposit at Monte delle Fate, Latium. The ranges in ?? 13C and ??18O of calcite are -1.3 to 3.4 and 9.5 to 17.3, respectively. ??D values of water extracted from fluid inclusions are -49 to -39 for calcite and -41 to -34 for fluorite. Fluid inclusion filling temperatures (225??-240??C) and salinites (3.75) are nearly the same for both fluorite and sparry calcite. An elongated form of calcite, of minor abundance, precipitated at lower temperatures. The data indicate that (1) the CO2 involved in the mineralization was provided by the local marine limestones, (2) the waters were meteoric in origin and underwent an 18O shift of ??? 10 permil by exchange with marine country rocks, and (3) all geochemical features can be explained by the action of two hydrothermal fluids. Hot brines recently discovered in the Cesano geothermal area, 30 km to the east, have temperatures and some chemical characteristics similar to the hydrothermal fluids at Monte delle Fate. ?? 1980 Springer-Verlag.

Mineralium Deposita↗

USGS2002

Annual emissions of mercury to the atmosphere from natural sources in Nevada and California

The impact of natural source emissions on atmospheric mercury concentrations and the biogeochemical cycle of mercury is not known. To begin to assess this impact, mercury emissions to the atmosphere were scaled up for three areas naturally enriched in mercury: the Steamboat Springs geothermal area, Nevada, the New Idria mercury mining district, California, and the Medicine Lake volcano, California. Data used to scale up area emissions included mercury fluxes, measured in-situ using field flux chambers, from undisturbed and disturbed geologic substrates, and relationships between mercury emissions and geologic rock types, soil mercury concentrations, and surface heat flux. At select locations mercury fluxes were measured for 24 h and the data were used to adjust fluxes measured at different times of the day to give an average daily flux. This adjustment minimized daily temporal variability, which is observed for mercury flux because of light and temperature effects. Area emissions were scaled spatially and temporally with GIS software. Measured fluxes ranged from 0.3 to approximately 50 ng m-2 h-1 at undisturbed sites devoid of mercury mineralization, and to greater than 10,000 ng m-2 h-1 from substrates that were in areas of mercury mining. Area-averaged fluxes calculated for bare soil at Steamboat Springs, New Idria, and Medicine Lake of 181, 9.2, and 2 ng m-2 h-1, respectively, are greater than fluxes previously ascribed to natural non-point sources, indicating that these sources may be more significant contributors of mercury to the atmosphere than previously realized.

Environmental Geology↗

USGS2006

Evidence of regional subsidence and associated interior wetland loss induced by hydrocarbon production, Gulf Coast region, USA

Analysis of remote images, elevation surveys, stratigraphic cross-sections, and hydrocarbon production data demonstrates that extensive areas of wetland loss in the northern Gulf Coast region of the United States were associated with large-volume fluid production from mature petroleum fields. Interior wetland losses at many sites in coastal Louisiana and Texas are attributed largely to accelerated land subsidence and fault reactivation induced by decreased reservoir pressures as a result of rapid or prolonged extraction of gas, oil, and associated brines. Evidence that moderately-deep hydrocarbon production has induced land-surface subsidence and reactivated faults that intersect the surface include: (1) close temporal and spatial correlation of fluid production with surficial changes including rapid subsidence of wetland sediments near producing fields, (2) measurable offsets of shallow strata across the zones of wetland loss, (3) large reductions in subsurface pressures where subsidence rates are high, (4) coincidence of orientation and direction of displacement between surface fault traces and faults that bound the reservoirs, and (5) accelerated subsidence rates near producing fields compared to subsidence rates in surrounding areas or compared to geological rates of subsidence. Based on historical trends, subsidence rates in the Gulf Coast region near producing fields most likely will decrease in the future because most petroleum fields are nearly depleted. Alternatively, continued extraction of conventional energy resources as well as potential production of alternative energy resources (geopressured-geothermal fluids) in the Gulf Coast region could increase subsidence and land losses and also contribute to inundation of areas of higher elevation. ?? Springer-Verlag 2006.

Environmental Geology↗

USGS1999

Gas hydrate accumulation at the Hakon Mosby Mud Volcano

Gas hydrate (GH) accumulation is characterized and modeled for the Hakon Mosby mud volcano, ca. 1.5 km across, located on the Norway-Barents-Svalbard margin. Pore water chemical and isotopic results based on shallow sediment cores as well as geothermal and geomorphological data suggest that the GH accumulation is of a concentric pattern controlled by and formed essentially from the ascending mud volcano fluid. The gas hydrate content of sediment peaks at 25% by volume, averaging about 1.2% throughout the accumulation. The amount of hydrate methane is estimated at ca. 108 m3 STP, which could account for about 1-10% of the gas that has escaped from the volcano since its origin.

Geo-Marine Letters↗

USGS2009

Accessory mineral U–Th–Pb ages and 40Ar/39Ar eruption chronology, and their bearing on rhyolitic magma evolution in the Pleistocene Coso volcanic field, California

We determined Ar/Ar eruption ages of eight extrusions from the Pleistocene Coso volcanic field, a long-lived series of small volume rhyolitic domes in eastern California. Combined with ion-microprobe dating of crystal ages of zircon and allanite from these lavas and from granophyre geothermal well cuttings, we were able to track the range of magma-production rates over the past 650 ka at Coso. In ≤230 ka rhyolites we find no evidence of protracted magma residence or recycled zircon (or allanite) from Pleistocene predecessors. A significant subset of zircon in the ~85 ka rhyolites yielded ages between ~100 and 200 Ma, requiring that generation of at least some rhyolites involves material from Mesozoic basement. Similar zircon xenocrysts are found in an ~200 ka granophyre. The new age constraints imply that magma evolution at Coso can occur rapidly as demonstrated by significant changes in rhyolite composition over short time intervals (≤10’s to 100’s ka). In conjunction with radioisotopic age constraints from other young silicic volcanic fields, dating of Coso rhyolites highlights the fact that at least some (and often the more voluminous) rhyolites are produced relatively rapidly, but that many small-volume rhyolites likely represent separation from long-lived mushy magma bodies.

California↗

USGS2005

Seismic signature of a phreatic explosion: Hydrofracturing damage at Karthala volcano, Grande Comore Island, Indian Ocean

Karthala volcano is a basaltic shield volcano with an active hydrothermal system that forms the southern two-thirds of the Grande Comore Island, off the east coat of Africa, northwest of Madagascar. Since the start of volcano monitoring by the local volcano observatory in 1988, the July 11th, 1991 phreatic eruption was the first volcanic event seismically recorded on this volcano, and a rare example of a monitored basaltic shield. From 1991 to 1995 the VT locations, 0.5<Ml<4.3, show a crack shaped pattern (3 km long, 1 km wide) within the summit caldera extending at depth from -2 km to +2 km relative to sea level. This N-S elongated pattern coincides with the direction of the regional maximum horizontal stress as deduced from regional focal mechanism solutions. This brittle signature of the damage associated with the 1991 phreatic eruption is a typical pattern of the seismicity induced by controlled fluid injections such as those applied at geothermal fields, in oil and gas recovery, or for stress measurements. It suggests the 1991 phreatic eruption was driven by hydraulic fracturing induced by forced fluid flow. We propose that the extremely high LP and VT seismicity rates, relative to other effusive volcanoes, during the climax of the 1991 phreatic explosion, are due to the activation of the whole hydrothermal system, as roughly sized by the distribution of VT hypocenters. The seismicity rate in 1995 was still higher than the pre-eruption seismicity rate, and disagrees with the time pattern of thermo-elastic stress readjustment induced by single magma intrusions at basaltic volcanoes. We propose that it corresponds to the still ongoing relaxation of pressure heterogeneity within the hydrothermal system as suggested by the few LP events that still occurred in 1995. ?? Springer-Verlag 2005.

Bulletin of Volcanology↗

USGS2012

Recent explosive eruptions and volcano hazards at Soputan volcano—a basalt stratovolcano in north Sulawesi, Indonesia

Soputan is a high-alumina basalt stratovolcano located in the active North Sulawesi-Sangihe Islands magmatic arc. Although immediately adjacent to the still geothermally active Quaternary Tondono Caldera, Soputan’s magmas are geochemically distinct from those of the caldera and from other magmas in the arc. Unusual for a basalt volcano, Soputan produces summit lava domes and explosive eruptions with high-altitude ash plumes and pyroclastic flows—eight explosive eruptions during the period 2003–2011. Our field observations, remote sensing, gas emission, seismic, and petrologic analyses indicate that Soputan is an open-vent-type volcano that taps basalt magma derived from the arc-mantle wedge, accumulated and fractionated in a deep-crustal reservoir and transported slowly or staged at shallow levels prior to eruption. A combination of high phenocryst content, extensive microlite crystallization and separation of a gas phase at shallow levels results in a highly viscous basalt magma and explosive eruptive style. The open-vent structure and frequent eruptions indicate that Soputan will likely erupt again in the next decade, perhaps repeatedly. Explosive eruptions in the Volcano Explosivity Index (VEI) 2–3 range and lava dome growth are most probable, with a small chance of larger VEI 4 eruptions. A rapid ramp up in seismicity preceding the recent eruptions suggests that future eruptions may have no more than a few days of seismic warning. Risk to population in the region is currently greatest for villages located on the southern and western flanks of the volcano where flow deposits are directed by topography. In addition, Soputan’s explosive eruptions produce high-altitude ash clouds that pose a risk to air traffic in the region.

Bulletin of Volcanology↗

USGS2021

Near-real-time volcanic cloud monitoring: Insights into global explosive volcanic eruptive activity through analysis of Volcanic Ash Advisories

Understanding the location, intensity, and likely duration of volcanic hazards is key to reducing risk from volcanic eruptions. Here, we use a novel near-real-time dataset comprising Volcanic Ash Advisories (VAAs) issued over 10 years to investigate global rates and durations of explosive volcanic activity. The VAAs were collected from the nine Volcanic Ash Advisory Centres (VAACs) worldwide. Information extracted allowed analysis of the frequency and type of explosive behaviour, including analysis of key eruption source parameters (ESPs) such as volcanic cloud height and duration. The results reflect changes in the VAA reporting process, data sources, and volcanic activity through time. The data show an increase in the number of VAAs issued since 2015 that cannot be directly correlated to an increase in volcanic activity. Instead, many represent increased observations, including improved capability to detect low- to mid-level volcanic clouds (FL101–FL200, 3–6 km asl), by higher temporal, spatial, and spectral resolution satellite sensors. Comparison of ESP data extracted from the VAAs with the Mastin et al. (J Volcanol Geotherm Res 186:10–21, 2009a ) database shows that traditional assumptions used in the classification of volcanoes could be much simplified for operational use. The analysis highlights the VAA data as an exceptional resource documenting global volcanic activity on timescales that complement more widely used eruption datasets.

Bulletin of Volcanology↗

USGS1998

Seismic image of the Mount Spurr magmatic system

The three-dimenisonal P-wave velocity structure of Moutn Spurr is determined to depths of 10km by tomographic inversion of 3754 first-arriving P-wave times from local earthquakes recorded by a permanent network of 11 seismographs. Results show a prominent low-velocity zone extending from the surface to 3-4km below sea level beneath the southeastern flank of Crater Peak, spatially coincident with a geothermal system. P-wave velocities in this low-velocity zone are approximately 20% slower than those in the shallow crystalline basement rocks. Beneath Crater Peak an approximately 3km-wide zone of relative low velocities correlates with a near-vertical band of seismicity, suggestive of a magmatic conduit. No large low-velocity zone indicative of a magma chamber occurs within the upper 10km of the crust. These observations are consistent with petrologic and geochemical studies suggesting that Crater Peak magmas originate in the lower crust or upper mantle and have s short residence time in the shallow crust. Earthquakes relocaetd using the three-dimensional velocity structure correlate well with surface geology and other geophysical observations; thus, they provide additional constraints on the kinematics of the Mount Spurr magmatic system.

Bulletin of Volcanology↗

USGS2009

Oxygen, hydrogen, and helium isotopes for investigating groundwater systems of the Cape Verde Islands, West Africa

Stable isotopes (?? 18 O, ?? 2 H), tritium ( 3 H), and helium isotopes ( 3 He, 4 He) were used for evaluating groundwater recharge sources, flow paths, and residence times of three watersheds in the Cape Verde Islands (West Africa). Stable isotopes indicate the predominance of high-elevation precipitation that undergoes little evaporation prior to groundwater recharge. In contrast to other active oceanic hotspots, environmental tracers show that deep geothermal circulation does not strongly affect groundwater. Low tritium concentrations at seven groundwater sites indicate groundwater residence times of more than 50 years. Higher tritium values at other sites suggest some recent recharge. High 4 He and 3 He/ 4 He ratios precluded 3 H/ 3 He dating at six sites. These high 3 He/ 4 He ratios (R/R a values of up to 8.3) are consistent with reported mantle derived helium of oceanic island basalts in Cape Verde and provided end-member constraints for improved dating at seven other locations. Tritium and 3 H/ 3 He dating shows that S??o Nicolau Island's Ribeira Faj?? Basin has groundwater residence times of more than 50 years, whereas Fogo Island's Mosteiros Basin and Santo Ant??o Island's Ribeira Paul Basin contain a mixture of young and old groundwater. Young ages at selected sites within these two basins indicate local recharge and potential groundwater susceptibility to surface contamination and/or salt-water intrusion. ?? Springer-Verlag 2009.

Hydrogeology Journal↗

USGS2014

Assessing the magnitude and timing of anthropogenic warming of a shallow aquifer: example from Virginia Beach, USA

Groundwater temperature measurements in a shallow coastal aquifer in Virginia Beach, Virginia, USA, suggest groundwater warming of +4.1 °C relative to deeper geothermal gradients. Observed warming is related to timing and depth of influence of two potential thermal drivers—atmospheric temperature increases and urbanization. Results indicate that up to 30 % of groundwater warming at the water table can be attributed to atmospheric warming while up to 70 % of warming can be attributed to urbanization. Groundwater temperature readings to 30-m depth correlate positively with percentage of impervious cover and negatively with percentage of tree canopy cover; thus, these two land-use metrics explain up to 70 % of warming at the water table. Analytical and numerical modeling results indicate that an average vertical groundwater temperature profile for the study area, constructed from repeat measurement at 11 locations over 15 months, is consistent with the timing of land-use change over the past century in Virginia Beach. The magnitude of human-induced warming at the water table (+4.1 °C) is twice the current seasonal temperature variation, indicating the potential for ecological impacts on wetlands and estuaries receiving groundwater discharge from shallow aquifers.

Virginia↗

USGS2007

Subaqueous geology and a filling model for Crater Lake, Oregon

Results of a detailed bathymetric survey of Crater Lake conducted in 2000, combined with previous results of submersible and dredge sampling, form the basis for a geologic map of the lake floor and a model for the filling of Crater Lake with water. The most prominent landforms beneath the surface of Crater Lake are andesite volcanoes that were active as the lake was filling with water, following caldera collapse during the climactic eruption of Mount Mazama 7700 cal. yr B.P. The Wizard Island volcano is the largest and probably was active longest, ceasing eruptions when the lake was 80 m lower than present. East of Wizard Island is the central platform volcano and related lava flow fields on the caldera floor. Merriam Cone is a symmetrical andesitic volcano that apparently was constructed subaqueously during the same period as the Wizard Island and central platform volcanoes. The youngest postcaldera volcanic feature is a small rhyodacite dome on the east flank of the Wizard Island edifice that dates from 4800 cal. yr B.P. The bathymetry also yields information on bedrock outcrops and talus/debris slopes of the caldera walls. Gravity flows transport sediment from wall sources to the deep basins of the lake. Several debris-avalanche deposits, containing blocks up to 280 m long, are present on the caldera floor and occur below major embayments in the caldera walls. Geothermal phenomena on the lake floor are bacterial mats, pools of solute-rich warm water, and fossil subaqueous hot spring deposits. Lake level is maintained by a balance between precipitation and inflow versus evaporation and leakage. High-resolution bathymetry reveals a series of up to nine drowned beaches in the upper 30 m of the lake that we propose reflect stillstands subsequent to filling of Crater Lake. A prominent wave-cut platform between 4 m depth and present lake level that commonly is up to 40 m wide suggests that the surface of Crater Lake has been at this elevation for a very long time. Lake level apparently is limited by leakage through a permeable layer in the northeast caldera wall. The deepest drowned beach approximately corresponds to the base of the permeable layer. Among a group of lake filling models, our preferred one is constrained by the drowned beaches, the permeable layer in the caldera wall, and paleoclimatic data. We used a precipitation rate 70% of modern as a limiting case. Satisfactory models require leakage to be proportional to elevation and the best fit model has a linear combination of 45% leakage proportional to elevation and 55% of leakage proportional to elevation above the base of the permeable layer. At modern precipitation rates, the lake would have taken 420 yr to fill, or a maximum of 740 yr if precipitation was 70% of the modern value. The filling model provides a chronology for prehistoric passage zones on postcaldera volcanoes that ceased erupting before the lake was filled.

Hydrobiologia↗

USGS2011

Kansas Energy Sources: A Geological Review

Kansas produces both conventional energy (oil, gas, and coal) and nonconventional (coalbed gas, wind, hydropower, nuclear, geothermal, solar, and biofuels) and ranks the 22nd in state energy production in the U.S. Nonrenewable conventional petroleum is the most important energy source with nonrenewable, nonconventional coalbed methane gas becoming increasingly important. Many stratigraphic units produce oil and/or gas somewhere in the state with the exception of the Salina Basin in north-central Kansas. Coalbed methane is produced from shallow wells drilled into the thin coal units in southeastern Kansas. At present, only two surface coal mines are active in southeastern Kansas. Although Kansas has been a major exporter of energy in the past (it ranked first in oil production in 1916), now, it is an energy importer.

Kansas↗

USGS2002

Study on 3-D velocity structure of crust and upper mantle in Sichuan-yunnan region, China

Based on the first arrival P and S data of 4 625 regional earthquakes recorded at 174 stations dispersed in the Yunnan and Sichuan Provinces, the 3-D velocity structure of crust and upper mantle in the region is determined, incorporating with previous deep geophysical data. In the upper crust, a positive anomaly velocity zone exists in the Sichuan basin, whereas a negative anomaly velocity zone exists in the western Sichuan plateau. The boundary between the positive and negative anomaly zones is the Longmenshan fault zone. The images of lower crust and upper mantle in the Longmenshan fault, Xianshuihe fault, Honghe fault and others appear the characteristic of tectonic boundary, indicating that the faults litely penetrate the Moho discontinuity. The negative velocity anomalies at the depth of 50 km in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with the temperature and composition variations in the upper mantle. The overall features of the crustal and the upper mantle structures in the Sichuan-Yunnan region are the lower average velocity in both crust and uppermost mantle, the large crustal thickness variations, and the existence of high conductivity layer in the crust or/and upper mantle, and higher geothermal value. All these features are closely related to the collision between the Indian and the Asian plates. The crustal velocity in the Sichuan-Yunnan rhombic block generally shows normal.value or positive anomaly, while the negative anomaly exists in the area along the large strike-slip faults as the block boundary. It is conducive to the crustal block side-pressing out along the faults. In the major seismic zones, the seismicity is relative to the negative anomaly velocity. Most strong earthquakes occurred in the upper-mid crust with positive anomaly or normal velocity, where the negative anomaly zone generally exists below.

Acta Seismologica Sinica↗

USGS2015

Geochemical evolution of groundwater in the Mud Lake area, eastern Idaho, USA

Groundwater with elevated dissolved-solids concentrations—containing large concentrations of chloride, sodium, sulfate, and calcium—is present in the Mud Lake area of Eastern Idaho. The source of these solutes is unknown; however, an understanding of the geochemical sources and processes controlling their presence in groundwater in the Mud Lake area is needed to better understand the geochemical sources and processes controlling the water quality of groundwater at the Idaho National Laboratory. The geochemical sources and processes controlling the water quality of groundwater in the Mud Lake area were determined by investigating the geology, hydrology, land use, and groundwater geochemistry in the Mud Lake area, proposing sources for solutes, and testing the proposed sources through geochemical modeling with PHREEQC. Modeling indicated that sources of water to the eastern Snake River Plain aquifer were groundwater from the Beaverhead Mountains and the Camas Creek drainage basin; surface water from Medicine Lodge and Camas Creeks, Mud Lake, and irrigation water; and upward flow of geothermal water from beneath the aquifer. Mixing of groundwater with surface water or other groundwater occurred throughout the aquifer. Carbonate reactions, silicate weathering, and dissolution of evaporite minerals and fertilizer explain most of the changes in chemistry in the aquifer. Redox reactions, cation exchange, and evaporation were locally important. The source of large concentrations of chloride, sodium, sulfate, and calcium was evaporite deposits in the unsaturated zone associated with Pleistocene Lake Terreton. Large amounts of chloride, sodium, sulfate, and calcium are added to groundwater from irrigation water infiltrating through lake bed sediments containing evaporite deposits and the resultant dissolution of gypsum, halite, sylvite, and bischofite.

Idaho↗

USGS1991

A statistical approach to the interpretation of aliphatic hydrocarbon distributions in marine sediments

Q-mode factor analysis was used to quantitate the distribution of the major aliphatic hydrocarbon (n-alkanes, pristane, phytane) systems in sediments from a variety of marine environments. The compositions of the pure end members of the systems were obtained from factor scores and the distribution of the systems within each sample was obtained from factor loadings. All the data, from the diverse environments sampled (estuarine (San Francisco Bay), fresh-water (San Francisco Peninsula), polar-marine (Antarctica) and geothermal-marine (Gorda Ridge) sediments), were reduced to three major systems: a terrestrial system (mostly high molecular weight aliphatics with odd-numbered-carbon predominance), a mature system (mostly low molecular weight aliphatics without predominance) and a system containing mostly high molecular weight aliphatics with even-numbered-carbon predominance. With this statistical approach, it is possible to assign the percentage contribution from various sources to the observed distribution of aliphatic hydrocarbons in each sediment sample. ?? 1991.

Chemical Geology↗