Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Algorithms”

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

Methods for rapidly estimating velocity precision from GNSS time series in the presence of temporal correlation: A new method and comparison of existing methods

Time series of position estimates from Global Navigational Satellite System (GNSS) are used to measure the velocities of points on the surface of the Earth. Along with the velocity estimates, a measure of the precision is needed to assess the quality of the velocity measurement. Here, I evaluate rate uncertainties provided by four different methods that have been applied to geodetic time series. The most rigorous approach uses a data covariance that incorporates a variety of noise processes relevant to geodetic time series but is computationally demanding. Two other approaches are efficient algorithms and are used widely, but both can provide less rigorous estimates of the rate uncertainty. I propose and evaluate a fourth method, which provides estimates of rate uncertainty closer to the rigorous approach but is significantly less computationally demanding. I have evaluated all three methods against the more rigorous method using both simulations and time series from 190 GNSS sites. For data best characterized as having a flicker type noise process, one of the widely used methods overestimates the uncertainty by up to a factor of 2, while the other widely used method underestimates the uncertainty by less than a factor of 2. For a random-walk process, both methods underestimate the rate uncertainty by a factor of 3 to 5.

Journal of Geophysical Research↗
USGS2020

Some experiments in extreme-value statistical modeling of magnetic superstorm intensities

In support of projects for forecasting and mitigating the deleterious effects of extreme space-weather storms, an examination is made of the intensities of magnetic superstorms recorded in the Dst index time series (1957-2016). Modified peak-over-threshold and solar-cycle, block-maximum sampling of the Dst time series are performed to obtain compi-lations of storm-maximum −Dstm intensity values. Lognormal, upper-limit lognormal, generalized Pareto, and generalized extreme-value model distributions are fitted to the−Dstm data using a maximum-likelihood algorithm. All four candidate models provide good representations of the data. Comparisons of the statistical significance and good-ness of fits of the various models gives no clear indication as to which model is best. The statistical models are used to extrapolate to extreme-value intensities, such as would be expected (on average) to occur once per century. An upper-limit lognormal fit to peak-over-threshold −Dstm data above a superstorm threshold of 283 nT gives a 100-year ex-trapolated intensity of 542 nT and a 68% confidence interval (obtained by bootstrap re-sampling) of [466, 583] nT. An upper-limit lognormal fit to solar-cycle, block-maximum−DstBM data gives a 9-solar-cycle (approximately 100-year) extrapolated intensity of 553 nT. The Dst data are found to be insufficient for providing usefully accurate esti-mates of a statistically theoretical upper limit for magnetic storm intensity. Secular change in storm intensities is noted, as is a need for improved estimates of pre-1957 magnetic storm intensities.

Space Weather↗
USGS2021

Extreme-event magnetic storm probabilities derived from rank statistics of historical Dst intensities for solar cycles 14-24

A compilation is made of the largest and second-largest magnetic-storm-maximum intensities, − Dst 1 and − Dst 2 , for solar cycles 14–24 (1902–2016) by sampling Oulu Dcx for cycles 19–24, using published − Dst m values for 4 intense storms in cycles 14, 15, and 18 (1903, 1909, 1921, 1946), and calculating 15 new storm-maximum − Dst m values (reported here) for cycles 14–18. Three different models are fitted to the cycle-ranked − Dst 1 and − Dst 2 values using a maximum-likelihood algorithm: A Gumbel model, an unconstrained Generalized-Extreme-Value model, and a Weibull model constrained to have a physically justified maximum storm intensity of − Dst m = 2500 nT. All three models are good descriptions of the data. Since the best model is not clearly revealed with standard statistical tests, inference is precluded of the source process giving rise to storm-maximum − Dst m values. Of the three candidate models, the constrained Weibull gives the lowest superstorm occurrence probabilities. Using the compiled data and the constrained Weibull model, a once-per-century storm intensity is estimated to be − Dst 1 = 663 nT, with a bootstrap 68% confidence interval of [497, 694] nT. Similarly, the probability that a future storm will have an intensity exceeding that of the March 1989 superstorm, − Dst m > 565 nT, is 0.246 per cycle with a 68% confidence interval of [0.140, 0.311] per cycle. Noting (possibly slight) ambiguity in the rankings of storm intensities, using the same methods, but storms more intense than those identified for cycles 14–16, would yield a higher once-per-century intensity and a higher probability for a − Dst m > 565 nT storm.

Space Weather↗
USGS2021

Estimating and forecasting time-varying groundwater recharge in fractured rock: A state-space formulation with preferential and diffuse flow to the water table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

