The effects of 180 years of aging on the physical and seismic properties of partially saturated sands

Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(6), (2021): e2020JB021341, https://doi.org/10.1029/2020JB021341.Constraining how the physical properties and seismic responses of recently deposited sands change with time is important for understanding earthquake site response, subsurface fluid flow, and early stages of lithification. Currently, however, there is no detailed (cm-scale) assessment of how sand's physical properties and associated seismic velocities evolve over the first two centuries after deposition. Here, we integrate sedimentation rates with seismic velocity and sediment physical properties data to assess how the vadose zone sands at Port Royal Beach, Jamaica, change within 180 years after deposition. We show that compressional and shear wave velocities increase with sediment age, whereas porosity, grain size, sorting, mineralogy, and cementation fraction remain relatively unchanged during the same period. Rock physics models (constrained by the measured physical properties) predict constant seismic velocities at all sites regardless of sediment age, though misfits between modeled and observed velocities increase with sediment age. We explain these misfits by proposing that shallow sands undergo microstructural grain reorganization that leads to a more uniform distribution of grain contact forces with time. Our results imply that beach sands undergo a previously undocumented lithification process that occurs before compaction.The Society of Exploration Geophysicists Geoscientists without Borders Grant and the Institute for Earth, Science, and Man at Southern Methodist University partially supported this work

Widespread gas hydrate instability on the upper U.S. Beaufort margin

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 119 (2014): 8594–8609, doi:10.1002/2014JB011290.The most climate-sensitive methane hydrate deposits occur on upper continental slopes at depths close to the minimum pressure and maximum temperature for gas hydrate stability. At these water depths, small perturbations in intermediate ocean water temperatures can lead to gas hydrate dissociation. The Arctic Ocean has experienced more dramatic warming than lower latitudes, but observational data have not been used to study the interplay between upper slope gas hydrates and warming ocean waters. Here we use (a) legacy seismic data that constrain upper slope gas hydrate distributions on the U.S. Beaufort Sea margin, (b) Alaskan North Slope borehole data and offshore thermal gradients determined from gas hydrate stability zone thickness to infer regional heat flow, and (c) 1088 direct measurements to characterize multidecadal intermediate ocean warming in the U.S. Beaufort Sea. Combining these data with a three-dimensional thermal model shows that the observed gas hydrate stability zone is too deep by 100 to 250 m. The disparity can be partially attributed to several processes, but the most important is the reequilibration (thinning) of gas hydrates in response to significant (~0.5°C at 2σ certainty) warming of intermediate ocean temperatures over 39 years in a depth range that brackets the upper slope extent of the gas hydrate stability zone. Even in the absence of additional ocean warming, 0.44 to 2.2 Gt of methane could be released from reequilibrating gas hydrates into the sediments underlying an area of ~5–7.5 × 103 km2 on the U.S. Beaufort Sea upper slope during the next century.This work was supported by the U.S. Department of Energy (DOE), grant DE-FE0010180 to SMU and a USGS-DOE interagency agreement DE-FE0005806.2015-06-0

Triggering mechanism and tsunamogenic potential of the Cape Fear Slide complex, U.S. Atlantic margin

Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 8 (2007): Q12008, doi:10.1029/2007GC001722.Analysis of new multibeam bathymetry data and seismic Chirp data acquired over the Cape Fear Slide complex on the U.S. Atlantic margin suggests that at least 5 major submarine slides have likely occurred there within the past 30,000 years, indicating that repetitive, large-scale mass wasting and associated tsunamis may be more common in this area than previously believed. Gas hydrate deposits and associated free gas as well as salt tectonics have been implicated in previous studies as triggers for the major Cape Fear slide events. Analysis of the interaction of the gas hydrate phase boundary and the various generations of slides indicates that only the most landward slide likely intersected the phase boundary and inferred high gas pressures below it. For much of the region, we believe that displacement along a newly recognized normal fault led to upward migration of salt, oversteepening of slopes, and repeated slope failures. Using new constraints on slide morphology, we develop the first tsunami model for the Cape Fear Slide complex. Our results indicate that if the most seaward Cape Fear slide event occurred today, it could produce waves in excess of 2 m at the present-day 100 m bathymetric contour.Acquisition of new data was funded by NOAA Ocean Exploration grant NA03OAR4600100 to C.R., and we thank the National Science Foundation for contributing to transit costs for the ship

