Geology Reports⌕ Search

Geology topics

T. S. Collett

Publications and source records attributed to T. S. Collett.

At least 19 recordsLinked to original sources

Anisotropic models to account for large borehole washouts to estimate gas hydrate saturations in the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II Alaminos 21 B well

Through the use of 3-D seismic amplitude mapping, several gashydrate prospects were identified in the Alaminos Canyon (AC) area of the Gulf of Mexico. Two locations were drilled as part of the Gulf of MexicoGasHydrate Joint Industry Project Leg II (JIP Leg II) in May of 2009 and a comprehensive set of logging-while-drilling (LWD) logs were acquired at each well site. LWD logs indicated that resistivity in the range of ~2 ohm-m and P-wave velocity in the range of ~1.9 km/s were measured in the target sand interval between 515 and 645 feet below sea floor. These values were slightly elevated relative to those measured in the sediment above and below the target sand. However, the initial well log analysis was inconclusive regarding the presence of gashydrate in the logged sand interval, mainly because largewashouts caused by drilling in the target interval degraded confidence in the well log measurements. To assess gashydratesaturations in the sedimentary section drilled in the Alaminos Canyon 21B (AC21-B) well, a method of compensating for the effect of washouts on the resistivity and acoustic velocities was developed. The proposed method models the washed-out portion of the borehole as a vertical layer filled with sea water (drilling fluid) and the apparent anisotropic resistivity and velocities caused by a vertical layer are used to correct the measured log values. By incorporating the conventional marine seismic data into the well log analysis, the average gashydratesaturation in the target sand section in the AC21-Bwell can be constrained to the range of 8–28%, with 20% being our best estimate.

Gulf Of Mexico;Alaminos Canyon↗

Pore- and fracture-filling gas hydrate reservoirs in the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II Green Canyon 955 H well

High-quality logging-while-drilling (LWD) downhole logs were acquired in seven wells drilled during the Gulf of MexicoGasHydrateJointIndustryProjectLegII in the spring of 2009. Well logs obtained in one of the wells, the GreenCanyon Block 955Hwell (GC955-H), indicate that a 27.4-m thick zone at the depth of 428 m below sea floor (mbsf; 1404 feet below sea floor (fbsf)) contains gashydrate within sand with average gashydrate saturations estimated at 60% from the compressional-wave (P-wave) velocity and 65% (locally more than 80%) from resistivity logs if the gashydrate is assumed to be uniformly distributed in this mostly sand-rich section. Similar analysis, however, of log data from a shallow clay-rich interval between 183 and 366 mbsf (600 and 1200 fbsf) yielded average gashydrate saturations of about 20% from the resistivity log (locally 50-60%) and negligible amounts of gashydrate from the P-wave velocity logs. Differences in saturations estimated between resistivity and P-wave velocities within the upper clay-rich interval are caused by the nature of the gashydrate occurrences. In the case of the shallow clay-rich interval, gashydrate fills vertical (or high angle) fractures in rather than fillingpore space in sands. In this study, isotropic and anisotropic resistivity and velocity models are used to analyze the occurrence of gashydrate within both the clay-rich and sand dominated gas-hydrate-bearing reservoirs in the GC955-Hwell.

Gulf Of Mexico↗

Examination of core samples from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Effects of retrieval and preservation

Collecting and preserving undamaged core samples containing gas hydrates from depth is difficult because of the pressure and temperature changes encountered upon retrieval. Hydrate-bearing core samples were collected at the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well in February 2007. Coring was performed while using a custom oil-based drilling mud, and the cores were retrieved by a wireline. The samples were characterized and subsampled at the surface under ambient winter arctic conditions. Samples thought to be hydrate bearing were preserved either by immersion in liquid nitrogen (LN), or by storage under methane pressure at ambient arctic conditions, and later depressurized and immersed in LN. Eleven core samples from hydrate-bearing zones were scanned using x-ray computed tomography to examine core structure and homogeneity. Features observed include radial fractures, spalling-type fractures, and reduced density near the periphery. These features were induced during sample collection, handling, and preservation. Isotopic analysis of the methane from hydrate in an initially LN-preserved core and a pressure-preserved core indicate that secondary hydrate formation occurred throughout the pressurized core, whereas none occurred in the LN-preserved core, however no hydrate was found near the periphery of the LN-preserved core. To replicate some aspects of the preservation methods, natural and laboratory-made saturated porous media samples were frozen in a variety of ways, with radial fractures observed in some LN-frozen sands, and needle-like ice crystals forming in slowly frozen clay-rich sediments. Suggestions for hydrate-bearing core preservation are presented.

