Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Offshore”

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 1,297 records · Page 72Linked to original sources

Coseismic slip distribution of the February 27, 2010 Mw 8.9 Maule, Chile earthquake

[1] Static offsets produced by the February 27, 2010 Mw = 8.8 Maule, Chile earthquake as measured by GPS and InSAR constrain coseismic slip along a section of the Andean megathrust of dimensions 650 km (in length) × 180 km (in width). GPS data have been collected from both campaign and continuous sites sampling both the near-field and far field. ALOS/PALSAR data from several ascending and descending tracks constrain the near-field crustal deformation. Inversions of the geodetic data for distributed slip on the megathrust reveal a pronounced slip maximum of order 15 m at ∼15–25 km depth on the megathrust offshore Lloca, indicating that seismic slip was greatest north of the epicenter of the bilaterally propagating rupture. A secondary slip maximum appears at depth ∼25 km on the megathrust just west of Concepción. Coseismic slip is negligible below 35 km depth. Estimates of the seismic moment based on different datasets and modeling approaches vary from 1.8 to 2.6 × 1022 N m. Our study is the first to model the static displacement field using a layered spherical Earth model, allowing us to incorporate both near-field and far-field static displacements in a consistent manner. The obtained seismic moment of 1.97 × 1022 N m, corresponding to a moment magnitude of 8.8, is similar to that obtained by previous seismic and geodetic inversions.

Geophysical Research Letters↗

Geodetic slip model of the 2011 M9.0 Tohoku earthquake

The three-dimensional crustal displacement field as sampled by GPS is used to determine the coseismic slip of the 2011 M9.0 Tohoku Earthquake. We employ a spherically layered Earth structure and use a combination of onland GPS, out to ∼4000 km from the rupture, and offshore GPS, which samples the high-slip region on the interplate boundary along the Japan trench. Inversion of the displacement field for dip slip, assuming an interplate boundary of variable dip and striking 195°, yields a compact slip maximum of about 33 m located 200 km east of Sendai. The geodetic moment is 4.06 × 10 22 N m, corresponding to M w = 9.0. The area of maximum slip is concentrated at a depth of about 10 km, is updip of the rupture areas of the M ≳ 7 Miyagi-oki earthquakes of 1933, 1936, 1937, and 1978, and roughly coincides with the rupture area of the M7.1 1981 Miyagi-oki earthquake. The overlap of the 2011 slip area with several preceding ruptures suggests that the same asperities may rupture repeatedly with M ≳ 7 events within several decades of one another.

Geophysical Research Letters↗

Holocene aridification of India

Spanning a latitudinal range typical for deserts, the Indian peninsula is fertile instead and sustains over a billion people through monsoonal rains. Despite the strong link between climate and society, our knowledge of the long‐term monsoon variability is incomplete over the Indian subcontinent. Here we reconstruct the Holocene paleoclimate in the core monsoon zone (CMZ) of the Indian peninsula using a sediment core recovered offshore from the mouth of Godavari River. Carbon isotopes of sedimentary leaf waxes provide an integrated and regionally extensive record of the flora in the CMZ and document a gradual increase in aridity‐adapted vegetation from ∼4,000 until 1,700 years ago followed by the persistence of aridity‐adapted plants after that. The oxygen isotopic composition of planktonic foraminifer Globigerinoides ruber detects unprecedented high salinity events in the Bay of Bengal over the last 3,000 years, and especially after 1,700 years ago, which suggest that the CMZ aridification intensified in the late Holocene through a series of sub‐millennial dry episodes. Cultural changes occurred across the Indian subcontinent as the climate became more arid after ∼4,000 years. Sedentary agriculture took hold in the drying central and south India, while the urban Harappan civilization collapsed in the already arid Indus basin. The establishment of a more variable hydroclimate over the last ca. 1,700 years may have led to the rapid proliferation of water‐conservation technology in south India.

Geophysical Research Letters↗

