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 379 records · Page 21Linked to original sources
USGS2011

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server—for example, Excel®, Visual Basic®, Python, or MATLAB". PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations.

Computers & Geosciences↗
USGS2015

The data quality analyzer: a quality control program for seismic data

The U.S. Geological Survey's Albuquerque Seismological Laboratory (ASL) has several initiatives underway to enhance and track the quality of data produced from ASL seismic stations and to improve communication about data problems to the user community. The Data Quality Analyzer (DQA) is one such development and is designed to characterize seismic station data quality in a quantitative and automated manner. The DQA consists of a metric calculator, a PostgreSQL database, and a Web interface: The metric calculator, SEEDscan, is a Java application that reads and processes miniSEED data and generates metrics based on a configuration file. SEEDscan compares hashes of metadata and data to detect changes in either and performs subsequent recalculations as needed. This ensures that the metric values are up to date and accurate. SEEDscan can be run as a scheduled task or on demand. The PostgreSQL database acts as a central hub where metric values and limited station descriptions are stored at the channel level with one-day granularity. The Web interface dynamically loads station data from the database and allows the user to make requests for time periods of interest, review specific networks and stations, plot metrics as a function of time, and adjust the contribution of various metrics to the overall quality grade of the station. The quantification of data quality is based on the evaluation of various metrics (e.g., timing quality, daily noise levels relative to long-term noise models, and comparisons between broadband data and event synthetics). Users may select which metrics contribute to the assessment and those metrics are aggregated into a “grade” for each station. The DQA is being actively used for station diagnostics and evaluation based on the completed metrics (availability, gap count, timing quality, deviation from a global noise model, deviation from a station noise model, coherence between co-located sensors, and comparison between broadband data and synthetics for earthquakes) on stations in the Global Seismographic Network and Advanced National Seismic System.

Computers & Geosciences↗
USGS2007

Barite-forming environments along a rifted continental margin, Southern California Borderland

The Southern California Continental Borderland (SCCB) is part of the broad San Andreas transform-fault plate boundary that consists of a series of fault-bounded, petroleum-generating basins. The SCCB has high heat flow and geothermal gradients produced by thinned continental crust and Neogene volcanism. Barite deposits in the SCCB occur along faults. Barite samples from two sea-cliff sites and four offshore sites in the SCCB were analyzed for mineralogy, chemical (54 elements) and isotopic (S, Sr) compositions, and petrography. Barite from Palos Verdes (PV) Peninsula sea-cliff outcrops is hosted by the Miocene Monterey Formation and underlying basalt; carbonate rocks from those outcrops were analyzed for C, O, and Sr isotopes and the basalt for S isotopes. Cold-seep barite from Monterey Bay, California was analyzed for comparison. SCCB offshore samples occur at water depths from about 500 to 1800 m. Those barites vary significantly in texture and occurrence, from friable, highly porous actively growing seafloor mounds to dense, brecciated, vein barite. This latter type of barite contrasts with cold-seep barite in being much more coarse grained, forms thick veins in places, and completely replaced rock clasts in breccia. The barite samples range from 94 to 99 wt% BaSO4, with low trace-element contents, except for high Sr, Zr, Br, U, and Hg concentrations compared to their crustal abundances. δ34S for SCCB offshore barites range from 21.6‰ to 67.4‰, and for PV barite from 62‰ to 70‰. Pyrite from PV sea-cliff basalt and sedimentary rocks that host the barites averages 7.8‰ and 2.2‰, respectively. Two offshore barite samples have δ34S values (21.6‰, 22.1‰) close to that of modern seawater sulfate, whereas all other samples are enriched to strongly enriched in 34S. 87Sr/86Sr ratios for the barites vary over a narrow range of 0.70830–0.70856 and are much lower than that of modern seawater and also lower than the middle Miocene seawater ratio, the time of deposition of the host rocks for the PV samples. This indicates that the mineralizing fluids were not unaltered seawater. We develop a model in which the barites precipitated both below the sediment–water interface and at the seafloor from low-temperature fluids that circulated along faults. The isotopic, chemical, and textural data indicate that the barites formed by several processes. Mesozoic and Cenozoic basement rocks (crystalline and overlying sedimentary rocks), Quaternary basin fill, turbidite fans, and seawater provided various elements for the barites in different environments. The fluids had a deep-seated source and were hydrothermal in the deeper parts of the system for all the barite types, including the seafloor cold-seep deposit, based on Sr isotopes and regional geothermal gradients. These deep-seated fluids mixed with other fluids as they ascended, including seawater at and near the seafloor. The high δ34S values may have resulted from extreme Rayleigh fractionation during bacterially mediated (semi)closed-system sulfate reduction, probably driven by the flux of methane- and hydrocarbon-bearing fluids through basement rocks and basin sediments. Early diagenetic dolomite and pyrite in the host Monterey Formation in the PV Headland also formed in a zone of sulfate reduction, but within sediment buried only centimeters to a few meters and with a predominantly seawater source for the sulfur. Dolomite occurring with vein barite in the PV Headland formed at temperatures in the range of 40–90 °C. The cold-seep barites have δ34S values near that of modern seawater, although still somewhat fractionated. The barites that precipitated below the sediment–water interface have higher δ34S values, suggesting that the fluids were relatively reduced with molar dissolved barium in excess of dissolved sulfate. Those samples were exposed at the seafloor by uplift along faults and are composed predominantly of massive, brecciated, and vein barite.