Alaska↗

Gas geochemistry of the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: implications for gas hydrate exploration in the Arctic

Gases were analyzed from well cuttings, core, gas hydrate, and formation tests at the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well, drilled within the Milne Point Unit, Alaska North Slope. The well penetrated a portion of the Eileen gas hydrate deposit, which overlies the more deeply buried Prudhoe Bay, Milne Point, West Sak, and Kuparuk River oil fields. Gas sources in the upper 200 m are predominantly from microbial sources (C 1 isotopic compositions ranging from −86.4 to −80.6‰). The C 1 isotopic composition becomes progressively enriched from 200 m to the top of the gas hydrate-bearing sands at 600 m. The tested gas hydrates occur in two primary intervals, units D and C, between 614.0 m and 664.7 m, containing a total of 29.3 m of gas hydrate-bearing sands. The hydrocarbon gases in cuttings and core samples from 604 to 914 m are composed of methane with very little ethane. The isotopic composition of the methane carbon ranges from −50.1 to −43.9‰ with several outliers, generally decreasing with depth. Gas samples collected by the Modular Formation Dynamics Testing (MDT) tool in the hydrate-bearing units were similarly composed mainly of methane, with up to 284 ppm ethane. The methane isotopic composition ranged from −48.2 to −48.0‰ in the C sand and from −48.4 to −46.6‰ in the D sand. Methane hydrogen isotopic composition ranged from −238 to −230‰, with slightly more depleted values in the deeper C sand. These results are consistent with the concept that the Eileen gas hydrates contain a mixture of deep-sourced, microbially biodegraded thermogenic gas, with lesser amounts of thermogenic oil-associated gas, and coal gas. Thermal gases are likely sourced from existing oil and gas accumulations that have migrated up-dip and/or up-fault and formed gas hydrate in response to climate cooling with permafrost formation.

Alaska↗

Seismic imaging of a fractured gas hydrate system in the Krishna-Godavari Basin offshore India

Gas hydrate was discovered in the Krishna-Godavari (KG) Basin during the India National Gas Hydrate Program (NGHP) Expedition 1 at Site NGHP-01-10 within a fractured clay-dominated sedimentary system. Logging-while-drilling (LWD), coring, and wire-line logging confirmed gas hydrate dominantly in fractures at four borehole sites spanning a 500m transect. Three-dimensional (3D) seismic data were subsequently used to image the fractured system and explain the occurrence of gas hydrate associated with the fractures. A system of two fault-sets was identified, part of a typical passive margin tectonic setting. The LWD-derived fracture network at Hole NGHP-01-10A is to some extent seen in the seismic data and was mapped using seismic coherency attributes. The fractured system around Site NGHP-01-10 extends over a triangular-shaped area of ~2.5 km2 defined using seismic attributes of the seafloor reflection, as well as " seismic sweetness" at the base of the gas hydrate occurrence zone. The triangular shaped area is also showing a polygonal (nearly hexagonal) fault pattern, distinct from other more rectangular fault patterns observed in the study area. The occurrence of gas hydrate at Site NGHP-01-10 is the result of a specific combination of tectonic fault orientations and the abundance of free gas migration from a deeper gas source. The triangular-shaped area of enriched gas hydrate occurrence is bound by two faults acting as migration conduits. Additionally, the fault-associated sediment deformation provides a possible migration pathway for the free gas from the deeper gas source into the gas hydrate stability zone. It is proposed that there are additional locations in the KG Basin with possible gas hydrate accumulation of similar tectonic conditions, and one such location was identified from the 3D seismic data ~6 km NW of Site NGHP-01-10. ?? 2010.

Marine and Petroleum Geology↗

Gas hydrate saturations estimated from fractured reservoir at Site NGHP-01-10, Krishna-Godavari Basin, India

