SEARCH · USGS Science

Results for “Earth Observer”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,279 records · Page 71Linked to original sources

USGS2003

Transient rheology of the uppermost mantle beneath the Mojave Desert, California

Geodetic data indicate that the M7.1 Hector Mine, California, earthquake was followed by a brief period (a few weeks) of rapid deformation preceding a prolonged phase of slower deformation. We find that the signal contained in continuous and campaign global positioning system data for 2.5 years after the earthquake may be explained with a transient rheology. Quantitative modeling of these data with allowance for transient (linear biviscous) rheology in the lower crust and upper mantle demonstrates that transient rheology in the upper mantle is dominant, its material properties being typified by two characteristic relaxation times ???0.07 and ???2 years. The inferred mantle rheology is a Jeffreys solid in which the transient and steady-state shear moduli are equal. Consideration of a simpler viscoelastic model with a linear univiscous rheology (2 fewer parameters than a biviscous model) shows that it consistently underpredicts the amplitude of the first ???3 months signal, and allowance for a biviscous rheology is significant at the 99.0% confidence level. Another alternative model - deep postseismic afterslip beneath the coseismic rupture - predicts a vertical velocity pattern opposite to the observed pattern at all time periods considered. Despite its plausibility, the advocated biviscous rheology model is non-unique and should be regarded as a viable alternative to the non-linear mantle rheology model for governing postseismic flow beneath the Mojave Desert. Published by Elsevier B.V.

Earth and Planetary Science Letters↗

USGS2009

Microseismicity at the North Anatolian Fault in the Sea of Marmara offshore Istanbul, NW Turkey

The North Anatolian Fault Zone (NAFZ) below the Sea of Marmara forms a “seismic gap” where a major earthquake is expected to occur in the near future. This segment of the fault lies between the 1912 Ganos and 1999 İzmit ruptures and is the only NAFZ segment that has not ruptured since 1766. To monitor the microseismic activity at the main fault branch offshore of Istanbul below the Çınarcık Basin, a permanent seismic array (PIRES) was installed on the two outermost Prince Islands, Yassiada and Sivriada, at a few kilometers distance to the fault. In addition, a temporary network of ocean bottom seismometers was deployed throughout the Çınarcık Basin. Slowness vectors are determined combining waveform cross correlation and P wave polarization. We jointly invert azimuth and traveltime observations for hypocenter determination and apply a bootstrap resampling technique to quantify the location precision. We observe seismicity rates of 20 events per month for M < 2.5 along the basin. The spatial distribution of hypocenters suggests that the two major fault branches bounding the depocenter below the Çınarcık Basin merge to one single master fault below ∼17 km depth. On the basis of a cross-correlation technique we group closely spaced earthquakes and determine composite focal mechanisms implementing recordings of surrounding permanent land stations. Fault plane solutions have a predominant right-lateral strike-slip mechanism, indicating that normal faulting along this part of the NAFZ plays a minor role. Toward the west we observe increasing components of thrust faulting. This supports the model of NW trending, dextral strike-slip motion along the northern and main branch of the NAFZ below the eastern Sea of Marmara.

Istanbul↗

USGS2012

Geodetic analysis of reservoir depletion at the Geyser steam field in northern California

Reservoir depletion at the Geysers from 1974 to 1977 is evident in measured changes in gravity, surface strain, and pore pressure drainage. The drainage area increased about 20%, the maximum gravity decrease was about −120 μGal, and the maximum elevation change was about 6 cm during this period. Since the net mass withdrawal is known, it may be combined with the gravity change to estimate a drainage volume. The maximum drainage volume is 25 km 3 . Because the depth of the caprock and extent of the pore pressure drainage is known, this volume limit implies that no significant mass withdrawal occurs below a depth of 4 km. The ratios of surface elevation changes to horizontal contraction coinciding with the drainage area imply an equant drainage geometry, assuming that fluid production produces negative dilatation. Using the same cylindrical geometry used to model the gravity, negative dilatational strain rates of 4 to 5 × 10 −5 /yr were found to produce the observed surface displacements. The likelihood of boiling in the system to produce steam combined with apparent large bulk moduli (from seismic velocities) and small pore pressure declines suggest that most of the reservoir contraction is due to cooling.