Water Resources Research↗
USGS2021

Robust earthquake early warning at a fraction of the cost: ASTUTI Costa Rica

We show that a fixed smartphone network can provide robust Earthquake Early Warning for at least two orders of magnitude less cost than scientific-grade networks. Our software and cloud-based data architecture that we have constructed for the Alerta Sismica Temprana Utilizando Teléfonos Inteligentes (ASTUTI; Earthquake Early Warning Utilizing Smartphones) network in Costa Rica is easily scaled and exported. Implementation comprises provisioning and installing modern smartphones in judicious locations. Stand-up time for regionally operational networks can be on the order of days. We evaluated a non-parametric ground-motion detection and alerting strategy that would alert the entire Costa Rican population of any event with a ground motion detection threshold of 0.55–0.65 %g at four neighboring stations. During a 6-month evaluation period ASTUTI detected and alerted on five of 13 earthquakes with M w 4.8–5.3 that caused felt Modified Mercalli Intensity shaking levels of 4.3–6. The system did not produce any false alerts and the undetected events did not produce wide-spread or significant felt shaking. System latencies were less than or similar to scientific-grade latencies. Alerts for all five detected events would have reached the capital city, San Jose, before strong S -wave shaking. This would have afforded time for Drop Cover Hold On actions by most residents. Two of the five alerts were triggered by P -waves suggesting that smartphone-based networks could approach the fastest theoretical EEW performance, especially with future expected improvements in smartphone sensors and processing algorithms.

AGU Advances↗
USGS2021

Evidence for latent crustal fluid injection transients in southern California from long-duration earthquake swarms

Earthquake swarms are manifestations of aseismic driving processes deep in the crust. We examine the spatiotemporal distribution of aseismic processes in Southern California using a 12-years catalog of swarms derived with deep learning algorithms. In a core portion of the plate boundary region, which is not associated with elevated heat flow, we identify 92 long-duration swarms ranging from 6 months to 7 years that constitute 26.4% of the total seismicity. We find that 53% of the swarms exhibit ultra-slow diffusive patterns with propagating backfronts, consistent with expectations for natural fluid injection processes. The chronology of the swarms indicates that the aseismic driving processes were active at all times during 2008–2020. The observations challenge common views about the nature of swarms, which would characterize any one of these sequences as anomalous. The regional prevalence of these sequences suggests that transient fluid injection processes play a key role in crustal fluid transport.

California↗
USGS2021

Clays are not created equal: How clay mineral type affects soil parameterization

Clay minerals dominate the soil colloidal fraction and its specific surface area. Differences among clay mineral types significantly influence their effects on soil hydrological and mechanical behavior. Presently, the soil clay content is used to parameterize soil hydraulic and mechanical properties (SHMP) for land surface models while disregarding the type of clay mineral. This undifferentiated use of clay leads to inconsistent parameterization, particularly between tropical and temperate soils, as shown herein. We capitalize on recent global maps of clay minerals that exhibit strong climatic and spatial segregation of active and inactive clays to consider spatially resolved clay mineral types in SHMP estimation. Clay mineral-informed pedotransfer functions and machine learning algorithms trained with datasets including different clay types and soil structure formation processes improve SHMP representation regionally with broad implications for hydrological and geomechanical Earth surface processes.

Geophysical Research Letters↗
USGS2022

Application of recursive estimation to heat tracing for groundwater/surface-water exchange

We present and demonstrate a recursive-estimation framework to infer groundwater/surface-water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat-transport problem as a state-space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time-varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch-Tung-Striebel Smoother (ERTSS). Whereas the EKF is suited to real-time (“online”) applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near-real time or batch-processing (“offline”) applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field-experimental data and are shown to be efficient and rapid for the estimation of time-varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time-varying groundwater/surface-water exchange.

Water Resources Research↗
USGS2022

Remote sensing of visible dye concentrations during a tracer experiment on a large, turbid river

