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 973 records · Page 54Linked to original sources

Spatial-temporal variation of low-frequency earthquake bursts near Parkfield, California

Tectonic tremor (TT) and low-frequency earthquakes (LFEs) have been found in the deeper crust of various tectonic environments globally in the last decade. The spatial-temporal behaviour of LFEs provides insight into deep fault zone processes. In this study, we examine recurrence times from a 12-yr catalogue of 88 LFE families with ∼730 000 LFEs in the vicinity of the Parkfield section of the San Andreas Fault (SAF) in central California. We apply an automatic burst detection algorithm to the LFE recurrence times to identify the clustering behaviour of LFEs (LFE bursts) in each family. We find that the burst behaviours in the northern and southern LFE groups differ. Generally, the northern group has longer burst duration but fewer LFEs per burst, while the southern group has shorter burst duration but more LFEs per burst. The southern group LFE bursts are generally more correlated than the northern group, suggesting more coherent deep fault slip and relatively simpler deep fault structure beneath the locked section of SAF. We also found that the 2004 Parkfield earthquake clearly increased the number of LFEs per burst and average burst duration for both the northern and the southern groups, with a relatively larger effect on the northern group. This could be due to the weakness of northern part of the fault, or the northwesterly rupture direction of the Parkfield earthquake.

California↗

Using a modified time-reverse imaging technique to locate low-frequency earthquakes on the San Andreas Fault near Cholame, California

We present a new method to locate low-frequency earthquakes (LFEs) within tectonic tremor episodes based on time-reverse imaging techniques. The modified time-reverse imaging technique presented here is the first method that locates individual LFEs within tremor episodes within 5 km uncertainty without relying on high-amplitude P-wave arrivals and that produces similar hypocentral locations to methods that locate events by stacking hundreds of LFEs without having to assume event co-location. In contrast to classic time-reverse imaging algorithms, we implement a modification to the method that searches for phase coherence over a short time period rather than identifying the maximum amplitude of a superpositioned wavefield. The method is independent of amplitude and can help constrain event origin time. The method uses individual LFE origin times, but does not rely on a priori information on LFE templates and families.We apply the method to locate 34 individual LFEs within tremor episodes that occur between 2010 and 2011 on the San Andreas Fault, near Cholame, California. Individual LFE location accuracies range from 2.6 to 5 km horizontally and 4.8 km vertically. Other methods that have been able to locate individual LFEs with accuracy of less than 5 km have mainly used large-amplitude events where a P-phase arrival can be identified. The method described here has the potential to locate a larger number of individual low-amplitude events with only the S-phase arrival. Location accuracy is controlled by the velocity model resolution and the wavelength of the dominant energy of the signal. Location results are also dependent on the number of stations used and are negligibly correlated with other factors such as the maximum gap in azimuthal coverage, source–station distance and signal-to-noise ratio.

California↗

Implications of the earthquake cycle for inferring fault locking on the Cascadia megathrust

GPS velocity fields in the Western US have been interpreted with various physical models of the lithosphere-asthenosphere system: (1) time-independent block models; (2) time-dependent viscoelastic-cycle models, where deformation is driven by viscoelastic relaxation of the lower crust and upper mantle from past faulting events; (3) viscoelastic block models, a time-dependent variation of the block model. All three models are generally driven by a combination of loading on locked faults and (aseismic) fault creep. Here we construct viscoelastic block models and viscoelastic-cycle models for the Western US, focusing on the Pacific Northwest and the earthquake cycle on the Cascadia megathrust. In the viscoelastic block model, the western US is divided into blocks selected from an initial set of 137 microplates using the method of Total Variation Regularization, allowing potential trade-offs between faulting and megathrust coupling to be determined algorithmically from GPS observations. Fault geometry, slip rate, and locking rates (i.e. the locking fraction times the long term slip rate) are estimated simultaneously within the TVR block model. For a range of mantle asthenosphere viscosity (4.4 × 10 18 to 3.6 × 10 20 Pa s) we find that fault locking on the megathrust is concentrated in the uppermost 20 km in depth, and a locking rate contour line of 30 mm yr −1 extends deepest beneath the Olympic Peninsula, characteristics similar to previous time-independent block model results. These results are corroborated by viscoelastic-cycle modelling. The average locking rate required to fit the GPS velocity field depends on mantle viscosity, being higher the lower the viscosity. Moreover, for viscosity ≲ 10 20 Pa s, the amount of inferred locking is higher than that obtained using a time-independent block model. This suggests that time-dependent models for a range of admissible viscosity structures could refine our knowledge of the locking distribution and its epistemic uncertainty.

