Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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

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

At least 775 records · Page 43Linked to original sources
USGS2023

Modeling the dynamic penetration depth of post-1950s water in unconfined aquifers using environmental tracers: Central Valley, California

The penetration depth of post-1950s recharge (D-1950) in aquifers is a marker that is frequently used to identify groundwater that is susceptible to anthropogenic contamination. Here, we compute D-1950 values at wells, interpolate them in space, and project them across time to map the moving front of modern recharge in four dimensions in the Central Valley aquifer system, California, USA. Tracers of groundwater age (tritium, carbon-14, noble gases, sulfur hexafluoride, and chlorofluorocarbons) were collected at 650 wells spatially distributed throughout the Central Valley and were fit to a lumped-parameter model that assumes a logarithmic age-depth profile in the aquifer. For samples where tritium was present (>0.3 tritium units), the model was used to predict D-1950 at wells screened above or across the modern-premodern interface (n = 484). Wells with samples where tritium was absent (≤0.3 tritium units) were used to define the depth beyond which groundwater is completely premodern (n = 166). Predicted D-1950 values were below the depth of screen bottoms for wells where groundwater is completely modern, and above the depth of screen tops for wells where groundwater is completely premodern. The interpolated surface of D-1950 is dynamic, less prone to extreme values, and produces maps with lower interpolation errors due to a higher spatial density of wells than maps based on the depth of premodern groundwater. Between 2005 and 2025, D-1950 is expected to deepen by 11 and 12 m in the northern and southern parts of the Central Valley, respectively. Areas where D-1950 increases rapidly are likely to see increases in nitrate and other anthropogenic contaminants associated with the downward moving front of modern water.

California↗
USGS2023

Incorporating temperature into seepage loss estimates for a large unlined irrigation canal

Quantifying seepage losses from unlined irrigation canals is necessary to improve water use and conservation. The use of heat as a tracer is widely used in quantifying seepage rates across the sediment–water interface. In this study, field observations and two-dimensional numerical models were used to simulate seepage losses during the 2018 and 2019 irrigation season in the Truckee Canal system. Nineteen transects were instrumented with temperature probes and stage recording devices for inverse modeling to derive seepage flux and volumetric losses over the 39 km length of canal. The numerical models for each transect were calibrated and validated using the two-year dataset. Soil zones and observation data were used in each numerical model to help guide calibration of vertical and lateral heat and fluid fluxes. Model simulations were used to derive multivariable regression equations that consider stage, temperature, and hydraulic gradient. The results demonstrate the value of long-term datasets that illustrate the seasonality of groundwater levels, siltation, stage, and temperature on seepage rates. Seepage rates estimated by the numerical models range from 0.16 to 4.6 m 3 /d m −1 . Total annual volumetric losses estimated for 2018 and 2019 were 1.6 × 10 -2 to 1.2 × 10 -2 km 3 , respectively. The seepage losses estimated by this study account for 32 % to 41 % of the inflow volumes. Regression models were able to reproduce seepage time-series simulated by the numerical models reasonably well. In arid environments, water diverted into irrigation canals may be influenced by seasonal variations in temperature sufficient to influence the water accounting of conveyed surface flows.

Nevada↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2025

Impact of gas/liquid phase change of CO2 during injection for sequestration

CO 2 sequestration in deep saline formations is an effective and important process to control the rapid rise in CO 2 emissions. The process of injecting CO 2 requires reliable predictions of the stress in the formation and the fluid pressure distributions – particularly since monitoring of the CO 2 migration is difficult – to mitigate leakage, prevent induced seismicity, and analyze wellbore stability. A key aspect of CO 2 is the gas–liquid phase transition at the temperatures and pressures of relevance to leakage and sequestration, which has been recognized as being critical for accurate predictions but has been challenging to model without ad hoc empiricisms. This paper presents a robust multiphase thermodynamics-based poromechanics model to capture the complex phase transition behavior of CO 2 and predict the stress and pressure distribution under super- and sub- critical conditions during the injection process. A finite element implementation of the model is applied to analyze the behavior of a multiphase porous system with CO 2 as it displaces the fluid brine phase. We find that if CO 2 undergoes a phase transition in the geologic reservoir, the spatial variation of the density is significantly affected, and the migration mobility of CO 2 decreases in the reservoir. A key feature of our approach is that we do not a priori assume the location of the CO 2 gas/liquid interface – or even if it occurs at all – but rather, this is a prediction of the model, along with the spatial variation of the phase of CO 2 and the change of the saturation profile due to the phase change.