During the Indian National Gas Hydrate Program Expedition 01 (NGHP-Ol), one of the richest marine gas hydrate accumulations was discovered at Site NGHP-01-10 in the Krishna-Godavari Basin. The occurrence of concentrated gas hydrate at this site is primarily controlled by the presence of fractures. Assuming the resistivity of gas hydratebearing sediments is isotropic, th?? conventional Archie analysis using the logging while drilling resistivity log yields gas hydrate saturations greater than 50% (as high as ???80%) of the pore space for the depth interval between ???25 and ???160 m below seafloor. On the other hand, gas hydrate saturations estimated from pressure cores from nearby wells were less than ???26% of the pore space. Although intrasite variability may contribute to the difference, the primary cause of the saturation difference is attributed to the anisotropic nature of the reservoir due to gas hydrate in high-angle fractures. Archie's law can be used to estimate gas hydrate saturations in anisotropic reservoir, with additional information such as elastic velocities to constrain Archie cementation parameters m and the saturation exponent n. Theory indicates that m and n depend on the direction of the measurement relative to fracture orientation, as well as depending on gas hydrate saturation. By using higher values of m and n in the resistivity analysis for fractured reservoirs, the difference between saturation estimates is significantly reduced, although a sizable difference remains. To better understand the nature of fractured reservoirs, wireline P and S wave velocities were also incorporated into the analysis.

Journal of Geophysical Research B: Solid Earth↗

Preliminary report on the commercial viability of gas production from natural gas hydrates

Economic studies on simulated gas hydrate reservoirs have been compiled to estimate the price of natural gas that may lead to economically viable production from the most promising gas hydrate accumulations. As a first estimate, $CDN 2005 12/Mscf is the lowest gas price that would allow economically viable production from gas hydrates in the absence of associated free gas, while an underlying gas deposit will reduce the viability price estimate to $CDN 2005 7.50/Mscf. Results from a recent analysis of the simulated production of natural gas from marine hydrate deposits are also considered in this report; on an IROR basis, it is $US 2008 3.50-4.00/Mscf more expensive to produce marine hydrates than conventional marine gas assuming the existence of sufficiently large marine hydrate accumulations. While these prices represent the best available estimates, the economic evaluation of a specific project is highly dependent on the producibility of the target zone, the amount of gas in place, the associated geologic and depositional environment, existing pipeline infrastructure, and local tariffs and taxes. ?? 2009 Elsevier B.V.

Energy Economics↗

Unique problems associated with seismic analysis of partially gas-saturated unconsolidated sediments

Gas hydrate stability conditions restrict the occurrence of gas hydrate to unconsolidated and high water-content sediments at shallow depths. Because of these host sediments properties, seismic and well log data acquired for the detection of free gas and associated gas hydrate-bearing sediments often require nonconventional analysis. For example, a conventional method of identifying free gas using the compressional/shear-wave velocity (V p /V s ) ratio at the logging frequency will not work, unless the free-gas saturations are more than about 40%. The P-wave velocity dispersion of partially gas-saturated sediments causes a problem in interpreting well log velocities and seismic data. Using the White, J.E. [1975. Computed seismic speeds and attenuation in rocks with partial gas saturation. Geophysics 40, 224-232] model for partially gas-saturated sediments, the difference between well log and seismic velocities can be reconciled. The inclusion of P-wave velocity dispersion in interpreting well log data is, therefore, essential to identify free gas and to tie surface seismic data to synthetic seismograms.

Marine and Petroleum Geology↗

Toward production from gas hydrates: Current status, assessment of resources, and simulation-based evaluation of technology and potential

Gas hydrates (GHs) are a vast energy resource with global distribution in the permafrost and in the oceans. Even if conservative estimates are considered and only a small fraction is recoverable, the sheer size of the resource is so large that it demands evaluation as a potential energy source. In this review paper, we discuss the distribution of natural GH accumulations, the status of the primary international research and development (R&D) programs, and the remaining science and technological challenges facing the commercialization of production. After a brief examination of GH accumulations that are well characterized and appear to be models for future development and gas production, we analyze the role of numerical simulation in the assessment of the hydrate-production potential, identify the data needs for reliable predictions, evaluate the status of knowledge with regard to these needs, discuss knowledge gaps and their impact, and reach the conclusion that the numerical-simulation capabilities are quite advanced and that the related gaps either are not significant or are being addressed. We review the current body of literature relevant to potential productivity from different types of GH deposits and determine that there are consistent indications of a large production potential at high rates across long periods from a wide variety of hydrate deposits. Finally, we identify (a) features, conditions, geology and techniques that are desirable in potential production targets; (b) methods to maximize production; and (c) some of the conditions and characteristics that render certain GH deposits undesirable for production. Copyright ?? 2009 Society of Petroleum Engineers.

