USGS Science⌕ Search

SEARCH · USGS Science

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 955 records · Page 53Linked to original sources

A random forest approach for bounded outcome variables

Random forests have become an established tool for classication and regres- sion, in particular in high-dimensional settings and in the presence of non-additive predictor-response relationships. For bounded outcome variables restricted to the unit interval, however, classical modeling approaches based on mean squared error loss may severely suer as they do not account for heteroscedasticity in the data. To address this issue, we propose a random forest approach for relating a beta dis- tributed outcome to a set of explanatory variables. Our approach explicitly makes use of the likelihood function of the beta distribution for the selection of splits dur- ing the tree-building procedure. In each iteration of the tree-building algorithm it chooses one explanatory variable in combination with a split point that maximizes the log-likelihood function of the beta distribution with the parameter estimates de- rived from the nodes of the currently built tree. Results of several simulation studies and an application using data from the U.S.A. National Lakes Assessment Survey demonstrate the properties and usefulness of the method, in particular when com- pared to random forest approaches based on mean squared error loss and parametric regression models.

Journal of Computational and Graphical Statistics↗

Multiple agricultural cropland products of South Asia developed using Landsat-8 30 m and MODIS 250 m data using machine learning on the Google Earth Engine (GEE) cloud and spectral matching techniques (SMTs) in support of food and water security

Cropland products are of great importance in water and food security assessments, especially in South Asia, which is home to nearly 2 billion people and 230 million hectares of net cropland area. In South Asia, croplands account for about 90% of all human water use. Cropland extent, cropping intensity, crop watering methods, and crop types are important factors that have a bearing on the quantity, quality, and location of production. Currently, cropland products are produced using mainly coarse-resolution (250–1000 m) remote sensing data. As multiple cropland products are needed to address food and water security challenges, our study was aimed at producing three distinct products that would be useful overall in South Asia. The first of these, Product 1, was meant to assess irrigated versus rainfed croplands in South Asia using Landsat 30 m data on the Google Earth Engine (GEE) platform. The second, Product 2, was tailored for major crop types using Moderate Resolution Imaging Spectroradiometer (MODIS) 250 m data. The third, Product 3, was designed for cropping intensity (single, double, and triple cropping) using MODIS 250 m data. For the kharif season (the main cropping season in South Asia, Jun–Oct), 10 major crops (5 irrigated crops: rice, soybean, maize, sugarcane, cotton; and 5 rainfed crops: pulses, rice, sorghum, millet, groundnut) were mapped. For the rabi season (post-rainy season, Nov–Feb), five major crops (three irrigated crops: rice, wheat, maize; and two rainfed crops: chickpea, pulses) were mapped. The irrigated versus rainfed 30 m product showed an overall accuracy of 79.8% with the irrigated cropland class providing a producer’s accuracy of 79% and the rainfed cropland class 74%. The overall accuracy demonstrated by the cropping intensity product was 85.3% with the producer’s accuracies of 88%, 85%, and 67% for single, double, and triple cropping, respectively. Crop types were mapped to accuracy levels ranging from 72% to 97%. A comparison of the crop-type area statistics with national statistics explained 63–98% variability. The study produced multiple-cropland products that are crucial for food and water security assessments, modeling, mapping, and monitoring using multiple-satellite sensor big-data, and Random Forest (RF) machine learning algorithms by coding, processing, and computing on the GEE cloud.

GIScience & Remote Sensing↗

Riverbank stability assessment of lower Osage River during hydro-peak flow events: The lower Osage River case (Missouri, USA)

The fluctuation of water level downstream from dams due to hydropower flow releases can negatively affect riverbank stability. This research aims to investigate riverbank instability resulting from the outflow variation of hydropower plants, using Bagnell Dam and the lower Osage River (Missouri, USA) as the basis of analysis. The effects of the water releases from the Bagnell Dam were investigated by computing a series of safety factors in relation to outflow events for 78 cross sections along the 130-km stretch of the lower Osage River using the BSTEM algorithm integrated into the HEC-RAS model package. The results showed that the rate of change in the flow magnitude at each cross section impacts the calculated safety factor. The flow release fluctuations will change the value of the safety factor along the river over time, as the safety factor value closely follows the variation in the flow rate. Also, the results show that the cross sections closer to Bagnell Dam are subjected to be more unstable than those farther downstream of Osage River.