Journal of the Mechanics and Physics of Solids↗
USGS2006

Geodetic observations and modeling of magmatic inflation at the Three Sisters volcanic center, central Oregon Cascade Range, USA

Tumescence at the Three Sisters volcanic center began sometime between summer 1996 and summer 1998 and was discovered in April 2001 using interferometric synthetic aperture radar (InSAR). Swelling is centered about 5 km west of the summit of South Sister, a composite basaltic-andesite to rhyolite volcano that last erupted between 2200 and 2000 yr ago, and it affects an area ∼20 km in diameter within the Three Sisters Wilderness. Yearly InSAR observations show that the average maximum displacement rate was 3–5 cm/yr through summer 2001, and the velocity of a continuous GPS station within the deforming area was essentially constant from June 2001 to June 2004. The background level of seismic activity has been low, suggesting that temperatures in the source region are high enough or the strain rate has been low enough to favor plastic deformation over brittle failure. A swarm of about 300 small earthquakes ( M max = 1.9) in the northeast quadrant of the deforming area on March 23–26, 2004, was the first notable seismicity in the area for at least two decades. The U.S. Geological Survey (USGS) established tilt-leveling and EDM networks at South Sister in 1985–1986, resurveyed them in 2001, the latter with GPS, and extended them to cover more of the deforming area. The 2001 tilt-leveling results are consistent with the inference drawn from InSAR that the current deformation episode did not start before 1996, i.e., the amount of deformation during 1995–2001 from InSAR fully accounts for the net tilt at South Sister during 1985–2001 from tilt-leveling. Subsequent InSAR, GPS, and leveling observations constrain the source location, geometry, and inflation rate as a function of time. A best-fit source model derived from simultaneous inversion of all three datasets is a dipping sill located 6.5 ± 2.5 km below the surface with a volume increase of 5.0 × 10 6 ± 1.5 × 10 6 m 3 /yr (95% confidence limits). The most likely cause of tumescence is a pulse of basaltic magma intruding the upper crust along the brittle–ductile interface — a process that must occur episodically beneath the Cascade Range but in the past would have escaped detection in the absence of unusual seismicity. We speculate that such intrusive episodes last from days to years and are separated by quiescent periods of decades to centuries. The likelihood that the current episode at Three Sisters will culminate in an eruption is judged to be low, but the impact of an eruption could be great. The USGS has updated its volcano hazards assessment for the Three Sisters region, notified appropriate agencies and the public, and is helping to prepare an emergency coordination and communication plan.

Central Oregon Cascade Range, Three Sisters volcan↗
USGS2023

Estimating proximity effects to wildfire fuels treatments on house prices in Cibola National Forest, New Mexico, USA

Forested landscapes in the Western United States are subject to growing size and severity of wildfires, in part due to historical management strategies focusing on wildfire suppression. Forest restoration treatments and fuels reductions, including thinning and prescribed burning, can reduce the frequency and intensity of wildfires. Extensive restoration and fuels treatment efforts are underway across many areas in the Southwestern United States, including New Mexico. The tradeoff between amenity values provided by forested landscapes and the wildfire risk associated with forested landscapes is becoming increasingly important to understand as development in the wildland-urban interface increases. Understanding how house proximity, relative to forest restoration or fuels treatments, is capitalized into home sale prices can provide useful information about how individuals value forested landscapes that have been altered to reduce wildfire risk or severity. We use a Hedonic Property Model to estimate the average treatment effect of proximity to forest restoration or fuel treatments in New Mexico, United States. We use matching methods to estimate the average treatment effect of proximity to forest restoration. We find that proximity to the forest has a positive amenity value; however, proximity to recent forest restoration or fuel treatments results in a decrease in house sale prices. We combine the results of our two models and calculate that homes not within one kilometer of a treated forest and within one kilometer of Cibola National Forest sell for an average $73,626 premium. The average premium drops to $22,996 for homes within one kilometer of a forest that has been recently treated and within one kilometer of Cibola National Forest.