California↗
USGS2020

History and evolution of seepage meters for quantifying flow between groundwater and surface water: Part 1 – Freshwater settings

More than 75 years after its introduction, the seepage meter remains the only device for directly quantifying exchange across the sediment-water interface between groundwater and surface water. This device, first presented in the literature in the 1940s, has been in a state of near-constant improvement and design change, necessitating a review of the history and evolution of the device and a description of current best-measurement practices. Part 1 of this two-part review documents the evolution of seepage meters deployed in freshwater settings, including a listing of suggestions for best-measurement and deployment practices. Part 2 covers the same scope for seepage meters deployed in marine settings. Traditional seepage meters isolate a portion of the sediment bed; seepage commonly is determined by routing the volume of flow across that isolated interface to or from a submerged measurement bag over a known time interval. The time-integrated volume is then divided by the bed area covered by the meter to obtain a seepage flux expressed in distance per time. Both the instrument and the measurement are deceptively simple, leading some early users to question the viability of the measurement. Numerous sources of error have been identified and addressed over the decades, resulting in large improvements in measurement consistency and accuracy. Duration of each measurement depends on the seepage rate and can vary from minutes to days, leading to the erroneous and yet common assumption that seepage is relatively stable over time. Designs that replace the measurement bag with a flowmeter eliminate bag-related errors and provide much finer temporal resolution. Resulting data indicate seepage is highly variable in many settings and responds to numerous sub-daily processes, including evapotranspiration, rainfall, seiches and waves. Combining direct measurements from seepage meters with other measurements, such as vertical hydraulic gradients and vertical temperature profiles, provides far better understanding of the processes that control exchange between groundwater and surface water.

Earth-Science Reviews↗
USGS2014

Triggered aseismic slip adjacent to the 6 February 2013 Mw 8.0 Santa Cruz Islands megathrust earthquake

Aseismic or slow slip events have been observed in many subduction zones, but whether they affect the occurrence of earthquakes or result from stress changes caused by nearby events is unclear. In an area lacking direct geodetic observations, inferences can be made from seismological studies of co-seismic slip, associated stress changes and the spatiotemporal nature of aftershocks. These observations indicate that the February 2013 Mw 8.0 Santa Cruz Islands earthquake may have triggered slow or aseismic slip on an adjacent section of the subduction thrust over the following hours to days. This aseismic event was equivalent to Mw 7.6, significantly larger than any earthquakes in the aftershock sequence. The aseismic slip was situated within the seismogenic portion of the subduction interface, and must have occurred to the south of the main seismic slip and most aftershocks in order to promote right-lateral faulting in the upper plate, the dominant deformation style of the aftershock sequence. This plate boundary segment can support either stable sliding (aseismic) or stick-slip (seismic) deformation in response to different driving conditions. The complete lack of aftershocks on the thrust interface implies this pair of megathrust slip episodes (seismic and aseismic) released a substantial portion of the stored strain on the northernmost section of the Vanuatu subduction zone.