Journal of Geophysical Research Solid Earth↗

USGS2022

Diverse tsunamigenesis triggered by the Hunga Tonga-Hunga Ha’apai eruption

On the evening of 15 January 2022, the Hunga Tonga-Hunga Ha’apai volcano 1 unleashed a violent underwater eruption, blanketing the surrounding land masses in ash and debris . The eruption generated tsunamis observed around the world. An event of this type last occurred in 1883 during the eruption of Krakatau , and thus we have the first observations of a tsunami from a large emergent volcanic eruption captured with modern instrumentation. Here we show that the explosive eruption generated waves through multiple mechanisms, including: (1) air–sea coupling with the initial and powerful shock wave radiating out from the explosion in the immediate vicinity of the eruption; (2) collapse of the water cavity created by the underwater explosion; and (3) air–sea coupling with the air-pressure pulse that circled the Earth several times, leading to a global tsunami. In the near field, tsunami impacts are strongly controlled by the water-cavity source whereas the far-field tsunami, which was unusually persistent, can be largely described by the air-pressure pulse mechanism. Catastrophic damage in some harbours in the far field was averted by just tens of centimetres, implying that a modest sea level rise combined with a future, similar event would lead to a step-function increase in impacts on infrastructure. Piecing together the complexity of this event has broad implications for coastal hazards in similar geophysical settings, suggesting a currently neglected source of global tsunamis.

Nature↗

USGS2010

Upper mantle rheology from GRACE and GPS postseismic deformation after the 2004 Sumatra-Andaman earthquake

Mantle rheology is one of the essential, yet least understood, material properties of our planet, controlling the dynamic processes inside the Earth's mantle and the Earth's response to various forces. With the advent of GRACE satellite gravity, measurements of mass displacements associated with many processes are now available. In the case of mass displacements related to postseismic deformation, these data may provide new constraints on the mantle rheology. We consider the postseismic deformation due to the M w = 9.2 Sumatra 26 December 2004 and M w = 8.7 Nias 28 March 2005 earthquakes. Applying wavelet analyses to enhance those local signals in the GRACE time varying geoids up to September 2007, we detect a clear postseismic gravity signal. We supplement these gravity variations with GPS measurements of postseismic crustal displacements to constrain postseismic relaxation processes throughout the upper mantle. The observed GPS displacements and gravity variations are well explained by a model of viscoelastic relaxation plus a small amount of afterslip at the downdip extension of the coseismically ruptured fault planes. Our model uses a 60 km thick elastic layer above a viscoelastic asthenosphere with Burgers body rheology. The mantle below depth 220 km has a Maxwell rheology. Assuming a low transient viscosity in the 60–220 km depth range, the GRACE data are best explained by a constant steady state viscosity throughout the ductile portion of the upper mantle (e.g., 60–660 km). This suggests that the localization of relatively low viscosity in the asthenosphere is chiefly in the transient viscosity rather than the steady state viscosity. We find a 8.10 18 Pa s mantle viscosity in the 220–660 km depth range. This may indicate a transient response of the upper mantle to the high amount of stress released by the earthquakes. To fit the remaining misfit to the GRACE data, larger at the smaller spatial scales, cumulative afterslip of about 75 cm at depth should be added over the period spanned by the GRACE models. It produces only small crustal displacements. Our results confirm that satellite gravity data are an essential complement to ground geodetic and geophysical networks in order to understand the seismic cycle and the Earth's inner structure.

Andaman Sea, Sumatra↗

USGS2010

Southern high latitude dune fields on Mars: Morphology, aeolian inactivity, and climate change