New Mexico↗
USGS2008

Eddy correlation measurements of submarine groundwater discharge

This paper presents a new, non-invasive means of quantifying groundwater discharge into marine waters using an eddy correlation approach. The method takes advantage of the fact that, in virtually all aquatic environments, the dominant mode of vertical transport near the sediment–water interface is turbulent mixing. The technique thus relies on measuring simultaneously the fluctuating vertical velocity using an acoustic Doppler velocimeter and the fluctuating salinity and/or temperature using rapid-response conductivity and/or temperature sensors. The measurements are typically done at a height of 5–15 cm above the sediment surface, at a frequency of 16 to 64 Hz, and for a period of 15 to 60 min. If the groundwater salinity and/or temperature differ from that of the water column, the groundwater specific discharge (cm d− 1) can be quantified from either a heat or salt balance. Groundwater discharge was estimated with this new approach in Salt Pond, a small estuary on Cape Cod (MA, USA). Estimates agreed well with previous estimates of discharge measured using seepage meters and 222Rn as a tracer. The eddy correlation technique has several desirable characteristics: 1) discharge is quantified under in-situ hydrodynamic conditions; 2) salinity and temperature can serve as two semi-independent tracers of discharge; 3) discharge can be quantified at high temporal resolution, and 4) long-term records of discharge may be possible, due to the low power requirements of the instrumentation.

Massachussetts↗
USGS2017

Perspectives on chemical oceanography in the 21st century: Participants of the COME ABOARD Meeting examine aspects of the field in the context of 40 years of DISCO

The questions that chemical oceanographers prioritize over the coming decades, and the methods we use to address these questions, will define our field's contribution to 21st century science. In recognition of this, the U.S. National Science Foundation and National Oceanic and Atmospheric Administration galvanized a community effort (the Chemical Oceanography MEeting: A BOttom-up Approach to Research Directions, or COME ABOARD) to synthesize bottom-up perspectives on selected areas of research in Chemical Oceanography. Representing only a small subset of the community, COME ABOARD participants did not attempt to identify targeted research directions for the field. Instead, we focused on how best to foster diverse research in Chemical Oceanography, placing emphasis on the following themes: strengthening our core chemical skillset; expanding our tools through collaboration with chemists, engineers, and computer scientists; considering new roles for large programs; enhancing interface research through interdisciplinary collaboration; and expanding ocean literacy by engaging with the public. For each theme, COME ABOARD participants reflected on the present state of Chemical Oceanography, where the community hopes to go and why, and actionable pathways to get there. A unifying concept among the discussions was that dissimilar funding structures and metrics of success may be required to accommodate the various levels of readiness and stages of knowledge development found throughout our community. In addition to the science, participants of the concurrent Dissertations Symposium in Chemical Oceanography (DISCO) XXV, a meeting of recent and forthcoming Ph.D. graduates in Chemical Oceanography, provided perspectives on how our field could show leadership in addressing long-standing diversity and early-career challenges that are pervasive throughout science. Here we summarize the COME ABOARD Meeting discussions, providing a synthesis of reflections and perspectives on the field.

Marine Chemistry↗
USGS2012

The FOBIMO (FOraminiferal BIo-MOnitoring) initiative—Towards a standardised protocol for soft-bottom benthic foraminiferal monitoring studies