Geophysical Journal International↗

SSR_pipeline: a bioinformatic infrastructure for identifying microsatellites from paired-end Illumina high-throughput DNA sequencing data

SSR_pipeline is a flexible set of programs designed to efficiently identify simple sequence repeats (e.g., microsatellites) from paired-end high-throughput Illumina DNA sequencing data. The program suite contains 3 analysis modules along with a fourth control module that can automate analyses of large volumes of data. The modules are used to 1) identify the subset of paired-end sequences that pass Illumina quality standards, 2) align paired-end reads into a single composite DNA sequence, and 3) identify sequences that possess microsatellites (both simple and compound) conforming to user-specified parameters. The microsatellite search algorithm is extremely efficient, and we have used it to identify repeats with motifs from 2 to 25bp in length. Each of the 3 analysis modules can also be used independently to provide greater flexibility or to work with FASTQ or FASTA files generated from other sequencing platforms (Roche 454, Ion Torrent, etc.). We demonstrate use of the program with data from the brine fly Ephydra packardi (Diptera: Ephydridae) and provide empirical timing benchmarks to illustrate program performance on a common desktop computer environment. We further show that the Illumina platform is capable of identifying large numbers of microsatellites, even when using unenriched sample libraries and a very small percentage of the sequencing capacity from a single DNA sequencing run. All modules from SSR_pipeline are implemented in the Python programming language and can therefore be used from nearly any computer operating system (Linux, Macintosh, and Windows).

Journal of Heredity↗

Pasture and diurnal temperature are key predictors of regional Plains Spotted Skunk (Spilogale interrupta) distribution

The Plains Spotted Skunk ( Spilogale interrupta ) is a small carnivore native to central North America that has experienced significant population reductions, and there is a lack of information about the species that could inform conservation. Our study aimed to address knowledge gaps about the distribution and habitat associations of the species in South Dakota using species distribution modeling. We used species location data collected from state natural resource managers, trappers, and members of online social media groups dedicated to hunting and wildlife conservation; environmental predictors; and 6 predictive modeling algorithms (i.e., artificial neural networks, artificial classification tree analysis, generalized boosting models, maximum entropy, multivariate adaptive regression splines, and random forests) to develop climate and landcover ensemble distribution models. The most important climate and landcover predictors were mean temperature diurnal range (i.e., average monthly differences between daily high and low temperatures) and proportion of area classified as pasture. Ensemble model concordance identified approximately 31,300 km 2 of potential Plains Spotted Skunk habitat primarily in eastern South Dakota and between the watersheds of the Missouri and James rivers. Our results offer insights that can guide conservation and inform effective management strategies for conserving Plains Spotted Skunk populations in the northern Great Plains. The promotion of low-intensity agricultural practices such as maintaining pastures, farm buildings, fences rows, and the management of woodland encroachment may improve habitat suitability and facilitate the recovery of plains spotted skunks in the region.

South Dakota↗

Satellite imagery can predict bird species occupancy and inform multispecies management in pine savannas