In a study area spanning the martian surface poleward of 50?? S., 1190 dune fields have been identified, mapped, and categorized based on dune field morphology. Dune fields in the study area span ??? 116400km2, leading to a global dune field coverage estimate of ???904000km2, far less than that found on Earth. Based on distinct morphological features, the dune fields were grouped into six different classes that vary in interpreted aeolian activity level from potentially active to relatively inactive and eroding. The six dune field classes occur in specific latitude zones, with a sequence of reduced activity and degradation progressing poleward. In particular, the first signs of stabilization appear at ???60?? S., which broadly corresponds to the edge of high concentrations of water-equivalent hydrogen content (observed by the Neutron Spectrometer) that have been interpreted as ground ice. This near-surface ground ice likely acts to reduce sand availability in the present climate state on Mars, stabilizing high latitude dunes and allowing erosional processes to change their morphology. As a result, climatic changes in the content of near-surface ground ice are likely to influence the level of dune activity. Spatial variation of dune field classes with longitude is significant, suggesting that local conditions play a major role in determining dune field activity level. Dune fields on the south polar layered terrain, for example, appear either potentially active or inactive, indicating that at least two generations of dune building have occurred on this surface. Many dune fields show signs of degradation mixed with crisp-brinked dunes, also suggesting that more than one generation of dune building has occurred since they originally formed. Dune fields superposed on early and late Amazonian surfaces provide potential upper age limits of ???100My on the south polar layered deposits and ???3Ga elsewhere at high latitudes. No craters are present on any identifiable dune fields, which can provide a lower age limit through crater counting: assuming all relatively stabilized dune fields represent a single noncontiguous surface of uniform age, their estimated crater retention age is <???10000years. An average-sized uncratered dune field (94km2) has a crater retention age <???8My. This apparent youth suggests that present-day climate conditions are responsible for the observed degradation and reduced level of aeolian activity. A lack of observed transport pathways and the absence of large dune fields in the largest basins (Hellas and Argyre Planitiae) are consistent with the previously proposed idea that dune sands are not typically transported far from their source regions on Mars. ?? 2009 Elsevier B.V.

Geomorphology↗

USGS2007

U/Th series radionuclides as coastal groundwater tracers

The study of coastal groundwater has recently surfaced as an active interdisciplinary area of research, driven foremost by its importance as a poorly quantified pathway for subsurface material transport into coastal ecosystems. Key issue in coastal groundwater research include a complete geochemical characterization of the groundwater(s); quantification of the kinetics of subsurface transport, including rock-water interactions; determination of groundwater ages; tracing of groundwater discharge into coastal waters using radiochemical fingerprints; and an assessment of the potential ecological impact of such subsurface flow to a reviving water body. For such applications, the isotopic systemics of select naturally occurring radionucludes in the U/Th series has proven to be particularly useful. These radionuclides (e.g., U, Th, Ram and Rn) are ubiquitous in all groundwaters ad are represented by several isotopes with widely different half-lives and chemistries (Figure 1). As a result, varied biogeochemical processes occurring over a broad range of time scales can be studied. In source rock, most U/Th series isotopes in secular equilibrium; that is, the rate of decay of a daughter isotope is equal to that of it radiogenic parent, and so will have equal activities (in this context, the specific activity is simply a measure of the amount of radioactivity per unit amount). In contrast, these nuclides exhibit strong fractionations within the surrounding groundwaters because of their respective physiochemical differences. Disequilibria in U/Th series radionuclides can thus be used to identify distinct water masses, quantify release rates from source rocks, assess groundwater migration rates, and assess groundwater discharge rates in coastal waters., Large isotopic variations also have the potential for providing precise fingerprints for groundwaters from specific aquifers and have been explored as a means for calculating groundwater ages and estuarine water mass transit times. The highly fractionated nature of U/Th series nuclides in groundwater is illustrated by the range in some measured activities. highest activities are typically observed for 222 Rn, reflecting the inert nature of this noble gas. Groundwater 222 Rn (t 1/2 =3.8) activities are thus controlled only by rapid in situ decay (Table 1) and production within host rocks, without the added complications of reversible removal via absorption or precipitation. Uranium, which is soluble as U(VI) in oxidizing waters, is present in intermediate activities in groundwaters that are moderated by redox-initiated removal onto aquifer rocks. The alkaline earth Ra and, to a greater extent, the less soluble actinide Th are readily removed from groundwater by water -- rock interactions and so are strongly depleted. Both of these elements have very short-lived as well as longer-lived isotopes, and so isotopes compositions reflect processes over a range of time scales. Many studies have evaluated and behavior of select radionuclides in groundwater and surface water systems. Recent advances in high-=precision mass spectrometry have opted new possibilities for more subtle interpretations in select long-lived U/Th series isotopes, such as U, Ra, Pa, and Th. However, these techniques have yet to be fully developed, ahns as a consequence, such data remain largely scarce and underutilized. Although many different approaches have been developed to study radionucluide behavior in groundwater, all are based on principles of radioactive production and decay and knowledge of source terms from weathering and recoil processes, as well as removal terms from the interaction with aquifer host rock surface by sorption and precipitation. This review is structured to present first a brief description of the background, driving forces, scales, and ecological significance of submarine groundwater discharge. Following this, a description of the geochemistry and behavior of select radionuclides in groundwater will be presented, and their application to tracing submarine groundwater discharge will be discussed.

