SEARCH · Search USGS

Results for “Developments in Earth Surface Processes”

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 235 records · Page 13Linked to original sources

USGS2019

Holocene thermokarst lake dynamics in northern Interior Alaska: The interplay of climate, fire, and subsurface hydrology

The current state of permafrost in Alaska and meaningful expectations for its future evolution are informed by long-term perspectives of previous permafrost degradation. Thermokarst processes in permafrost landscapes often lead to widespread lake formation and the spatial and temporal evolution of thermokarst lake landscapes reflects the combined effects of climate, ground conditions, vegetation, and fire. This study provides detailed analyses of thermokarst lake sediments of Holocene age from the southern loess uplands of the Yukon Flats; including bathymetry and sediment core analyses across a water depth transect. The sediment core results, dated by radiocarbon and 210Pb, indicate the onset of finely laminated lacustrine sedimentation between ~10,000 and 9,000 cal yr BP following basin development through inferred thermokarst processes. Thermokarst expansion to modern shoreline configurations continued until ~5000 cal yr BP, which may have been influenced by increased fire. Between ~5000 and 2000 cal yr BP, the preservation of fine laminations at intermediate and deep-water depths indicate higher lake levels than present. At that time, the lake likely overflowed into an over-deepened gully system that is no longer occupied by perennial streams. By ~2000 cal yr BP, massive sedimentation at intermediate water depths indicates that lake levels lowered, which is interpreted to reflect a response to drier conditions based on correspondence with Yukon Flats regional fire and local paleoclimate reconstructions. Consideration of additional contributing mechanisms include the possible influence of catastrophic lake drainages on downgradient base flow levels that may have enhanced subsurface water loss, although this mechanism is untested. The overall consistency between the millennial lake level trends documented here with regional paleoclimate trends indicates that after lakes formed, their size and depth has likely been affected directly by North Pacific atmospheric circulation changes and indirectly through evolution of permafrost, ground ice and sub-surface hydrology. As the first detailed study of Holocene thermokarst basin expansion, stabilization and subsequent climate-driven lake level variations in a loess upland, results provide a framework for future investigations of paleoclimatic signals from similar lake systems that characterize large regions of Alaska and Siberia.

Alaska↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗

USGS2024

Subsurface redox interactions regulate ebullitive methane flux in heterogeneous Mississippi River deltaic wetland

As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive-transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient-driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

Journal of Advances in Modeling Earth Systems↗

USGS2016

Volcano dome dynamics at Mount St. Helens: Deformation and intermittent subsidence monitored by seismicity and camera imagery pixel offsets

The surface deformation field measured at volcanic domes provides insights into the effects of magmatic processes, gravity- and gas-driven processes, and the development and distribution of internal dome structures. Here we study short-term dome deformation associated with earthquakes at Mount St. Helens, recorded by a permanent optical camera and seismic monitoring network. We use Digital Image Correlation (DIC) to compute the displacement field between successive images and compare the results to the occurrence and characteristics of seismic events during a 6 week period of dome growth in 2006. The results reveal that dome growth at Mount St. Helens was repeatedly interrupted by short-term meter-scale downward displacements at the dome surface, which were associated in time with low-frequency, large-magnitude seismic events followed by a tremor-like signal. The tremor was only recorded by the seismic stations closest to the dome. We find a correlation between the magnitudes of the camera-derived displacements and the spectral amplitudes of the associated tremor. We use the DIC results from two cameras and a high-resolution topographic model to derive full 3-D displacement maps, which reveals internal dome structures and the effect of the seismic activity on daily surface velocities. We postulate that the tremor is recording the gravity-driven response of the upper dome due to mechanical collapse or depressurization and fault-controlled slumping. Our results highlight the different scales and structural expressions during growth and disintegration of lava domes and the relationships between seismic and deformation signals.

Washington↗

USGS2012

Displacement fields from point cloud data: Application of particle imaging velocimetry to landslide geodesy