Missouri↗

Classifying physiographic regimes on terrain and hydrologic factors for adaptive generalization of stream networks

Automated generalization software must accommodate multi-scale representations of hydrographic networks across a variety of geographic landscapes, because scale-related hydrography differences are known to vary in different physical conditions. While generalization algorithms have been tailored to specific regions and landscape conditions by several researchers in recent years, the selection and characterization of regional conditions have not been formally defined nor statistically validated. This paper undertakes a systematic classification of landscape types in the conterminous United States to spatially subset the country into workable units, in preparation for systematic tailoring of generalization workflows that preserve hydrographic characteristics. The classification is based upon elevation, standard deviation of elevation, slope, runoff, drainage and bedrock density, soil and bedrock permeability, area of inland surface water, infiltration-excess of overland flow, and a base flow index. A seven class solution shows low misclassification rates except in areas of high landscape diversity such as the Appalachians, Rocky Mountains, and Western coastal regions.

International Journal of Cartography↗

Remote sensing of tracer dye concentrations to support dispersion studies in river channels

In river channels the flow field influences the dispersion of biota, contaminants, and other suspended or dissolved materials. Insight on patterns and rates of dispersion can be gained by injecting a pulse of visible dye and observing spatial and temporal variations in dye concentration as the pulse moves downstream. We evaluated the potential of passive optical remote sensing to enhance such tracer experiments by providing spatially distributed concentration information. During tests performed in both an experimental flume facility and a large natural channel, we made field measurements of Rhodamine WT dye concentration and above-water spectral reflectance. At Korea's River Experiment Center, a small unmanned aircraft system (sUAS) was used to acquire hyperspectral images of a sinuous outdoor flume. On the Kootenai River in northern Idaho, USA, field spectra were collected from a boat and hyperspectral image data and high resolution aerial photographs were obtained from manned aircraft. We modified an Optimal Band Ratio Analysis (OBRA) algorithm to identify wavelength combinations that yielded strong correlations between a spectrally based quantity X and dye concentration C. For both the flume and field tests, we obtained very strong (R^2 from 0.94 to 0.99) relationships between X and C across a broad range of visible wavelengths. On the Kootenai, we found that X vs. C relations derived from field spectra could be applied to airborne hyperspectral images and that dye concentrations could be estimated nearly as reliably from relatively simple three-band images as from hyperspectral data. These results imply that remote sensing could become a powerful tool for mapping dye patterns. Such a capability would advance our understanding of dispersion processes by enabling more rigorous testing of numerical flow models.

Journal of Ecohydraulics↗

Real-time control of the robotic lunar observatory telescope

The US Geological Survey operates an automated observatory dedicated to the radiometry of the Moon with the objective of developing a multispectral, spatially resolved photometric model of the Moon to be used in the calibration of Earth-orbiting spacecraft. Interference filters are used with two imaging instruments to observe the Moon in 32 passbands from 350-2500 nm. Three computers control the telescope mount and instruments with a fourth computer acting as a master system to control all observation activities. Real-time control software has been written to operate the instrumentation and to automate the observing process. The observing software algorithms use information including the positions of objects in the sky, the phase of the Moon, and the times of evening and morning twilight to decide how to observe program objects. The observatory has been operating in a routine mode since late 1995 and is expected to continue through at least 2002 without significant modifications.

Publications of the Astronomical Society of the Pa↗

Re-oligotrophy in the Upper Mississippi River, USA, occurred in just a few years