Chemical Reviews↗

USGS2023

Investigating permafrost carbon dynamics in Alaska with artificial intelligence

Positive feedbacks between permafrost degradation and the release of soil carbon into the atmosphere impact land–atmosphere interactions, disrupt the global carbon cycle, and accelerate climate change. The widespread distribution of thawing permafrost is causing a cascade of geophysical and biochemical disturbances with global impacts. Currently, few earth system models account for permafrost carbon feedback (PCF) mechanisms. This research study integrates artificial intelligence (AI) tools and information derived from field-scale surveys across the tundra and boreal landscapes in Alaska. We identify and interpret the permafrost carbon cycling links and feedback sensitivities with GeoCryoAI, a hybridized multimodal deep learning (DL) architecture of stacked convolutionally layered, memory-encoded recurrent neural networks (NN). This framework integrates in-situ measurements and flux tower observations for teacher forcing and model training. Preliminary experiments to quantify, validate, and forecast permafrost degradation and carbon efflux across Alaska demonstrate the fidelity of this data-driven architecture. More specifically, GeoCryoAI logs the ecological memory and effectively learns covariate dynamics while demonstrating an aptitude to simulate and forecast PCF dynamics—active layer thickness (ALT), carbon dioxide flux (CO 2 ), and methane flux (CH 4 )—with high precision and minimal loss (i.e. ALT RMSE : 1.327 cm [1969–2022]; CO 2 RMSE : 0.697 µ molCO 2 m −2 s −1 [2003–2021]; CH 4 RMSE : 0.715 nmolCH 4 m −2 s −1 [2011–2022]). ALT variability is a sensitive harbinger of change, a unique signal characterizing the PCF, and our model is the first characterization of these dynamics across space and time.

Alaska↗

USGS2023

First observation of the ground-state electron-capture of 40K

Potassium-40 is a widespread, naturally occurring isotope whose radioactivity impacts estimated geological ages spanning billions of years, nuclear structure theory, and subatomic rare-event searches—including those for dark matter and neutrinoless double-beta decay. The decays of this long-lived isotope must be precisely known for its use as a geochronometer, and to account for its presence in low-background experiments. There are several known decay modes for potassium-40, but a predicted electron-capture decay directly to the ground state of argon-40 has never been observed. The existence of this decay mode impacts several fields, while theoretical predictions span an order of magnitude. Here we report on the first, successful observation of this rare decay mode, obtained by the KDK (potassium decay) Collaboration using a novel combination of a low-threshold x-ray detector surrounded by a tonne-scale, high-efficiency γ -ray tagger at Oak Ridge National Laboratory. A blinded analysis reveals a distinctly nonzero ratio of intensities of ground-state electron-captures ( I EC 0 ) over excited-state ones ( I EC ∗ ) of I EC 0 / I EC ∗ = 0.0095 stat ± 0.0022 sys ± 0.0010 (68% CL), with the null hypothesis rejected at 4 σ [Stukel et al. , Phys. Rev. Lett. 131 , 052503 (2023) ]. In terms of branching ratio, this unambiguous signal yields I EC 0 = 0.098 % stat ± 0.023 % sys ± 0.010 % , roughly half of the commonly used prediction. This first observation of a third-forbidden unique electron capture improves our understanding of low-energy backgrounds in dark-matter searches and has implications for nuclear-structure calculations. For example, a shell-model based theoretical estimate for the neutrinoless double-beta decay half-life of calcium-48 is increased by a factor of 7 + 3 − 2 . Our nonzero measurement shifts geochronological ages by up to a percent; implications are illustrated for Earth and solar system chronologies.