The European Community Marine Strategy Framework Directive (MSFD) was established to provide guidelines for monitoring the quality of marine ecosystems. Monitoring the status of marine environments is traditionally based on macrofauna surveys, for which standardised methods have been established. Benthic foraminifera are also good indicators of environmental status because of their fast turnover rates, high degree of specialisation, and the preservation of dead assemblages in the fossil record. In spite of the growing interest in foraminiferal bio-monitoring during the last decades, no standardised methodology has been proposed until today. The aim of the FOraminiferal BIo-MOnitoring (FOBIMO) expert workshop, held in June 2011 at Fribourg, Switzerland, which assembled 37 scientists from 24 research groups and 13 countries, was to develop a suite of standard methods. This paper presents the main outcome of the workshop, a list of motivated recommendations with respect to sampling devices, sample storage, treatment, faunal analysis and documentation. Our recommendations fulfil the criteria imposed both by scientific rigour and by the practical limitations of routine studies. Hence, our aim is to standardise methodologies used in bio-monitoring only and not to limit the use of different methods in pure scientific studies. Unless otherwise stated, all recommendations concern living (stained) benthic foraminiferal assemblages. We have chosen to propose two types of recommendations. Mandatory recommendations have to be followed if a study wants to qualify as sound and compatible to the norms. The most important of these recommendations are the interval from 0 to 1 cm below the sediment surface has to be sampled, and an interface corer or box corer that keeps the sediment surface intact is to be used for offshore surveys. A grab sampler must not be deployed in soft sediments. Three replicate samples are to be taken and analysed separately. Samples are to be washed on a 63-μm screen, and the living benthic foraminiferal fauna of the > 125 μm fraction is to be analysed. Splits are to be picked and counted entirely, and all counted foraminifera from at least one replicate per station have to be stored in micropalaeontological slides. Census data, supplementary laboratory data and microslides have to be archived. Advisory recommendations are to sample in autumn, to have a sample size of 50 cm 2 or a tube of 8 cm inner diameter, to use > 70% ethanol as a preservative, rose Bengal at a concentration of 2 grams per litre for staining, and a staining time of at least 14 days. The split size should be defined by a target value of 300 specimens, heavy liquid separation should be avoided, and the 63–125 μm fraction or deeper sediment levels may be considered in some environments. We are convinced that the application of this protocol by a large number of scientists is a necessary first step to a general acceptance of benthic foraminifera as a reliable tool in bio-monitoring studies.

Marine Micropaleontology↗
USGS2018

Benthic foraminifera from the Carnarvon Ramp reveal variability in Leeuwin Current activity (Western Australia) since the Pliocene

Benthic foraminiferal assemblages from a ~300 m deep core from an outer carbonate-ramp site off Western Australia (International Ocean Discovery Program Core U1460A) were examined to reconstruct the paleoceanographic evolution of the Carnarvon Ramp and the warm surficial Leeuwin Current (LC) for the last 3.54 Ma. Of the identified 179 benthic foraminiferal species, occurrences of the 15 most abundant taxa were assessed using Q-mode Cluster Analysis and Non-Metric Dimensional Scaling. Diversity, equitability, planktonic/benthic index, microhabitat preference, and sedimentary parameters such as lithology and sponge spicule content were analyzed to gather information about past intermediate- and surface-water circulation. Relative abundances of infaunal and epifaunal species were applied to indicate changes in organic-matter supply and oxygenation at the sea floor. Influence of upwelling was recognized by a high infaunal species ratio, with dominance by Uvigerina peregrina , Lagena annellatrachia and Trifarina bradyi . Epifaunal species such as Hanzawaia nipponica and Hyalinea florenceae gradually became more abundant around 1.14 Ma, indicating increased ventilation and establishment of the paleo-LC. A more substantial change was initiated by 0.91 Ma as marked by key species Spirorutilus carinatus and Rotorbinella sp. , together with increased faunal diversity, benthic foraminiferal accumulation rates, and evidence for suspension feeding sponges. With the LC flow suppressing upwelling, and better ventilated waters entering the shelf, the environment favored epifaunal agglutinates, rotalids, and miliolids, while buliminids decreased. Under high-flow conditions of the LC, sponge spicules and skeletal carbonate production reached an optimum at ~0.6 Ma before returning to modern conditions. Supported by these observations, we propose the following paleoceanographic evolution of the Carnarvon Ramp: During the late Pliocene to mid Pleistocene (3.54–0.91 Ma) conditions of deep-water upwelling from the Western Australian Current and Indian Ocean Gyre indicate the absence of the capping LC on the outer carbonate ramp. A transitional phase started in the mid Pleistocene (1.14–0.61 Ma). The paleo-LC triggered gradual oxygenation at the sediment-water interface, which coincided with an increase in carbonate sedimentation rates, and waning sea-surface productivity. During a third phase, mid Pleistocene to present (0.91–0 Ma), the LC’s intensity and flow rates peaked at ~0.6 Ma. Benthic foraminiferal accumulation rates reached a high, then decreased to present-day rates. For short periods, sea-surface productivity was moderately enhanced, likely due to fluctuating LC persistence or landward shift during glacial maxima.