Understanding dispersion in rivers is critical for numerous applications, such as characterizing larval drift for endangered fish species and responding to spills of hazardous materials. Injecting a visible dye into the river can yield insight on dispersion processes, but conventional field instrumentation yields limited data on variations in dye concentration over time at a few, fixed points. Remote sensing can provide more detailed, spatially distributed information on the dye's motion, but this approach has only been tested in clear-flowing streams. The purpose of this study was to assess the potential of remote sensing to facilitate tracer studies in more turbid rivers. To pursue this objective, we injected Rhodamine WT dye into the Missouri River and collected field spectra from a boat, videos from a small unoccupied aircraft system (sUAS), and orthophotos from an airplane. Applying an optimal band ratio analysis (OBRA) algorithm to the field spectra revealed strong correlations ( R 2 = 0.936) between a spectrally based quantity and in situ concentration measurements. OBRA also performed well for broadband RGB (red, green, blue) images extracted from the sUAS-based videos; the resulting concentration maps were used to produce animations that captured movement of the dye pulse. Spectral mixture analysis of repeat orthophoto coverage yielded relative concentration estimates that provided a synoptic perspective on dispersion of the dye throughout the entire 13.8 km reach over the full 2.5-hr duration of the experiment. The results of this study demonstrate the potential to remotely sense tracer dye concentrations in large, highly turbid rivers.

Missouri↗
USGS2022

Tree rings reveal unmatched 2nd century drought in the Colorado River Basin

The ongoing 22 year drought in the Upper Colorado River Basin (UCRB) has been extremely severe, even in the context of the longest available tree-ring reconstruction of annual flow at Lees Ferry, Arizona, dating back to 762 CE. While many southwestern drought assessments have been limited to the past 1200 years, longer paleorecords of moisture variability do exist for the UCRB. Here, gridded drought-atlas data in the UCRB domain along with naturalized streamflow data from the instrumental period (1906–2021) are used in a K nearest neighbor (KNN) nonparametric algorithm to develop a streamflow reconstruction for the Lees Ferry gage starting in 1 CE. The reconstruction reveals a 2nd century drought unmatched in severity by the current drought or by well-documented medieval period droughts in the UCRB. Although data are sparse, analysis of individual long tree ring records and other paleoclimatic data also support the occurrence of an exceptional 2nd century drought.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2023

Identification of seismo-volcanic regimes at Whakaari/White Island (New Zealand) via systematic tuning of an unsupervised classifier

We present an algorithm based on Self-Organizing Maps (SOM) and k-means clustering to recognize patterns in a continuous 12.5-year tremor time series recorded at Whakaari/White Island volcano, New Zealand (hereafter referred to as Whakaari). The approach is extendable to a variety of volcanic settings through systematic tuning of the classifier. Hyperparameters are evaluated by statistical means, yielding a combination of “ideal” SOM parameters for the given data set. Extending from this, we applied a Kernel Density Estimation approach to automatically detect changes within the observed seismicity. We categorize the Whakaari seismic time series into regimes representing distinct volcano-seismic states during recent unrest episodes at Whakaari (2012/2013, 2016, and 2019). There is a clear separation in classification results between background regimes and those representing elevated levels of unrest. Onset of unrest is detected by the classifier 6 weeks before the August 2012 eruption, and ca. 3.5 months before the December 2019 eruption, respectively. Regime changes are corroborated by changes in commonly monitored tremor proxies as well as with reported volcanic activity. The regimes are hypothesized to represent diverse mechanisms including: system pressurization and depressurization, degassing, and elevated surface activity. Labeling these regimes improves visualization of the 2012/2013 and 2019 unrest and eruptive episodes. The pre-eruptive 2016 unrest showed a contrasting shape and nature of seismic regimes, suggesting differing onset and driving processes. The 2016 episode is proposed to result from rapid destabilization of the shallow hydrothermal system, while rising magmatic gases from new injections of magma better explain the 2012/2013 and 2019 episodes.

Whakaari, White Island↗
USGS2023

Moving Aircraft River Velocimetry (MARV): Framework and proof-of-concept on the Tanana River

Information on velocity fields in rivers is critical for designing infrastructure, modeling contaminant transport, and assessing habitat. Although non-contact approaches to measuring flow velocity are well established, these methods assume a stationary imaging platform. This study eliminates this constraint by introducing a framework for moving aircraft river velocimetry (MARV). The workflow takes as input images acquired from an airplane and involves orthorectification, frame overlap analysis, image enhancement, particle image velocimetry (PIV), and aggregation of the resulting velocity vectors onto a prediction grid. We also use new metrics to quantify the agreement between image-derived and field-measured velocity vectors in terms of both orientation and magnitude. The potential of MARV was evaluated using data from two Alaskan rivers: a large, highly turbid channel and its smaller, clearer tributary. Sediment boil vortices on the mainstem provided natural features trackable via PIV and estimated velocities corresponded closely with field measurements ( R 2 up to 0.911). We compared an exhaustive approach that evaluates overlap for all frame combinations to a simpler rolling window implementation and found that the more efficient algorithm did not compromise accuracy. Sensitivity analysis suggested that the method was robust to window parameterization. Comparing PIV output from different flying heights and imaging systems indicated that larger pixels led to higher accuracy and that a more advanced dual-camera system provided superior performance. Results from the tributary were less encouraging, presumably due to a lack of trackable features in visible images. Testing across a range of rivers is needed to assess the generality of MARV.