Santa Cruz Islands↗
USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗
USGS2008

Weathering of the Rio Blanco quartz diorite, Luquillo Mountains, Puerto Rico: Coupling oxidation, dissolution, and fracturing

In the mountainous Rio Icacos watershed in northeastern Puerto Rico, quartz diorite bedrock weathers spheroidally, producing a 0.2-2 m thick zone of partially weathered rock layers (???2.5 cm thickness each) called rindlets, which form concentric layers around corestones. Spheroidal fracturing has been modeled to occur when a weathering reaction with a positive ??V of reaction builds up elastic strain energy. The rates of spheroidal fracturing and saprolite formation are therefore controlled by the rate of the weathering reaction. Chemical, petrographic, and spectroscopic evidence demonstrates that biotite oxidation is the most likely fracture-inducing reaction. This reaction occurs with an expansion in d (0 0 1) from 10.0 to 10.5 A??, forming 'altered biotite'. Progressive biotite oxidation across the rindlet zone was inferred from thin sections and gradients in K and Fe(II). Using the gradient in Fe(II) and constraints based on cosmogenic age dates, we calculated a biotite oxidation reaction rate of 8.2 ?? 10-14 mol biotite m-2 s-1. Biotite oxidation was documented within the bedrock corestone by synchrotron X-ray microprobe fluorescence imaging and XANES. X-ray microprobe images of Fe(II) and Fe(III) at 2 ??m resolution revealed that oxidized zones within individual biotite crystals are the first evidence of alteration of the otherwise unaltered corestone. Fluids entering along fractures lead to the dissolution of plagioclase within the rindlet zone. Within 7 cm surrounding the rindlet-saprolite interface, hornblende dissolves to completion at a rate of 6.3 ?? 10-13 mol hornblende m-2 s-1: the fastest reported rate of hornblende weathering in the field. This rate is consistent with laboratory-derived hornblende dissolution rates. By revealing the coupling of these mineral weathering reactions to fracturing and porosity formation we are able to describe the process by which the quartz diorite bedrock disaggregates and forms saprolite. In the corestone, biotite oxidation induces spheroidal fracturing, facilitating the influx of fluids that react with other minerals, dissolving plagioclase and chlorite, creating additional porosity, and eventually dissolving hornblende and precipitating secondary minerals. The thickness of the resultant saprolite is maintained at steady state by a positive feedback between the denudation rate and the weathering advance rate driven by the concentration of pore water O2 at the bedrock-saprolite interface. ?? 2008 Elsevier Ltd. All rights reserved.

Geochimica et Cosmochimica Acta↗
USGS2017

Lithological influences on contemporary and long-term regolith weathering at the Luquillo Critical Zone Observatory