Western Australia↗
USGS2008

Natural gas geochemistry of sediments drilled on the 2005 Gulf of Mexico JIP cruise

In April and May 2005, cores were acquired and sub-sampled for gases in lease blocks Atwater Valley 13 and 14 and Keathley Canyon 151 during deep subseafloor drilling conducted as part of the JIP study of gas hydrates in the northern Gulf of Mexico. Sample types included sediment headspace gas, free gas derived from sediment gas exsolution, and gas exsolution from controlled degassing of pressurized cores. The gases measured both onboard and in shore-based labs were nitrogen, oxygen, hydrogen sulfide, carbon dioxide, and the hydrocarbons methane through hexane. The presence of seafloor mounds, seismic anomalies, a shallow sulfate-methane interface, and similar gas compositions and isotopic compositions near the seafloor and at depth suggest an upward flux of methane at both sites. Sediment gases at the Atwater Valley sites, where seafloor mounds and adjacent sediments were cored, strongly suggest a microbial source of methane, with very little thermogenic gas input. Sediment gas from all cores contained from about 96 to 99.9% methane, with the balance composed primarily of carbon dioxide. Methane to ethane ratios were greater than 1000, and often over 10,000. Gases from cores at Keathley Canyon were similar to those at Atwater Valley, however, deeper cores from Keathley Canyon contained more ethane, propane, and butane suggesting mixing with minor concentrations thermogenic gas. The isotopic composition of methane, ethane, and carbon dioxide were measured, and ??13C values range from -84.3 to -71.5???, -65.2 to -46.8???, and -23.5 to -3.0???, respectively, all consistent with microbial gas sources, early diagenesis of organic matter and perhaps biodegradation of petroleum. The presence of deep microbial gas at these sites here and elsewhere highlights a potentially significant, predominantly microbial gas source in the northern Gulf of Mexico.

Marine and Petroleum Geology↗
USGS2020

Potential freshening impacts on fines migration and pore-throat clogging during gas hydrate production: 2-D micromodel study with Diatomaceous UBGH2 sediments

The methane gas hydrate stored in natural sediments is considered a potential gas resource. Countries such as China, India, Japan, and Korea are interested in commercializing this resource, and offshore field pilot tests for gas production have been conducted using depressurization methods to destabilize gas hydrate and facilitate the migration of methane to the production well. However, fine-grained sediments (fines), which are present even in coarse-grained, gas hydrate-bearing sediments, can be resuspended in the production fluid, subsequently clogging pore throats in the formation and reducing the overall production efficiency. We conducted laboratory tests to evaluate the suspension and clogging potential of fines collected from the Ulleung Basin, East Sea, Korea during the 2010 Ulleung Basin Gas Hydrate Expedition 2 (UBGH2). Experimental results reveal that diatoms are prevalent in the sediment and largely control the suspension and clogging behavior. Fluid flow experiments in 2D micromodels show clogging occurs even when injecting the minimum sediment concentration (0.1wt% in the fluid) through micromodels with pore-throat widths at the high end of the anticipated range for UBGH2 gas hydrate-bearing sands (100µm). Mobile gas/fluid interfaces forming during gas hydrate dissociation accentuate clogging by concentrating and mobilizing fines. Sedimentation tests show pore-water freshening during dissociation is not anticipated to change the potential for diatoms to become entrained in the pore water flow, even for the observed gas hydrate saturations of ~80%. Muscovite and illite are also significant components of the tested sediment, however, and pore-water freshening increases their potential for resuspension and clogging. Overall, the resuspension and clogging potential of these fine sediments should increase as gas hydrate dissociation progresses in the thin, gas hydrate-bearing sands investigated in the Ulleung Basin.

Marine and Petroleum Geology↗
USGS2025

