Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth in Space”

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 685 records · Page 38Linked to original sources
USGS2020

Calibrating geosynchronous and polar orbiting satellites: Sharing best practices

Earth remote sensing optical satellite systems are often divided into two categories—geosynchronous and sun-synchronous. Geosynchronous systems essentially rotate with the Earth and continuously observe the same region of the Earth. Sun-synchronous systems are generally in a polar orbit and view differing regions of the Earth at the same local time. Although similar in instrument design, there are enough differences in these two types of missions that often the calibration of the instruments can be substantially different. Thus, respective calibration teams develop independent methods and do not interact regularly or often. Yet, there are numerous areas of overlap and much to learn from one another. To address this issue, a panel of experts from both types of systems was convened to discover common areas of concern, areas where improvements can be made, and recommendations for the future. As a result of the panelist’s efforts, a set of eight recommendations were developed. Those that are related to improvements of current technologies include maintaining sun-synchronous orbits (not allowing orbital decay), standardization of spectral bandpasses, and expanded use of well-developed calibration techniques such as deep convective clouds, pseudo invariant calibration sites, and lunar methodologies. New techniques for expanded calibration capability include using geosynchronous instruments as transfer radiometers, continued development of ground-based prelaunch calibration technologies, expansion of RadCalNet, and development of space-based calibration radiometer systems.

Remote Sensing↗
USGS2009

Experimental Advanced Airborne Research Lidar (EAARL) Data Processing Manual

The Experimental Advanced Airborne Research Lidar (EAARL) is an example of a Light Detection and Ranging (Lidar) system that utilizes a blue-green wavelength (532 nanometers) to determine the distance to an object. The distance is determined by recording the travel time of a transmitted pulse at the speed of light (fig. 1). This system uses raster laser scanning with full-waveform (multi-peak) resolving capabilities to measure submerged topography and adjacent coastal land elevations simultaneously (Nayegandhi and others, 2009). This document reviews procedures for the post-processing of EAARL data using the custom-built Airborne Lidar Processing System (ALPS). ALPS software was developed in an open-source programming environment operated on a Linux platform. It has the ability to combine the laser return backscatter digitized at 1-nanosecond intervals with aircraft positioning information. This solution enables the exploration and processing of the EAARL data in an interactive or batch mode. ALPS also includes modules for the creation of bare earth, canopy-top, and submerged topography Digital Elevation Models (DEMs). The EAARL system uses an Earth-centered coordinate and reference system that removes the necessity to reference submerged topography data relative to water level or tide gages (Nayegandhi and others, 2006). The EAARL system can be mounted in an array of small twin-engine aircraft that operate at 300 meters above ground level (AGL) at a speed of 60 meters per second (117 knots). While other systems strive to maximize operational depth limits, EAARL has a narrow transmit beam and receiver field of view (1.5 to 2 milliradians), which improves the depth-measurement accuracy in shallow, clear water but limits the maximum depth to about 1.5 Secchi disk depth (~20 meters) in clear water. The laser transmitter [Continuum EPO-5000 yttrium aluminum garnet (YAG)] produces up to 5,000 short-duration (1.2 nanosecond), low-power (70 microjoules) pulses each second. Each pulse is focused into an illumination area that has a radius of about 20 centimeters on the ground. The pulse-repetition frequency of the EAARL transmitter varies along each across-track scan to produce equal cross-track sample spacing and near uniform density (Nayegandhi and others, 2006). Targets can have varying physical and optical characteristics that cause extreme fluctuations in laser backscatter complexity and signal strength. To accommodate this dynamic range, EAARL has the real-time ability to detect, capture, and automatically adapt to each laser return backscatter. The backscattered energy is collected by an array of four high-speed waveform digitizers connected to an array of four sub-nanosecond photodetectors. Each of the four photodetectors receives a finite range of the returning laser backscatter photons. The most sensitive channel receives 90% of the photons, the least sensitive receives 0.9%, and the middle channel receives 9% (Wright and Brock, 2002). The fourth channel is available for detection but is not currently being utilized. All four channels are digitized simultaneously into 65,536 samples for every laser pulse. Receiver optics consists of a 15-centimeter-diameter dielectric-coated Newtonian telescope, a computer-driven raster scanning mirror oscillating at 12.5 hertz (25 rasters per second), and an array of sub-nanosecond photodetectors. The signal emitted by the pulsed laser transmitter is amplified as backscatter by the optical telescope receiver. The photomultiplier tube (PMT) then converts the optical energy into electrical impulses (Nayegandhi and others, 2006). In addition to the full-waveform resolving laser, the EAARL sensor suite includes a down-looking 70-centimeter-resolution Red-Green-Blue (RGB) digital network camera, a high-resolution color infrared (CIR) multispectral camera (14-centimeter-resolution), two precision dual-frequency kinematic carrier-phase global positioning system (GPS) receivers, and an

