Search USGS⌕ Search

SEARCH · Search USGS

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 937 records · Page 52Linked to original sources

Parameter estimation method and updating of regional prediction equations for ungaged sites in the desert region of California

The U.S. Geological Survey (USGS) is currently updating at-site flood frequency estimates for USGS streamflow-gaging stations in the desert region of California. The at-site flood-frequency analysis is complicated by short record lengths (less than 20 years is common) and numerous zero flows/low outliers at many sites. Estimates of the three parameters (mean, standard deviation, and skew) required for fitting the log Pearson Type 3 (LP3) distribution are likely to be highly unreliable based on the limited and heavily censored at-site data. In a generalization of the recommendations in Bulletin 17B, a regional analysis was used to develop regional estimates of all three parameters (mean, standard deviation, and skew) of the LP3 distribution. A regional skew value of zero from a previously published report was used with a new estimated mean squared error (MSE) of 0.20. A weighted least squares (WLS) regression method was used to develop both a regional standard deviation and a mean model based on annual peak-discharge data for 33 USGS stations throughout California’s desert region. At-site standard deviation and mean values were determined by using an expected moments algorithm (EMA) method for fitting the LP3 distribution to the logarithms of annual peak-discharge data. Additionally, a multiple Grubbs-Beck (MGB) test, a generalization of the test recommended in Bulletin 17B, was used for detecting multiple potentially influential low outliers in a flood series. The WLS regression found that no basin characteristics could explain the variability of standard deviation. Consequently, a constant regional standard deviation model was selected, resulting in a log-space value of 0.91 with a MSE of 0.03 log units. Yet drainage area was found to be statistically significant at explaining the site-to-site variability in mean. The linear WLS regional mean model based on drainage area had a Pseudo- 2 R of 51 percent and a MSE of 0.32 log units. The regional parameter estimates were then used to develop a set of equations for estimating flows with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities for ungaged basins. The final equations are functions of drainage area.Average standard errors of prediction for these regression equations range from 214.2 to 856.2 percent.

California↗

Fractal analysis of multiscale spatial autocorrelation among point data

The analysis of spatial autocorrelation among point-data quadrats is a well-developed technique that has made limited but intriguing use of the multiscale aspects of pattern. In this paper are presented theoretical and algorithmic approaches to the analysis of aggregations of quadrats at or above a given density, in which these sets are treated as multifractal regions whose fractal dimension, D, may vary with phenomenon intensity, scale, and location. The technique is illustrated with Matui's quadrat house-count data, which yield measurements consistent with a nonautocorrelated simulated Poisson process but not with an orthogonal unit-step random walk. The paper concludes with a discussion of the implications of such analysis for multiscale geographic analysis systems. -Author

Environment and Planning A↗

Linking decomposition rates of soil organic amendments to their chemical composition

The stock of organic carbon contained within a soil represents the balance between inputs and losses. Inputs are defined by the ability of vegetation to capture and retain carbon dioxide, effects that management practices have on the proportion of captured carbon that is added to soil and the application organic amendments. The proportion of organic amendment carbon retained is defined by its rate of mineralisation. In this study, the rate of carbon mineralisation from 85 different potential soil organic amendments (composts, manures, plant residues and biosolids) was quantified under controlled environmental conditions over a 547 day incubation period. The composition of each organic amendment was quantified using nuclear magnetic resonance and mid- and near-infrared spectroscopies. Cumulative mineralisation of organic carbon from the amendments was fitted to a two-pool exponential model. Multivariate chemometric algorithms were derived to allow the size of the fast and slow cycling pools of carbon to be predicted from the acquired spectroscopic data. However, the fast and slow decomposition rate constants could not be predicted suggesting that prediction of the residence time of organic amendment carbon in soil would likely require additional information related to soil type, environmental conditions, and management practices in use at the site of application.

Soil Research↗

Cost-effective fuel treatment planning: A theoretical justification and case-study

Modelling the spatial prioritisation of fuel treatments and their net effect on values at risk is an important area for applied work as economic damages from wildfire continue to grow. We model and demonstrate a cost-effective fuel treatment planning algorithm using two ecosystem services as benefits for which fuel treatments are prioritised. We create a surface of expected fuel treatment costs to incorporate the heterogeneity in factors affecting the revenue and costs of fuel treatments, and then prioritise treatments based on a cost-effectiveness ratio to maximise the averted loss of ecosystem services from fire. We compare treatment scenarios that employ cost-effectiveness with those that do not, and use common tools and models in a case study of the Sisters Ranger District on the Deschutes National Forest in central Oregon, USA. Using cost-effectiveness not only increases the expected averted losses from fuel treatments, but it also allows a larger area to be treated for the same cost, simply by incorporating costs and cost-effectiveness into the prioritisation routine. These results have considerable implications for policymakers and land managers trying to minimise risk. Incorporating costs into the spatial planning of treatments could allow more effective outcomes without increasing fuel treatment budgets.