Avian influenza spillover into poultry: Environmental influences and biosecurity protections

With the continued spread of highly pathogenic avian influenza (HPAI), understanding the complex dynamics of virus transfer at the wild – agriculture interface is paramount. Spillover events (i.e., virus transfer from wild birds into poultry) are related to proximity to infected wild bird populations and environmental conditions. By accounting for such dynamics, we can take a combined approach to assess the impacts of biosecurity measures implemented at poultry farms while simultaneously accounting for their local risk levels. We implemented a Bayesian joint-likelihood logistic regression for the Continental U.S. comparing models of spatiotemporal risk according to land use, weather, and predicted waterfowl distributions followed by integrating a farm-level case-control questionnaire dataset focused on identifying trends in HPAI spillover risk associated with a farm's biosecurity practices. We found that estimates of waterfowl abundance, along with mean precipitation and temperature during winter, were most correlated with spatiotemporal HPAI risk. Additionally, we identified multiple biosecurity practices associated with reduced risk to HPAI, where the strongest relationships were related to litter decontamination treatments, vehicle wash stations, and avoiding shared dead-bird disposal sites with other farms. This model broadly guides surveillance of HPAI in wild and domestic populations, identifying when and where we are most likely to see increased instances of the virus while also providing insights into how poultry farms can better protect themselves from risk.

contiguous United States↗
USGS2006

Earthquake lights and the stress-activation of positive hole charge carriers in rocks

Earthquake-related luminous phenomena (also known as earthquake lights) may arise from (1) the stress-activation of positive hole (p-hole) charge carriers in igneous rocks and (2) the accumulation of high charge carrier concentrations at asperities in the crust where the stress rates increase very rapidly as an earthquake approaches. It is proposed that, when a critical charge carrier concentration is reached, the p-holes form a degenerated solid state plasma that can break out of the confined rock volume and propagate as a rapidly expanding charge cloud. Upon reaching the surface the charge cloud causes dielectric breakdown at the air-rock interface, i.e. corona discharges, accompanied by the emission of light and high frequency electromagnetic radiation. ?? 2006 Elsevier Ltd. All rights reserved.

Physics and Chemistry of the Earth↗
USGS2009

Geomechanical response of permafrost-associated hydrate deposits to depressurization-induced gas production

In this simulation study, we analyzed the geomechanical response during depressurization production from two known hydrate-bearing permafrost deposits: the Mallik (Northwest Territories, Canada) deposit and Mount Elbert (Alaska, USA) deposit. Gas was produced from these deposits at constant pressure using horizontal wells placed at the top of a hydrate layer (HL), located at a depth of about 900??m at the Mallik site and 600??m at the Mount Elbert site. The simulation results show that general thermodynamic and geomechanical responses are similar for the two sites, but with substantially higher production and more intensive geomechanical responses at the deeper Mallik deposit. The depressurization-induced dissociation begins at the well bore and then spreads laterally, mainly along the top of the HL. The depressurization results in an increased shear stress within the body of the receding hydrate and causes a vertical compaction of the reservoir. However, its effects are partially mitigated by the relatively stiff permafrost overburden, and compaction of the HL is limited to less than 0.4%. The increased shear stress may lead to shear failure in the hydrate-free zone bounded by the HL overburden and the downward-receding upper dissociation interface. This zone undergoes complete hydrate dissociation, and the cohesive strength of the sediment is low. We determined that the likelihood of shear failure depends on the initial stress state as well as on the geomechanical properties of the reservoir. The Poisson's ratio of the hydrate-bearing formation is a particularly important parameter that determines whether the evolution of the reservoir stresses will increase or decrease the likelihood of shear failure.

Journal of Petroleum Science and Engineering↗
USGS2011

A working environment for digital planetary data processing and mapping using ISIS and GRASS GIS