Open-File Report↗
USGS2012

Formation of left-lateral fractures within the Summit Ridge Shear Zone, 1989 Loma Prieta, California, Earthquake

The 1989 Loma Prieta, California, earthquake is characterized by the lack of major, throughgoing, coseismic, right-lateral faulting along strands of the San Andreas fault zone in the epicentral area. Instead, throughout the Summit Ridge area there are zones of tension cracks and left-lateral fracture zones oriented about N45°W, that is, roughly parallel to the San Andreas fault in this area. The left-lateral fractures zones are enigmatic because their left-lateral slip is opposite to the right-lateral sense of the relative motion between the Pacific and North American plates. We suggest that the enigmatic fractures can be understood if we assume that coseismic deformation was by right-lateral shear across a broad zone, about 0.5 km wide and 4 km long, beneath Summit Ridge. The shear zone has an orientation of about N70° W and is oblique to the San Andreas fault zone, which has a trend of N45° to 50°W, so that right-lateral shearing, along with some dilation, can account for the orientation of the fracture zones. Using an analog, kinematic model of the Summit Ridge shear zone and the observable geometric parameters (a shear zone about 0.5 km wide, an acute angle of 25° between the tension cracks and the shear-zone walls, left-lateral slip of 5 cm, and spacing of about 100 m in the tension cracks) we estimate that the amount of right-lateral shift across the Summit Ridge shear zone was on the order of 1.4 m. This is comparable to the amount of slip for coseismic faulting at depth, 1.9 to 2.4 m, as computed by several investigators. Thus contrary to most previous reports on the Loma Prieta earthquake, which assert that coseismic, right-lateral ground rupture was restricted to considerable (greater than 4 km) depths in the epicentral area, we find that nearly all the right-lateral offset is represented at the ground surface by the Summit Ridge shear zone.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2012

Analysis of two-color geodimeter measurements of deformation within the Long Valley caldera: June 1983 to October 1985

Line length changes from several baselines in a trilateration network within the Long Valley caldera clearly define a decrease in strain rate from June 1983 through October 1985. The data consist of more than 1600 length measurements on 23 baselines using a two-color geodimeter, which has a precision of 0.2 ppm of the line length. Initial measurements made during the summer of 1983 show extension rates were as high as 6 ppm/yr on some baselines. The rates subsequently decreased steadily, to less than 1 ppm/yr during the summer of 1985. Ih contrast, the strain rates on a few other baselines were less than 0.2 ppm/yr averaged over the 2.3-year interval. A model is constructed using these observations as well as yearly observations of a trilateration network within and near the caldera. The model contains two points of inflation located at 5 and 10 km beneath the resurgent dome of the caldera plus dextral slip on a fault plane within the south moat within an elastic half-space. To mimic the observed time dependence in the two-color geodimeter data, we conclude that the rates of both slip and inflation have steadily decreased over the survey interval. Because the two-color geodimeter network is closer to all three sources than the single-color geodimeter network, these data are important for defining the factor of 2–3 rate decrease of slip and inflation. Comparing our estimates of slip and inflation with similar estimates for the year prior to July 1983, we conclude that the rate decrease is probably a result of postseismic and postinflation relaxation following the earthquake swarm in the south moat in January 1983.

Journal of Geophysical Research Solid Earth↗
USGS2012

Strain accumulation across the central Nevada seismic zone, 1973–1994