Oregon↗

Mapping fire regime ecoregions in California

The fire regime is a central framing concept in wildfire science and ecology and describes how a range of wildfire characteristics vary geographically over time. Understanding and mapping fire regimes is important for guiding appropriate management and risk reduction strategies and for informing research on drivers of global change and altered fire patterns. Most efforts to spatially delineate fire regimes have been conducted by identifying natural groupings of fire parameters based on available historical fire data. This can result in classes with similar fire characteristics but wide differences in ecosystem types. We took a different approach and defined fire regime ecoregions for California to better align with ecosystem types, without using fire as part of the definition. We used an unsupervised classification algorithm to segregate the state into spatial clusters based on distinctive biophysical and anthropogenic attributes that drive fire regimes – and then used historical fire data to evaluate the ecoregions. The fire regime ecoregion map corresponded well with the major land cover types of the state and provided clear separation of historical patterns in fire frequency and size, with lower variability in fire severity. This methodology could be used for mapping fire regimes in other regions with limited historical fire data or forecasting future fire regimes based on expected changes in biophysical characteristics.

California↗

Fuel models and fire potential from satellite and surface observations

A national 1-km resolution fire danger fuel model map was derived through use of previously mapped land cover classes and ecoregions, and extensive ground sample data, then refined through review by fire managers familiar with various portions of the U.S. The fuel model map will be used in the next generation fire danger rating system for the U.S., but it also made possible immediate development of a satellite and ground based fire potential index map. The inputs and algorithm of the fire potential index are presented, along with a case study of the correlation between the fire potential index and fire occurrence in California and Nevada. Application of the fire potential index in the Mediterranean ecosystems of Spain, Chile, and Mexico will be tested.

International Journal of Wildland Fire↗

A direct-gradient multivariate index of biotic condition

Multimetric indexes constructed by summing metric scores have been criticized despite many of their merits. A leading criticism is the potential for investigator bias involved in metric selection and scoring. Often there is a large number of competing metrics equally well correlated with environmental stressors, requiring a judgment call by the investigator to select the most suitable metrics to include in the index and how to score them. Data-driven procedures for multimetric index formulation published during the last decade have reduced this limitation, yet apprehension remains. Multivariate approaches that select metrics with statistical algorithms may reduce the level of investigator bias and alleviate a weakness of multimetric indexes. We investigated the suitability of a direct-gradient multivariate procedure to derive an index of biotic condition for fish assemblages in oxbow lakes in the Lower Mississippi Alluvial Valley. Although this multivariate procedure also requires that the investigator identify a set of suitable metrics potentially associated with a set of environmental stressors, it is different from multimetric procedures because it limits investigator judgment in selecting a subset of biotic metrics to include in the index and because it produces metric weights suitable for computation of index scores. The procedure, applied to a sample of 35 competing biotic metrics measured at 50 oxbow lakes distributed over a wide geographical region in the Lower Mississippi Alluvial Valley, selected 11 metrics that adequately indexed the biotic condition of five test lakes. Because the multivariate index includes only metrics that explain the maximum variability in the stressor variables rather than a balanced set of metrics chosen to reflect various fish assemblage attributes, it is fundamentally different from multimetric indexes of biotic integrity with advantages and disadvantages. As such, it provides an alternative to multimetric procedures.

Arkansas, Mississippi↗

Making Recursive Bayesian inference accessible

Bayesian models provide recursive inference naturally because they can formally reconcile new data and existing scientific information. However, popular use of Bayesian methods often avoids priors that are based on exact posterior distributions resulting from former studies. Two existing Recursive Bayesian methods are: Prior- and Proposal-Recursive Bayes. Prior-Recursive Bayes uses Bayesian updating, fitting models to partitions of data sequentially, and provides a way to accommodate new data as they become available using the posterior from the previous stage as the prior in the new stage based on the latest data. ProposalRecursive Bayes is intended for use with hierarchical Bayesian models and uses a set of transient priors in first stage independent analyses of the data partitions. The second stage of Proposal-Recursive Bayes uses the posteriors from the first stage as proposals in an MCMC algorithm to fit the full model. We combine Prior- and Proposal-Recursive concepts to fit any Bayesian model, and often with computational improvements. We demonstrate our method with two case studies. Our approach has implications for big data, streaming data, and optimal adaptive design situations.