Physical Review C↗

USGS1996

Application of ground-penetrating radar to investigation of near-surface fault properties in the San Francisco Bay region

In many geologic environments, ground-penetrating radar (GPR) provides high-resolution images of near-surface Earth structure. GPR data collection is nondestructive and very economical. The scale of features detected by GPR lies between those imaged by high-resolution seismic reflection surveys and those exposed in trenches and is therefore potentially complementary to traditional techniques for fault location and mapping. Sixty-two GPR profiles were collected at 12 sites in the San Francisco Bay region. Results show that GPR data correlate with large-scale features in existing trench observations, can be used to locate faults where they are buried or where their positions are not well known, and can identify previously unknown fault segments. The best data acquired were on a profile across the San Andreas fault, traversing Pleistocene terrace deposits south of Olema in Marin County; this profile shows a complicated multi-branched fault system from the ground surface down to about 40 m, the maximum depth for which data were recorded.

California↗

USGS2022

Medium spatial resolution mapping of global land cover and land cover change across multiple decades from Landsat

Land cover maps are essential for characterizing the biophysical properties of the Earth’s land areas. Because land cover information synthesizes a rich array of information related to both the ecological condition of land areas and their exploitation by humans, they are widely used for basic and applied research that requires information related to land surface properties (e.g., terrestrial carbon models, water balance models, weather, and climate models) and are core inputs to models and analyses used by natural resource scientists and land managers. As the Earth’s global population has grown over the last several decades rates of land cover change have increased dramatically, with enormous impacts on ecosystem services (e.g., biodiversity, water supply, carbon sequestration, etc.). Hence, accurate information related to land cover is essential for both managing natural resources and for understanding society’s ecological, biophysical, and resource management footprint. To address the need for high-quality land cover information we are using the global record of Landsat observations to compile annual maps of global land cover from 2001 to 2020 at 30 m spatial resolution. To create these maps we use features derived from time series of Landsat imagery in combination with ancillary geospatial data and a large database of training sites to classify land cover at annual time step. The algorithm that we apply uses temporal segmentation to identify periods with stable land cover that are separated by breakpoints in the time series. Here we provide an overview of the methods and data sets we are using to create global maps of land cover. We describe the algorithms used to create these maps and the core land cover data sets that we are creating through this effort, and we summarize our approach to accuracy assessment. We also present a synthesis of early results and discuss the strengths and weaknesses of our early map products and the challenges that we have encountered in creating global land cover data sets from Landsat. Initial accuracy assessment for North America shows good overall accuracy (77.0 ± 2.0% correctly classified) and 79.8% agreement with the European Space Agency (ESA) WorldCover product. The land cover mapping results we report provide the foundation for robust, repeatable, and accurate mapping of global land cover and land cover change across multiple decades at 30 m spatial resolution from Landsat.

Frontiers in Remote Sensing↗

USGS2017

Reviews and syntheses: Field data to benchmark the carbon cycle models for tropical forests