Ecological systems can undergo large changes and regime shifts that are either catastrophic, neutral, or desirable. Rivers worldwide have recently undergone desirable regime shifts related to re-oligotrophy, which is a notable and ongoing reduction in concentrations of total suspended solids (TSS), total N, total P, or phytoplankton. For example, the Upper Mississippi River, USA, has experienced major water-quality changes in multiple river reaches in recent decades. In this study, we sought to understand the timing and magnitude of re-oligotrophy in the Mississippi River over a 20-y period. We used 2 topological data analysis algorithms to address hypotheses related to the following questions: What were the order and timing of water-quality changes? What was the time period over which the major changes occurred? What was the magnitude of water-quality change before and after change points (i.e., specific years when water-quality conditions transitioned abruptly to new states)? We examined 6 water-quality state variables that defined the ecological regime for the Upper Mississippi River. In one river reach, we found that strong reductions in phytoplankton/chlorophyll a had occurred first (2008), followed by total P (2013), and last in TSS (2014). In a downriver reach, we found notable reductions for chlorophyll a (2007) but substantial increases in TSS (2013). In both reaches, the water-quality changes trended over ≥15 y, but the largest changes and a likely regime shift occurred in just 6 y. The timing (2007–2014) and range (~6 y) of water-quality changes were similar between the 2 river reaches, but the directionality of the regime shift indicated re-oligotrophy for the upstream reach and water-quality degradation for the downstream reach. Topological methods applied to long-term datasets can aid our understanding of re-oligotrophication and degradation processes and may help resource managers restore desirable regimes.

Freshwater Science↗

3D near-to-surface conductivity reconstruction by inversion of VETEM data using the distorted Born iterative method

Three-dimensional (3D) subsurface imaging by using inversion of data obtained from the very early time electromagnetic system (VETEM) was discussed. The study was carried out by using the distorted Born iterative method to match the internal nonlinear property of the 3D inversion problem. The forward solver was based on the total-current formulation bi-conjugate gradient-fast Fourier transform (BCCG-FFT). It was found that the selection of regularization parameter follow a heuristic rule as used in the Levenberg-Marquardt algorithm so that the iteration is stable.

Inverse Problems↗

Limiting the effects of earthquakes on gravitational-wave interferometers

Ground-based gravitational wave interferometers such as the Laser Interferometer Gravitational-wave Observatory (LIGO) are susceptible to ground shaking from high-magnitude teleseismic events, which can interrupt their operation in science mode and significantly reduce their duty cycle. It can take several hours for a detector to stabilize enough to return to its nominal state for scientific observations. The down time can be reduced if advance warning of impending shaking is received and the impact is suppressed in the isolation system with the goal of maintaining stable operation even at the expense of increased instrumental noise. Here, we describe an early warning system for modern gravitational-wave observatories. The system relies on near real-time earthquake alerts provided by the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA). Preliminary low latency hypocenter and magnitude information is generally available in 5 to 20 min of a significant earthquake depending on its magnitude and location. The alerts are used to estimate arrival times and ground velocities at the gravitational-wave detectors. In general, 90% of the predictions for ground-motion amplitude are within a factor of 5 of measured values. The error in both arrival time and ground-motion prediction introduced by using preliminary, rather than final, hypocenter and magnitude information is minimal. By using a machine learning algorithm, we develop a prediction model that calculates the probability that a given earthquake will prevent a detector from taking data. Our initial results indicate that by using detector control configuration changes, we could prevent interruption of operation from 40 to 100 earthquake events in a 6-month time-period.

Classical and Quantum Gravity↗

Predicting flood damage probability across the conterminous United States

Floods are the leading cause of natural disaster damages in the United States, with billions of dollars incurred every year in the form of government payouts, property damages, and agricultural losses. The Federal Emergency Management Agency oversees the delineation of floodplains to mitigate damages, but disparities exist between locations designated as high risk and where flood damages occur due to land use and climate changes and incomplete floodplain mapping. We harnessed publicly available geospatial datasets and random forest algorithms to analyze the spatial distribution and underlying drivers of flood damage probability caused by excessive rainfall and overflowing water bodies across the conterminous United States. From this, we produced the first spatially complete map of flood damage probability for the nation, along with spatially explicit standard errors for four selected cities. We trained models using the locations of historical reported flood damage events ( n = 71,434) and a suite of geospatial predictors (e.g., flood severity, climate, socio-economic exposure, topographic variables, soil properties, and hydrologic characteristics). We developed independent models for each hydrologic unit code level 2 watershed and generated a flood damage probability for each 100-m pixel. Our model classified damage or no damage with an average area under the curve accuracy of 0.75; however, model performance varied by environmental conditions, with certain land cover classes (e.g., forest) resulting in higher error rates than others (e.g., wetlands). Our results identified flood damage probability hotspots across multiple spatial and regional scales, with high probabilities common in both inland and coastal regions. The highest flood damage probabilities tended to be in areas of low elevation, in close proximity to streams, with extreme precipitation, and with high urban road density. Given rapid environmental changes, our study demonstrates an efficient approach for updating flood damage probability estimates across the nation.