Lithologic differences give rise to the differential weatherability of the Earth’s surface and globally variable silicate weathering fluxes, which provide an important negative feedback on climate over geologic timescales. To isolate the influence of lithology on weathering rates and mechanisms, we compare two nearby catchments in the Luquillo Critical Zone Observatory in Puerto Rico, which have similar climate history, relief and vegetation, but differ in bedrock lithology. Regolith and pore water samples with depth were collected from two ridgetops and at three sites along a slope transect in the volcaniclastic Bisley catchment and compared to existing data from the granitic Río Icacos catchment. The depth variations of solid-state and pore water chemistry and quantitative mineralogy were used to calculate mass transfer (tau) and weathering solute profiles, which in turn were used to determine weathering mechanisms and to estimate weathering rates. Regolith formed on both lithologies is highly leached of most labile elements, although Mg and K are less depleted in the granitic than in the volcaniclastic profiles, reflecting residual biotite in the granitic regolith not present in the volcaniclastics. Profiles of both lithologies that terminate at bedrock corestones are less weathered at depth, near the rock-regolith interfaces. Mg fluxes in the volcaniclastics derive primarily from dissolution of chlorite near the rock-regolith interface and from dissolution of illite and secondary phases in the upper regolith, whereas in the granitic profile, Mg and K fluxes derive from biotite dissolution. Long-term mineral dissolution rates and weathering fluxes were determined by integrating mass losses over the thickness of solid-state weathering fronts, and are therefore averages over the timescale of regolith development. Resulting long-term dissolution rates for minerals in the volcaniclastic regolith include chlorite: 8.9 × 10 −14 mol m −2 s −1 , illite: 2.1 × 10 −14 mol m −2 s −1 and kaolinite: 4.0 × 10 −14 mol m −2 s −1 . Long-term weathering fluxes are several orders of magnitude lower in the granitic regolith than in the volcaniclastic, despite higher abundances of several elements in the granitic regolith. Contemporary weathering fluxes were determined from net (rain-corrected) solute profiles and thus represent rates over the residence time of water in the regolith. Contemporary weathering fluxes within the granitic regolith are similar to the long-term fluxes. In contrast, the long-term fluxes are faster than the contemporary fluxes in the volcaniclastic regolith. Contemporary fluxes in the granitic regolith are generally also slightly faster than in the volcaniclastic. The differences in weathering fluxes over space and time between these two watersheds indicate significant lithologic control of chemical weathering mechanisms and rates.

Bisley watersheds↗
USGS2006

Ge/Si and 87Sr/86Sr tracers of weathering reactions and hydrologic pathways in a tropical granitoid system

Ge/Si and 87Sr/86Sr data from primary and secondary minerals, soil waters, and stream waters in a tropical granitoid catchment quantitatively reflect mineral alteration reactions that occur at different levels within the bedrock-saprolite-soil zone. Near the bedrock-saprolite interface, plagioclase to kaolinite reaction yields low Ge/Si and 87Sr/86Sr. Higher in the regolith column, biotite weathering and kaolinite dissolution drive Ge/Si and 87Sr/86Sr to high values. Data from streams at base flow sample the bedrock-saprolite interface zone, while at high discharge solutes are derived from upper saprolite-soil zone. Coupled Ge/Si and 87Sr/86Sr can be effective tools for quantifying the importance of specific weathering reactions, and for geochemical hydrograph separation. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Geochemical Exploration↗
USGS2009

Shallow radar (SHARAD) sounding observations of the Medusae Fossae Formation, Mars

The SHARAD (shallow radar) sounding radar on the Mars Reconnaissance Orbiter detects subsurface reflections in the eastern and western parts of the Medusae Fossae Formation (MFF). The radar waves penetrate up to 580 m of the MFF and detect clear subsurface interfaces in two locations: west MFF between 150 and 155° E and east MFF between 209 and 213° E. Analysis of SHARAD radargrams suggests that the real part of the permittivity is ∼3.0, which falls within the range of permittivity values inferred from MARSIS data for thicker parts of the MFF. The SHARAD data cannot uniquely determine the composition of the MFF material, but the low permittivity implies that the upper few hundred meters of the MFF material has a high porosity. One possibility is that the MFF is comprised of low-density welded or interlocked pyroclastic deposits that are capable of sustaining the steep-sided yardangs and ridges seen in imagery. The SHARAD surface echo power across the MFF is low relative to typical martian plains, and completely disappears in parts of the east MFF that correspond to the radar-dark Stealth region. These areas are extremely rough at centimeter to meter scales, and the lack of echo power is most likely due to a combination of surface roughness and a low near-surface permittivity that reduces the echo strength from any locally flat regions. There is also no radar evidence for internal layering in any of the SHARAD data for the MFF, despite the fact that tens-of-meters scale layering is apparent in infrared and visible wavelength images of nearby areas. These interfaces may not be detected in SHARAD data if their permittivity contrasts are low, or if the layers are discontinuous. The lack of closely spaced internal radar reflectors suggests that the MFF is not an equatorial analog to the current martian polar deposits, which show clear evidence of multiple internal layers in SHARAD data.

Icarus↗
USGS2004

Field investigation into unsaturated flow and transport in a fault: Model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model. ?? 2004 Published by Elsevier B.V.

Journal of Contaminant Hydrology↗
USGS2018