Holocene faulting in the Bellingham forearc basin: Upper-plate deformation at the northern end of the Cascadia subduction zone

The northern Cascadia forearc takes up most of the strain transmitted northward via the Oregon Coast block from the northward-migrating Sierra Nevada block. The north-south contractional strain in the forearc manifests in upper-plate faults active during the Holocene, the northern-most components of which are faults within the Bellingham Basin. The Bellingham Basin is the northern of four basins of the actively deforming northern Cascadia forearc. A set of Holocene faults, Drayton Harbor, Birch Bay, and Sandy Point faults, occur within the Bellingham Basin and can be traced from onshore to offshore using a combination of aeromagnetic lineaments, paleoseismic investigations and scarps identified using LiDAR imagery. With the recognition of such Holocene faults, the northernmost margin of the actively deforming Cascadia forearc extends 60 km north of the previously recognized limit of Holocene forearc deformation. Although to date no Holocene faults are recognized at the northern boundary of the Bellingham Basin, which is 15 km north of the international border, there is no compelling tectonic reason to expect that Holocene faults are limited to south of the international border.

Washington↗

Effects of smectite to illite transformation on the frictional strength and sliding stability of intact marine mudstones

At subduction zones, earthquake nucleation and coseismic slip occur only within a limited depth range, known as the “seismogenic zone”. One leading hypothesis for the upper aseismic-seismic transition is that transformation of smectite to illite at ∼100–150°C triggers a change from rate-strengthening frictional behavior that allows only stable sliding, to rate weakening behavior considered a prerequisite for unstable slip. Previous studies on powdered gouges have shown that changes in clay mineralogy alone are unlikely to control this transition, but associated fabric and cementation developed during diagenesis remain possible candidates. We conducted shearing experiments designed specifically to evaluate this hypothesis, by using intact wafers of mudstone from Ocean Drilling Program Site 1174, offshore SW Japan, which have undergone progressive smectite transformation in situ. We sheared specimens along a sawcut in a triaxial configuration, oriented parallel to bedding, at normal stresses of ∼20–150 MPa and a pore pressure of 1 MPa. During shearing, we conducted velocity-stepping tests to measure the friction rate parameter (a-b). Friction coefficient ranges from 0.28–0.40 and values of (a-b) are uniformly positive; both are independent of clay transformation progress. Our work represents the most direct and comprehensive test of the clay transformation hypothesis to date, and suggests that neither illitization, nor accompanying fabric development and cementation, trigger a transition to unstable frictional behavior. We suggest that strain localization, in combination with precipitation of calcite and quartz, is a viable alternative that is consistent with both field observations and recent conceptual models of a heterogeneous seismogenic zone.

Philippine Sea Plate;Eurasian Plate↗

Minimum distribution of subsea ice-bearing permafrost on the US Beaufort Sea continental shelf

Starting in Late Pleistocene time (~19 ka), sea level rise inundated coastal zones worldwide. On some parts of the present-day circum-Arctic continental shelf, this led to flooding and thawing of formerly subaerial permafrost and probable dissociation of associated gas hydrates. Relict permafrost has never been systematically mapped along the 700-km-long U.S. Beaufort Sea continental shelf and is often assumed to extend to ~120 m water depth, the approximate amount of sea level rise since the Late Pleistocene. Here, 5,000 km of multichannel seismic (MCS) data acquired between 1977 and 1992 were examined for high-velocity (>2.3 km s −1 ) refractions consistent with ice-bearing, coarse-grained sediments. Permafrost refractions were identified along <5% of the tracklines at depths of ~5 to 470 m below the seafloor. The resulting map reveals the minimum extent of subsea ice-bearing permafrost, which does not extend seaward of 30 km offshore or beyond the 20 m isobath.

Geophysical Research Letters↗

Mechanisms of wave‐driven water level variability on reef‐fringed coastlines