American Statistician↗

Development and application of a mark-recapture model incorporating predicted sex and transitory behaviour

We developed an extension of Cormack-Jolly-Seber models to handle a complex mark-recapture problem in which (a) the sex of birds cannot be determined prior to first moult, but can be predicted on the basis of body measurements, and (b) a significant portion of captured birds appear to be transients (i.e. are captured once but leave the area or otherwise become ' untrappable'). We applied this methodology to a data set of 4184 serins (Serinus serinus) trapped in northeastern Spain during 1985-96, in order to investigate age-, sex-, and time-specific variation in survival rates. Using this approach, we were able to successfully incorporate the majority of ringings of serins. Had we eliminated birds not previously captured (as has been advocated to avoid the problem of transience) we would have reduced our sample sizes by >2000 releases. In addition, we were able to include 1610 releases of birds of unknown (but predicted) sex; these data contributed to the precision of our estimates and the power of statistical tests. We discuss problems with data structure, encoding of the algorithms to compute parameter estimates, model selection, identifiability of parameters, and goodness-of-fit, and make recommendations for the design and analysis of future studies facing similar problems.

Bird Study↗

Risk assessment and screening for potentially invasive fishes

Preventing the introduction of potentially invasive species is becoming more important as this worldwide problem continues to grow. The ability to predict the identity or range of potential invaders could influence regulatory decisions and help to optimally allocate resources to deal with ongoing invasions. One screening tool presented here, using species life history and environmental tolerances to identify potential invaders similar to past invaders, can be used to predict potential invading species. Another screening tool, genetic algorithms, can be used to predict the potential range of an invading species. Use and further development of tools such as these, that are quantitative and relatively transparent, would give managers and other decision makers more information for making better-informed decisions.

New Zealand Journal of Marine and Freshwater Resea↗

A case study of green tree frog population size estimation by repeated capture-mark-recapture method with individual tagging: A parametric bootstrap method vs. Jolly-Seber method

This paper deals with estimation of a green tree frog population in an urban setting using repeated capture–mark–recapture (CMR) method over several weeks with an individual tagging system which gives rise to a complicated generalization of the hypergeometric distribution. Based on the maximum likelihood estimation, a parametric bootstrap approach is adopted to obtain interval estimates of the weekly population size which is the main objective of our work. The method is computation-based; and programming intensive to implement the algorithm for re-sampling. This method can be applied to estimate the population size of any species based on repeated CMR method at multiple time points. Further, it has been pointed out that the well-known Jolly–Seber method, which is based on some strong assumptions, produces either unrealistic estimates, or may have situations where its assumptions are not valid for our observed data set.

Journal of Statistical Computation and Simulation↗

Water vapour correction of the daily 1 km AVHRR global land dataset: Part I validation and use of the Water Vapour input field

An atmospheric correction algorithm developed for the 1 km Advanced Very High Resolution Radiometer (AVHRR) global land dataset was modified to include a near real-time total column water vapour data input field to account for the natural variability of atmospheric water vapour. The real-time data input field used for this study is the Television and Infrared Observational Satellite (TIROS) Operational Vertical Sounder (TOVS) Pathfinder A global total column water vapour dataset. It was validated prior to its use in the AVHRR atmospheric correction process using two North American AVHRR scenes, namely 13 June and 28 November 1996. The validation results are consistent with those reported by others and entail a comparison between TOVS, radiosonde, experimental sounding, microwave radiometer, and data from a hand-held sunphotometer. The use of this data layer as input to the AVHRR atmospheric correction process is discussed.

International Journal of Remote Sensing↗

The use of multi-temporal Landsat Normalized Difference Vegetation Index (NDVI) data for mapping fuels in Yosemite National Park, USA

The objective of this study was to test the applicability of using Normalized Difference Vegetation Index (NDVI) values derived from a temporal sequence of six Landsat Thematic Mapper (TM) scenes to map fuel models for Yosemite National Park, USA. An unsupervised classification algorithm was used to define 30 unique spectral-temporal classes of NDVI values. A combination of graphical, statistical and visual techniques was used to characterize the 30 classes and identify those that responded similarly and could be combined into fuel models. The final classification of fuel models included six different types: short annual and perennial grasses, tall perennial grasses, medium brush and evergreen hardwoods, short-needled conifers with no heavy fuels, long-needled conifers and deciduous hardwoods, and short-needled conifers with a component of heavy fuels. The NDVI, when analysed over a season of phenologically distinct periods along with ancillary data, can elicit information necessary to distinguish fuel model types. Fuels information derived from remote sensors has proven to be useful for initial classification of fuels and has been applied to fire management situations on the ground.