For more accurate projections of both the global carbon (C) cycle and the changing climate, a critical current need is to improve the representation of tropical forests in Earth system models. Tropical forests exchange more C, energy, and water with the atmosphere than any other class of land ecosystems. Further, tropical-forest C cycling is likely responding to the rapid global warming, intensifying water stress, and increasing atmospheric CO 2 levels. Projections of the future C balance of the tropics vary widely among global models. A current effort of the modeling community, the ILAMB (International Land Model Benchmarking) project, is to compile robust observations that can be used to improve the accuracy and realism of the land models for all major biomes. Our goal with this paper is to identify field observations of tropical-forest ecosystem C stocks and fluxes, and of their long-term trends and climatic and CO 2 sensitivities, that can serve this effort. We propose criteria for reference-level field data from this biome and present a set of documented examples from old-growth lowland tropical forests. We offer these as a starting point towards the goal of a regularly updated consensus set of benchmark field observations of C cycling in tropical forests.

Biogeosciences↗

USGS2014

Time-averaged discharge rate of subaerial lava at Kīlauea Volcano, Hawai‘i, measured from TanDEM-X interferometry: Implications for magma supply and storage during 2011-2013

Differencing digital elevation models (DEMs) derived from TerraSAR add-on for Digital Elevation Measurements (TanDEM-X) synthetic aperture radar imagery provides a measurement of elevation change over time. On the East Rift Zone (EZR) of Kīlauea Volcano, Hawai‘i, the effusion of lava causes changes in topography. When these elevation changes are summed over the area of an active lava flow, it is possible to quantify the volume of lava emplaced at the surface during the time spanned by the TanDEM-X data—a parameter that can be difficult to measure across the entirety of an ~100 km 2 lava flow field using ground-based techniques or optical remote sensing data. Based on the differences between multiple TanDEM-X-derived DEMs collected days to weeks apart, the mean dense-rock equivalent time-averaged discharge rate of lava at Kīlauea between mid-2011 and mid-2013 was approximately 2 m 3 /s, which is about half the long-term average rate over the course of Kīlauea's 1983–present ERZ eruption. This result implies that there was an increase in the proportion of lava stored versus erupted, a decrease in the rate of magma supply to the volcano, or some combination of both during this time period. In addition to constraining the time-averaged discharge rate of lava and the rates of magma supply and storage, topographic change maps derived from space-based TanDEM-X data provide insights into the four-dimensional evolution of Kīlauea's ERZ lava flow field. TanDEM-X data are a valuable complement to other space-, air-, and ground-based observations of eruptive activity at Kīlauea and offer great promise at locations around the world for aiding with monitoring not just volcanic eruptions but any hazardous activity that results in surface change, including landslides, floods, earthquakes, and other natural and anthropogenic processes.

Hawai'i↗

USGS2011

Bedform response to flow variability

Laboratory observations and computational results for the response of bedform fields to rapid variations in discharge are compared and discussed. The simple case considered here begins with a relatively low discharge over a flat bed on which bedforms are initiated, followed by a short high-flow period with double the original discharge, during which the morphology of the bedforms adjusts, followed in turn by a relatively long period of the original low discharge. For the grain size and hydraulic conditions selected, the Froude number remains subcritical during the experiment, and sediment moves predominantly as bedload. Observations show rapid development of quasi-two-dimensional bedforms during the initial period of low flow with increasing wavelength and height over the initial low-flow period. When the flow increases, the bedforms rapidly increase in wavelength and height, as expected from other empirical results. When the flow decreases back to the original discharge, the height of the bedforms quickly decreases in response, but the wavelength decreases much more slowly. Computational results using an unsteady two-dimensional flow model coupled to a disequilibrium bedload transport model for the same conditions simulate the formation and initial growth of the bedforms fairly accurately and also predict an increase in dimensions during the high-flow period. However, the computational model predicts a much slower rate of wavelength increase, and also performs less accurately during the final low-flow period, where the wavelength remains essentially constant, rather than decreasing. In addition, the numerical results show less variability in bedform wavelength and height than the measured values; the bedform shape is also somewhat different. Based on observations, these discrepancies may result from the simplified model for sediment particle step lengths used in the computational approach. Experiments show that the particle step length varies spatially and temporally over the bedforms during the evolution process. Assuming a constant value for the step length neglects the role of flow alterations in the bedload sediment-transport process, which appears to result in predicted bedform wavelength changes smaller than those observed. However, observations also suggest that three-dimensional effects play at least some role in the decrease of bedform wavelength, so incorporating better models for particle hop lengths alone may not be sufficient to improve model predictions.