Wave‐driven water level variability (and runup at the shoreline) is a significant cause of coastal flooding induced by storms. Wave runup is challenging to predict, particularly along tropical coral reef‐fringed coastlines due to the steep bathymetric profiles and large bottom roughness generated by reef organisms, which can violate assumptions in conventional models applied to open sandy coastlines. To investigate the mechanisms of wave‐driven water level variability on a reef‐fringed coastline, we performed a set of laboratory flume experiments on an along‐shore uniform bathymetric profile with and without bottom roughness. Wave setup and waves at frequencies lower than the incident sea‐swell forcing (infragravity waves) were found to be the dominant components of runup. These infragravity waves were positively correlated with offshore wave groups, signifying they were generated in the surf zone by the oscillation of the breakpoint. On the reef flat and at the shoreline, the low‐frequency waves formed a standing wave pattern with energy concentrated at the natural frequencies of the reef flat, indicating resonant amplification. Roughness elements used in the flume to mimic large reef bottom roughness reduced low frequency motions on the reef flat and reduced wave run up by 30% on average, compared to the runs over a smooth bed. These results provide insight into sea‐swell and infragravity wave transformation and wave setup dynamics on steep‐sloped coastlines, and the effect that future losses of reef bottom roughness may have on coastal flooding along reef‐fringed coasts.

Journal of Geophysical Research C: Oceans↗

Field observations of alongshore runup variability under dissipative conditions in presence of a shoreline sandwave