Multispecies management can contribute to meeting growing challenges of preserving biodiversity, yet current game and threatened species management often focuses on individual species. Satellite imagery available at high spatial and temporal resolution provides a potential tool to overcome the challenge posed by multispecies management of linking patterns of habitat use among species. We sought to determine whether satellite imagery could be used to describe patterns of species occupancy and inform multispecies management in pine savannas in Georgia, USA. We conducted point-count surveys at 7 sites in 2022 for 3 bird species: Colinus virginianus (Northern Bobwhite), Dryobates borealis (Red-cockaded Woodpecker), and Peucaea aestivalis (Bachman’s Sparrow). We built single-season occupancy models comparing a set of models using covariates collected from field vegetation surveys and another set using covariates extracted from Sentinel-2 satellite imagery. We then used a multi-objective optimization algorithm to identify quasi-optimal management solutions (i.e., sets of covariate values from top satellite imagery metric models). We found that models created using satellite imagery performed well at predicting occupancy of all 3 species as measured by the area under the receiver operating characteristic curve (AUC > 0.8) and had higher AUC scores than field-derived habitat covariate-based models. We found combinations of metrics that could result in high rates of predicted probability of occupancy for all species (within 86% of highest possible occupancy probability), but these combinations did not exist at any of the sites. Our results demonstrate that (1) satellite imagery can allow users to build reliable occupancy models without intensive field-based vegetation surveys; and (2) C. virginianus , D. borealis, and P. aestivalis in pine savanna ecosystems could be simultaneously managed through more frequent burning, changes in canopy cover or by producing suitable heterogeneity of habitats after identifying an appropriate scale of management.

Georgia↗

Cyanobacterial bloom occurrence and emergency department visits for asthma or wheeze, Wisconsin, 2017–2019

Background: Cyanobacterial harmful algal blooms (cyanoHABs) pose risks to human and animal health. Methods: We investigated the relationship between cyanoHABs and asthma or wheeze-related emergency department (ED) visits near three Wisconsin cities (Green Bay, Madison, and Oshkosh) during 2017–2019. CyanoHAB exposure was approximated using the Cyanobacterial Assessment Network remotely sensed satellite indicator of cyanobacterial biomass, a chlorophyl algorithm (Chl BS ) aggregated by water-adjacent ZIP Code Tabulation Areas (ZCTA), and distance weighted from the nearest waterbody. Weekly counts of ED visits for asthma or wheeze were aggregated by ZCTA. Poisson generalized linear models estimated the association between the weekly number of ED visits and weekly Chl BS, adjusting for maximum temperature, dewpoint, fine particulate matter (PM 2.5 ), month, and correlation within ZCTA. Results: During 2017–2019, 7,057 ED visits for asthma or wheeze occurred in the study area (42 ZCTAs). Peaks in Chl BS occurred between June and October, with higher values in Lake Winnebago and Lake Mendota compared to Green Bay. Chl BS was not associated with ED visits for asthma or wheeze (adjusted rate ratio = 1.00, 95% confidence interval = 0.99, 1.00), and the presence of onshore winds did not change this result. Monthly aggregations of ED visits and Chl BS showed a monotonic trend between increasing Chl BS and ED visits during July–September. Conclusion: This study demonstrates the utility of remote sensing data in environmental health research. Future studies could explore individual-level exposure and outcomes to refine health risks associated with cyanoHABs.

Wisconsin↗

Ungulate personality and the human shield contribute to long-distance migration loss

Long-distance ungulate migrations are declining and past research has focused on preserving migration paths where habitat fragmentation and loss disrupts movement corridors. However, changing residency-migration tradeoffs are the stronger driver of long-distance migration loss in some populations. The human shield effect relative to predation risk and anthropogenic food resources likely shapes these tradeoffs, but individual animals also vary in their propensity to tolerate proximity to humans and developed areas. We investigated how personality relative to human-habituation affects migration behavior. We categorized elk as bold or shy based on use of anthropogenic food resources identified through a clustering algorithm applied to GPS collar data. Bold elk were 4 times more likely to select wintering areas close to human activity and migrated 60% shorter distances compared to shy elk. As a result, elk wintering grounds were spatially structured such that conflict- and disease-prone individuals selected areas adjacent to human activity. Our results suggest that bold personality traits act as a precursor to human-habituation, which permits bold elk to reap the forage and predation rewards that occur in suburban landscapes. A multi-pronged approach beyond just maintaining habitat corridors may be necessary to conserve long-distance migrations for species that can become human-habituated.