Earth Surface Processes and Landforms↗

USGS2012

Calving seismicity from iceberg-sea surface interactions

Iceberg calving is known to release substantial seismic energy, but little is known about the specific mechanisms that produce calving icequakes. At Yahtse Glacier, a tidewater glacier on the Gulf of Alaska, we draw upon a local network of seismometers and focus on 80 hours of concurrent, direct observation of the terminus to show that calving is the dominant source of seismicity. To elucidate seismogenic mechanisms, we synchronized video and seismograms to reveal that the majority of seismic energy is produced during iceberg interactions with the sea surface. Icequake peak amplitudes coincide with the emergence of high velocity jets of water and ice from the fjord after the complete submergence of falling icebergs below sea level. These icequakes have dominant frequencies between 1 and 3 Hz. Detachment of an iceberg from the terminus produces comparatively weak seismic waves at frequencies between 5 and 20 Hz. Our observations allow us to suggest that the most powerful sources of calving icequakes at Yahtse Glacier include iceberg-sea surface impact, deceleration under the influence of drag and buoyancy, and cavitation. Numerical simulations of seismogenesis during iceberg-sea surface interactions support our observational evidence. Our new understanding of iceberg-sea surface interactions allows us to reattribute the sources of calving seismicity identified in earlier studies and offer guidance for the future use of seismology in monitoring iceberg calving.

Journal of Geophysical Research F: Earth Surface↗

USGS2006

Submarine landslides in the Santa Barbara Channel as potential tsunami sources

Recent investigations using the Monterey Bay Aquarium Research Institutes (MBARI) Remotely Operated Vehicles (ROVs) "Ventana" and "Tiburon" and interpretation of MBARI's EM 300 30 kHz multibeam bathymetric data show that the northern flank of the Santa Barbara Basin has experienced massive slope failures. Of particular concern is the large (130 km2) Goleta landslide complex located off Coal Oil Point near the town of Goleta, that measures 14.6-km long extending from a depth of 90 m to nearly 574 m deep and is 10.5 km wide. We estimate that approximately 1.75 km3 has been displaced by this slide during the Holocene. This feature is a complex compound submarine landslide that contains both surfical slump blocks and mud flows in three distinct segments. Each segment is composed of a distinct head scarp, down-dropped head block and a slide debris lobe. The debris lobes exhibit hummocky topography in the central areas that appear to result from compression during down slope movement. The toes of the western and eastern lobes are well defined in the multibeam image, whereas the toe of the central lobe is less distinct. Continuous seismic reflection profiles show that many buried slide debris lobes exist and comparison of the deformed reflectors with ODP Drill Site 149, Hole 893 suggest that at least 200 000 years of failure have occurred in the area (Fisher et al., 2005a). Based on our interpretation of the multibeam bathymetry and seismic reflection profiles we modeled the potential tsunami that may have been produced from one of the three surfical lobes of the Goleta slide. This model shows that a 10 m high wave could have run ashore along the cliffs of the Goleta shoreline. Several other smaller (2 km2 and 4 km2) slides are located on the northern flank of the Santa Barbara Basin, both to the west and east of Goleta slide and on the Concepcion fan along the western flank of the basin. One slide, named the Gaviota slide, is 3.8 km2, 2.6 km long and 1.7 km wide. A distinct narrow scar extends from near the eastern head wall of this slide for over 2 km eastward toward the Goleta slide and may represent either an incipient failure or a remnant of a previous failure. Push cores collected within the main head scar of this slide consisted of hydrogen sulfide bearing mud, possibly suggesting active fluid seepage and a vibra-core penetrated ???50 cm of recent sediment overlying colluvium or landslide debris confirming the age of ???300 years as proposed by Lee et al. (2004). However, no seeps or indications of recent movement were observed during our ROV investigation within this narrow head scar indicating that seafloor in the scar is draped with mud. ?? 2006 Author(s). This work is licensed under a Creative Commons License.