Water Resources Research↗
USGS2023

Predicting daily river chlorophyll concentrations at a continental scale

Eutrophication is one of the largest threats to aquatic ecosystems and chlorophyll a measurements are relevant indicators of trophic state and algal abundance. Many studies have modeled chlorophyll a in rivers but model development and testing has largely occurred at individual sites which hampers creating generalized models capable of making broad-scale predictions. To address this gap, we compiled a large data set of chlorophyll a concentrations matched to other water quality, meteorological, and reach characteristic data for a diverse set of 82 streams and rivers across the United States. We used this data set and extreme gradient boosting, a tree-based machine learning algorithm, to predict daily chlorophyll a concentrations. Furthermore, we tested several practical considerations of broad-scale models, such as making predictions at sites not included in model training or the utility of in situ water quality data versus universally available remotely estimated model inputs. Predictions were very strongly correlated to observations when compared against a randomly withheld subset of days; however, the model had lower accuracy when applied to completely novel sites withheld from model training. Turbidity and total nitrogen were the two most important variables for predicting chlorophyll a . Although in situ variables improved modeled estimates and were identified as more important during model interpretation, using only remote inputs still resulted in highly correlated predictions with small bias. Testing a model across many sites allowed for identification of common variables relevant to chlorophyll a and highlighted several challenges for applying data-driven models to new sites or at larger spatial scales.

Water Resources Research↗
USGS2023

Exploring the influence of input feature space on CNN-based geomorphic feature extraction from digital terrain data

Many studies of Earth surface processes and landscape evolution rely on having accurate and extensive data sets of surficial geologic units and landforms. Automated extraction of geomorphic features using deep learning provides an objective way to consistently map landforms over large spatial extents. However, there is no consensus on the optimal input feature space for such analyses. We explore the impact of input feature space for extracting geomorphic features from land surface parameters (LSPs) derived from digital terrain models (DTMs) using convolutional neural network (CNN)-based semantic segmentation deep learning. We compare four input feature space configurations: (a) a three-layer composite consisting of a topographic position index (TPI) calculated using a 50 m radius circular window, square root of topographic slope, and TPI calculated using an annulus with a 2 m inner radius and 10 m outer radius, (b) a single illuminating position hillshade, (c) a multidirectional hillshade, and (d) a slopeshade. We test each feature space input using three deep learning algorithms and four use cases: two with natural features and two with anthropogenic features. The three-layer composite generally provided lower overall losses for the training samples, a higher F1-score for the withheld validation data, and better performance for generalizing to withheld testing data from a new geographic extent. Results suggest that CNN-based deep learning for mapping geomorphic features or landforms from LSPs is sensitive to input feature space. Given the large number of LSPs that can be derived from DTM data and the variety of geomorphic mapping tasks that can be undertaken using CNN-based methods, we argue that additional research focused on feature space considerations is needed and suggest future research directions. We also suggest that the three-layer composite implemented here can offer better performance in comparison to using hillshades or other common terrain visualization surfaces and is, thus, worth considering for different mapping and feature extraction tasks.

Earth and Space Science↗
USGS2024

Climate change will impact surface water extents and dynamics across the central United States

Climate change is projected to impact river, lake, and wetland hydrology, with global implications for the condition and productivity of aquatic ecosystems. We integrated Sentinel-1 and Sentinel-2 based algorithms to track monthly surface water extent (2017–2021) for 32 sites across the central United States (U.S.). Median surface water extent was highly variable across sites, ranging from 3.9% to 45.1% of a site. To account for landscape-based differences (e.g., water storage capacity, land use) in the response of surface water extents to meteorological conditions, individual statistical models were developed for each site. Future changes to climate were defined as the difference between 2006–2025 and 2061–2080 using MACA-CMIP5 (MACAv2-METDATA) Global Circulation Models. Time series of climate change adjusted surface water extents were projected. Annually, 19 of the 32 sites under RCP4.5 and 22 of the 32 sites under RCP8.5 were projected to show an average decline in surface water extent, with drying most consistent across the southeast central, southwest central, and midwest central U.S. Projected declines under surface water dry conditions at these sites suggest greater impacts of drought events are likely in the future. Projected changes were seasonally variable, with the greatest decline in surface water extent expected in summer and fall seasons. In contrast, many north central sites showed a projected increase in surface water in most seasons, relative to the 2017–2021 period, likely attributable to projected increases in winter and spring precipitation exceeding increases in projected temperature.