Since the beginning of planetary exploration, mapping has been fundamental to summarize observations returned by scientific missions. Sensor-based mapping has been used to highlight specific features from the planetary surfaces by means of processing. Interpretative mapping makes use of instrumental observations to produce thematic maps that summarize observations of actual data into a specific theme. Geologic maps, for example, are thematic interpretative maps that focus on the representation of materials and processes and their relative timing. The advancements in technology of the last 30 years have allowed us to develop specialized systems where the mapping process can be made entirely in the digital domain. The spread of networked computers on a global scale allowed the rapid propagation of software and digital data such that every researcher can now access digital mapping facilities on his desktop. The efforts to maintain planetary missions data accessible to the scientific community have led to the creation of standardized digital archives that facilitate the access to different datasets by software capable of processing these data from the raw level to the map projected one. Geographic Information Systems (GIS) have been developed to optimize the storage, the analysis, and the retrieval of spatially referenced Earth based environmental geodata; since the last decade these computer programs have become popular among the planetary science community, and recent mission data start to be distributed in formats compatible with these systems. Among all the systems developed for the analysis of planetary and spatially referenced data, we have created a working environment combining two software suites that have similar characteristics in their modular design, their development history, their policy of distribution and their support system. The first, the Integrated Software for Imagers and Spectrometers (ISIS) developed by the United States Geological Survey, represents the state of the art for processing planetary remote sensing data, from the raw unprocessed state to the map projected product. The second, the Geographic Resources Analysis Support System (GRASS) is a Geographic Information System developed by an international team of developers, and one of the core projects promoted by the Open Source Geospatial Foundation (OSGeo). We have worked on enabling the combined use of these software systems throughout the set-up of a common user interface, the unification of the cartographic reference system nomenclature and the minimization of data conversion. Both software packages are distributed with free open source licenses, as well as the source code, scripts and configuration files hereafter presented. In this paper we describe our work done to merge these working environments into a common one, where the user benefits from functionalities of both systems without the need to switch or transfer data from one software suite to the other one. Thereafter we provide an example of its usage in the handling of planetary data and the crafting of a digital geologic map. ?? 2010 Elsevier Ltd. All rights reserved.

Conference Paper↗
USGS2025

Pysochron: A Python-based solution for calculating cosmogenic 26Al/10Be isochron burial ages

Cosmogenic 26 Al/ 10 Be isochron burial dating is a powerful tool for dating sediment burial over the past several million years. By measuring in-situ 26 Al and 10 Be in a suite of samples from the same depth in a buried deposit, it is possible to quantify the inventory of cosmogenic nuclides produced after burial, date the burial of shallow sediments, identify sediment reworking, and calculate paleo-erosion rates. While this approach has been used to date materials around the world for over a decade, few published codes exist for performing 26 Al/ 10 Be isochron calculations. The isochron calculation options that are available typically rely on numerous files and libraries, rendering modification and troubleshooting difficult. Moreover, the widespread use of proprietary programming languages – and their associated addon packages – can place an additional financial burden on an already costly endeavor. Pysochron ( https://code.usgs.gov/recon/pysochron ) provides a solution to these issues. In its base form, it exists as a single script that can be easily modified, upgraded, and shared. Because it was developed in an open-source environment, all required computational packages are available free of charge. A user-friendly interface allows rapid modification of calculation parameters, and an automated commentary on isochron results provides insights and recommendations. Pysochron has been validated with 40 published cosmogenic 26 Al/ 10 Be burial isochrons around the world, with burial ages ranging from ∼5 Ma to ∼180 ka. As such, it is a promising option for members of the cosmogenic nuclide community seeking a straightforward, cost-effective, and flexible solution to isochron burial dating challenges.

Quaternary Geochronology↗
USGS2010

Mass spectrometric U-series dating of the Chaoxian hominin site at Yinshan, eastern China

The fossils of Chaoxian hominin, widely accepted as representing archaic Homo sapiens in eastern China, were recovered from the middle or slightly higher levels of Layer 2 deposits of a collapsed cave at Yinshan, Anhui Province. Results of mass spectrometric U-series dating of intercalated speleothem calcites are presented. Based mainly on four broadly coeval calcite samples, the hominin fossils should be bracketed in the range of 310–360 ka or somewhat older. These ages are much older than the previous estimate at 160–200 ka based on the U-series dating of fossil teeth and bones, and may be cited as supporting evidence for an earlier H. erectus –archaic H. sapiens interface in China.

Yinshan↗