International Journal of Remote Sensing↗

A multi-scale segmentation approach to filling gaps in Landsat ETM+ SLC-off images

On 31 May 2003, the Landsat Enhanced Thematic Plus (ETM+) Scan Line Corrector (SLC) failed, causing the scanning pattern to exhibit wedge-shaped scan-to-scan gaps. We developed a method that uses coincident spectral data to fill the image gaps. This method uses a multi-scale segment model, derived from a previous Landsat SLC-on image (image acquired prior to the SLC failure), to guide the spectral interpolation across the gaps in SLC-off images (images acquired after the SLC failure). This paper describes the process used to generate the segment model, provides details of the gap-fill algorithm used in deriving the segment-based gap-fill product, and presents the results of the gap-fill process applied to grassland, cropland, and forest landscapes. Our results indicate this product will be useful for a wide variety of applications, including regional-scale studies, general land cover mapping (e.g. forest, urban, and grass), crop-specific mapping and monitoring, and visual assessments. Applications that need to be cautious when using pixels in the gap areas include any applications that require per-pixel accuracy, such as urban characterization or impervious surface mapping, applications that use texture to characterize landscape features, and applications that require accurate measurements of small or narrow landscape features such as roads, farmsteads, and riparian areas.

International Journal of Remote Sensing↗

Small-footprint, waveform-resolving lidar estimation of submerged and sub-canopy topography in coastal environments

The experimental advanced airborne research lidar (EAARL) is an airborne lidar instrument designed to map near-shore submerged topography and adjacent land elevations simultaneously. This study evaluated data acquired by the EAARL system in February 2003 and March 2004 along the margins of Tampa Bay, Florida, USA, to map bare-earth elevations under a variety of vegetation types and submerged topography in shallow, turbid water conditions. A spatial filtering algorithm, known as the iterative random consensus filter (IRCF), was used to extract ground elevations from a point cloud of processed last-surface EAARL returns. Filtered data were compared with acoustic and field measurements acquired in shallow submerged (0-2.5 m water depth) and sub-canopy environments. Root mean square elevation errors (RMSEs) ranged from 10-14 cm for submerged topography to 16-20 cm for sub-canopy topography under a variety of vegetation communities. The effect of lidar sampling angles and global positioning system (GPS) satellite configuration on accuracy was investigated. Results show high RMSEs for data acquired during periods of poor satellite configuration and at large sampling angles along the edges of the lidar scan. The results presented in this study confirm the cross-environment capability of a green-wavelength, waveform-resolving lidar system, making it an ideal tool for mapping coastal environments.

International Journal of Remote Sensing↗

Synergistic use of optical and InSAR data for urban impervious surface mapping: A case study in Hong Kong

A wide range of urban ecosystem studies, including urban hydrology, urban climate, land use planning and watershed resource management, require accurate and up‐to‐date geospatial data of urban impervious surfaces. In this study, the potential of the synergistic use of optical and InSAR data in urban impervious surface mapping at the sub‐pixel level was investigated. A case study in Hong Kong was conducted for this purpose by applying a classification and regression tree (CART) algorithm to SPOT 5 multispectral imagery and ERS‐2 SAR data. Validated by reference data derived from high‐resolution colour‐infrared (CIR) aerial photographs, our results show that the addition of InSAR feature information can improve the estimation of impervious surface percentage (ISP) in comparison with using SPOT imagery alone. The improvement is especially notable in separating urban impervious surface from the vacant land/bare ground, which has been a difficult task in ISP modelling with optical remote sensing data. In addition, the results demonstrate the potential to map urban impervious surface by using InSAR data alone. This allows frequent monitoring of world's cities located in cloud‐prone and rainy areas.

International Journal of Remote Sensing↗

Global irrigated area map (GIAM), derived from remote sensing, for the end of the last millennium