Acquiring spatially continuous ground-surface displacement fields from Terrestrial Laser Scanners (TLS) will allow better understanding of the physical processes governing landslide motion at detailed spatial and temporal scales. Problems arise, however, when estimating continuous displacement fields from TLS point-clouds because reflecting points from sequential scans of moving ground are not defined uniquely, thus repeat TLS surveys typically do not track individual reflectors. Here, we implemented the cross-correlation-based Particle Image Velocimetry (PIV) method to derive a surface deformation field using TLS point-cloud data. We estimated associated errors using the shape of the cross-correlation function and tested the method's performance with synthetic displacements applied to a TLS point cloud. We applied the method to the toe of the episodically active Cleveland Corral Landslide in northern California using TLS data acquired in June 2005–January 2007 and January–May 2010. Estimated displacements ranged from decimeters to several meters and they agreed well with independent measurements at better than 9% root mean squared (RMS) error. For each of the time periods, the method provided a smooth, nearly continuous displacement field that coincides with independently mapped boundaries of the slide and permits further kinematic and mechanical inference. For the 2010 data set, for instance, the PIV-derived displacement field identified a diffuse zone of displacement that preceded by over a month the development of a new lateral shear zone. Additionally, the upslope and downslope displacement gradients delineated by the dense PIV field elucidated the non-rigid behavior of the slide.

California↗

USGS2023

Steady-state forms of channel profiles shaped by debris flow and fluvial processes

Debris flows regularly traverse bedrock channels that dissect steep landscapes, but our understanding of bedrock erosion by debris flows and their impact on steepland morphology is still rudimentary. Quantitative models of steep bedrock channel networks are based on geomorphic transport laws designed to represent erosion by water-dominated flows. To quantify the impact of debris flow erosion on steep channel network form, it is first necessary to develop methods to estimate spatial variations in bulk debris flow properties (e.g., flow depth, velocity) throughout the channel network that can be integrated into landscape evolution models. Here, we propose and evaluate two methods to estimate spatial variations in bulk debris flow properties along the length of a channel profile. We incorporate both methods into a model designed to simulate the evolution of longitudinal channel profiles that evolve in response to debris flow and fluvial processes. To explore this model framework, we propose a general family of debris flow erosion laws where erosion rate is a function of debris flow depth and channel slope. Model results indicate that erosion by debris flows can explain the occurrence of a scaling break in the slope–area curve at low-drainage areas and that upper-network channel morphology may be useful for inferring catchment-averaged erosion rates in quasi-steady landscapes. Validating specific forms of a debris flow incision law, however, would require more detailed model–data comparisons in specific landscapes where input parameters and channel morphometry can be better constrained. Results improve our ability to interpret topographic signals within steep channel networks and identify observational targets critical for constraining a debris flow incision law.

Earth Surface Dynamics↗

USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗

USGS2000

A test of a climatic index of dune mobility using measurements from the southwestern United States

The climatic index of dune mobility developed by Lancaster (1988) has been applied to a variety of different environments. The index is, however, untested and unverified. We tested the index by comparison of values of the dune mobility index calculated from climate data with rates of sand transport measured at three stations in Arizona and New Mexico over the period 1985 to 1997. Our results show that changes in measured rates of sand transport closely parallel temporal changes in the dune mobility index. The mobility index is, however, a relatively poor predictor of the magnitude of actual sand transport on a year-to-year basis. This discrepancy is probably due to the fact that sand transport rates at these sites are strongly influenced by vegetation cover, the state of which may lag changes in annual precipitation. There is, however, a good relation between the mean annual mobility index and mean annual rates of sand transport. This indicates that the dune mobility index is a valid predictor of the long-term state of the aeolian system and can be used confidently for the purposes for which it was originally intended.

Arizona, New Mexico↗

USGS2005

Macropolygon morphology, development, and classification on North Panamint and Eureka playas, Death Valley National Park CA

Panamint and Eureka playas, both located within Death Valley National Park, exhibit a host of surficial features including fissures, pits, mounds, and plant-covered ridges, representing topographic highs and lows that vary up to 2 m of relief from the playa surface. Aerial photographs reveal that these linear strands often converge to form polygons, ranging in length from several meters to nearly a kilometer. These features stand out in generally dark contrast to the brighter intervening expanse of flat, plant-free, desiccated mud of the typical playa surface. Ground-truth mapping of playa features with differential GPS (Global Positioning System) was conducted in 1999 (North Panamint Valley) and 2002 (Eureka Valley). High-resolution digital maps reveal that both playas possess macropolygons of similar scale and geometry, and that fissures may be categorized into one of two genetic groups: (1) shore-parallel or playa-interior desiccation and shrinkage; and (2) tectonic-induced cracks. Early investigations of these features in Eureka Valley concluded that their origin may have been related to agricultural activity by paleo-Indian communities. Although human artifacts are abundant at each locale, there is no evidence to support the inference that surface features reported on Eureka Playa are anthropogenic in origin. Our assumptions into the genesis of polygons on playas is based on our fortuitous experience of witnessing a fissure in the process of formation on Panamint Playa after a flash flood (May 1999); our observations revealed a paradox that saturation of the upper playa crusts contributes to the establishment of some desiccation features. Follow-up visits to the same feature over 2 yrs' time are a foundation for insight into the evolution and possible longevity of these features. ?? 2005 Elsevier B.V. All rights reserved.