Environmental Research Letters↗

Variability in terrestrial characteristics and erosion rates on the Alaskan Beaufort Sea coast

Arctic coastal environments are eroding and rapidly changing. A lack of pan-Arctic observations limits our ability to understand controls on coastal erosion rates across the entire Arctic region. Here, we capitalize on an abundance of geospatial and remotely sensed data, in addition to model output, from the North Slope of Alaska to identify relationships between historical erosion rates and landscape characteristics to guide future modeling and observational efforts across the Arctic. Using existing datasets from the Alaska Beaufort Sea coast and a hierarchical clustering algorithm, we developed a set of 16 coastal typologies that captures the defining characteristics of environments susceptible to coastal erosion. Relationships between landscape characteristics and historical erosion rates show that no single variable alone is a good predictor of erosion rates. Variability in erosion rate decreases with increasing coastal elevation, but erosion rate magnitudes are highest for intermediate elevations. Areas along the Alaskan Beaufort Sea coast (ABSC) protected by barrier islands showed a three times lower erosion rate on average, suggesting that barrier islands are critical to maintaining mainland shore position. Finally, typologies with the highest erosion rates are not broadly representative of the ABSC and are generally associated with low elevation, north- to northeast-facing shorelines, a peaty pebbly silty lithology, and glaciomarine deposits with high ice content. All else being equal, warmer permafrost is also associated with higher erosion rates, suggesting that warming permafrost temperatures may contribute to higher future erosion rates on permafrost coasts. The suite of typologies can be used to guide future modeling and observational efforts by quantifying the distribution of coastlines with specific landscape characteristics and erosion rates.

Alaska↗

Maps of active layer thickness in northern Alaska by upscaling P-band polarimetric synthetic aperture radar retrievals

Extensive, detailed information on the spatial distribution of active layer thickness (ALT) in northern Alaska and how it evolves over time could greatly aid efforts to assess the effects of climate change on the region and also help to quantify greenhouse gas emissions generated due to permafrost thaw. For this reason, we have been developing high-resolution maps of ALT throughout northern Alaska. The maps are produced by upscaling from high-resolution swaths of estimated ALT retrieved from airborne P-band synthetic aperture radar (SAR) images collected for three different years. The upscaling was accomplished by using hundreds of thousands of randomly selected samples from the SAR-derived swaths of ALT to train a machine learning regression algorithm supported by numerous spatial data layers. In order to validate the maps, thousands of randomly selected samples of SAR-derived ALT were excluded from the training in order to serve as validation pixels; error performance calculations relative to these samples yielded root-mean-square errors (RMSEs) of 7.5–9.1 cm, with bias errors of magnitude under 0.1 cm. The maps were also compared to ALT measurements collected at a number of in situ test sites; error performance relative to the site measurements yielded RMSEs of approximately 11–12 cm and bias of 2.7–6.5 cm. These data are being used to investigate regional patterns and underlying physical controls affecting permafrost degradation in the tundra biome.

Environmental Research Letters↗

Some considerations on the use of ecological models to predict species' geographic distributions

Peterson (2001) used Genetic Algorithm for Rule-set Prediction (GARP) models to predict distribution patterns from Breeding Bird Survey (BBS) data. Evaluations of these models should consider inherent limitations of BBS data: (1) BBS methods may not sample species and habitats equally; (2) using BBS data for both model development and testing may overlook poor fit of some models; and (3) BBS data may not provide the desired spatial resolution or capture temporal changes in species distributions. The predictive value of GARP models requires additional study, especially comparisons with distribution patterns from independent data sets. When employed at appropriate temporal and geographic scales, GARP models show considerable promise for conservation biology applications but provide limited inferences concerning processes responsible for the observed patterns.