Five trilateration networks extending for 280 km along the central Nevada seismic zone (1915 Pleasant Valley, M = 7.3; 1954 Dixie Valley, M = 6.8; 1954 Stillwater, M = 6.8; 1954 Rainbow Mountain, M = 6.6; 1954 Fairview Peak, M = 7.1; and 1932 Cedar Mountain, M = 7.2) have been surveyed 6 times since 1973 to determine deformation along the zone. Within the precision of measurement the deformation appears uniform along the zone and is described by the principal strain rates 0.036±0.008 μstrain/yr N60°W±3° and −0.031±0.008 μstrain/yr N30°E±3°, extension reckoned positive. The observed strain rates are consistent with simple, right‐lateral, tensor shear at the rate of 0.033 μstrain/yr across a shear zone striking N15°W. This central Nevada shear zone appears to be the northward continuation of the eastern California shear zone. The orientation of the strike‐slip and normal‐slip ruptures within the central Nevada seismic zone are consistent with principal stress axes parallel to the measured principal strain rate axes. Space‐based geodetic measurements (very long baseline interferometry) indicate that the relative motion accommodated across the Basin and Range province west of Ely, Nevada, is about 9.1±1.5 mm/yr N16°W±8° (Dixon et al., 1995.) Notice that the right‐lateral shear zone postulated to explain deformation in the central Nevada seismic zone is properly oriented to accommodate that relative motion. However, a 135‐km effective width of the shear zone would be required to accommodate all of the 9.1 mm/yr relative motion at the strain rates observed in the Nevada seismic zone; only about 3 mm/yr of that relative motion is accommodated within the span of the trilateration networks.

California, Nevada↗
USGS2021

The weight of cities: Urbanization effects on Earth’s subsurface

Across the world, people increasingly choose to live in cities. By 2050, 70% of Earth's population will live in large urban areas. Upon considering a large city, questions arise such as, how much does that weigh? What are its effects on the landscape? Does it cause measurable subsidence? Here I calculate the weight of San Francisco Bay region urbanization, where 7.75 million people live at, or near the coast. It is difficult to account for everything that is in a city. I assume that most of the weight is buildings and their contents, which allows the use of base outline and height data to approximate their mass, which is cumulatively 1.6·10 12 kg. I build a series of finite element models to study effects of pressure exerted by the weight distribution. Within the elastic realm, I look at compression, flexure, isostatic compensation, stress change, dilatation, and fluid flow changes. Within the nonlinear realm I show example calculations of primary and secondary settlement of soils under load. The combined modeled subsidence from building loads is at least 5–80 mm, with the largest contributions coming from nonlinear settlement and creep in soils. A general result is closing of pore space and redirection of pore fluids. While the calculated subsidence of the Bay Area is relatively small compared with other sources of elevation change such as pumping and recharge of aquifers, all sources of subsidence are concerning given an expected 200–300 mm sea level rise at San Francisco by the year 2050.

AGU Advances↗
USGS2023

Landslides triggered by the 2016–2017 storm season, eastern San Francisco Bay region, California

Discussion The winter rainy season of 2016–2017 brought abundant rainfall to the State of California and to the San Francisco Bay region. In January and February of 2017, intense rainfall from strong winter storms saturated soils in the region and triggered thousands of shallow landslides. The highest concentration of these landslides was in the eastern part of the bay region, where landslides in the hills east of the Cities of Richmond, Berkeley, Oakland, Hayward, and Fremont (see main map for locations) damaged homes, displaced a major electrical transmission-line tower, and blocked several heavily traveled roadways. This map shows 8,928 landslides manually mapped from rectified high-resolution (0.25-meter [m]) satellite imagery (March 11, 2017, from Google Earth) in a three-dimensional geographic information system (GIS) framework. The map area encompasses approximately 1,050 square kilometers (km 2 ), bounded by the Carquinez Strait and San Francisco Bay to the north and west, respectively, and extending to the Interstate Highway 680 corridor to the south and east. Individual landslides were mapped as polygons, but for ease of display, they are shown here as points denoting the highest elevation of each landslide headscarp. The greatest calculated landslide concentration (measured as the total number of landslides per unit area) exceeded 80 landslides per 0.25 km 2 in the hills east of the City of Berkeley. This zone is illustrated on the map in red, a color which denotes areas of greater than or equal to 30 landslides per 0.25 km 2 . The source area and deposit polygon datasets are available as a data release (Corbett and Collins, 2023). Complementary field investigations at more than 150 landslides of those displayed here indicate that most mapped landslides are shallow (less than 1 m deep) debris slides and debris flows. Although most landslides occurred in undeveloped areas (for example, parks and open space), about 1,400 landslides, or 16 percent of the total number, affected structures or infrastructure (including paved and dirt roads) in some way. The inset figures (figs. 1–5) show photographs taken in January and February 2017 that depict some areas of damage that occurred during and after the storms; the locations of where these photographs were taken are shown on the map. In the San Francisco Bay region, landslides usually take place each year during the winter rainy season, but widespread, intense rainfall events such as those whose effects are depicted here typically occur as a result of considerably above-average precipitation conditions. Although the landslides of early 2017 are not the most damaging to affect the region historically, they are still a potent reminder of the potential for landslide hazards present in the San Francisco Bay region.