BioRxiv↗

Detection of buried targets using a new enhanced very early time electromagnetic (VETEM) prototype system

In this paper, numerical simulations of a new enhanced very early time electromagnetic (VETEM) prototype system are presented, where a horizontal transmitting loop and two horizontal receiving loops are used to detect buried targets, in which three loops share the same axis and the transmitter is located at the center of receivers. In the new VETEM system, the difference of signals from two receivers is taken to eliminate strong direct-signals from the transmitter and background clutter and furthermore to obtain a better SNR for buried targets. Because strong coupling exists between the transmitter and receivers, accurate analysis of the three-loop antenna system is required, for which a loop-tree basis function method has been utilized to overcome the low-frequency breakdown problem. In the analysis of scattering problem from buried targets, a conjugate gradient (CG) method with fast Fourier transform (FFT) is applied to solve the electric field integral equation. However, the convergence of such CG-FFT algorithm is extremely slow at very low frequencies. In order to increase the convergence rate, a frequency-hopping approach has been used. Finally, the primary, coupling, reflected, and scattered magnetic fields are evaluated at receiving loops to calculate the output electric current. Numerous simulation results are given to interpret the new VETEM system. Comparing with other single-transmitter-receiver systems, the new VETEM has better SNR and ability to reduce the clutter.

IEEE Transactions on Geoscience and Remote Sensing↗

Landsat-4/5 Band 6 relative radiometry

Relative radiometric responses for the thematic mapper (TM) band 6 data from Landsat-4 and Landsat-5 were analyzed, and an algorithm has been developed that significantly reduces the striping in Band 6 images due to detector mismatch. The TM internal calibration system as originally designed includes a DC restore circuit, which acts as a feedback system designed to keep detector bias at a constant value. There is a strong indication that the DC restore circuitry implemented in Band 6 does not function as it had been designed to. It operates as designed only during a portion of the calibration interval and not at all during acquisition of scene data. This renders the data acquired during the calibration shutter interval period virtually useless for correction of the individual responses of the four detectors in Band 6. It was observed and statistically quantified that the relative response of each of the detectors to the band average is stable over the dynamic range and throughout the lifetime of the instrument. This allows an alternate approach to relative radiometric correction of TM Band 6 images

IEEE Transactions on Geoscience and Remote Sensing↗

United states national land cover data base development 1992-2001 and beyond

An accurate, up-to-date and spatially-explicate national land cover database is required for monitoring the status and trends of the nation's terrestrial ecosystem, and for managing and conserving land resources at the national scale. With all the challenges and resources required to develop such a database, an innovative and scientifically sound planning must be in place and a partnership be formed among users from government agencies, research institutes and private sectors. In this paper, we summarize major scientific and technical issues regarding the development of the NLCD 1992 and 2001. Experiences and lessons learned from the project are documented with regard to project design, technical approaches, accuracy assessment strategy, and projecti imiplementation.Future improvements in developing next generation NLCD beyond 2001 are suggested, including: 1) enhanced satellite data preprocessing in correction of atmospheric and adjacency effect and the topographic normalization; 2) improved classification accuracy through comprehensive and consistent training data and new algorithm development; 3) multi-resolution and multi-temporal database targeting major land cover changes and land cover database updates; 4) enriched database contents by including additional biophysical parameters and/or more detailed land cover classes through synergizing multi-sensor, multi-temporal, and multi-spectral satellite data and ancillary data, and 5) transform the NLCD project into a national land cover monitoring program. ?? 2008 IEEE.