Natural Hazards and Earth System Sciences↗

USGS2008

Isotope geochemistry of mercury in source rocks, mineral deposits and spring deposits of the California Coast Ranges, USA

We present here the first study of the isotopic composition of mercury in rocks, ore deposits, and active spring deposits from the California Coast Ranges, a part of Earth's crust with unusually extensive evidence of mercury mobility and enrichment. The Franciscan Complex and Great Valley Sequence, which form the bedrock in the California Coast Ranges, are intruded and overlain by Tertiary volcanic rocks including the Clear Lake Volcanic Sequence. These rocks contain two types of mercury deposits, hot-spring deposits that form at shallow depths (< 300??m) and silica-carbonate deposits that extend to depths of 1000??m. Active springs and geothermal areas continue to precipitate Hg and Au and are modern analogues to the fossil hydrothermal systems preserved in the ore deposits. The Franciscan Complex and Great Valley Sequence contain clastic sedimentary rocks with higher concentrations of mercury than volcanic rocks of the Clear Lake Volcanic Field. Mean mercury isotopic compositions (??202Hg) for all three rock units are similar, although the range of values in Franciscan Complex rocks is greater than in either Great Valley or Clear Lake rocks. Hot spring and silica-carbonate mercury deposits have similar average mercury isotopic compositions that are indistinguishable from averages for the three rock units, although ??202Hg values for the mercury deposits have a greater variance than the country rocks. Precipitates from spring and geothermal waters in the area have similarly large variance and a mean ??202Hg value that is significantly lower than the ore deposits and rocks. These observations indicate that there is little or no isotopic fractionation (< ?? 0.5???) during release of mercury from its source rocks into hydrothermal solutions. Isotopic fractionation does appear to take place during transport and concentration of mercury in deposits, however, especially in their uppermost parts. Boiling of hydrothermal fluids, separation of a mercury-bearing CO2 vapor or reduction and volatilization of Hg(0) in the near-surface environment are likely the most important processes causing the observed Hg isotope fractionation. This should result in the release of mercury with low ??202Hg values into the atmosphere from the top of these hydrothermal systems. Estimates of mass balance suggest that residual Hg reservoirs are not measurably enriched in heavy Hg isotopes as a result of this process because only a small amount of Hg (< 4%) leaves actively ore-forming systems. ?? 2008 Elsevier B.V. All rights reserved.

Earth and Planetary Science Letters↗

USGS2002

Origins of high-frequency scattered waves near PKKP from large aperture seismic array data

This article identifies the likely origin of 1-Hz scattered waves in the vicinity of PKKP by comparing measurements of slowness and onset time to ray-theoretical predictions. The measurements are obtained from slant stacks of Large Aperture Seismic Array (LASA) data from 36 earthquakes and six explosions in the range 30??-116??. Three types of scattered waves explain the main features seen in the stacks, including: P scattered to PKP near the Earth's surface (P.PKP), PKKP scattered near its core-mantle-boundary (CMB) reflection point (PK.KP), and SKKP scattered near its CMB reflection point (SK.KP). The LASA stacks image the amplitude and slowness variations of the scattered waves with time. They also show where these waves can be detected and where they are free from contaminating arrivals. SK.KP waves rise above the noise approximately 100 sec before the onset time of the main SKKP arrival near 113??. Observations of PK.KP span 30??-100??. However, at distances greater than 50?? they suffer from P.PKP contamination. At distances less than 40?? the PK.KP last for about 280 sec. This is approximately 130 sec longer than the maximum ray-theoretical prediction for waves scattered at the CMB, indicating a possible combination of near-surface scattering and contributions from the overlying mantle.

Bulletin of the Seismological Society of America↗