California↗
USGS2018

Full-range, solar-reflected hyperspectral microscopy to support earth remote sensing research

Over the past 20 years, hyperspectral microscopy has grown into a robust field of analysis for a number of applications. The visible to near-infrared (VNIR; 400 to 1000 nm) region of the spectrum has demonstrated utility for the characterization of healthy and diseased tissue and of biomolecular indicators at the cellular level. Here, we describe the development of a hyperspectral imaging (HSI) microscope that is aimed at material characterization to complement traditional stand-off, earth remote sensing with hyperspectral sensors. We combine commercial off the shelf technology to build an HSI microscope to collect spectral data with illumination provided by a tunable laser. Hyperspectral imaging microscopy (HIM) facilitates detailed examination of target materials at the subcentimeter spatial scale. The custom-built, laser illumination HSI microscope covers the NIR to shortwave infrared (NIR/SWIR; 900 to 2500 nm) solar-reflected spectral range. It is combined with a separate VNIR sensor (400 to 900 nm) that utilizes quartz–tungsten–halogen lamps for illumination. The combined sensors provide a means to collect <10,000 s of spectra in the full VNIR/SWIR spectral range from both pure substances and precisely engineered linear and nonlinear mixtures. The large abundance of spectra allows for a more detailed understanding of the variability and multivariate probability distributions of spectral signatures. This additional information aids in understanding the variability observed in ground truth spectra collected from portable spectrometers, and it greatly enhances sample description and metadata content. In addition, HIM data cubes can serve as proxies, as “microscenes,” for systems engineering applications such as trade studies for HSI acquired by air- and space-borne sensors.

Journal of Applied Remote Sensing↗
USGS1999

Gravel sediment routing from widespread, low-intensity landscape disturbance, Current River basin, Missouri

During the last 160 years, land-use changes in the Ozarks have had the potential to cause widespread, low-intensity delivery of excess amounts of gravel-sized sediment to stream channels. Previous studies have indicated that this excess gravel bedload is moving in wave-like forms through Ozarks drainage basins. The longitudinal, areal distribution of gravel bars along 160 km of the Current River, Missouri, was evaluated to determine the relative effects of valley-scale controls, tributary basin characteristics, and lagged sediment transport in creating areas of gravel accumulations. The longitudinal distribution of gravel-bar area shows a broad scale wave-like form with increases in gravel-bar area weakly associated with tributary junctions. Secondary peaks of gravel area with 1·8–4·1 km spacing (disturbance reaches) are superimposed on the broad form. Variations in valley width explain some, but not all, of the short-spacing variation in gravel-bar area. Among variables describing tributary drainage basin morphometry, present-day land use and geologic characteristics, only drainage area and road density relate even weakly to gravel-bar areal inventories. A simple, channel network-based sediment routing model shows that many of the features of the observed longitudinal gravel distribution can be replicated by uniform transport of sediment from widespread disturbances through a channel network. These results indicate that lagged sediment transport may have a dominant effect on the synoptic spatial distribution of gravel in Ozarks streams; present-day land uses are only weakly associated with present-day gravel inventories; and valley-scale characteristics have secondary controls on gravel accumulations in disturbance reaches.

Earth Surface Processes and Landforms↗
USGS2005

Structure and regional significance of the Late Permian(?) Sierra Nevada - Death Valley thrust system, east-central California

An imbricate system of north-trending, east-directed thrust faults of late Early Permian to middle Early Triassic (most likely Late Permian) age forms a belt in east-central California extending from the Mount Morrison roof pendant in the eastern Sierra Nevada to Death Valley. Six major thrust faults typically with a spacing of 15-20 km, original dips probably of 25-35??, and stratigraphic throws of 2-5 km compose this structural belt, which we call the Sierra Nevada-Death Valley thrust system. These thrusts presumably merge into a de??collement at depth, perhaps at the contact with crystalline basement, the position of which is unknown. We interpret the deformation that produced these thrusts to have been related to the initiation of convergent plate motion along a southeast-trending continental margin segment probably formed by Pennsylvanian transform truncation. This deformation apparently represents a period of tectonic transition to full-scale convergence and arc magmatism along the continental margin beginning in the Late Triassic in central California. ?? 2005 Elsevier B.V. All rights reserved.