Conference Paper↗

Study on a pattern classification method of soil quality based on simplified learning sample dataset

Based on the massive soil information in current soil quality grade evaluation, this paper constructed an intelligent classification approach of soil quality grade depending on classical sampling techniques and disordered multiclassification Logistic regression model. As a case study to determine the learning sample capacity under certain confidence level and estimation accuracy, and use c-means algorithm to automatically extract the simplified learning sample dataset from the cultivated soil quality grade evaluation database for the study area, Long chuan county in Guangdong province, a disordered Logistic classifier model was then built and the calculation analysis steps of soil quality grade intelligent classification were given. The result indicated that the soil quality grade can be effectively learned and predicted by the extracted simplified dataset through this method, which changed the traditional method for soil quality grade evaluation. ?? 2011 IEEE.

Conference Paper↗

SBAS-InSAR analysis of surface deformation at Mauna Loa and Kilauea volcanoes in Hawaii

We investigate the deformation of Mauna Loa and Kllauea volcanoes, Hawai'i, by exploiting the advanced differential Synthetic Aperture Radar Interferometry (InSAR) technique referred to as the Small BAseline Subset (SBAS) algorithm. In particular, we present time series of line-of-sight (LOS) displacements derived from SAR data acquired by the ASAR instrument, on board the ENVISAT satellite, from the ascending (track 93) and descending (track 429) orbits between 2003 and 2008. For each coherent pixel of the radar images we compute time-dependent surface displacements as well as the average LOS deformation rate. Our results quantify, in space and time, the complex deformation of Mauna Loa and Kllauea volcanoes. The derived InSAR measurements are compared to continuous GPS data to asses the quality of the SBAS-InSAR products.

Conference Paper↗

Single receiver target localization in mobile marine acoustic telemetry

Many marine ecological studies track animal movement and migration using acoustic telemetry, in which animals are implanted with acoustic transmitters broadcasting a unique ID code. Receivers detect and decode these transmissions and provide an estimate of the animal’s location. As transmissions can be detected at distances of hundreds of meters to even several kilometers, localization performed using binary detection information is very coarse. Fine-scale localization can be performed using time difference of arrival systems, but such systems are costlier and require a greater number of receivers. We present a Bayesian method for target localization using only binary detection information collected by a single receiver mounted on a mobile platform. The approach reduces error in the location estimate by considering the distance-dependent probability of detection (detection efficiency curve), which must first be characterized. Stationary receivers are used to infer failed detections. This localization method is then used to inform waypoint selection for the mobile platform, using the second nondominated sorting genetic algorithm considering the tradeoff between the travel distance of the platform and localization error. Simulation results examine the effects of detection efficiency curve characteristics, number of transmissions, and travel distance on the localization error. The proposed method is further supported by experimental results from field trials conducted with an autonomous surface vehicle in Lake Huron.

IEEE Journal of Oceanic Engineering↗

Description of the U.S. Geological Survey Geo Data Portal data integration framework

The U.S. Geological Survey has developed an open-standard data integration framework for working efficiently and effectively with large collections of climate and other geoscience data. A web interface accesses catalog datasets to find data services. Data resources can then be rendered for mapping and dataset metadata are derived directly from these web services. Algorithm configuration and information needed to retrieve data for processing are passed to a server where all large-volume data access and manipulation takes place. The data integration strategy described here was implemented by leveraging existing free and open source software. Details of the software used are omitted; rather, emphasis is placed on how open-standard web services and data encodings can be used in an architecture that integrates common geographic and atmospheric data.

IEEE Journal of Selected Topics in Applied Earth O↗

Remote sensing with simulated unmanned aircraft imagery for precision agriculture applications