Research note: Mapping spatial patterns in sewer age, material, and proximity to surface waterways to infer sewer leakage hotspots

Identifying areas where deteriorating sewer infrastructure is in close proximity to surface waterways is needed to map likely connections between sewers and streams. We present a method to estimate sewer installation year and deterioration status using historical maps of the sewer network, parcel-scale property assessment data, and pipe material. Areas where streams were likely buried into the sewer system were mapped by intersecting the historical stream network derived from a 10-m resolution digital elevation model with sewer pipe locations. Potential sewer leakage hotspots were mapped by identifying where aging sewer pipes are in close proximity (50-m) to surface waterways. Results from Pittsburgh, Pennsylvania (USA), indicated 41% of the historical stream length was lost or buried and the potential interface between sewers and streams is great. The co-location of aging sewer infrastructure (>75 years old) near stream channels suggests that 42% of existing streams are located in areas with a high potential for sewer leakage if sewer infrastructure fails. Mapping the sewer-stream interface provides an approach to better understand areas were failing sewers may contribute a disproportional amount of nutrients and other pathogens to surface waterways.

Pennsylvania↗
USGS2019

2D micromodel study of clogging behavior of fine-grained particles associated with gas hydrate production in NGHP-02 gas hydrate reservoir sediments

Fine-grained particles (fines) commonly coexist with coarse-grained sediments that host gas hydrate. These fines can be mobilized by liquid and gas flow during gas hydrate production. Once mobilized, fines can clog pore throats and reduce reservoir permeability. Even where particle sizes are smaller than pore-throat sizes, clogs can form due to clusters of fines. For certain types of fines, particularly swelling clays, cluster sizes depend on pore-fluid chemistry, which changes as pore-fluid freshens during gas hydrate dissociation. Fines can also be concentrated by a moving gas/liquid interface, increasing the chances of pore-throat clogging regardless of fines type. To test the relative significance of these clogging mechanisms, 2D micromodel experiments have been conducted with different pore-throat widths (20, 40, 60 and 100 m), single-phase pore-fluids (deionized water and 2M-sodium-chloride solution), and moving gas/liquid interfaces on specimens from Sites NGHP-02-09 and NGHP-02-16 (NGHP-02: National Gas Hydrate Program Expedition 02) as well as a selection of pure fines (silica silt, mica, calcium carbonate, diatoms, kaolin, and bentonite). Clogging depended on the ratio of particle-to-pore throat size, and also on pore-fluid chemistry because the pore-fluid chemistry changes effectively increased or decreased the fines cluster size relative to the pore-throat width. These interactions can be predicted based on the fines electrical sensitivity (defined by Jang and Santamarina, 2016). The fine-grained sediment component (grain size < 75 m) from the primary gas hydrate reservoir layers at Sites NGHP-02-09 and -16 show clogging via blocking or size exclusion (sieving) due to the large particles. Clogs also formed due to bridging or blocking by clusters of the smaller particles. Clogging generally occurred for pore-water sediment concentrations so low (0.2% by mass or less), that it was difficult to resolve the enhanced clogging in the presence of the gas/liquid meniscus.

Bay of Bengal↗
USGS2019

Numerical simulations of depressurization-induced gas production from an interbedded turbidite gas hydrate-bearing sedimentary section in the offshore India: Site NGHP-02-16 (Area-B)