Conference Paper↗

Integrated analysis of well logs and seismic data to estimate gas hydrate concentrations at Keathley Canyon, Gulf of Mexico

Accurately detecting and quantifying gas hydrate or free gas in sediments from seismic data require downhole well-log data to calibrate the physical properties of the gas hydrate-/free gas-bearing sediments. As part of the Gulf of Mexico Joint Industry Program, a series of wells were either cored or drilled in the Gulf of Mexico to characterize the physical properties of gas hydrate-bearing sediments, to calibrate geophysical estimates, and to evaluate source and transport mechanisms for gas within the gas hydrates. Downhole acoustic logs were used sparingly in this study because of degraded log quality due to adverse wellbore conditions. However, reliable logging while drilling (LWD) electrical resistivity and porosity logs were obtained. To tie the well-log information to the available 3-D seismic data in this area, a velocity log was calculated from the available resistivity log at the Keathley Canyon 151-2 well, because the acoustic log or vertical seismic data acquired at the nearby Keathley Canyon 151-3 well were either of poor quality or had limited depth coverage. Based on the gas hydrate saturations estimated from the LWD resistivity log, the modified Biot-Gassmann theory was used to generate synthetic acoustic log and a synthetic seismogram was generated with a fairly good agreement with a seismic profile crossing the well site. Based on the well-log information, a faintly defined bottom-simulating reflection (BSR) in this area is interpreted as a reflection representing gas hydrate-bearing sediments with about 15% saturation overlying partially gas-saturated sediments with 3% saturation. Gas hydrate saturations over 30-40% are estimated from the resistivity log in two distinct intervals at 220-230 and 264-300 m below the sea floor, but gas hydrate was not physically recovered in cores. It is speculated that the poor recovery of cores and gas hydrate morphology are responsible for the lack of physical gas hydrate recovery.

Marine and Petroleum Geology↗

Geologic framework of the 2005 Keathley Canyon gas hydrate research well, northern Gulf of Mexico

The Keathley Canyon sites drilled in 2005 by the Chevron Joint Industry Project are located along the southeastern edge of an intraslope minibasin (Casey basin) in the northern Gulf of Mexico at 1335 m water depth. Around the drill sites, a grid of 2D high-resolution multichannel seismic data designed to image depths down to at least 1000 m sub-bottom reveals 7 unconformities and disconformities that, with the seafloor, bound 7 identifiable seismic stratigraphic units. A major disconformity in the middle of the units stands out for its angular baselapping geometry. From these data, three episodes of sedimentary deposition and deformation are inferred. The oldest episode consists of fine-grained muds deposited during a period of relative stability in the basin (units e, f, and g). Both the BSR and inferred gas hydrate occur within these older units. The gas hydrate occurs in near-vertical fractures. A second episode (units c and d) involved large vertical displacements associated with infilling and ponding of sediment. This second interval corresponds to deposition of intercalated fine and coarse-grained material that was recovered in the drill hole that penetrated the thin edges of the regionally much thicker units. The final episode of deposition (units a and b) occurred during more subdued vertical motions. Hemipelagic drape (unit a) characterizes the modern seafloor. The present-day Casey basin is mostly filled. Its sill is part of a subsiding graben structure that is only 10-20 m shallower than the deepest point in the basin, indicating that gravity-driven transport would mostly bypass the basin. Contemporary faulting along the basin margins has selectively reactivated an older group of faults. The intercalated sand and mud deposits of units c and d are tentatively correlated with Late Pleistocene deposition derived from the western shelf-edge delta/depocenter of the Mississippi River, which was probably most active from 320 ka to 70 ka [Winker, C.D., Booth, J., 2000. Sedimentary dynamics of the salt-dominated continental slope, Gulf of Mexico: integration of observations from the seafloor, near-surface, and deep subsurface. In: Proceedings of the GCSSEPM Foundation 20th Annual Research Conference, Deep-water Reservoirs of the World, pp. 1059-1086]. The presence of sand within the gas hydrate stability zone (in units c and d) is not sufficient to concentrate gas hydrate even though dispersed gas hydrate occurs deeper in the fractured mud/clay-rich sections of units e and f.