An important application of unmanned aircraft systems (UAS) may be remote-sensing for precision agriculture, because of its ability to acquire images with very small pixel sizes from low altitude flights. The objective of this study was to compare information obtained from two different pixel sizes, one about a meter (the size of a small vegetation plot) and one about a millimeter. Cereal rye (Secale cereale) was planted at the Beltsville Agricultural Research Center for a winter cover crop with fall and spring fertilizer applications, which produced differences in biomass and leaf chlorophyll content. UAS imagery was simulated by placing a Fuji IS-Pro UVIR digital camera at 3-m height looking nadir. An external UV-IR cut filter was used to acquire true-color images; an external red cut filter was used to obtain color-infrared-like images with bands at near-infrared, green, and blue wavelengths. Plot-scale Green Normalized Difference Vegetation Index was correlated with dry aboveground biomass ( ${mbi {r}} = 0.58$ ), whereas the Triangular Greenness Index (TGI) was not correlated with chlorophyll content. We used the SamplePoint program to select 100 pixels systematically; we visually identified the cover type and acquired the digital numbers. The number of rye pixels in each image was better correlated with biomass ( ${mbi {r}} = 0.73$ ), and the average TGI from only leaf pixels was negatively correlated with chlorophyll content ( ${mbi {r}} = -0.72$ ). Thus, better information for crop requirements may be obtained using very small pixel sizes, but new algorithms based on computer vision are needed for analysis. It may not be necessary to geospatially register large numbers of photographs with very small pixel sizes. Instead, images could be analyzed as single plots along field transects.

IEEE Journal of Selected Topics in Applied Earth O↗

TreeMAC: Localized TDMA MAC protocol for real-time high-data-rate sensor networks

Earlier sensor network MAC protocols focus on energy conservation in low-duty cycle applications, while some recent applications involve real-time high-data-rate signals. This motivates us to design an innovative localized TDMA MAC protocol to achieve high throughput and low congestion in data collection sensor networks, besides energy conservation. TreeMAC divides a time cycle into frames and frame into slots. Parent determines children's frame assigmnent based on their relative bandwidth demand, and each node calculates its own slot assignment based on its hop-count to the sink. This innovative 2-dimensional frame-slot assignment algorithm has the following nice theory properties. Firstly, given any node, at any time slot, there is at most one active sender in its neighborhood (includ ing itself). Secondly, the packet scheduling with TreelMAC is bufferless, which therefore minimizes the probability of network congestion. Thirdly, the data throughput to gateway is at least 1/3 of the optimum assuming reliable links. Our experiments on a 24 node test bed demonstrate that TreeMAC protocol significantly improves network throughput and energy efficiency, by comparing to the TinyOS's default CSMA MAC protocol and a recent TDMA MAC protocol Funneling-MAC[8]. ?? 2009 IEEE.

Conference Paper↗

Landsat-5 bumper-mode geometric correction

The Landsat-5 Thematic Mapper (TM) scan mirror was switched from its primary operating mode to a backup mode in early 2002 in order to overcome internal synchronization problems arising from long-term wear of the scan mirror mechanism. The backup bumper mode of operation removes the constraints on scan start and stop angles enforced in the primary scan angle monitor operating mode, requiring additional geometric calibration effort to monitor the active scan angles. It also eliminates scan timing telemetry used to correct the TM scan geometry. These differences require changes to the geometric correction algorithms used to process TM data. A mathematical model of the scan mirror's behavior when operating in bumper mode was developed. This model includes a set of key timing parameters that characterize the time-varying behavior of the scan mirror bumpers. To simplify the implementation of the bumper-mode model, the bumper timing parameters were recast in terms of the calibration and telemetry data items used to process normal TM imagery. The resulting geometric performance, evaluated over 18 months of bumper-mode operations, though slightly reduced from that achievable in the primary operating mode, is still within the Landsat specifications when the data are processed with the most up-to-date calibration parameters.

IEEE Transactions on Geoscience and Remote Sensing↗