The relationship between eruptive activity, flank collapse, and sea level at volcanic islands: A long-term (>1 Ma) record offshore Montserrat, Lesser Antilles

Hole U1395B, drilled southeast of Montserrat during Integrated Ocean Drilling Program Expedition 340, provides a long (>1 Ma) and detailed record of eruptive and mass-wasting events (>130 discrete events). This record can be used to explore the temporal evolution in volcanic activity and landslides at an arc volcano. Analysis of tephra fall and volcaniclastic turbidite deposits in the drill cores reveals three heightened periods of volcanic activity on the island of Montserrat (?930 ka to ?900 ka, ?810 ka to ?760 ka, and ?190 ka to ?120 ka) that coincide with periods of increased volcano instability and mass-wasting. The youngest of these periods marks the peak in activity at the Soufrière Hills volcano. The largest flank collapse of this volcano (?130 ka) occurred towards the end of this period, and two younger landslides also occurred during a period of relatively elevated volcanism. These three landslides represent the only large (>0.3 km3) flank collapses of the Soufrière Hills edifice, and their timing also coincides with periods of rapid sea-level rise (>5 m/ka). Available age data from other island arc volcanoes suggests a general correlation between the timing of large landslides and periods of rapid sea-level rise, but this is not observed for volcanoes in intra-plate ocean settings. We thus infer that rapid sea-level rise may modulate the timing of collapse at island arc volcanoes, but not in larger ocean-island settings

Permeability and pressure measurements in Lesser Antilles submarine slides: Evidence for pressure-driven slow-slip failure

Recent studies hypothesize that some submarine slides fail via pressure-driven slow-slip deformation. To test this hypothesis, this study derives pore pressures in failed and adjacent unfailed deep marine sediments by integrating rock physics models, physical property measurements on recovered sediment core, and wireline logs. Two drill sites (U1394 and U1399) drilled through interpreted slide debris; a third (U1395) drilled into normal marine sediment. Near-hydrostatic fluid pressure exists in sediments at site U1395. In contrast, results at both sites U1394 and U1399 indicate elevated pore fluid pressures in some sediment. We suggest that high pore pressure at the base of a submarine slide deposit at site U1394 results from slide shearing. High pore pressure exists throughout much of site U1399, and Mohr circle analysis suggests that only slight changes in the stress regime will trigger motion. Consolidation tests and permeability measurements indicate moderately low (~10⁻¹⁶–10⁻¹⁷ m²) permeability and overconsolidation in fine-grained slide debris, implying that these sediments act as seals. Three mechanisms, in isolation or in combination, may produce the observed elevated pore fluid pressures at site U1399: (1) rapid sedimentation, (2) lateral fluid flow, and (3) shearing that causes sediments to contract, increasing pore pressure. Our preferred hypothesis is this third mechanism because it explains both elevated fluid pressure and sediment overconsolidation without requiring high sedimentation rates. Our combined analysis of subsurface pore pressures, drilling data, and regional seismic images indicates that slope failure offshore Martinique is perhaps an ongoing, creep-like process where small stress changes trigger motion

Detecting hydrate and fluid flow from bottom simulating reflector depth anomalies

Methane hydrates, ice-like compounds that consist of water and methane, represent a potentially enormous unconventional methane resource that may play a critical role in climate change and ocean acidification; however, it remains unclear how much hydrate exists. Here, using a newly developed three-dimensional (3-D) thermal technique, we reveal a novel method for detecting and quantifying methane hydrate. The analysis reveals where fluids migrate in three dimensions across a continental margin and is used to quantify hydrate with meter-scale horizontal resolution. Our study, located at Hydrate Ridge, offshore Oregon (United States), suggests that heat flow and hydrate concentrations are coupled and that 3-D thermal analysis can be used to constrain hydrate and fluid flow in 3-D seismic data. Hydrate estimates using this technique are consistent with 1-D drilling results, but reveal large, previously unrecognized swaths of hydrate-rich sediments that have gone undetected due to spatially limited drilling and sampling techniques used in past studies. The 3-D analysis suggests that previous hydrate estimates based on drilling at this site are low by a factor of approximately three