Gulf of Mexico↗

Using open hole and cased-hole resistivity logs to monitor gas hydrate dissociation during a thermal test in the mallik 5L-38 research well, Mackenzie Delta, Canada

Gas hydrates, which are naturally occurring ice-like combinations of gas and water, have the potential to provide vast amounts of natural gas from the world's oceans and polar regions. However, producing gas economically from hydrates entails major technical challenges. Proposed recovery methods such as dissociating or melting gas hydrates by heating or depressurization are currently being tested. One such test was conducted in northern Canada by the partners in the Mallik 2002 Gas Hydrate Production Research Well Program. This paper describes how resistivity logs were used to determine the size of the annular region of gas hydrate dissociation that occurred around the wellbore during the thermal test in the Mallik 5L-38 well. An open-hole logging suite, run prior to the thermal test, included array induction, array laterolog, nuclear magnetic resonance and 1.1-GHz electromagnetic propagation logs. The reservoir saturation tool was run both before and after the thermal test to monitor formation changes. A cased-hole formation resistivity log was run after the test.Baseline resistivity values in each formation layer (Rt) were established from the deep laterolog data. The resistivity in the region of gas hydrate dissociation near the wellbore (Rxo) was determined from electromagnetic propagation and reservoir saturation tool measurements. The radius of hydrate dissociation as a function of depth was then determined by means of iterative forward modeling of cased-hole formation resistivity tool response. The solution was obtained by varying the modeled dissociation radius until the modeled log overlaid the field log. Pretest gas hydrate production computer simulations had predicted that dissociation would take place at a uniform radius over the 13-ft test interval. However, the post-test resistivity modeling showed that this was not the case. The resistivity-derived dissociation radius was greatest near the outlet of the pipe that circulated hot water in the wellbore, where the highest temperatures were recorded. The radius was smallest near the center of the test interval, where a conglomerate section with low values of porosity and permeability inhibited dissociation. The free gas volume calculated from the resistivity-derived dissociation radii yielded a value within 20 per cent of surface gauge measurements. These results show that the inversion of resistivity measurements holds promise for use in future gas hydrate monitoring. ?? 2008 Society of Petrophysicists and Well Log Analysts. All rights reserved.

Conference Paper↗

Toward production from gas hydrates: Current status, assessment of resources, and simulation-based evaluation of technology and potential

Gas hydrates are a vast energy resource with global distribution in the permafrost and in the oceans. Even if conservative estimates are considered and only a small fraction is recoverable, the sheer size of the resource is so large that it demands evaluation as a potential energy source. In this review paper, we discuss the distribution of natural gas hydrate accumulations, the status of the primary international R&D programs, and the remaining science and technological challenges facing commercialization of production. After a brief examination of gas hydrate accumulations that are well characterized and appear to be models for future development and gas production, we analyze the role of numerical simulation in the assessment of the hydrate production potential, identify the data needs for reliable predictions, evaluate the status of knowledge with regard to these needs, discuss knowledge gaps and their impact, and reach the conclusion that the numerical simulation capabilities are quite advanced and that the related gaps are either not significant or are being addressed. We review the current body of literature relevant to potential productivity from different types of gas hydrate deposits, and determine that there are consistent indications of a large production potential at high rates over long periods from a wide variety of hydrate deposits. Finally, we identify (a) features, conditions, geology and techniques that are desirable in potential production targets, (b) methods to maximize production, and (c) some of the conditions and characteristics that render certain gas hydrate deposits undesirable for production. Copyright 2008, Society of Petroleum Engineers.

Conference Paper↗

Estimates of in situ gas hydrate concentration from resistivity monitoring of gas hydrate bearing sediments during temperature equilibration