Condor↗

Assessment of a claimed ultra-low frequency electromagnetic (ULFEM) earthquake precursor

The rate of occurrence of anomalous ultra-low frequency electromagnetic (ULFEM) pulses has been claimed to have increased days to weeks prior to the M5.4 2007 and M4.0 2010 Alum Rock earthquakes. We re-examine the previously reported ultra-low frequency (ULF: 0.01–10 Hz) magnetic data recorded at a QuakeFinder site located 9 km from the earthquake hypocentre, and compare to data from a nearby Stanford-USGS site located 42 km from the hypocentre, to analyse the characteristics of the pulses and assess their origin. Using pulse definitions and pulse-counting algorithms analogous to those previously reported, we corroborate the increase in pulse counts before the 2007 Alum Rock earthquake at the QuakeFinder station, but we note that the number of pulses depends on chosen temporal and amplitude detection thresholds. These thresholds are arbitrary because we lack a clear physical model or basis for their selection. We do not see the same increase in pulse counts before the 2010 Alum Rock earthquake at the QuakeFinder or Stanford-USGS stations. In addition, the majority of pulses in the QuakeFinder data and Stanford-USGS data do not match temporally, indicating the pulses lack a common origin and are not from lightning or solar-driven ionospheric/magnetospheric disturbances. Our assessment of the temporal distribution of pulse counts shows pulse counts increase during peak human activity hours, suggesting these pulses result from local cultural noise and are not tectonic in origin. The many unknowns about the character and even existence of precursory earthquake pulses means that standard numerical and statistical tests cannot easily be applied. Yet here we show that exhaustive investigation of many different aspects of ULFEM signals can be used to properly characterize their origin.

California↗

Fast rupture of the 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California dynamically triggers seismicity in California

In the Gulf of California, Mexico, the relative motion across the North America-Pacific boundary is accommodated by a series of marine transform faults and spreading centers. About 40 M>6 earthquakes have occurred in the region since 1960. On 3 August 2009, an Mw 6.9 earthquake occurred near Canal de Ballenas in the region. The earthquake was a strike-slip event with a shallow hypocenter that is likely close to the seafloor. In contrast to an adjacent M7 earthquake, this earthquake triggered a ground-motion-based earthquake early warning algorithm being tested in southern California (∼600 km away). This observation suggests that the abnormally large ground motions and dynamic strains observed for this earthquake relate to its rupture properties. To investigate this possibility, we image the rupture process and resolve the slip distribution of the event using a P-wave back-projection approach and a teleseismic, finite-fault inversion method. Results from these two independent analyses indicate a relatively simple, unilateral rupture propagation directed along-strike in the northward direction. However, the average rupture speed is estimated around 4 km/s, suggesting a possible supershear rupture. The supershear speed is also supported by a Rayleigh wave Mach cone analysis, although uncertainties in local velocity structure preclude a definitive conclusion. The Canal de Ballenas earthquake dynamically triggered seismicity at multiple sites in California, with triggering response characteristics varying from location-to-location. For instance, some of the triggered earthquakes in California occurred up to 24 hours later, suggesting that nonlinear triggering mechanisms likely have modulated their occurrence.

Baja California↗

Converted-wave reverse time migration imaging in subduction zone settings