Heat Flow on the U.S. Beaufort Margin, Arctic Ocean: Implications for Ocean Warming, Methane Hydrate Stability, and Regional Tectonics

Abstract Results from the first focused heat flow study on the U.S. Beaufort Margin provide insight into decadal‐scale Arctic Ocean temperature change and raise new questions regarding Beaufort Margin evolution. This study measured heat flow using a 3.5‐m Lister probe at 103 sites oriented along four north‐south transects perpendicular to the ~700‐km long U.S. Beaufort Margin. The new heat flow measurements, corrected both for seasonal ocean temperature fluctuations and bathymetric effects, reveal low average heat flow values (~35 mW/m2) at seafloor depths of 300–900 m below sea level (mbsl) and anomalously high (~80 mW/m2) values at seafloor depths of >1,000 mbsl, near the predicted continent‐ocean transition. Anomalously low heat flow values measured on the upper margin are consistent with previous studies suggesting decadal‐scale ocean temperature warming to ~500 mbsl. Our results, however, indicate this ocean warming likely extends to depths as great at 900 mbsl—400 m deeper than previous studies suggest—implying widespread, ongoing, methane hydrate destabilization across much of the U.S. Beaufort Margin. The cause of the anomalously high heat flow values observed at seafloor depths >1,000 at the continent‐ocean transition is unclear. We suggest three candidate processes: (1) higher heat production and lower thermal conductivity on the margin edge due to the thickest sedimentary cover at the ocean‐continent transition, (2) seaward migrating subsurface advection, and (3) possible fault‐reactivation at the northern boundary of the Alaskan Microplate

Ten years of measurements and modeling of soil temperature changes and their effects on permafrost in Northwestern Alaska

AbstractMultiple studies demonstrate Northwest Alaska and the Alaskan North Slope are warming. Melting permafrost causes surface destabilization and ecological changes. Here, we use thermistors permanently installed in 1996 in a borehole in northwestern Alaska to study past, present, and future ground and subsurface temperature change, and from this, forecast future permafrost degradation in the region. We measure and model Ground Surface Temperature (GST) warming trends for a 10year period using equilibrium Temperature-Depth (TD) measurements from borehole T96-012, located near the Red Dog Mine in northwestern Alaska—part of the Arctic ecosystem where a continuous permafrost layer exists. Temperature measurements from 1996 to 2006 indicate the subsurface has clearly warmed at depths shallower than 70m. Seasonal climate effects are visible in the data to a depth of 30m based on a visible sinusoidal pattern in the TD plots that correlate with season patterns. Using numerical models constrained by thermal conductivity and temperature measurements at the site, we show that steady warming at depths of ~30 to 70m is most likely the direct result of longer term (decadal-scale) surface warming. The analysis indicates the GST in the region is warming at ~0.44±0.05°C/decade, a value consistent with Surface Air Temperature (SAT) warming of ~1.0±0.8°C/decade observed at Red Dog Mine, but with much lower uncertainty. The high annual variability in the SAT signal produces significant uncertainty in SAT trends. The high annual variability is filtered out of the GST signal by the low thermal diffusivity of the subsurface. Comparison of our results to recent permafrost monitoring studies suggests changes in latitude in the polar regions significantly impacts warming rates. North Slope average GST warming is ~0.9±0.5°C/decade, double our observations at RDM, but within error. The RDM warming rate is within the warming variation observed in eastern Alaska, 0.36–0.71°C/decade, which suggests changes in longitude produce a smaller impact but have warming variability likely related to ecosystem, elevation, microclimates, etc. changes. We also forward model future warming by assuming a 1D diffusive heat flow model and incorporating latent heat effects for permafrost melting. Our analysis indicates ~1 to 4m of loss at the upper permafrost boundary, a ~145±100% increase in the active layer thickness by 2055. If warming continues at a constant rate of ~0.44±0.05°C/decade, we estimate the 125m thick zone of permafrost at this site will completely melt by ~2150. Permafrost is expected to melt by ~2200, ~2110, or ~2080, if the rate of warming is altered to 0.25, 0.90, or 2.0°C/decade, respectively, as an array of different climate models suggest. Since our model assumes no advection of heat (a more efficient heat transport mechanism), and no accelerated warming, our current prediction of complete permafrost loss by 2150 may overestimate the residence time of permafrost in this region of Northwest Alaska