Earth's Future↗
USGS2024

Incorporating intensity distance attenuation into PLUM ground-motion-based earthquake early warning in the United States: The APPLES configuration

We develop Attenuated ProPagation of Local Earthquake Shaking (APPLES), a new configuration for the United States West Coast version of the Propagation of Local Undamped Motion (PLUM) earthquake early warning (EEW) algorithm that incorporates attenuation into its ground-motion prediction procedures. Under APPLES, instead of using a fixed radius to forward-predict observed peak ground shaking to the area surrounding a seismic station, the forward-predicted intensity at a location depends on the distance from the station using an intensity prediction relationship. We conduct conceptual tests of maximum intensity distribution predictions in APPLES and PLUM using a catalog of ShakeMaps to confirm that the attenuation relationship in APPLES is appropriately modeling shaking distributions for West Coast earthquakes. Then, we run APPLES and PLUM in simulated real-time tests to determine warning time performance. Finally, we compare real-time alert behavior during the 2022 M 6.4 Ferndale, California, earthquake and other recent events. We find that APPLES presents two potential improvements to PLUM by reducing over-alerting during smaller magnitude earthquakes and by increasing warning times in some locations during larger earthquakes. APPLES can produce missed and late alerts in locations that experience shaking intensities close to the level used to issue alerts, so preferred alerting strategies with APPLES would use alert thresholds that are lower than the intensities targeted for EEW alerts. We find alerts using APPLES are also similar to those for the source-based approaches currently used in the ShakeAlert EEW system, which will make APPLES easier to integrate into the system.

California, Oregon, Washington↗
USGS2023

Probabilistic source classification of large tephra producing eruptions using supervised machine learning: An example from the Alaska-Aleutian arc

Alaska contains over 130 volcanoes and volcanic fields that have been active within the last 2 million years. Of these, roughly 90 have erupted during the Holocene, with many characterized by at least one large explosive eruption. These large tephra-producing eruptions (LTPEs) generate orders of magnitude more erupted material than a “typical” arc explosive eruption and distribute ash thousands of kilometers from their source. Because LTPEs occur infrequently, and the proximal explosive deposit record in Alaska is generally limited to the Holocene, we require a method that links distal deposits to a source volcano where the correlative proximal deposits from that eruption are no longer preserved. We present a model that accurately and confidently identifies LTPE volcanic sources in the Alaska-Aleutian arc using only in situ geochemistry. The model is a voting ensemble classifier comprised of six conceptually different machine learning algorithms trained on proximal tephra deposits that have had their source positively identified. We show that incompatible trace element ratios (e.g., Nb/U, Th/La, Rb/Sm) help produce a feature space that contains significantly more variance than one produced by major element concentrations, ultimately creating a model that can achieve high accuracy, precision, and recall on predicted volcanic sources, regardless of the perceived 2D data distribution (i.e., bimodal, uniform, normal) or composition (i.e., andesite, trachyte, rhyolite) of that source. Finally, we apply our model to unidentified distal marine tephra deposits in the region to better understand explosive volcanism in the Alaska-Aleutian arc, specifically its pre-Holocene spatiotemporal distribution.

Alaska↗
USGS2023

VogCast: A framework for modeling volcanic air pollution and its application to the 2022 eruption of Mauna Loa Volcano, Hawai'i

Volcanic activity and the associated gas emissions into the atmosphere often result in adverse air quality conditions and present a hazard to human health and the environment. Building on a decade-long effort to provide operational surface sulfur dioxide and sulfate aerosol forecasts for the State of Hawai'i, we present an air quality modeling framework called VogCast. VogCast is designed to simplify ensemble air quality prediction on a regional scale by linking together multiple state-of-the-art models of meteorology, emissions, and dispersion. The framework is open-source and introduces a new dynamic plume-rise algorithm for distributing pollutants vertically. Using radar and satellite data, we demonstrate that VogCast reasonably captured the mean injection height, the location, and the general envelope of the vog plume during Mauna Loa's 2022 eruption. The results suggest that during the 12-day eruption period model performance varied between days with trade and non-trade wind conditions. Our findings also highlight the importance of sulfur dioxide emission rate and vent parameter inputs for improving forecast accuracy. The broad goal of this work is to better our understanding of vog dispersion and improve air quality prediction for impacted communities.

Hawaii↗