Video measurements of runup were collected at low tide along several profiles covering an alongshore distance of 500 m. The morphology displayed a complex shape with a shoreline sandwave in the lower beach face of about 250 m long mirrored in the inner sandbar. Wave conditions were stationary and moderate (offshore height of 2 m and peak period of nearly 13 s) but yet dissipative. Runup energy was dominated by infragravity frequencies. Alongshore variations in runup (by a factor up to 3) observed both in the incident and infragravity bands were much higher than reported previously (e.g., Guedes et al., 2012, https://doi.org/10.1016/j.csr.2012.08.022 ; Ruggiero et al., 2004, https://doi.org/10.1029/2003JC002160 ) while the alongshore variations in other environmental parameters (e.g., foreshore beach slope) appear to be much lower. Our data suggest that the beach morphology in the inner surf zone plays a crucial role by inducing rapid and significant modification in the incident wave pattern and the alongshore coherence length scales were consistent with the typical alongshore length scale of the morphology.

Journal of Geophysical Research C: Oceans↗

How often can Earthquake Early Warning systems alert sites with high intensity ground motion?

Although numerous Earthquake Early Warning (EEW) algorithms have been developed we still lack a detailed understanding of how often and under what circumstances useful ground motion alerts can be provided to end-users. Here we analyze the alerting performance of the PLUM, EPIC and FinDer algorithms by running them retrospectively on the seismic strong motion data of the 219 earthquakes in Japan since 1996 that exceeded Modified Mercalli Intensity (MMI) of 4.5 on at least 10 sites (Mw 4.5-9.1). Our analysis suggests that, irrespective of the algorithm, EEW end-users should be prepared that EEW can often but not always provide useful ground motion alerts. A majority of sites with moderate-strong ground motion (MMI 5-6) can generally get at least a few seconds of warning time from all algorithms. If such shaking is caused by a shallow crustal event, around 50% of such sites receive alerts with warning times >5 s. Many sites with severe-extreme ground motion (MMI >=8) can be alerted successfully in the case of very large offshore earthquakes, but less than 20% can be alerted ahead of time if such shaking is caused by a shallow crustal event. Our results provide detailed quantitative insight into the expected alerting performance for EEW algorithms under realistic conditions. The main caveat is that the largest shallow crustal event in our data set has Mw7.0, i.e. the data set does not contain very large strike slip events.

Journal of Geophysical Research↗

Introduction to special issue on gas hydrate in porous media: Linking laboratory and field‐scale phenomena

The proliferation of drilling expeditions focused on characterizing natural gas hydrate as a potential energy resource has spawned widespread interest in gas hydrate reservoir properties and associated porous media phenomena. Between 2017 and 2019, a Special Section of this journal compiled contributed papers elucidating interactions between gas hydrate and sediment based on laboratory, numerical modeling, and field studies. Motivated mostly by field observations in the northern Gulf of Mexico and offshore Japan, several papers focus on the mechanisms for gas hydrate formation and accumulation, particularly with vapor phase gas, not dissolved gas, as the precursor to hydrate. These studies rely on numerical modeling or laboratory experiments using sediment packs or benchtop micromodels. A second focus of the Special Section is the role of fines in inhibiting production of gas from methane hydrate, controlling the distribution of hydrate at a pore scale, and influencing the bulk behavior of seafloor sediments. Other papers fill knowledge gaps related to the physical properties of hydrate-bearing sediments and advance new approaches in coupled thermal-mechanical modeling of these sediments during hydrate dissociation. Finally, one study addresses the long-standing question about the fate of methane hydrate at the molecular level when CO2 is injected into natural reservoirs under hydrate-forming conditions.

Journal of Geophysical Research B: Solid Earth↗

Structural controls on slope failure within the western Santa Barbara Channel based on 2D and 3D seismic imaging

The Santa Barbara Channel, offshore California, contains several submarine landslides and ample evidence for incipient failure. This region hosts active thrust and reverse faults that accommodate several mm/yr of convergence, yet the relationships between tectonic deformation and slope failure remain unclear. We present 3‐D and 2‐D multichannel seismic reflection (MCS) data sets, multibeam bathymetry, and chronostratigraphic constraints to investigate the controls on slope failure. Splay faulting along the North Channel Deformation Trend (NCDT) coincides with a distinct zone of compressional uplift and onlapping of steeply dipping Quaternary strata. The NCDT is spatially correlated with seafloor fissures, and 3‐D seismic analyses reveal an intricate system of en echelon reverse faults that offset sediments younger than ~25 ka. Localized uplift zones are located between faults, one of which underlies the Gaviota landslide headscarp. We observe a direct relationship between slope failure and along‐strike variations in the tectonostratigraphic framework. Based on geophysical properties at Ocean Drilling Program (ODP) Site 893, we predict a trend in compaction and porosity reduction in the basin that drives pore fluids up‐dip, toward the zone of onlap above the NCDT, thus reducing slope stability. This interplay between tectonic, sedimentary, and fluid‐flow processes along the NCDT has created a confluence of preconditioning factors, with Gaviota and Goleta landslides being distinguished from the surrounding slopes by their position above the NCDT. The distribution of seafloor fissures suggests sections of the slope remain unstable and are prone to future landsliding. These results provide insights into the processes and 3‐D feedbacks that lead to slope instability along other convergent margins.

California↗

Commentary: The role of geodetic algorithms for earthquake early warning in Cascadia

The ShakeAlert earthquake early warning (EEW) system issues public alerts in California and will soon extend to Oregon and Washington. The Cascadia subduction zone presents significant new challenges and opportunities for EEW. Initial publications suggested that EEW algorithms based on Global Navigation Satellite System (GNSS) data could provide improved warning for intraslab events and dramatically improved warning for offshore megathrust events, both of which contribute significantly to hazard in Cascadia. We find that some expectations in these publications were unrealistic, and we demonstrate that in general geodetic algorithms would not produce timely warnings for intraslab events nor warning times of two minutes or more for severe shaking from megathrust earthquakes. Nonetheless, lessons from recent earthquakes in Japan and California, for which alerts from seismic algorithms suffered from magnitude saturation and high data latencies, demonstrate the urgent need for rigorous testing of geodetic EEW as a potential complement to seismic EEW.

California, Oregon, Washington↗

Shallow basin structure and attenuation are key to predicting long shaking duration in Los Angeles Basin

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three-dimensional simulations using the CVM-S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM-H do not match the timing of the initial arrivals as well as CVM-S4.26.M01, and the strong late arrivals in the CVM-H simulation travel with an apparent velocity slower than observed. A higher-quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off-azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long-period ground motions in basins.

California↗

Constraints on the geometry of the subducted Gorda Plate with converted phases generated by local earthquakes

The largest slip in great megathrust earthquakes often occurs in the 10–30 km depth range, yet seismic imaging of the material properties in this region has proven difficult. We utilize a dense onshore-offshore passive seismic dataset from the southernmost Cascadia subduction zone where seismicity in the mantle of the subducted Gorda Plate produces S -to- P and P -to- S conversions generated within a few km of the plate interface. These conversions typically occur in the 10–20 km depth range at either the top or bottom of a ∼5 km thick layer with a high Vp/Vs that we infer to be primarily the subducted crust. We use their arrival times and amplitudes to infer the location of the top and bottom of the subducted crust as well as the velocity contrasts across these discontinuities. Comparing with both the Slab1.0 and the updated Slab2 interface models, the Slab2 model is generally consistent with the converted phases, while the Slab1.0 model is 1–2 km deeper in the 2–20 km depth range and ∼6–8 km too deep in the 10–20 km depth range between 40.25°N and 40.4°N. Comparing the amplitudes of the converted phases to synthetics for simplified velocity structures, the amplitude of the converted phases requires models containing a ∼5 km thick zone with at least a ∼10%–20% reduction in S wave velocity. Thus, the plate boundary is likely contained within or at the top of this low velocity zone, which potentially indicates a significant porosity and fluid content within the seismogenic zone.

California↗

The contribution of currents, sea-swell waves, and infragravity waves to suspended-sediment transport across a coral reef-lagoon system.

Coral reefs generate substantial volumes of carbonate sediment, which is redistributed throughout the reef‐lagoon system. However, there is little understanding of the specific processes that transport this sediment produced on the outer portions of coral reefs throughout a reef‐lagoon system. Furthermore, the separate contributions of currents, sea‐swell waves, and infragravity waves to transport, which are all strongly influenced by the presence of a reef, is not fully understood. Here, we show that in reef‐lagoon systems most suspended sediment is transported close to the seabed and can, at times, be suspended higher in the water column during oscillatory flow transitions (i.e., near slack flow) at sea‐swell wave frequencies, and during the peak onshore oscillatory velocity phase at infragravity wave frequencies. While these wave frequencies contribute to the transport of suspended sediment offshore and onshore, respectively, the net flux is small. Mean currents are the primary transport mechanism and responsible for almost 2 orders of magnitude more suspended‐sediment flux than sea‐swell and infragravity waves. Whilst waves may not be the primary mechanism for the transport of sediment, our results suggest they are an important driver of sediment suspension from the seabed, as well as contributing to the partitioning of sediment grain sizes from the reef to the shoreline. As the ocean wave climate changes, sea level rises, and the composition of reef benthic communities change, the relative importance of mean currents, sea‐swell waves, and infragravity waves is likely to change, and this will affect how sediment is redistributed throughout reef‐lagoon systems.

Ningaloo Reef↗

Spatio-temporal analysis of hypoxia in the Central Basin of Lake Erie of North America

We develop a spatio-temporal geostatistical interpolation framework to estimate hypoxia extent (dissolved oxygen [DO] concentrations below 2 mg/L) with data from a network of DO loggers. The framework uses empirical orthogonal functions and Bayesian kriging to identify the spatially varying temporal pattern and estimate the distribution of hypoxia, including estimation uncertainty. A prototype web application is also developed in R. The framework is applied to analyze spatio-temporal dynamics of DO in the central basin of Lake Erie in North America using data sampled from a logger network placed on the lake bottom during the summers of 2014, 2015, and 2016. Cross-validation results demonstrate that the framework is capable of capturing the dynamic nature of bottom hypoxia over offshore areas, but nearshore areas have poor interpolation performance due to the impacts of complex physical processes such as seiche events. The findings showed that in the central basin, hypoxia started to emerge in early August of 2014, while in 2015 and 2016 hypoxia began in July. The peak hypoxia extent occurred in late September 2014, mid-August 2015, and early September 2016. The prediction error of the overall spatial extent of hypoxia was as large as 25% of the interpolation area based on current logger deployment. Based on the cross-validation and interpolation error, we suggest placing more loggers in nearshore areas to reduce prediction error near the margins of the hypoxic zone.

central Lake Erie↗

Isolating detrital and diagenetic signals in magnetic susceptibility records from methane-bearing marine sediments

Volume-dependent magnetic susceptibility (κ) is commonly used for paleoenvironmental reconstructions in both terrestrial and marine sedimentary environments where it reflects a mixed signal between primary deposition and secondary diagenesis. In the marine environment, κ is strongly influenced by the abundance of ferrimagnetic minerals regulated by sediment transport processes. Post-depositional alteration by H 2 S, however, can dissolve titanomagnetite, releasing reactive Fe that promotes pyritization and subsequently decreases κ. Here, we provide a new approach for isolating the detrital signal in κ and identifying intervals of diagenetic alteration of κ driven by organoclastic sulfate reduction (OSR) and the anaerobic oxidation of methane (AOM) in methane-bearing marine sediments offshore India. Using the correlation of a heavy mineral proxy from X-ray fluorescence data (Zr/Rb) and κ in unaltered sediments, we predict the primary detrital κ signal and identify intervals of decreased κ, which correspond to increased total sulfur content. Our approach is a rapid, high-resolution method that can identify overprinted κ resulting from pyritization of titanomagnetite due to H 2 S production in marine sediments. In addition, total organic carbon, total sulfur, and authigenic carbonate δ 13 C measurements indicate that both OSR and AOM can drive the observed κ loss, but AOM drives the greatest decreases in κ. Overall, our approach can enhance paleoenvironmental reconstructions and provide insight into paleo-positions of the sulfate-methane transition zone, past enhancements of OSR or paleo-methane seepage, and the role of detrital iron oxide minerals on the marine sediment sulfur sink, with consequences influencing the development of chemosynthetic biological communities at methane seeps.

Geochemistry, Geophysics, Geosystems↗

Earthquake magnitude distributions on northern Caribbean faults from combinatorial optimization models

On-fault earthquake magnitude distributions are calculated for northern Caribbean faults using estimates of fault slip and regional seismicity parameters. Integer programming, a combinatorial optimization method, is used to determine the optimal spatial arrangement of earthquakes sampled from a truncated Gutenberg-Richter distribution that minimizes the global misfit in slip rates on a complex fault system. Slip rates and their uncertainty on major faults are derived from a previously published GPS block model for the region, with fault traces determined from offshore geophysical mapping and previously published onshore studies. The optimal spatial arrangement of the sampled earthquakes is compared with the 500-year history of earthquake observations. Rupture segmentation of the subduction interface along the Hispaniola-Puerto Rico Trench (PRT) fault and seismic coupling on the PRT fault appear to exert the primary control over this spatial arrangement. Introducing a rupture barrier for the Hispaniola-PRT fault northwest of Mona Passage, based on geophysical and seismicity observations, and assigning a low slip rate of 2 mm/yr on the PRT fault are most consistent with historical earthquakes in the region. The addition of low slip-rate secondary faults as well as segmentation of the Hispaniola and Septentrional strike-slip fault improves the consistency with historical seismicity. An important observation from the modeling is that varying the slip rate on the PRT fault and different segmentation scenarios result in significant changes to the optimal magnitude distribution on faults farther away. In general, optimal on-fault magnitude distributions are more complex and inter-dependent than is typically assumed in probabilistic seismic hazard analysis and probabilistic tsunami hazard analysis.

JGR Solid Earth↗