As part of Ocean Drilling Program Leg 204 at southern Hydrate Ridge off Oregon we have monitored changes in sediment electrical resistivity during controlled gas hydrate dissociation experiments. Two cores were used, each filled with gas hydrate bearing sediments (predominantly mud/silty mud). One core was from Site 1249 (1249F-9H3), 42.1 m below seafloor (mbsf) and the other from Site 1248 (1248C-4X1), 28.8 mbsf. At Site 1247, a third experiment was conducted on a core without gas hydrate (1247B-2H1, 3.6 mbsf). First, the cores were imaged using an infra-red (IR) camera upon recovery to map the gas hydrate occurrence through dissociation cooling. Over a period of several hours, successive runs on the multi-sensor track (includes sensors for P-wave velocity, resistivity, magnetic susceptibility and gamma-ray density) were carried out complemented by X-ray imaging on core 1249F-9H3. After complete equilibration to room temperature (17-18??C) and complete gas hydrate dissociation, the final measurement of electrical resistivity was used to calculate pore-water resistivity and salinities. The calculated pore-water freshening after dissociation is equivalent to a gas hydrate concentration in situ of 35-70% along core 1249F-9H3 and 20-35% for core 1248C-4X1 assuming seawater salinity of in situ pore fluid. Detailed analysis of the IR scan, X-ray images and split-core photographs showed the hydrate mainly occurred disseminated throughout the core. Additionally, in core 1249F-9H3, a single hydrate filled vein, approximately 10 cm long and dipping at about 65??, was identified. Analyses of the logging-while-drilling (LWD) resistivity data revealed a structural dip of 40-80?? in the interval between 40 and 44 mbsf. We further analyzed all resistivity data measured on the recovered core during Leg 204. Generally poor data quality due to gas cracks allowed analyses to be carried out only at selected intervals at Sites 1244, 1245, 1246, 1247, 1248, 1249, and 1252. With a few exceptions, data from these intervals yield low to no gas hydrate concentration, which corresponds to estimates from downhole resistivity logs. However, since the gas cracking may be the result of gas hydrate dissociation, this is a biased sampling. Cores that had contained some gas hydrate may have been excluded. ?? 2005 Elsevier B.V. All rights reserved.

Marine Geology↗

Results at Mallik highlight progress in gas hydrate energy resource research and development

The recent studies that project the role of gas hydrates in the future energy resource management are reviewed. Researchers have long speculated that gas hydrates could eventually be a commercial resource for the future. A Joint Industry Project led by ChevronTexaco and the US Department of Energy is designed to characterize gas hydrates in the Gulf of Mexico. Countries including Japan, canada, and India have established large gas hydrate research and development projects, while China, Korea and Mexico are investigating the viability of forming government-sponsored gas hydrate research projects.

Conference Paper↗

Three-dimensional distribution of gas hydrate beneath southern Hydrate Ridge: Constraints from ODP Leg 204

Large uncertainties about the energy resource potential and role in global climate change of gas hydrates result from uncertainty about how much hydrate is contained in marine sediments. During Leg 204 of the Ocean Drilling Program (ODP) to the accretionary complex of the Cascadia subduction zone, we sampled the gas hydrate stability zone (GHSZ) from the seafloor to its base in contrasting geological settings defined by a 3D seismic survey. By integrating results from different methods, including several new techniques developed for Leg 204, we overcome the problem of spatial under-sampling inherent in robust methods traditionally used for estimating the hydrate content of cores and obtain a high-resolution, quantitative estimate of the total amount and spatial variability of gas hydrate in this structural system. We conclude that high gas hydrate content (30-40% of pore space or 20-26% of total volume) is restricted to the upper tens of meters below the seafloor near the summit of the structure, where vigorous fluid venting occurs. Elsewhere, the average gas hydrate content of the sediments in the gas hydrate stability zone is generally <2% of the pore space, although this estimate may increase by a factor of 2 when patchy zones of locally higher gas hydrate content are included in the calculation. These patchy zones are structurally and stratigraphically controlled, contain up to 20% hydrate in the pore space when averaged over zones ???10 m thick, and may occur in up to ???20% of the region imaged by 3D seismic data. This heterogeneous gas hydrate distribution is an important constraint on models of gas hydrate formation in marine sediments and the response of the sediments to tectonic and environmental change. ?? 2004 Published by Elsevier B.V.

Earth and Planetary Science Letters↗