A Global Irrigated Area Map (GIAM) has been produced for the end of the last millennium using multiple satellite sensor, secondary, Google Earth and groundtruth data. The data included: (a) Advanced Very High Resolution Radiometer (AVHRR) 3-band and Normalized Difference Vegetation Index (NDVI) 10 km monthly time-series for 1997-1999, (b) Syste me pour l'Observation de la Terre Vegetation (SPOT VGT) NDVI 1 km monthly time series for 1999, (c) East Anglia University Climate Research Unit (CRU) rainfall 50km monthly time series for 1961-2000, (d) Global 30 Arc-Second Elevation Data Set (GTOPO30) 1 km digital elevation data of the World, (e) Japanese Earth Resources Satellite-1 Synthetic Aperture Radar (JERS-1 SAR) data for the rain forests during two seasons in 1996 and (f) University of Maryland Global Tree Cover 1 km data for 1992-1993. A single mega-file data-cube (MFDC) of the World with 159 layers, akin to hyperspectral data, was composed by re-sampling different data types into a common 1 km resolution. The MFDC was segmented based on elevation, temperature and precipitation zones. Classification was performed on the segments. Quantitative spectral matching techniques (SMTs) used in hyperspectral data analysis were adopted to group class spectra derived from unsupervised classification and match them with ideal or target spectra. A rigorous class identification and labelling process involved the use of: (a) space-time spiral curve (ST-SC) plots, (b) brightness-greenness-wetness (BGW) plots, (c) time series NDVI plots, (d) Google Earth very-high-resolution imagery (VHRI) 'zoom-in views' in over 11 000 locations, (e) groundtruth data broadly sourced from the degree confluence project (3 864 sample locations) and from the GIAM project (1 790 sample locations), (f) high-resolution Landsat-ETM+ Geocover 150m mosaic of the World and (g) secondary data (e.g. national and global land use and land cover data). Mixed classes were resolved based on decision tree algorithms and spatial modelling, and when that did not work, the problem class was used to mask and re-classify the MDFC, and the class identification and labelling protocol repeated. The sub-pixel area (SPA) calculations were performed by multiplying full-pixel areas (FPAs) with irrigated area fractions (IAFs) for every class. A 28 class GIAM was produced and the area statistics reported as: (a) annualized irrigated areas (AIAs), which consider intensity of irrigation (i.e. sum of irrigated areas from different seasons in a year plus continuous year-round irrigation or gross irrigated areas), and (b) total area available for irrigation (TAAI), which does not consider intensity of irrigation (i.e. irrigated areas at any given point of time plus the areas left fallow but 'equipped for irrigation' at the same point of time or net irrigated areas). The AIA of the World at the end of the last millennium was 467 million hectares (Mha), which is sum of the non-overlapping areas of: (a) 252 Mha from season one, (b) 174 Mha from season two and (c) 41 Mha from continuous year-round crops. The TAAI at the end of the last millennium was 399 Mha. The distribution of irrigated areas is highly skewed amongst continents and countries. Asia accounts for 79% (370 Mha) of all AIAs, followed by Europe (7%) and North America (7%). Three continents, South America (4%), Africa (2%) and Australia (1%), have a very low proportion of the global irrigation. The GIAM had an accuracy of 79-91%, with errors of omission not exceeding 21%, and the errors of commission not exceeding 23%. The GIAM statistics were also compared with: (a) the United Nations Food and Agricultural Organization (FAO) and University of Frankfurt (UF) derived irrigated areas and (b) national census data for India. The relationships and causes of differences are discussed in detail. The GIAM products are made available through a web portal (http://www.iwmigiam.org). ?? 2009 Taylor & Francis.

International Journal of Remote Sensing↗

An initial validation of Landsat 5 and 7 derived surface water temperature for U.S. lakes, reservoirs, and estuaries

The United States Harmful Algal Bloom and Hypoxia Research Control Act of 2014 identified the need for forecasting and monitoring harmful algal blooms (HAB) in lakes, reservoirs, and estuaries across the nation. Temperature is a driver in HAB forecasting models that affects both HAB growth rates and toxin production. Therefore, temperature data derived from the U.S. Geological Survey Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus thermal band products were validated across 35 lakes and reservoirs, and 24 estuaries. In situ data from the Water Quality Portal (WQP) were used for validation. The WQP serves data collected by state, federal, and tribal groups. Discrete in situ temperature data included measurements at 11,910 U.S. lakes and reservoirs from 1980 through 2015. Landsat temperature measurements could include 170,240 lakes and reservoirs once an operational product is achieved. The Landsat-derived temperature mean absolute error was 1.34°C in lake pixels >180 m from land, 4.89°C at the land-water boundary, and 1.11°C in estuaries based on comparison against discrete surface in situ measurements. This is the first study to quantify Landsat resolvable U.S. lakes and reservoirs, and large-scale validation of an operational satellite provisional temperature climate data record algorithm. Due to the high performance of open water pixels, Landsat satellite data may supplement traditional in situ sampling by providing data for most U.S. lakes, reservoirs, and estuaries over consistent seasonal intervals (even with cloud cover) for an extended period of record of more than 35 years.

International Journal of Remote Sensing↗