The recent National Gas Hydrate Program Expedition 02 (NGHP-02) identified the existence of gas hydrate-bearing sand reservoirs at a number of sites in the offshore of India including Site NGHP-02-16 in Area-B of the Krishna-Godavari Basin. The architecture of that gas hydrate accumulation is characterized by thin, gas hydrate-bearing, high quality sand layers interbedded with mud layers within a turbidite interval. The lowest gas hydrate-bearing layer contacting a thinly-interbedded saline aquifer designates the base of the gas hydrate stability zone (BGHSZ). The proximity of the BGHSZ and the average temperature around 20oC make the reservoir a favorable target for hydrate destabilization by means of the depressurization method. The results of the reservoir simulations indicate high gas production potential from this marine gas hydrate deposit with manageable concomitant water production using a well completion design that hydraulically isolates layers with water-saturated sands. Using a detailed geological input model, the predicted cumulative gas rates reach 3.1 x 104 m3/day (1.1 mmscf/day) after 90 days of continuous depressurization and demonstrate sustained production rates of 3.0 x 104 m3/day (1.0 mmscf/day) after 5 years of production. The interbedded nature of this gas hydrate occurrence promotes the development of horizontal dissociation interfaces between gas hydrate-bearing sand and mud layers. As a result, non-uniform gas production along the horizontal interfaces becomes a primary determinant of reservoir performance. Simulation cases have been executed to determine the impact of the uncertainty in in situ reservoir permeability and the manner in which intrinsic permeability dynamically changes during dissociation in response to the imposed effective stress increase. The cases where effective permeability is initially low and compaction is highly sensitive to stress result in the least favorable production predictions.

Journal of Marine and Petroleum Geology↗
USGS2004

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.

Tectonophysics↗
USGS2026

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗
USGS2001

Reconstructing the rise of recent coastal anoxia; molybdenum in Chesapeake Bay sediments

Sporadic, direct observations over a 50 yr period inadequately characterize the history of seasonal hypoxia and anoxia in Chesapeake Bay, alarge estuary threatened by eutrophication. Here, we undertake a reconstruction of 20th century oxygen depletion in this estuary using Mo concentrations in 210Pb-dated sediments; Cu concentrations are used to control for anthropogenic influences. Cores from the central channel display mild Mo enrichments above crustal backgrounds (up to 5 ??g/g) and strong Cu enrichments (up to 35 ??g/g). Temporally, Cu enrichment (mostly anthropogenic) began earlier and stabilized in the last two thirds of the 20th century. In contrast, Mo enrichment has grown during the last two thirds of the century. Molybdenum enrichment is mostly hydrogenic, except in a section of the channel that receives additional Mo from erosion of Early Miocene shore deposits. Two geochemical mechanisms promote Mo enrichment: Manganese refluxing concentrates dissolved MoO24- at the sediment-water interface and sulfide substitution into MoO24- produces thiomolybdates, which can be fixed by particles. The Mo enrichment mechanisms operate primarily during periods when bottom waters are anoxic and thiomolybdate formation can occur near the sediment-water interface. This implies a temporal coupling between water-column anoxia and Mo fixation even though fixation occurs only within sediments. The Mo enrichment profiles suggest that Chesapeake Bay has experienced growing O2 depletion since the first half of the 20th century, but especially after 1960. Copyright ?? 2001 Elsevier Science Ltd.

Geochimica et Cosmochimica Acta↗
USGS1998

New approach to analyzing soil-building systems

A new method of analyzing seismic response of soil-building systems is introduced. The method is based on the discrete-time formulation of wave propagation in layered media for vertically propagating plane shear waves. Buildings are modeled as an extension of the layered soil media by assuming that each story in the building is another layer. The seismic response is expressed in terms of wave travel times between the layers, and the wave reflection and transmission coefficients at layer interfaces. The calculation of the response is reduced to a pair of simple finite-difference equations for each layer, which are solved recursively starting from the bedrock. Compared with commonly used vibration formulation, the wave propagation formulation provides several advantages, including the ability to incorporate soil layers, simplicity of the calculations, improved accuracy in modeling the mass and damping, and better tools for system identification and damage detection.A new method of analyzing seismic response of soil-building systems is introduced. The method is based on the discrete-time formulation of wave propagation in layered media for vertically propagating plane shear waves. Buildings are modeled as an extension of the layered soil media by assuming that each story in the building is another layer. The seismic response is expressed in terms of wave travel times between the layers, and the wave reflection and transmission coefficients at layer interfaces. The calculation of the response is reduced to a pair of simple finite-difference equations for each layer, which are solved recursively starting from the bedrock. Compared with commonly used vibration formulation, the wave propagation formulation provides several advantages, including the ability to incorporate soil layers, simplicity of the calculations, improved accuracy in modeling the mass and damping, and better tools for system identification and damage detection.

Soil Dynamics and Earthquake Engineering↗