SEARCH · Geology Reports

Results for “Ecology and Evolution”

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.

1,657 records · Page 8Linked to original sources

USGS2025

Chemical and isotopic fractionation during melt inclusion formation

Melt inclusions are used to study the origin and evolution of magmas. The extent to which they represent equilibrium melt compositions, however, critically hinges on the ratio of crystal growth rate to diffusion rate in melts. If the rate of crystal growth is limited by the supply of nutrients and the dissipation of unwanted components, the trapped melt will be depleted in compatible elements and enriched in incompatible elements. Despite widespread recognition of the potential for melt inclusions to trap disequilibrium compositions, their identification remains a challenge. Here, we demonstrate how stable isotopes provide a solution to this problem. Melt inclusions in basaltic tephra from Kīlauea volcano, Hawai‘i, display Mg isotope fractionation up to 0.49‰ ( 26 Mg/ 24 Mg) relative to bulk rock. These observations indicate that kinetic processes, such as chemical diffusion, are at play, as equilibrium mineral-melt isotope fractionation at high temperatures is an order of magnitude smaller. We propose that the heavy Mg isotopic compositions of the melt inclusions were generated by the faster diffusion and thus preferential incorporation of 24 Mg by the growing olivine, leaving the slower diffusant ( 26 Mg) behind. The incompatible elements (e.g., Ca, Ti, and K) are predicted to display similarly large fractionations when boundary layer effects are significant. These findings show that stable isotopes can be a useful tool in identifying disequilibrium compositions in melt inclusions.

Geology↗

USGS2024

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed↗

USGS2026

The United States Magnetotelluric Array and the National Impedance Map

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

conterminous United States↗

USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗

USGS2024

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington↗

USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗

USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗

USGS2024

Characterization of change in tree cover state and condition over the conterminous United States

Variability in the effects of disturbances and extreme climate events can lead to changes in tree cover over time, including partial or complete loss, with diverse ecological consequences. It is therefore critical to identify in space and time the change processes that lead to tree cover change. Studies of change are often hampered by the lack of data capable of consistently detecting different types of change. Using the Landsat satellite record to create a long time-series of land cover and land cover change, the U.S. Geological Survey Land Change Monitoring Assessment and Projection (LCMAP) project has made an annual time series of land cover across the conterminous United States for the period 1985 to 2018. Multiple LCMAP products analyzed together with map validation reference plots provide a robust basis for understanding tree cover change. In LCMAP (Collection 1.2), annual change detection is based on harmonic model breaks calculated at each Landsat pixel from the Continuous Change Detection and Classification (CCDC) algorithm. The results showed that the majority of CCDC harmonic model breaks (signifying change) indicated partial tree cover loss (associated with management practices such as tree cover thinning) as compared to complete tree cover loss (associated with practices like clearcut harvest or fire disturbance). Substantially fewer occurrences of complete tree cover loss were associated with change in land cover state. The area of annual tree cover change increased after the late 1990s and stayed high for the rest of the study period. The reference data showed that tree harvest dominated across the conterminous United States. The majority of tree cover change occurred in evergreen forests. Large estimates of disturbance-related tree cover change indicated that tree cover loss may have previously been underreported due to omission of partial tree cover loss in prior studies. This has considerable implications for forest carbon accounting along with tracking ecosystem goods and services.

Forests↗

USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗

USGS2025

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts↗

USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗

USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗

USGS2023

Preliminary geologic map of the Southern Santa Rosa Mountains and Borrego Badlands, San Diego County, Southern California

This investigation delineates the geologic framework of an area of 75 square kilometers (km 2 ) located west of the Salton Sea in southern California (fig. 1, on sheet 1). The study area encompasses the south flank of the Santa Rosa Mountains and the eastern part of the Borrego Badlands (sheet 1). In this study area, regionally important stratigraphic and structural elements collectively inform the late Cenozoic geologic evolution of the Anza-Borrego sector of the Salton Trough province. Critical stratigraphic and structural elements in the map area include the following: The well exposed sequence of late Cenozoic, nonmarine sedimentary rocks that filled the Anza-Borrego subbasin (fig. 1) of the Salton Trough; A tectonic boundary that—in the southern Santa Rosa Mountains—separates the sedimentary strata from underlying crystalline rocks of Peninsular Ranges type. This tectonic boundary, named the West Salton Detachment Fault System by Axen and Fletcher (1998), is projected to underlie all late Cenozoic sedimentary strata in the Anza-Borrego subbasin of the Salton Trough; A variety of transpressional, transtensional, and strike-slip structures that have deformed the late Cenozoic sedimentary strata and collectively guided syntectonic and posttectonic depositional events within the Anza-Borrego subbasin of the Salton Trough; and The southeasternmost surface expression of the Clark Fault, a major strand of the dextral San Jacinto Fault Zone. Geologic mapping and analysis for this investigation focused on clarifying geologic relations among these four stratigraphic and structural aspects in the map area.

California↗

USGS2025

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming↗

USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗

USGS2024

California State Waters Map Series—Benthic habitat characterization in the region offshore Humboldt Bay, California

Coastal and Marine Ecological Classification Standard (CMECS) geoform, substrate, and biotic component geographic information system (GIS) products were developed for the California State Waters of northern California in the region offshore of Humboldt Bay. The study was motivated by interest in development of offshore wind-energy capacity and infrastructure in Federal waters offshore. This project, carried out by the U.S. Geological Survey (USGS), resulted in four data releases for individual map blocks that are part of the “California State Waters Map Series”: (1) Offshore of Arcata, (2) Offshore of Eureka, (3) Offshore of the Eel River, and (4) Offshore of Cape Mendocino. The study area consists of 436 square kilometers of multibeam echo sounder (MBES) data acquired by Fugro Pelagos, Inc., in 2007. Towed camera-sled video was acquired in 2009 and 2010 to supervise the classification of the MBES data into habitats, and single channel sparker data were collected to calculate sediment thickness above the transgressive unconformity. Using video observations of habitat as ground truth, derivatives of the MBES data were classified into 3 seafloor character types (hard-rugose, hard-flat, and soft-flat), 26 induration-slope-depth groups, and 15 geoforms. The study area substrate is predominantly soft-flat sediment (mud and fine sand) covering 73.6 percent of the area. Hard-flat substrate areas, predominantly coarse sediment in scour depressions, cover 5.4 percent of the study area. The hard-rugose substrate areas are primarily outcrops of layered sedimentary bedrock and constitute 20.9 percent of the study area. Fifteen geoforms were identified in the analysis. The predominant geoforms mirror the seafloor character results, shelf geoforms, rock outcrop geoforms, and scour depression geoforms. Rock and scour areas are restricted to the southern portion of the study area off Cape Mendocino where uplift has exposed bedrock. On the flat shelf area post-transgressive sediment varies in thickness from 1.7 meters (m) nearshore to 28.1 m offshore.

California↗

USGS2026

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗

USGS2024

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