Earth-Science Reviews↗

USGS2021

Projecting climate dependent coastal flood risk with a hybrid statistical dynamical model

Numerical models for tides, storm surge, and wave runup have demonstrated ability to accurately define spatially varying flood surfaces. However these models are typically too computationally expensive to dynamically simulate the full parameter space of future oceanographic, atmospheric, and hydrologic conditions that will constructively compound in the nearshore to cause both extreme event and nuisance flooding during the 21st century. A surrogate modeling framework of waves, winds, and tides is developed in this study to efficiently predict spatially varying nearshore and estuarine water levels contingent on any combination of offshore forcing conditions. The surrogate models are coupled with a time-dependent stochastic climate emulator that provides efficient downscaling for hypothetical iterations of offshore conditions. Together, the hybrid statistical-dynamical framework can assess present day and future coastal flood risk, including the chronological characteristics of individual flood and wave-induced dune overtopping events and their changes into the future. The framework is demonstrated at Naval Base Coronado in San Diego, CA, utilizing the regional Coastal Storm Modeling System (CoSMoS; composed of Delft3D and XBeach) as the dynamic simulator and Gaussian process regression as the surrogate modeling tool. Validation of the framework uses both in-situ tide gauge observations within San Diego Bay, and a nearshore cross-shore array deployment of pressure sensors in the open beach surf zone. The framework reveals the relative influence of large-scale climate variability on future coastal flood resilience metrics relevant to the management of an open coast artificial berm, as well as the stochastic nature of future total water levels.

California↗

USGS2016

Bayesian estimation of magma supply, storage, and eruption rates using a multiphysical volcano model: Kīlauea Volcano, 2000–2012

Estimating rates of magma supply to the world's volcanoes remains one of the most fundamental aims of volcanology. Yet, supply rates can be difficult to estimate even at well-monitored volcanoes, in part because observations are noisy and are usually considered independently rather than as part of a holistic system. In this work we demonstrate a technique for probabilistically estimating time-variable rates of magma supply to a volcano through probabilistic constraint on storage and eruption rates. This approach utilizes Bayesian joint inversion of diverse datasets using predictions from a multiphysical volcano model, and independent prior information derived from previous geophysical, geochemical, and geological studies. The solution to the inverse problem takes the form of a probability density function which takes into account uncertainties in observations and prior information, and which we sample using a Markov chain Monte Carlo algorithm. Applying the technique to Kīlauea Volcano, we develop a model which relates magma flow rates with deformation of the volcano's surface, sulfur dioxide emission rates, lava flow field volumes, and composition of the volcano's basaltic magma. This model accounts for effects and processes mostly neglected in previous supply rate estimates at Kīlauea, including magma compressibility, loss of sulfur to the hydrothermal system, and potential magma storage in the volcano's deep rift zones. We jointly invert data and prior information to estimate rates of supply, storage, and eruption during three recent quasi-steady-state periods at the volcano. Results shed new light on the time-variability of magma supply to Kīlauea, which we find to have increased by 35–100% between 2001 and 2006 (from 0.11–0.17 to 0.18–0.28 km 3 /yr), before subsequently decreasing to 0.08–0.12 km 3 /yr by 2012. Changes in supply rate directly impact hazard at the volcano, and were largely responsible for an increase in eruption rate of 60–150% between 2001 and 2006, and subsequent decline by as much as 60% by 2012. We also demonstrate the occurrence of temporal changes in the proportion of Kīlauea's magma supply that is stored versus erupted, with the supply “surge” in 2006 associated with increased accumulation of magma at the summit. Finally, we are able to place some constraints on sulfur concentrations in Kīlauea magma and the scrubbing of sulfur by the volcano's hydrothermal system. Multiphysical, Bayesian constraint on magma flow rates may be used to monitor evolving volcanic hazard not just at Kīlauea but at other volcanoes around the world.

Hawai'i↗

USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗

USGS2017

Zirconium and hafnium