We use a newly developed 2-D elastic reverse time migration (RTM) imaging algorithm based on the Helmholtz decomposition to test approaches for imaging the descending slab in subduction zone regions using local earthquake sources. Our elastic RTM method is designed to reconstruct incident and scattered wavefields at depth, isolate constituent P- and S- wave components via Helmholtz decomposition, and evaluate normalized imaging functions that leverage dominant P and S signals. This method allows us to target particular converted-wave scattering geometries, for example incident S to scattered P , which may be expected to have dominant signals in any given data set. The method is intended to be applied to dense seismic array observations that adequately capture both incident and converted wavefields. We draw a direct connection between our imaging functions and the first-order contrasts in shear wave material properties across seismic discontinuities. Through tests on synthetic data using either S → P or P → S conversions, we find that our technique can successfully recover the structure of a subducting slab using data from a dense wide-angle array of surface stations. We also calculate images with a small-aperture array to test the impact of array geometry on image resolution and interpretability. Our results show that our imaging technique is capable of imaging multiple seismic discontinuities at depth, even with a small number of earthquakes, but that limitations arise when a small aperture array is considered. In this case, the presence of artefacts makes it more difficult to determine the location of seismic discontinuities.

Geophysical Journal International↗

Seismotectonic framework of the 2010 February 27 M w 8.8 Maule, Chile earthquake sequence

After the 2010 M w 8.8 Maule earthquake, an international collaboration involving teams and instruments from Chile, the US, the UK, France and Germany established the International Maule Aftershock Deployment temporary network over the source region of the event to facilitate detailed, open-access studies of the aftershock sequence. Using data from the first 9-months of this deployment, we have analyzed the detailed spatial distribution of over 2500 well-recorded aftershocks. All earthquakes have been relocated using a hypocentral decomposition algorithm to study the details of and uncertainties in both their relative and absolute locations. We have computed regional moment tensor solutions for the largest of these events to produce a catalogue of 465 mechanisms, and have used all of these data to study the spatial distribution of the aftershock sequence with respect to the Chilean megathrust. We refine models of co-seismic slip distribution of the Maule earthquake, and show how small changes in fault geometries assumed in teleseismic finite fault modelling significantly improve fits to regional GPS data, implying that the accuracy of rapid teleseismic fault models can be substantially improved by consideration of existing fault geometry model databases. We interpret all of these data in an integrated seismotectonic framework for the Maule earthquake rupture and its aftershock sequence, and discuss the relationships between co-seismic rupture and aftershock distributions. While the majority of aftershocks are interplate thrust events located away from regions of maximum co-seismic slip, interesting clusters of aftershocks are identified in the lower plate at both ends of the main shock rupture, implying internal deformation of the slab in response to large slip on the plate boundary interface. We also perform Coulomb stress transfer calculations to compare aftershock locations and mechanisms to static stress changes following the Maule rupture. Without the incorporation of uncertainties in earthquake locations, just 55 per cent of aftershock nodal planes align with faults promoted towards failure by co-seismic slip. When epicentral uncertainties are considered (on the order of just ±2–3 km), 90 per cent of aftershocks are consistent with occurring along faults demonstrating positive stress transfer. These results imply large sensitivities of Coulomb stress transfer calculations to uncertainties in both earthquake locations and models of slip distributions, particularly when applied to aftershocks close to a heterogeneous fault rupture; such uncertainties should therefore be considered in similar studies used to argue for or against models of static stress triggering.

Maule↗

A comparison of methods to estimate seismic phase delays--Numerical examples for coda wave interferometry

Time-shift estimation between arrivals in two seismic traces before and after a velocity perturbation is a crucial step in many seismic methods. The accuracy of the estimated velocity perturbation location and amplitude depend on this time shift. Windowed cross correlation and trace stretching are two techniques commonly used to estimate local time shifts in seismic signals. In the work presented here, we implement Dynamic Time Warping (DTW) to estimate the warping function – a vector of local time shifts that globally minimizes the misfit between two seismic traces. We illustrate the differences of all three methods compared to one another using acoustic numerical experiments. We show that DTW is comparable to or better than the other two methods when the velocity perturbation is homogeneous and the signal-to-noise ratio is high. When the signal-to-noise ratio is low, we find that DTW and windowed cross correlation are more accurate than the stretching method. Finally, we show that the DTW algorithm has better time resolution when identifying small differences in the seismic traces for a model with an isolated velocity perturbation. These results impact current methods that utilize not only time shifts between (multiply) scattered waves, but also amplitude and decoherence measurements. DTW is a new tool that may find new applications in seismology and other geophysical methods (e.g., as a waveform inversion misfit function).

Geophysical Journal International↗