Earth-Science Reviews↗
USGS2019

Spatial and temporal variation in river corridor exchange across a 5th order mountain stream network

Although most field and modeling studies of river corridor exchange have been conducted a scales ranging from 10’s to 100’s of meters; results of these studies are used to predict their ecological and hydrological influences at the scale of river networks. Further complicating prediction, exchanges are expected to vary with hydrologic forcing and the local geomorphic setting. While we desire predictive power, we lack a complete spatiotemporal relationship relating discharge to the variation in geologic setting and hydrologic forcing that are expected across a river basin. Indeed, Wondzell’s (2011) conceptual model predicts systematic variation in river corridor exchange as a function of (1) variation in baseflow over time at a fixed location, (2) variation in discharge with location in the river network, and (3) local geomorphic setting. To test this conceptual model we conducted more than 60 solute tracer studies including a synoptic campaign in the 5th order river network of the H.J. Andrews Experimental Forest (Oregon, USA) and replicate-in-time experiments in four watersheds. We interpret the data using a series of metrics describing river corridor exchange and solute transport, testing for consistent direction and magnitude of relationships relating these metrics to discharge and local geomorphic setting. We confirmed systematic decrease in river corridor exchange space through the river networks, from headwaters to the larger mainstem. However, we did not find systematic variation with changes in discharge through time, nor with local geomorphic setting. While interpretation of our results is complicated by problems with the analytical methods, they are sufficiently robust for us to conclude that space-for-time and time-for-space substitutions are not appropriate in our study system. Finally, we suggest two strategies that will improve the interpretability of tracer test results and help the hyporheic community develop robust data sets that will enable comparisons across multiple sites and/or discharge conditions.

Hydrology and Earth System Sciences↗
USGS2019

Volcano deformation: Insights into magmatic systems

Volcano geodesy is the branch of geodetic science that deals with the changing shapes of volcanoes, whether large or small, deep-seated or surficial. Together with seismicity and volcanic gas flux, deformation of the ground surface can be a key indicator of subsurface conditions and processes at volcanoes—information that not only improves scientific understanding of magmatic systems but also is useful for assessing volcano hazards and mitigating their potential consequences. To take full advantage of such information requires detailed characterization of the deformation field in space and time. Currently, no single geodetic technique is capable of providing both the high spatial and temporal resolution required. However, important advances have been made recently by combining information from real-time in situ sensors such as continuous GPS, strainmeters, and tiltmeters with repeated campaign-style GPS, microgravity, interferometric synthetic aperture radar (InSAR), lidar, and photogrammetric observations. Continuous real-time data constrain the timing but not necessarily the spatial extent and pattern of deformation, whereas InSAR provides detailed spatial information but only at intervals of several days to weeks. Repeated microgravity surveys, when combined with independent measurements of surface height, are uniquely sensitive to changes in subsurface mass distribution and therefore can be used to distinguish among processes driven by magma, hydrous fluids, or gas. Simultaneous analysis of multiple geodetic datasets, especially when guided by information from global volcano databases and numerical simulations of volcanic processes and products, can provide a statistical basis for outcome prediction and serve as a guide for additional observations. In the foreseeable future, interactions among magmatic, tectonic, and hydrothermal systems will be monitored and modeled at regional scale in real time, enabling new insights into Earth’s subsurface environment.

Book chapter↗
USGS2005

Dating floodplain sediments using tree-ring response to burial