Zirconium and hafnium are corrosion-resistant metals that are widely used in the chemical and nuclear industries. Most zirconium is consumed in the form of the main ore mineral zircon (ZrSiO 4 , or as zirconium oxide or other zirconium chemicals. Zirconium and hafnium are both refractory lithophile elements that have nearly identical charge, ionic radii, and ionic potentials. As a result, their geochemical behavior is generally similar. Both elements are classified as incompatible because they have physical and crystallochemical properties that exclude them from the crystal lattices of most rock-forming minerals. Zircon and another, less common, ore mineral, baddeleyite (ZrO 2 ), form primarily as accessory minerals in igneous rocks. The presence and abundance of these ore minerals in igneous rocks are largely controlled by the element concentrations in the magma source and by the processes of melt generation and evolution. The world’s largest primary deposits of zirconium and hafnium are associated with alkaline igneous rocks, and, in one locality on the Kola Peninsula of Murmanskaya Oblast, Russia, baddeleyite is recovered as a byproduct of apatite and magnetite mining. Otherwise, there are few primary igneous deposits of zirconium- and hafnium-bearing minerals with economic value at present. The main ore deposits worldwide are heavy-mineral sands produced by the weathering and erosion of preexisting rocks and the concentration of zircon and other economically important heavy minerals, such as ilmenite and rutile (for titanium), chromite (for chromium), and monazite (for rare-earth elements) in sedimentary systems, particularly in coastal environments. In coastal deposits, heavy-mineral enrichment occurs where sediment is repeatedly reworked by wind, waves, currents, and tidal processes. The resulting heavy-mineral-sand deposits, called placers or paleoplacers, preferentially form at relatively low latitudes on passive continental margins and supply 100 percent of the world’s zircon. Zircon makes up a relatively small percentage of the economic heavy minerals in most deposits and is produced primarily as a byproduct of heavy-mineral-sand mining for titanium minerals. From 2003 to 2012, world zirconium mineral concentrates production increased by more than 40 percent, and Australia and South Africa were the leading producers. Global consumption of zirconium mineral concentrates generally increased during the same time period, largely as a result of increased demand in developing economies in Asia and the Middle East. Global demand weakened in 2012, causing a decrease in world production of zirconium mineral concentrates and delaying the development of several new mining projects. Global consumption is expected to increase in the future, however, as demand from the ceramics, chemicals, and metals industries increases (driven by renewed growth in developing economies) and demand for zirconium and hafnium metal increases (driven by the construction and operation of new nuclear powerplants). The behaviors of zirconium and hafnium in the environment are very similar to one another in that most zirconium- and hafnium-bearing minerals have limited solubility and reactivity. Anthropogenic sources of zirconium, and likely hafnium, are from industrial zirconium-containing byproducts and emissions from the processing of sponge zirconium, and exposure to the general population from these sources is small. Zirconium and hafnium are likely not essential to human health and generally are considered to be of low toxicity to humans. The main exposure risks are associated with industrial inhalation and dermal exposure. Because of the low solubility of zirconium and hafnium, ecological health concerns in the aquatic environment and in soils are minimal. Heavy-mineral-sand mining may lead to increased erosion rates when the mining is managed improperly. In addition, surface mining requires removal of the overlying organic soil layer and produces waste material that includes tailings and slimes. The soil removal and mining activity disturbs the surrounding ecosystem and alters the character of the landscape. Dry mineral separation processes create high amounts of airborne dust, whereas wet mineral separation processes do not. In operations that restore the landscape to pre-mining conditions, the volume of waste and the impact on the landscape may be relatively temporary.

Professional Paper↗

USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗

USGS1985

Description and preliminary testing of the CDSN Seismic Sensor Systems

The China Digital Seismograph Network (CDSN) is being designed and installed to provide the People's Republic of China with the facilities needed to create a national digital database for earthquake research. The CDSN, which is being developed jointly by the PRC State Seismological Bureau and the U.S. Geological Survey, will consist initially of nine digitally-recording seismograph stations, a data management system to be used for compiling network-day tapes, and a depot maintenance center. Data produced by the network will be shared with research scientists throughout the world. A national seismograph network must be designed to support a variety of research objectives. From this standpoint, the choices and tradeoffs involved in specifying signal bandwidth, resolution, and dynamic range are the most important decisions in system design. As in the case of the CDSN, these decisions are made during the selection and design of the seismic sensor system and encoder components. The purpose of this report is to describe the CDSN sensor systems, their important signal characteristics, and the results of preliminary tests that have been performed on the instruments. Four overlapping data bands will be recorded at each station: short period (SP), broadband (BB), long period (LP), and very long period (VLP). Amplitude response curves are illustrated in Figure I. Vertical and horizontal components will be recorded for each data band. The SP and LP channels will be recorded with sufficient sensitivities to resolve earth background noise at seismically quiet sites. The BB channels will have a lower sensitivity and are intended for broadband recording of moderate-to-large body-wave signals and for increasing the effective amplitude range in the short- and long-period bands. The VLP channel does not provide additional spectral coverage at long periods; its purpose is to make use of on-site filtration and decimation to reduce post processing requirements for VLP studies. Early plans also included a triaxial set of low-sensitivity accelerometers for recording strong signals from large local and regional earthquakes. The accelerometers are not being installed; however, they may be added in the future. The short-period signals will be derived from a three-component set of PRC-supplied Model DJ-I SP seismometers and US-supplied SP amplifiers. The seismometers will be installed in surface or shallow subsurface vaults, except at two of the stations where they will be installed in boreholes. The BB, LP, and VLP signals will be derived from Streckeisen STS-1 broadband sensor systems installed in vaults, except at one site where the LP signals only will-be derived from a KS-36000 borehole seismometer installed at a depth of 100 meters. Analog signals will be sampled and quantized by an analog-to-digital converter (ADC) that is part of the recording system. Sampling rates chosen for the CDSN are as follows: * SP 40 samples/second * BB 20 samples/second * LP 1 sample/second * VLP 6 samples/minute The ADC 16-bit data word format makes use of 14 bits to quantize the signal and 2 bits to specify an automatically ranged gain of 1, 8, 32, or 128. This will provide 84 dB of resolution and up to 42 dB of gain ranging for a total opera- ting range of 126 dB peak to peak. Magnetic tape cartridges, each having a capacity of 67 megabytes, will be used for recording the digital data. LP and VLP data will be recorded continu- ously. SP and BB data will be processed through an automatic signal detector of the type described by Murdock and Hutt (1983), and only detected events will be stored on tape. Detection parameters, such as turn-on sensitivity and mini- mum recording duration for the SP and BB channels, will be fully programmable and easily changed. One or more of the data channels may also be recorded on analog recorders. A CDSN recording system was not available at the time that the preliminary tests were performed on the CDSN sensor systems. This did not interfere with the principal goals of the testing which were to determine the best sensor installation techniques and to demonstrate the operational performance of the sensor systems, especially with regard to instrument noise and detection capability. When the CDSN recording system is available, additional tests will be performed to measure distortion levels in the full system and components. Distortion tests have been performed on the Streckeisen STS-1 broadband seismometers, the most critical component from the standpoint of linearity, and the results of the tests are reported by Wielandt and Streckeisen (1982). Calibration accuracy and stability are also important sensor system characteristics that cannot be defined without additional testing. These tests will be performed as the instruments are installed at the stations and periodically during operation.

Open-File Report↗

USGS2023

A model integrating satellite-derived shoreline observations for predicting fine-scale shoreline response to waves and sea-level rise across large coastal regions

Satellite-derived shoreline observations combined with dynamic shoreline models enable fine-scale predictions of coastal change across large spatiotemporal scales. Here, we present a satellite-data-assimilated, “littoral-cell”-based, ensemble Kalman-filter shoreline model to predict coastal change and uncertainty due to waves, sea-level rise (SLR), and other natural and anthropogenic processes. We apply the developed ensemble model to the entire California coastline (approximately 1,760 km), much of which is sparsely monitored with traditional survey methods (e.g., Lidar/GPS). Water-level-corrected, satellite-derived shoreline observations (obtained from the CoastSat toolbox) offer a nearly unbiased representation of in situ surveyed shorelines (e.g., mean sea-level elevation contours) at Ocean Beach, San Francisco. We demonstrate that model calibration with satellite observations during a 20-year hindcast period (1995–2015) provides nearly equivalent model forecast accuracy during a validation period (2015–2020) compared to model calibration with monthly in situ observations at Ocean Beach. When comparing model-predicted shoreline positions to satellite-derived observations, the model achieves an accuracy of <10 m RMSE for nearly half of the entire California coastline for the validation period. The calibrated/validated model is then applied for multi-decadal simulations of shoreline change due to projected wave and sea-level conditions, while holding the model parameters fixed. By 2100, the model estimates that 24%–75% of California's beaches may become completely eroded due to SLR scenarios of 1.0–3.0 m, respectively. The satellite-data-assimilated modeling system presented here is generally applicable to a variety of coastal settings around the world owing to the global coverage of satellite imagery.

JGR Earth Surface↗

USGS2020

Thresholds for post-wildfire debris flows: Insights from the Pinal Fire, Arizona, USA

Wildfire significantly alters the hydrologic properties of a burned area, leading to increases in overland flow, erosion, and the potential for runoff-generated debris flows. The initiation of debris flows in recently burned areas is well-characterized by rainfall intensity-duration (ID) thresholds. However, there is currently a paucity of data quantifying the rainfall intensities required to trigger post-wildfire debris flows, which limits our understanding of how and why rainfall ID thresholds vary in different climatic and geologic settings. In this study, we monitored debris-flow activity following the Pinal Fire in central Arizona, which differs from both a climatic and hydrogeomorphic perspective from other regions in the western U.S. where ID thresholds for post-wildfire debris flows are well-established, namely the Transverse Ranges of southern CA. Since the peak rainfall intensity within a rainstorm may exceed the rainfall intensity required to trigger a debris flow, the development of robust rainfall ID thresholds requires knowledge of the timing of debris flows within rainstorms. Existing post-wildfire debris-flow studies in Arizona only constrain the peak rainfall intensity within debris-flow-producing storms, which may far exceed the intensity that actually triggered the observed debris flow. In this study, we used pressure transducers within 5 burned drainage basins to constrain the timing of debris flows within rainstorms. Rainfall ID thresholds derived here from triggering rainfall intensities are, on average, 22 mm/h lower than ID thresholds derived under the assumption that the triggering intensity is equal to the maximum rainfall intensity recorded during a rainstorm. We then use a hydrologic model to demonstrate that the magnitude of the 15-minute rainfall ID threshold at the Pinal Fire site is associated with the rainfall intensity required to exceed a recently proposed dimensionless discharge threshold for debris-flow initiation. Model results further suggest that previously observed differences in regional ID thresholds between Arizona and the San Gabriel Mountains of southern CA may be attributed, in large part, to differences in the hydraulic properties of burned soils.

Arizona↗

USGS2009

The 2006-2007 Kuril Islands great earthquake sequence

[1] The southwestern half of a ∼500 km long seismic gap in the central Kuril Island arc subduction zone experienced two great earthquakes with extensive preshock and aftershock sequences in late 2006 to early 2007. The nature of seismic coupling in the gap had been uncertain due to the limited historical record of prior large events and the presence of distinctive upper plate, trench and outer rise structures relative to adjacent regions along the arc that have experienced repeated great interplate earthquakes in the last few centuries. The intraplate region seaward of the seismic gap had several shallow compressional events during the preceding decades (notably an M S 7.2 event on 16 March 1963), leading to speculation that the interplate fault was seismically coupled. This issue was partly resolved by failure of the shallow portion of the interplate megathrust in an M W = 8.3 thrust event on 15 November 2006. This event ruptured ∼250 km along the seismic gap, just northeast of the great 1963 Kuril Island (M w = 8.5) earthquake rupture zone. Within minutes of the thrust event, intense earthquake activity commenced beneath the outer wall of the trench seaward of the interplate rupture, with the larger events having normal-faulting mechanisms. An unusual double band of interplate and intraplate aftershocks developed. On 13 January 2007, an M W = 8.1 extensional earthquake ruptured within the Pacific plate beneath the seaward edge of the Kuril trench. This event is the third largest normal-faulting earthquake seaward of a subduction zone on record, and its rupture zone extended to at least 33 km depth and paralleled most of the length of the 2006 rupture. The 13 January 2007 event produced stronger shaking in Japan than the larger thrust event, as a consequence of higher short-period energy radiation from the source. The great event aftershock sequences were dominated by the expected faulting geometries; thrust faulting for the 2006 rupture zone, and normal faulting for the 2007 rupture zone. A large intraplate compressional event occurred on 15 January 2009 (M w = 7.4) near 45 km depth, below the rupture zone of the 2007 event and in the vicinity of the 16 March 1963 compressional event. The fault geometry, rupture process and slip distributions of the two great events are estimated using very broadband teleseismic body and surface wave observations. The occurrence of the thrust event in the shallowest portion of the interplate fault in a region with a paucity of large thrust events at greater depths suggests that the event removed most of the slip deficit on this portion of the interplate fault. This great earthquake doublet demonstrates the heightened seismic hazard posed by induced intraplate faulting following large interplate thrust events. Future seismic failure of the remainder of the seismic gap appears viable, with the northeastern region that has also experienced compressional activity seaward of the megathrust warranting particular attention.

Journal of Geophysical Research B: Solid Earth↗