Floodplain sediments can be dated precisely based on the change in anatomy of tree rings upon burial. When a stem of tamarisk (Tamarix ramosissima) or sandbar willow (Salix exigua) is buried, subsequent annual rings in the buried section resemble the rings of roots: rings become narrower, vessels within the rings become larger, and transitions between rings become less distinct. We combined observations of these changes with tree-ring counts to determine the year of deposition of sedimentary beds exposed in a 150-m-long trench across the floodplain of the Rio Puerco, a rapidly filling arroyo in New Mexico. This method reliably dated most beds thicker than about 30 cm to within a year of deposition. Floodplain aggradation rates varied dramatically through time and space. Sediment deposition was mostly limited to brief overbank flows occurring every few years. The most rapid deposition occurred on channel-margin levees, which migrated laterally during channel narrowing. At the decadal timescale, the cross-section-average sediment deposition rate was steady, but there was a shift in the spatial pattern of deposition in the 1980s. From 1936 to 1986, sediment deposition occurred by channel narrowing, with little change in elevation of the thalweg. After 1986 sediment deposition occurred by vertical aggradation. From 1936 to 2000 about 27 per cent of the arroyo cross-section filled with sediment. The rate of filling from 1962 to 2000 was 0-8 vertical m/decade or 85 m2/decade. Published in 2005 by John Wiley & Sons, Ltd.

Earth Surface Processes and Landforms↗
USGS2008

Investigation and hazard assessment of the 2003 and 2007 Staircase Falls rock falls, Yosemite National Park, California, USA

Since 1857 more than 600 rock falls, rock slides, debris slides, and debris flows have been documented in Yosemite National Park, with rock falls in Yosemite Valley representing the majority of the events. On 26 December 2003, a rock fall originating from west of Glacier Point sent approximately 200 m 3 of rock debris down a series of joint-controlled ledges to the floor of Yosemite Valley. The debris impacted talus near the base of Staircase Falls, producing fragments of flying rock that struck occupied cabins in Curry Village. Several years later on 9 June 2007, and again on 26 July 2007, smaller rock falls originated from the same source area. The 26 December 2003 event coincided with a severe winter storm and was likely triggered by precipitation and/or frost wedging, but the 9 June and 26 July 2007 events lack recognizable triggering mechanisms. We investigated the geologic and hydrologic factors contributing to the Staircase Falls rock falls, including bedrock lithology, weathering, joint spacing and orientations, and hydrologic processes affecting slope stability. We improved upon previous geomorphic assessment of rock-fall hazards, based on a shadow angle approach, by using STONE, a three-dimensional rock-fall simulation computer program. STONE produced simulated rock-fall runout patterns similar to the mapped extent of the 2003 and 2007 events, allowing us to simulate potential future rock falls from the Staircase Falls detachment area. Observations of recent rock falls, mapping of rock debris, and simulations of rock fall runouts beneath the Staircase Falls detachment area suggest that rock-fall hazard zones extend farther downslope than the extent previously defined by mapped surface talus deposits.

Natural Hazards and Earth System Sciences↗
USGS1982

The mobility of uranium and other elements during alteration of rhyolite ash to montmorillonite: A case study in the Troublesome Formation, Colorado, U.S.A.

An unusual occurrence of juxtaposed glassy and clay-altered ash was sampled to estimate the degree and type of element mobility during alteration of glass to montmorillonite. The results are particularly interesting in that major mobilization of uranium is indicated. Closely spaced samples of glassy and montmorillonitic ash were collected from the same 20-50 cm thick stratigraphic horizon in the Troublesome Formation (Miocene) of northwestern Colorado. Sharp contacts exist between glassy ash and underlying pink montmorillonite and indicate that water-saturated conditions were restricted to basal ash layers. Formation of montmorillonite instead of zeolites suggests that the water was not highly saline or alkaline. Isotopic and chemical analyses of glassy and clay-altered samples indicate the following: 1. (1) Montmorillonite has U concentrations which are only 10-15% of the concentrations in coexisting glass. Similarly depleted elements include Cs, Rb, Na and K. Much smaller depletions of these elements in some glassy samples serve as sensitive indicators of incipient alteration of glass to montmorillonite. 2. (2) Abundances of relatively insoluble elements such as Th, Ta, Hf and Al are slightly higher (5-50%) in clay-altered ash and serve as indicators of the maximum levels of enrichment in residual material. Greater enrichment of elements such as Ca, Mg, Sr, Sc, P, Cr and Co indicate structural incorporation, adsorption, or ion-exchange uptake by clay or secondary hydrous oxides of Fe and Mn. 3. (3) The rare-earth-element patterns and abundances in glass are sufficiently mimicked by detritus-free montmorillonite to document the compositional equivalency of the two. 4. (4) Radioactive equilibrium exists between 238U and its decay products 234U and 230Th. This documents minimal open-system mobility of U within the last ??? 0.3 Ma. ?? 1982.

Chemical Geology↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗
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↗